poisson equ with neumann boundary question# Computation - 科学计算
l*n
1 楼
写了个matlab的code,
但是结果老是有问题，问题出在边值上。我直接用的central difference,然后带入
方程，再用dct(discrete cosine transform)来解。我做了几个例子，但是都是边
值有问题。
比如那个u=xy的例子，边值如果是
for i=1:n+1
f(1,i)=f(1,i)-dy*(i-1)*dx;
f(i,1)=f(i,1)-dx*(i-1)*dx;
f(n+1,i)=f(n+1,i)+dx*(i-1)*dy;
f(i,n+1)=f(i,n+1)+dx*(i-1)*dy;
end
结果就是正确的。我不知道为什么没有2，因为时central difference.其他几个例
子这是边上的数值有问题。
那位大侠，有时间帮忙看一下，我始终看不出问题。谢谢
% -lapla(u)=f
clear all;
n=input('Please enter dimension:');
dx=1/n;
dy=dx;
x=0:dx:1;
y=0:dy:1;
%计算特征值
for i=1:n+1
v(i)=2*(
但是结果老是有问题，问题出在边值上。我直接用的central difference,然后带入
方程，再用dct(discrete cosine transform)来解。我做了几个例子，但是都是边
值有问题。
比如那个u=xy的例子，边值如果是
for i=1:n+1
f(1,i)=f(1,i)-dy*(i-1)*dx;
f(i,1)=f(i,1)-dx*(i-1)*dx;
f(n+1,i)=f(n+1,i)+dx*(i-1)*dy;
f(i,n+1)=f(i,n+1)+dx*(i-1)*dy;
end
结果就是正确的。我不知道为什么没有2，因为时central difference.其他几个例
子这是边上的数值有问题。
那位大侠，有时间帮忙看一下，我始终看不出问题。谢谢
% -lapla(u)=f
clear all;
n=input('Please enter dimension:');
dx=1/n;
dy=dx;
x=0:dx:1;
y=0:dy:1;
%计算特征值
for i=1:n+1
v(i)=2*(