高斯若爾當(dāng)消元法步驟 用matlab程序法編出高斯消元法?
用matlab程序法編出高斯消元法?函數(shù)[x,XA]=gaussxqbyorder(a,b)%Gaussian階消去n=size(a)n=n(1)forI=1:(n-1)forj=(I 1):NIF(
用matlab程序法編出高斯消元法?
函數(shù)[x,XA]=gaussxqbyorder(a,b)%Gaussian階消去n=size(a)n=n(1)forI=1:(n-1)forj=(I 1):NIF(a(I,I)==0)disp(”對(duì)角線元素為0!“)%防止對(duì)角線元素0returnendl=a(J,I)M=a(I,I)a(J,1:n)=a(J,1:n)-L*a(I,1:n)/M%消去方程B(J)=B(J)-L*B(I)/mendendx=solveuptriangle(a,B)%。函數(shù)XA=a%消去系數(shù)矩陣函數(shù)x=solveuptriangle(a,b)n=size(a)n=n(1)對(duì)于I=n:-1:1如果(I<N)s=a(I,(I 1):n*x1):n,1)elses=0endx(I,1)=(b(I)-s)/a(I,I)結(jié)束程序如上,讓我們自己計(jì)算。這個(gè)東西自己應(yīng)該試試編程參考:MATLAB常用算法