vettelkimi 发表于 2012-12-3 20:50

请教多项式求根遇到的问题

今天,求解一个多项式的根,但是将解出来的数值回带到原多项式后发现其数值不为0,下面为原程序:
x=,
y=roots(x),
得到的解为:24.7525   12.4607    -1.5868 + 1.9261i-1.5868 - 1.9261i   0.2783 + 1.8979i0.2783 - 1.8979i   0.2841 + 0.7020i   0.2841 - 0.7020i0.0538 + 0.5245i0.0538 - 0.5245i0.1250 + 0.4995i 0.1250 - 0.4995i
将24.7525回带:polyval(x,24.7525),得到的值为: -1.6360e+061
请大大们帮忙看下这是为什么。用maple也解了一下,发现答案是一样的。

happy 发表于 2012-12-4 09:38

你的多项系数太大了,应该是舍入误差造成的
参考http://www.mathworks.cn/support/solutions/en/data/1-15ROS/index.html?solution=1-15ROS

vettelkimi 发表于 2012-12-4 09:49

happy 发表于 2012-12-4 09:38 static/image/common/back.gif
你的多项系数太大了,应该是舍入误差造成的
参考http://www.mathworks.cn/support/solutions/en/data/1-15 ...

谢谢了,这个多项式是一个7维含参的行列式展开得到的,不知道有没有什么好的解法?

happy 发表于 2012-12-4 10:06

vettelkimi 发表于 2012-12-4 09:49 static/image/common/back.gif
谢谢了,这个多项式是一个7维含参的行列式展开得到的,不知道有没有什么好的解法?

找其他算法看看,比如二分法之类的

chengjinm 发表于 2012-12-4 17:40

换算法吧
页: [1]
查看完整版本: 请教多项式求根遇到的问题