이 재귀 프로그램을 코딩하는 데 어려움을 겪고 있으며 누구든지 도울 수 있는지 궁금합니다.
이 재귀 방정식을 코딩하고 싶습니다.
for k=1,2,...
beta(k)=k-sum_(i=0)^(k-1)Kchoosei*beta(i)*exp(-i(k-i))
R에서 수동 방식으로 수행했지만 함수에 넣고 싶습니다.
beta0<-0
beta1<-1-choose(1,0)*beta0*exp(-0*lambdaL*(1-0))
beta2<-2-choose(2,0)*beta0*exp(-0*lambdaL*(2-0))-choose(2,1)*beta1*exp(-1*lambdaL*(2-1))
beta3<-3-choose(3,0)*beta0*exp(-0*lambdaL*(3-0))-choose(3,1)*beta1*exp(-1*lambdaL*(3-1))-choose(3,2)*beta2*exp(-2*lambdaL*(3-2))
beta4<-4-choose(4,0)*beta0*exp(-0*lambdaL*(4-0))-choose(4,1)*beta1*exp(-1*lambdaL*(4-1))-choose(4,2)*beta2*exp(-2*lambdaL*(4-2))-choose(4,3)*beta3*exp(-3*lambdaL*(4-3))
합계에 대한 두 번째 루프를 정의 할 수 있습니다. 여기서 인덱싱 은 "인덱스 이동"으로 이어지는 1
것이 아니라로 시작합니다 0
.
beta = numeric()
beta[1] <- 0
for (k in 1:10){
beta[k+1] <- k
for (i in 0:(k-1))
beta[k+1] <- beta[k+1] - choose(k, i)*beta[i+1]*exp(-i*(k-i))
}
beta
# [1] 0.000000 1.000000 1.264241 2.080705 3.247551 4.528104 5.748673
# [8] 6.876234 7.941197 8.972749 9.987645
이 기사는 인터넷에서 수집됩니다. 재 인쇄 할 때 출처를 알려주십시오.
침해가 발생한 경우 연락 주시기 바랍니다[email protected] 삭제
몇 마디 만하겠습니다