[問題] matlab可執行的位數計算

作者: etu4315 (成大海龍王)   2015-10-18 04:53:44
各位高手神人大家好
小弟菜比八,為了這個作業弄到快天亮
google了許久仍然沒能找到相關資料
故來版上尋求神人解惑
小弟目前碰到的問題是要解一個很簡單的五次方函數
從一階微分到四階微分
使用中央差分法來計算數值解跟解析解的誤差
中央差分法的公式如圖所示
http://imgur.com/QNDZHOp
上為四階CDS 下為三階CDS
當小弟設定dx為10^-10次方時
代進去五次方函數時相當於10^-50次方
太小了以至於電腦完全省略error都跑出100%
想請問這是不是電腦的極限,抑或是matlab的極限
有方法可以解決嗎?求神人開示了QQ
作者: tn00364361 (小氵斿)   2015-10-18 12:26:00
MATLAB的machine epsilon是2.2e-16

Links booklink

Contact Us: admin [ a t ] ucptt.com