二維差分緊致差分(matlab求解泊松方程)(編輯修改稿)
2025-07-20 17:16 本頁面
【文章內(nèi)容簡介】 h=1/N。S=h^2。x=0:h:1。y=0:h:1。%%% Stiff matrix A=zeros((N1)^2,(N1)^2)。B=zeros((N1)^2,(N1)^2)。F=zeros((N1)^2,1)。for i=1 A(i,i)=10/(3*h^2)。 A(i,i+1)=2/(3*h^2)。 A(i,i+(N1))=2/(3*h^2)。 A(i,i+(N1)+1)=1/(6*h^2)。 B(i,i)=25/36。 B(i,i+1)=5/72。 B(i,i+(N1))=5/72。 B(i,i+(N1)+1)=1/144。endfor i=N1 A(i,i1)=2/(3*h^2)。 A(i,i)=10/(3*h^2)。 A(i,i+(N1)1)=1/(6*h^2)。 A(i,i+(N1))=2/(3*h^2)。 B(i,i1)=5/72。 B(i,i)=25/36。 B(i,i+(N1)1)=1/144。 B(i,i+(N1))=5/72。endfor i=(N2)*(N1)+1 A(i,i(N1))=2/(3*h^2)。 A(i,i(N1)+1)=1/(6*h^2)。 A(i,i)=10/(3*h^2)。 A(i,i+1)=2/(3*h^2)。 B(i,i(N1))=5/72。 B(i,i(N1)+1)=1/144。 B(i,i)=25/36。 B(i,i+1)=5/72。endfor i=(N1)^2 A(i,i(N1)1)=1/(6*h^2)。 A(i,i(N1))=2/(3*h^2)。 A(i,i1)=2/(3*h^2)。 A(i,i)=10/(3*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。endfor n=2:N2 i=(N2)*(N1)+n。 A(i,i(N1)1)=1/(6*h^2)。 A(i,i(N1))=2/(3*h^2)。
法律信息相關(guān)推薦
鄂ICP備17016276號-1