微分方程數(shù)值解法實(shí)驗(yàn)報(bào)告(已改無錯(cuò)字)
2022-08-18 17:34:13 本頁面
【正文】 for mm=1:N1 for nn=1:N1 uu(nn+1,mm+1)=u2((mm1)*(N1)+nn,1)。 endend%畫圖x=[0:1/N:1]。y=[0:1/N:1]。figure(1)mesh(x,y,uu)。title(39。單參數(shù)PR迭代差分解圖像（c1=c2=*(1/sin(pi/16))39。) %畫圖數(shù)值解共軛梯度法function[u_5,m_5]=gongetidufa(A,b,n,err)m_5=0。u_5_0=zeros(size(b,1),1)。 p0=bA*u_5_0。r0=p0。flag=1。while flag a0=((r0)39。*r0)/((p0)39。*(A*p0))。 u_5=u_5_0+a0*p0。 if norm(u_5u_5_0,inf)err flag=0。 end u_5_0=u_5。 r1=r0a0*A*p0。 b0=((r1)39。*r1)/((r0)39。*r0)。 p0=r1+b0*p0。 r0=r1。 m_5=m_5+1。enduu=zeros(n+1)。for mm=1:n1 for nn=1:n1 uu(nn+1,mm+1)=u_5((mm1)*(n1)+nn,1)。 endend%共軛梯度法迭代差分解圖像x=[0:1/n:1]。y=[0:1/n:1]。figure(5)mesh(x,y,uu)。title(39。共軛梯度法迭代差分解圖像39。) 預(yù)處理共軛梯度法function[u_6,m_6]=yuchuligongetidufa(A,b,n,err)m_6=0。u_6_0=zeros(size(b,1),1)。 % D=diag(A)。% D=diag(D)。% L=tril(A,1)。% L=L+D\D。% B=L*D*(L)39。B=*eye((n1)^2)p0=bA*u_6_0。r0=p0。s0=B\r0。flag=1。while flag a0=((r0)39。*s0)/((p0)39。*(A*p0))。 u_6=u_6_0+a0*p0。 if norm(u_6u_6_0,inf)err flag=0。 end u_6_0=u_6。 r1=r0a0*A*p0。 s1=B\r1。 b0=((r1)39。*s1)/((r0)39。*s0)。 p0=s1+b0*p0。 r0=r1。 s0=s1。 m_6=m_6+1。enduu=zeros(n+1)。for mm=1:n1 for nn=1:n1 uu(nn+1,mm+1)=u_6((mm1)*(n1)+nn,1)。 endend%預(yù)處理共軛梯度法迭代差分解圖像x=[0:1/n:1]。y=[0:1/n
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1