三线性方程组直接解法.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 zhica 0 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 zhica 0 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阶行列式为0，不能分解，A的秩如下： RA 2 h1 0 扩展：lu分解 function [l,u] lvfenjie a %对非奇异矩阵a作lu分解 [n,n] size a ;%判断被分解矩阵的维数 for j 1:n %给出U的第一行L的第一列 u 1,j a 1,j ; l j,1 a j,1 /u 1,1 ; end i 1; for k 2:n for j k:n for m 1: k-1 uu m l k,m *u m,j ;%计算公式中对应的l*u的和 end u k,j a k,j -sum uu ;%u的分解结果 for m 1: k-1 uu m l j,m *u m,k ; end l j,k a j,k -sum uu /u k,k ; end end clear a [1 2 3;1 12 7;4 5 6]; [l,u] lvfenjie a l 1.0000 0 0 1.0000 1.0000 0 4.0000 -0.3000 1.0000 u 1.0