二維差分緊致差分(matlab求解泊松方程)-wenkub
2023-07-08 17:16:40 本頁(yè)面
【正文】 =N1 A(i,i1)=1/h^2。%%% Stiff matrix A=zeros((N1)^2,(N1)^2)。h=1/N。S=h^2。for i=1 A(i,i)=4/h^2。 A(i,i)=4/h^2。 A(i,i+1)=1/h^2。endfor n=2:N2 i=(N2)*(N1)+n。 A(i,i+1)=1/h^2。 A(i,i+(N1))=1/h^2。 A(i,i+1)=1/h^2。 A(i,i1)=1/h^2。% A(i,i(N1))=1/h^2。% endfor m=1:N3 for n=2:N2i=m*(N1)+n。A(i,i+1)=1/h^2。F(i)=2*pi^2*sin(pi*n*h)*sin(pi*(m+1)*h)。for m=1:N1 U1(m*(N+1)+2:m*(N+1)+N)=U((m1)*(N1)+1:m*(N1))。 %L2 errorErrL2err%%%% 真解u=sin(pi*x)*sin(pi*y) %%%%%%% 方程 Laplace(u)=f %%%%%%%%%% f=2*pi^2*sin(pi*x)*sin(pi*y) %%%%%%%%%% pact difference code for elliptic equations with constant coefficient %%%%%clear%clcN=20。y=0:h:1。for i=1 A(i,i)=10/(3*h^2)。 B(i,i)=25/36。endfor i=N1 A(i,i1)=2/(3*h^2)。 B(i,i1)=5/72。endfor i=(N2)*(N1)+1 A(i,i(N1))=2/(3*h^2)。 B(i,i(N1))=5/72。endfor i=(N1)^2 A(i,i(N1)1)=1/(6*h^2)。 B(i,i(N1)1)=1/144。
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1