R程序做牛顿迭代和定点迭代，求帮忙解决问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

y=function(x){x^3-x-1}  #所求函数#
z=function(x){3*x^2-1}  #导函数#
#下面是牛顿迭代法#
Newton=function(f,z,a,b,TOL,N){
n=0
x=0
while(n<=N){   
    xx=x-y(x)/z(x)
    n=n+1   
     if(abs(xx-x)>TOL)  x=xx
     }
  list(boot=xx,fun=f(x),n)
}

r0=fpi(y,z,1,2,1e-5,1e+5)；r0
运行出现的问题是：错误于fpi(y, z, 1, 2, 1e-05, 1e+05) : 参数((1e+05)) 没有用


#定点迭代#
fpi=function(f,a,b,TOL,N0){n=0
while( abs(a-b)>TOL & n<=N0 ) {   
a=(a+b)/2   #此处让a变成初始值X0，这样做我是为了使上面的循环abs(a-b)方便#  
b=f(a)        #让b代表成新生成的x，即对应Xn的Xn+1#  
n=n+1  
}
list(boot=b,fun=f(x),n)
}
r0=fpi(f,1,2,1e-5,1e+5)r0
出现的问题是：
错误于while (abs(a - b) > TOL & n <= N0) { :   
需要TRUE/FALSE值的地方不可以用缺少值
#这里我知道是算出来的结果发散了，而不是收敛的，不知道要怎么改#
这个程序我还用二分法写过，二分法运行正常，希望论坛里有人能帮下忙。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：解决问题 R程序 求帮忙 function Newton 程序

如果要用迭代法估计参数或者最优化，不必用R那么麻烦，1stopt操作简单，效率很高，论坛可以下载1.5版本。

xgxhf 发表于 2012-3-19 19:18
如果要用迭代法估计参数或者最优化，不必用R那么麻烦，1stopt操作简单，效率很高，论坛可以下载1.5版本。

谢谢，我现在是学习R的过程，不是单纯的为了做题。

前面的Newton没有问题，后面那个函数里面b=f(a)是干嘛的？

qoiqpwqr 发表于 2012-3-19 20:49
前面的Newton没有问题，后面那个函数里面b=f(a)是干嘛的？

也就相当于迭代啊    while循环里的a b不再是上面赋值a b的意思了。说白了，也就是相当于令x0=a，b=x1=f(x0)=f(a)。。。然后再这样迭代下去

我也遇到了同样的问题，只是我求的极值点是3维的，函数有点复杂。我带入了几个点算，梯度向量比较大了，不知道是不是迭代到最后发散到很大无法计算了？总是报错：需要TRUE/FALSE值的地方不可以用缺少值。。。不知道你是否已经解决？

楼主你好，我也碰到类似的问题，在做迭代时候 while((eui[m])^2>M&n<N)报错“需要TRUE/FALSE值的地方不可以用缺少值”。你是怎么解决的？

第一个报错原因，fpi只有五个参数，而你给了六个，所以出错了

第二个程序的问题我不知道是什么原因。但是你把两个程序分开运行，注意定义最后的函数f，就没有问题了。

江梦海凝 发表于 2013-9-13 10:49
楼主你好，我也碰到类似的问题，在做迭代时候 while((eui[m])^2>M&n

我昨天遇到过类似的问题，反复查程序，是因为需要TRUE或FALSE的时候，我给的是空值。不知道你遇到的问题和我的类似否。R中好像没有检查数组下标过界，反正我昨天的程序出问题，是因为下标过界，而R说，是我没有给他老人家TRUE或FALSE的值，哈哈