二維差分緊致差分(matlab求解泊松方程)-資料下載頁(yè)
2025-06-23 17:16本頁(yè)面
【正文】 h^2)。 B(i,i(N1)1)=1/144。 B(i,i(N1))=5/72。 B(i,i1)=5/72。 B(i,i)=25/36。 B(i,i+(N1)1)=1/144。 B(i,i+(N1))=5/72。endfor m=1:N3 for n=2:N2i=m*(N1)+n。A(i,i(N1)1)=1/(6*h^2)。A(i,i(N1))=2/(3*h^2)。A(i,i(N1)+1)=1/(6*h^2)。A(i,i1)=2/(3*h^2)。A(i,i)=10/(3*h^2)。A(i,i+1)=2/(3*h^2)。A(i,i+(N1)1)=1/(6*h^2)。A(i,i+(N1))=2/(3*h^2)。A(i,i+(N1)+1)=1/(6*h^2)。B(i,i(N1)1)=1/144。B(i,i(N1))=5/72。B(i,i(N1)+1)=1/144。B(i,i1)=5/72。B(i,i)=25/36。B(i,i+1)=5/72。B(i,i+(N1)1)=1/144。B(i,i+(N1))=5/72。B(i,i+(N1)+1)=1/144。Endfor m=0:N2 for n=1:N1i=m*(N1)+n。F(i)=2*pi^2*sin(pi*n*h)*sin(pi*(m+1)*h)。 endendU=zeros((N1)^2,1)。U=A\(B*F)。U1=zeros((N+1)^2,1)。for m=1:N1 U1(m*(N+1)+2:m*(N+1)+N)=U((m1)*(N1)+1:m*(N1))。endplot(x,y,U1,’k’)Uexact=inline(‘sin(pi*xx(:,1)).sin(pi*xx(:,2))’,’xx’)。plot(x,y,Uexact,’r’)Err=max(abs(U1Uexact))。 %最大模L2err=norm(U1Uexact)*h。 %L2 errorErrL2err
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1