[討論] 用迴圈跑sigma 但矩陣有上限?
小弟利用輸入熱傳導解析解
想用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)
推
10/15 14:09, , 1F
10/15 14:09, 1F
→
10/15 16:25, , 2F
10/15 16:25, 2F
→
10/16 00:08, , 3F
10/16 00:08, 3F
→
10/16 00:09, , 4F
10/16 00:09, 4F
→
10/16 20:44, , 5F
10/16 20:44, 5F