| 查看: 2872 | 回复: 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完全不正确。是代码的哪一部分出了问题呢？ |
回复此楼» 猜你喜欢
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
|
本帖内容被屏蔽 |
2楼2018-01-09 22:46:15
20
投票:
