列主元高斯消去法和列主元三角分解法解線性方程-文庫吧資料
2025-06-30 22:39本頁面
【正文】 nd x(i)=(y(i)h)/A(i,i)。 for i=n1:1:1 h=0。 end y(i)=A(i,n+1)h。 for i=2:n h=0。 %根據(jù)公式求解，并把結(jié)果存在矩陣A中 end for i=r+1:nA(i,r)=(A(i,r)A(i,1:r1)*A(1:r1,r))/A(r,r)。 A(r+1,:)=p。 %按照方程求取u(r,i) if abs(u(r,r))=eps %如果u（r,r）小于e，則交換行 p=A(r,:)。 A(2:n,1)=A(2:n,1)/A(1,1)。 A(1,i)=A(1,i)。 A(1,:)=A(d,:)。 else end d=d+11。 %將A與b合并，得到增廣矩陣for r=1:n if r==1 for i=1:n [c d]=max(abs(A(:,1)))。)n=length(A)。) %顯示方程組Abdisp(39。enddisp(39。輸入方程有誤！39。)x=b。disp(39。 %回代求解for i=n1:1:1 b(i)=(b(i)sum(A(i,i+1:n)*b(i+1:n)))/A(i,i)。消元后所得到的上三角陣是39。 %第i行消元處理 b(i)=b(i)m(i,k)*b(k)。end %消元 for i=k+1:n m(i,k)=A(i,k)/A(k,k)。 b(k,:)=b(q,:)。 元所在行交換(包括b) A(q,:)=temp1。 elseif qktemp1=A(k,:)。)。 %q在A(k:n,k)中的行號轉(zhuǎn)換為在A中的行號 if abs(p)c disp(39。for k=1:n1 %找列主元 [p,q]=max(abs(A(k:n,k)))。39。原方程為AX=b：39。)
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1