[讨论] 高次方常系数多项式求解

楼主: askooen (凤梨)   2014-06-06 00:27:09
最近在运算高次方常系数多项式的解,先用syms建构
函数的最高次方到12次,而用solve指令解出来的12个解。
如果带回原多项式应该要几乎等于零才是。
但是将解代回去原多项式之后,最夸张是等于10^-4左右,最好在10^-13
其它大部分的解都是落在10^-11~10^-12之间。
但我想要的误差是希望代回去的值大约在10^-15~10^-16之间。
请问这有什么解决方法吗?或是有什么方法可以让他精度提高呢?
烦请各位回答了,谢谢。
作者: wxtab019 (天霜凝月)   2014-06-06 00:48:00
看你带回去的时候 小数点是不是被四舍五入掉了?
作者: alamabarry (.............)   2014-06-06 02:46:00
解析解代入数值很容易有误差值阿 所以要用数值解比较准直接用roots就好了
作者: vaakaa (我要欧趴)   2014-06-07 20:58:00
因为已经接近解了,再写一个解根的程式去"refine"这个根也行当然 要先确认函数计算本身的误差严不严重

Links booklink

Contact Us: admin [ a t ] ucptt.com