作者billiron (驚驚驚驚驚)
看板MATLAB
標題[討論] 用迴圈跑sigma 但矩陣有上限?
時間Sun Oct 14 22:25:18 2012
小弟利用輸入熱傳導解析解
想用sigma做傅立葉級數
程式如下
function exactcond = T( x,y )
a=1;
b=1;
for n=1:226
t(1,n)=(2/(a*sinh(n*pi*b/a)))*(a/(n*pi))
*sinh(n*pi*y./a).*sin(n*pi*x./a)*(1-(-1)^n);
end
exactcond=sum(t);
end
照理來說n越大會越接近解析解
但是當我輸入到n=226之後
出現的都是Nan
想請問是出了什麼問題??
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 1.162.248.46
※ 編輯: billiron 來自: 140.114.56.30 (10/15 13:13)
1F：推 sunev:分母爆掉了 10/15 14:09
2F：→ billiron:那想請問一下，我應該如何修正才好?? 因為精準度還不夠 10/15 16:25
3F：→ sunev:太大的值和太小的值都是行不通的，如果你確定最後的值 10/16 00:08
4F：→ sunev:沒那麼大也沒那麼小，那你得自己先把這些值給手算抵消掉 10/16 00:09
5F：→ billiron:感謝你!! 先放大再除回去好像就可以大一點點了 10/16 20:44