三、线性方组直接解法.docVIP

1、编写判断一个方程组是否有解的函数，并调用之。 function[RA,RB,n]=jiepb(A,b) B=[A b];n=length(b);RA=rank(A); RB=rank(B);zhica=RB-RA; if zhica0 disp(无解) return end if RA==RB if RA==n disp(唯一解) else disp(无穷解) end end A=[1 2 3 4;2 3 4 5;3 4 5 6;4 5 6 7]; b=[6 7 8 9]; [RA,RB,n]=jiepb(A,b) 无穷解 RA = 2 RB = 2 n = 4 A=[1 -1 1 -3;0 -1 -1 1;2 -2 -4 6;1 -2 -4 1]; b=[1 0 -1 -1]; [RA,RB,n]=jiepb(A,b) 唯一解 RA = 4 RB = 4 n = 4 2、编写用列主元消去法解线性方程组的函数，并调用之。 function[RA,RB,n,X]=gaus(A,b) B=[A b];n=length(b);RA=rank(A); RB=rank(B);zhica=RB-RA;X=zeros(n,1); if zhica0 disp(无解) return end if RA==RB if RA==n disp(唯一解) C=zeros(1,n+1); for p=1:n-1 for k=p+1:n l=B(k,p)/B(p,p); B(k,p:n+1)=B(k,p:n+1)-l*B(p,p:n+1); end end b=B(1:n,n+1);A=B(1:n,1:n);X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); end else disp(无穷多解) end end A=[1 -1 1 -3;0 -1 -1 1;2 -2 -4 6;1 -2 -4 1]; b=[1 0 -1 -1]; [RA,RB,n,X]=gaus(A,b) 唯一解 RA = 4 RB = 4 n = 4 X = 0 -0.5000 0.5000 0 3、编写判断矩阵A能否进行LU分解的函数，并调用之。 function h1=pdlufj(A) [n n]=size(A);RA=rank(A); if RA~=n disp(n阶行列式为0，不能分解，A的秩如下：),RA,h1=det(A); return end if RA==n for p=1:n h(p)=det(A(1:p,1:p)) end h1=h(1:n); for i=1:n if h(1,i)==0 disp(r阶主子式为0，不能分解，A的秩和各阶顺序主子式值如下：),RA,h1; return end end if h(1,i)~=0 disp(各阶主子式都不为0，能进行分解，A的秩和各阶顺序主子式值如下：),RA,h1; end end clear A=[1 2 3;1 12 7;4 5 6]; h1=pdlufj(A) h = 1 h = 1 10 h = 1 10 -48 各阶主子式都不为0，能进行分解，A的秩和各阶顺序主子式值如下： RA = 3 h1 = 1 10 -48 clear A=[1 2 3;1 2 7;4 5 6]; h1=pdlufj(A) h = 1 h = 1 0 h = 1 0 12 r阶主子式为0，不能分解，A的秩和各阶顺序主子式值如下： RA = 3 h1 = 1 0 12 clear A=[1 2 3;1 2 3;4 5 6]; h1=pdlufj(A) n阶行