作者askooen (刺刺的凤梨)
看板MATLAB
标题[讨论] 高次方常系数多项式求解
时间Fri Jun 6 00:27:09 2014
最近在运算高次方常系数多项式的解，先用syms建构
函数的最高次方到12次，而用solve指令解出来的12个解。
如果带回原多项式应该要几乎等於零才是。
但是将解代回去原多项式之後，最夸张是等於10^-4左右，最好在10^-13
其它大部分的解都是落在10^-11~10^-12之间。
但我想要的误差是希望代回去的值大约在10^-15~10^-16之间。
请问这有什麽解决方法吗？或是有什麽方法可以让他精度提高呢？
烦请各位回答了，谢谢。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 147.188.254.213
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1401985633.A.9BC.html
1F：→ wxtab019:看你带回去的时候 小数点是不是被四舍五入掉了? 06/06 00:48
2F：推 alamabarry:解析解代入数值很容易有误差值阿 所以要用数值解比较准 06/06 02:46
3F：→ alamabarry:直接用roots就好了 06/06 02:47
4F：→ vaakaa:因为已经接近解了,再写一个解根的程式去"refine"这个根也行 06/07 20:58
5F：→ vaakaa:当然 要先确认函数计算本身的误差严不严重 06/07 21:00
6F：推 tomichy:还不够接近根吗？ 有请 牛顿拉夫森先生出来射个两枪看看 07/03 14:24