下面迭代程序显示参数长度为0，怎么改

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

theta <- .5
iter=2
hattheta=c()
while ((hattheta[iter]-hattheta[iter-1])^2>10^(-2)){
        hattheta[1]<-theta
        hattheta[iter]=(159*hattheta[iter-1]+68)/(197*hattheta[iter-1]+144)
        print(hattheta)
        iter=iter+1

}

显示错误
Error in while ((hattheta[iter] - hattheta[iter - 1])^2 > 10^(-2)) { :
  参数长度为零

迭代公式为：

\[]\theta^{t+1}=\frac{159\theta^{t}+68}{197\theta^{t}+144}\]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Theta LaTeX Late atex The Error 程序

仔细体会一下代码的运行，你的while循环走不下去啊。

1. theta <- .5
   
2. iter=2
   
3. hattheta=c()
   
4. hattheta[1]<-theta
   
5. hattheta[iter]=(159*hattheta[iter-1]+68)/(197*hattheta[iter-1]+144)
   
6. # 前方判断，必须要有初始值，前两个，所以要先计算好
   
7. while ((hattheta[iter]-hattheta[iter-1])^2 >10^(-3)){
   
8. #这里，需要计算下一个值了
   
9.         hattheta[iter+1]=(159*hattheta[iter-1]+68)/(197*hattheta[iter-1]+144)
   
10.         print(hattheta[iter+1])
    
11. # 把循环体往前推进一个单位
    
12.         iter=iter+1
    
13. }

复制代码

jiangbeilu 发表于 2016-3-27 11:18
仔细体会一下代码的运行，你的while循环走不下去啊。

但这样出来的迭代结果有问题0.6082474

hattheta
[1] 0.5000000 0.6082474 0.6243211 0.6264889 0.6267773 0.6268156 0.6268207
[8] 0.6268214 0.6268215 0.6268215 0.6268215 0.6268215 0.6268215 0.6268215
[15] 0.6268215 0.6268215 0.6268215 0.6268215 0.6268215 0.6268215 0.6268215
>
这是我用for做20次循环做出来的可以看出稳定在0.628215