R语言：求解方程根公式uniroot

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求解方程根公式uniroot(f,interval,lower=min(interval),upper=max(interval),tol=,maxiter=,)

其中，f表示要求解方程，interval表示包含方程根的初始区间，lower表示interval中方程根初始区间的左端点，upper表示右端点，tol是想要达到的计算精度，maxiter是设置的最大迭代次数。

R：

> x<-rcauchy(1000,1)#表示的是Cauchy分布产生的1000个，每个值相距为1的数，见下面

> x

   [1]  2.116535e+00  1.820128e+00 -1.793299e+01  2.126767e+00  1.639780e+00

   [6] -1.611596e+01 -4.004089e+00  1.429597e+00  2.040831e+00 -2.248832e-01

  [11]  8.781683e-01  1.325409e+00  7.981114e-02  8.350597e-01  3.593475e-01

。。。。。。还有很多值，一共一千个。

> f<-function(p) sum((x-p)/(1+(x-p)^2))#Cauchy分布的方程，其中x的取值就是上面产生的随机数，目标是为了求P的值。

> out<-uniroot(f,c(0,5)) #利用这个函数就接P值，设置P值的区间为0到5

> out

$root    #求解出来的方程近似解

[1] 1.077677

$f.root  #在此近似解下的方程值

[1] -0.0002070086

$iter   #求解过程中，方程迭代的次数

[1] 6

$init.it

[1] NA

$estim.prec   #近似解与精确解的误差估计，即近似解与精确解之间误差的绝对值不超过下面的数值

[1] 6.103516e-05

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：root 解方程 R语言 Interval inter

谢谢指导

谢谢哦

请问uniroot函数里面c(,)这个上下界应该怎么确定呢？没改变一次结果都不一样，结果经常很差，函数值和0相差很大，很烦恼，求指教~