微分方程數(shù)值解法實驗報告-在線瀏覽
2024-08-31 17:34本頁面
【正文】 diag(A)。L=tril(A,1)。B1=(Dw*L)\((1w)*D+w*R2)。m_3=0。%初始迭代值flag=1。 if norm(u_3u_3_0,inf)err flag=0。 m_3=m_3+1。for mm=1:n1 for nn=1:n1 uu(nn+1,mm+1)=u_3((mm1)*(n1)+nn,1)。y=[0:1/n:1]。 title(39。) PR方法clc。pi=。for i=1:(N1)*(N1)1 L_1(i,i)=2。L_1(i+1,i)=1。%L2矩陣的定義L_2=zeros((N1)*(N1))。endL_2((N1)^2,(N1)^2)=2。L_2(i(N1),i)=1。for j=1:N1 for i=1:N1 b((j1)*(N1)+i)=(2*pi^2*exp(pi*(i/N+j/N))*sin(pi*(i/N+j/N)))/(N^2)。%PR格式的建立u=zeros((N1)^2,1)。endu1=inv(I+tao_opt*L_1)*((Itao_opt*L_2)*u+tao_opt*b)。count=1。 u=u2。 u2=inv(I+tao_opt*L_2)*((Itao_opt*L_1)*u1+tao_opt*b)。for j=1:N1 for i=1:N1 z((j1)*(N1)+i)=exp(pi*(i/N+j/N))*(sin(pi*i/N)*sin(pi*j/N))。for mm=1:N1 for nn=1:N1 uu(nn+1,mm+1)=u2((mm1)*(N1)+nn,1)。y=[0:1/N:1]。title(39。) %畫圖數(shù)值解共軛梯度法function[u_5,m_5]=gongetidufa(A,b,n,err)m_5=0。 p0=bA*u_5_0。flag=1。*r0)/((p0)39。 u_5=u_5_0+a0*p0。 end u_5_0=u_5。 b0=((r1)39。*r0)。 r0=r1。enduu=zeros(n+1)。 endend%共軛梯度法迭代差分解圖像x=[0:1/n:1]。figure(5)mesh(x,y,uu)。共軛梯度法迭代差分解圖像39。u_6_0=zeros(size(b,1),1)。% D=diag(D)。% L=L+D\D。B=*eye((n1)^2)p0=bA*u_6_0。s0=B\r0。while flag a0=((r0)39。*(A*p0))。 if norm(u_6u_6_0,inf)err flag=0。 r1=r0a0*A*p0。 b0=((r1)39。*s0)。 r0=r1。 m_6=m_6+1。for mm
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1