实验二向量组的线性相关性及线性方程组分析.ppt

实验二 二、齐次线性方程组解的结构 三、非齐次线性方程组解的结构 实验三 一、特征值与特征向量 四、习题 X=A\b X = 1.0000 -2.0000 -0.0000 0.5000 解法2：用Gramer法则求解 A=[1 -1 1 –2;2 0 -1 4;3 2 1 0;-1 2 -1 2]; b=[2;4;-1;-4]; D1=[b,A(:,2:4)]; D2=[A(:,1),b,A(:,3:4)]; D3=[A(:,1:2),b,A(:,4)]; D4=[A(:,1:3),b]; D=det(A); x1=det(D1)/D; x2=det(D2)/D; x3=det(D3)/D; x4=det(D4)/D; X=[x1;x2;x3;x4] X = 1.0000 -2.0000 0 0.5000 matlab中提供了一个求解线性代数方程组的函数linsolve(A,b)，其调用格式是： 解法3： linsolve(A,b) 该函数返回线性代数方程组 AX=b 的解。这里 A,b 允许是符号矩阵， A 必须行满秩（即rank(A)=m，ｍ是A的行数）。 A=[1 -1 1 -2;2 0 -1 4;3 2 1 0;-1 2 -1 2]; b=[2 4 -1 -4]; X=linsolve(A,b) X = [ 1] [ -2] [ 0] [ 1/2] 解1： A=[1 –1 1 -2;2 0 -1 4;3 2 1 0;-1 2 –1 2 ] ; b=[2 4 –1 -4]； det(A) ans = -2 x=inv（A）*b 或 x=A \ b x= 1 -2 0 1/2 x= 1 -2 0 1/2 A=[1 -1 1 -2;2 0 -1 4;3 2 1 0;-1 2 -1 2 ] ; b=[2 4 -1 -4]; D=det(A); n=size(A) D1=ones (n(1),1); for i=1:n(1) A1=A; A1(:,i)=b; D1(i)=det(A1); end x=D1/D 解2： x = 1.0000 -2.0000 0 0.5000 解： A=[1 2 3 4;2 2 1 1;2 4 6 8;4 4 2 2]; b=[1;3;2;6]; C=[A b]; [rank(A),rank(C)] ans = 2 2 B=null(A,r) x0=pinv(A)*b 或 x0= A\b x0 = 0.9542 0.7328 -0.0763 -0.2977 B = 2.0000 3.0000 -2.5000 -3.5000 1.0000 0 0 1.0000 答：非齐次线性方程组的通解为： 解：建立一个M-文件：line1 如下： A=[1 5 -1 -1;1 -2 1 3;3 8 -1 1;1 -9 3 7 ]; b=[-1;3;1;7]; B=[A b]; n=4; r1=rank(A); r2=rank(B); if (r1==r2r1==n) X=A\b %唯一解 else if (r1==r2r1n) C=null(A,r) %基础解系 x0=pinv(A)*b %特解 else fprintf(方程组无解) end end C = -0.4286 -1.8571 0.2857 0.5714 1.0000 0 0 1.0000 x0 = 0.3785 -0.0876 0.1873 0.7530 输入命令：line1 解：建立一个M-文件：line2 如下： A=[4 2 -1;3 -1 2;11 3 0 ]; b=[2;10;8