Question

基于有限差分法求解的二维谐振子的MATLAB程序如下。哪位大神能帮我做个注明啊，完全看不懂啊，，急

tic
clc
clear
L=20;
W=20;
N=20;
M=20;
hx=L/(2*N);
hy=W/(2*M);
S=zeros((2*M-1)*(2*N-1));
for m=1:2*M-1
D(m,m)=-1/(hy.^2);
end
for m=1:2*N-2
mx=(2*M-1)*(m-1)+1;
my=(2*M-1)*(m-1)+2*M-1;
nx=(2*M-1)*(m-1)+(2*M-1)+1;
ny=(2*M-1)*(m-1)+2*(2*M-1);
S(mx:my,nx:ny)=D;
S(nx:ny,mx:my)=D;
end
for n=1:2*N-1
C=zeros(2*M-1);
for m=1:2*M-1
C(m,m)=2*(1/hx^2+1/hy^2)+((n-M)*hx)^2+((m-N)*hy)^2;
end
for m=1:2*M-2
C(m,m+1)=-1/hx.^2;
C(m+1,m)=-1/hx.^2;
end
cx=(2*M-1)*(n-1)+1;
cy=(2*M-1)*(n-1)+2*M-1;
S(cx:cy,cx:cy)=C;
end
[V,E]=eig(S); Y=V(:,2); k=1;
for j=1:2*M-1;
for i=1:2*N-1;
ZZ(i,j)=Y(k);
k=k+1;
end
end
y=linspace(-W/2,W/2,2*N-1);
x=linspace(-L/2,L/2,2*M-1);
[x1,y1]=meshgrid(x,y);
surf(x1,y1,-ZZ);
shading interp;
zlabel('¦•');
toc

Answer 1

你是有限差分法不懂，还是matlab语言不懂？前者不懂可以去看数值计算方面的书籍，后者不懂得先去学matlab。