| 查看: 2892 | 回复: 1 | |||
[求助]
解泊松方程
|
|
想用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完全不正确。是代码的哪一部分出了问题呢？ |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
|
本帖内容被屏蔽 |
2楼2018-01-09 22:46:15
