二維差分緊致差分(matlab求解泊松方程)(參考版)
2025-06-26 17:16本頁面
【正文】 %L2 errorErrL2err。plot(x,y,Uexact,’r’)Err=max(abs(U1Uexact))。for m=1:N1 U1(m*(N+1)+2:m*(N+1)+N)=U((m1)*(N1)+1:m*(N1))。U=A\(B*F)。F(i)=2*pi^2*sin(pi*n*h)*sin(pi*(m+1)*h)。B(i,i+(N1)+1)=1/144。B(i,i+(N1)1)=1/144。B(i,i)=25/36。B(i,i(N1)+1)=1/144。B(i,i(N1)1)=1/144。A(i,i+(N1))=2/(3*h^2)。A(i,i+1)=2/(3*h^2)。A(i,i1)=2/(3*h^2)。A(i,i(N1))=2/(3*h^2)。endfor m=1:N3 for n=2:N2i=m*(N1)+n。 B(i,i+(N1)1)=1/144。 B(i,i1)=5/72。 B(i,i(N1)1)=1/144。 A(i,i+(N1)1)=1/(6*h^2)。 A(i,i1)=2/(3*h^2)。 A(i,i(N1)1)=1/(6*h^2)。 B(i,i+(N1)+1)=1/144。 B(i,i+1)=5/72。 B(i,i(N1)+1)=1/144。 A(i,i+(N1)+1)=1/(6*h^2)。 A(i,i+1)=2/(3*h^2)。 A(i,i(N1)+1)=1/(6*h^2)。endfor m=1:N3 i=m*(N1)+1。 B(i,i+(N1))=5/72。 B(i,i+1)=5/72。 B(i,i1)=5/72。 A(i,i+(N1))=2/(3*h^2)。 A(i,i+1)=2/(3*h^2)。 A(i,i1)=2/(3*h^2)。 B(i,i+1)=5/72。 B(i,i1)=5/72。 B(i,i(N1))=5/72。 A(i,i+1)=2/(3*h^2)
法律信息相關(guān)推薦
鄂ICP備17016276號-1