24小时热门版块排行榜    

查看: 2769  |  回复: 1

MattChen1316

新虫 (初入文坛)

[求助] 解泊松方程

想用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完全不正确。是代码的哪一部分出了问题呢?
回复此楼

» 猜你喜欢

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

FungEcho

禁虫 (正式写手)

本帖内容被屏蔽

2楼2018-01-09 22:46:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MattChen1316 的主题更新
信息提示
请填处理意见