24小时热门版块排行榜

返回列表

MattChen1316

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4.5
帖子: 1
在线: 57分钟
虫号: 2778480
注册: 2013-11-05
专业: 偏微分方程

[求助] 解泊松方程

想用2D傅里叶变换解个泊松方程Uxx+Uyy=f(x,y).
离散化以后得到

算法就三步： 1. 对f(x,y)做傅里叶变换得到f^(x,y)
2. 结果乘以公式里的系数
3. 做傅里叶逆变换得到数值解u(x,y)

n = 100;
h = 1/(n+1);

x=0:h:1;
y=0:h:1;

[j,k]=meshgrid(x,y);
f = 6*j+6*k;

lambda=2*(1-cos(2*(0:n+1)*pi/(n+1)));
lambda_rec=1./(eps+(lambda'*ones(1,n+2)+ones(n+2,1)*lambda));
fft2u = h*h*lambda_rec.*fft2(f);
u = ifft2(fft2u);
求得的u完全不正确。是代码的哪一部分出了问题呢？

回复此楼

» 猜你喜欢

1楼 2015-03-15 02:41:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FungEcho

禁虫 (正式写手)

本帖内容被屏蔽

2楼2018-01-09 22:46:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 MattChen1316 的主题更新

返回列表