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