基于matlab平臺的三種迭代法求解矩陣方程-資料下載頁
2025-03-30 04:01本頁面
【正文】 +1。 xx(:,k+1)=x。 if k=M disp(＇迭代次數太多可能不收斂?。?)。 break。 end n=k。 jingdu(k)=norm(r)。 end end 子函 31(共軛梯度法 ): function [x,xx,n,jingdu]=con_grad(A,b,eps,M) % 利用迭代方法求解矩陣方程 這里是共軛梯度迭代方法 % A 為系數矩陣 b 為右端向量 err 為精度大小 返回求解所得向量 x及迭代次數 % M 為最大迭代次數 jingdu 求解過程的精度 n 所需迭代次數 xx 存儲求解過程中每次迭代產生的解 x0=zeros(length(b),1)。 r0=bA*x0。 p0=r0。 % t0=r0＇ *r0/(r0＇ *A*r0)。 % x=x0+t0*r0。 % r=bA*x。 % xx(:,1)=x。 k=0。 x=x0。 r=r0。 p=p0。 while norm(r)eps x=x。 r=r。 p=p。 afa=r＇ *r/(p＇ *A*p)。 x1=x+afa*p。 r1=rafa*A*p。 beta=r1＇ *r1/(r＇ *r)。 p1=r1+beta*p。 x=x1。 r=r1。 p=p1。 k=k+1。 xx(:,k)=x。 if k=M disp(＇迭代次數太多可能不收斂?。?)。 break。 end n=k。 jingdu(k)=norm(r)。 end end
環(huán)評公示相關推薦
鄂ICP備17016276號-1