R求解方程

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位老师、大神。
                        小弟在本版块里浏览很久得知。新的“BB”包可以计算3次2元方程组。但好像不适用于小弟的数据。
                        经过变换，已将问题转化为求解一个一元四次方程，除了利用求根公式求得以外，R中有什么好的方法可以解算出
                        所有         值么？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：解方程 方程组 方程组

uniroot.all

danica33 发表于 2014-4-3 09:22
uniroot.all

感谢指导，已经下载rootsolve包并进行调试。但是发现，有个interval参数设置很有技巧啊。这个不是指搜索的区间范围么？为什么搜索的区间比较大的时候反而求不出所有的解了呢？

1. > uniroot.all(function(x) x^4-10*x^3+35*x^2-50*x+24,c(0,5))
   
2. [1] 1 2 3 4
   
3. > uniroot.all(function(x) x^4-10*x^3+35*x^2-50*x+24,c(-1,10))
   
4. [1] 0.9999999 2.0000189 3.0000053 3.9999991
   
5. > uniroot.all(function(x) x^4-10*x^3+35*x^2-50*x+24,c(-10,10))
   
6. [1] 1 2 3 4
   
7. > uniroot.all(function(x) x^4-10*x^3+35*x^2-50*x+24,c(-100,100))
   
8. [1] 2 4
   
9. >

复制代码

在实际使用时，应该怎么去设置这个范围呢？

一元高次方程用迭代就行了

cyoubo5 发表于 2014-4-3 12:06
感谢指导，已经下载rootsolve包并进行调试。但是发现，有个interval参数设置很有技巧啊。这个不是指搜索的 ...

你看一下帮助?uniroot.all

It is not guaranteed that all roots will be recovered.
This will depend on n, the number of subintervals in which the interval is divided.

区间长度较大，把参数n设置大一点，默认是100.

求证1加1 发表于 2014-4-3 13:22
一元高次方程用迭代就行了

不好意思哈。
这个方面我比较比较薄弱。学R和编程都不是很久。。。
您所说的迭代解是什么意思？或者R中应该怎么做呢？？

cyoubo5 发表于 2014-4-3 18:30
不好意思哈。
这个方面我比较比较薄弱。学R和编程都不是很久。。。
您所说的迭代解是什么意思？或者R中 ...

转换成一个优化问题 然后用optim()来做咯