[討論] 高次方常係數多項式求解

作者: 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