线性代数上机作业

  • 格式:doc
  • 大小:96.50 KB
  • 文档页数:19

1.A,B列向量组的一个最大无关组.2.A*B+B*A3.Ax=b,Bx=b的解4.A,B的行列式,逆,秩5.A-B,A+B,A*B的行列式,逆,秩>> a=[14 17 10 111 2 3 412 15 22 1721 74 14 70]a =14 17 10 111 2 3 412 15 22 1721 74 14 70>> rref(a)ans =1 0 0 00 1 0 00 0 1 00 0 0 1b =11 21 14 171 2 3 412 15 22 175 7 1 6>> a*b+b*aans =1046 2252 1197 2175205 460 228 4501300 2553 1601 24961038 1847 1071 1831 c =127410>> x=a\cx =-13.113110.847210.4454-9.4794>> x=b\cx =12.5267-7.55381.8938-0.2751>> rref(b)ans =1 0 0 00 1 0 00 0 1 00 0 0 1 >> det(a)ans =-10324>> det(b)ans =1403>> inv(a)ans =0.2367 0.8346 -0.1980 -0.0368-0.1281 -0.9845 0.1702 0.0351 -0.1083 -0.4669 0.1544 0.00620.0861 0.8838 -0.1514 -0.0130>> inv(b)ans =-0.1703 -0.3264 0.1390 0.30650.1696 -0.1561 -0.0798 -0.15040.0235 -0.0763 0.0478 -0.1511-0.0599 0.4669 -0.0306 0.1119 >> rank(a)ans =4>> rank(b)ans =4>> e=a-b,f=a+b,g=a*be =3 -4 -4 -60 0 0 00 0 0 016 67 13 64f =25 38 24 282 4 6 824 30 44 3426 81 15 76g =346 555 478 54269 98 90 100496 731 714 740823 1289 894 1311>> x1=det(e),x2=det(f),x3=det(g)x1 =x2 =-29508x3 =-14484572>> y1=inv(e),y2=inv(f),y3=inv(g)Warning: Matrix is singular to working precision.(Type "warning off MATLAB:singularMatrix" to suppress this warning.) y1 =Inf Inf Inf InfInf Inf Inf InfInf Inf Inf InfInf Inf Inf Infy2 =0.3637 1.0802 -0.3083 -0.1098-0.2116 -0.8926 0.2107 0.0777-0.1560 -0.5076 0.1643 0.03740.1319 0.6819 -0.1515 -0.0394y3 =0.0128 0.3852 -0.0468 -0.00830.0559 0.1996 -0.0497 -0.0103-0.0028 -0.0611 0.0126 -0.0013 -0.0610 -0.3964 0.0696 0.0169>> z1=rank(e),z2=rank(f),z3=rank(g)z1 =2z2 =4z3 =41-1. 4. (1-4)>> a=[4 2 -1-2 1 01 0 0]a =4 2 -1-2 1 01 0 0>> b=[0 31 -12 1]b =0 31 -12 1>> x=a*bx =0 91 -70 3>> b=[4 0 4 50 1 7 62 3 0 0]b =4 0 4 50 1 7 62 3 0 0 >> x=a*bx =20 6 16 206 107 6 >> a=[4 0 00 1 0-2 0 4]a =4 0 00 1 0-2 0 4>> b=[1 3 64 5 -14 7 0]b =1 3 64 5 -14 7 0>> c=[0 1 01 0 00 0 1]0 1 01 0 00 0 1 >> x=a*b*cx =12 4 245 4 -122 14 -12>> a=[1 1 31 1 01 0 0]a =1 1 31 1 01 0 0>> b=[1 1 31 1 01 0 0]b =1 1 31 1 01 0 0>> b=b'b =1 1 11 1 03 0 0>> x=a*bx =11 2 12 2 11 1 1 1-3. 1.2>> a=[4 1 -12 0 -43 3 -1]a =4 1 -12 0 -43 3 -1>> b=[125]b =125>> x=a\bx =-0.37501.8125-0.6875>> a=[4 1 12 2 31 2 0]a =4 1 12 2 31 2 0>> b=[150]b =15>> x=a\bx =-0.21050.10531.7368>> A=[0 1 01 0 00 0 4]A =0 1 01 0 00 0 4>> B=[0 1 32 1 34 5 2]B =0 1 32 1 34 5 2>> C=[4 3 00 1 00 0 4]C =4 3 00 1 00 0 4>> X=A*B,Y=B*A,Z=C*B,M=B*C,N=A*C X =2 1 30 1 316 20 8Y =1 0 121 2 125 4 8Z =6 7 212 1 316 20 8M =0 1 128 7 1216 17 8N =0 1 04 3 00 0 16>> P=C*AP =3 4 01 0 02-1. 2.(1-3)>> A=[4 0 2 0-1 0 3 00 2 0 -10 1 0 0]A =4 0 2 0-1 0 3 00 2 0 -10 1 0 0>> X=det(A)X =-14>> B=[0 0 0 40 0 4 30 4 3 24 3 2 4]B =0 0 0 40 0 4 30 4 3 24 3 2 4 >> Y=det(B)Y =256>> c=[4 2 -12 3 13 5 1]c =4 2 -12 3 1>> z=det(c)z =-72-4.1.(1,2,4)>> a=[1 3 0 02 1 0 00 0 6 50 0 8 1]a =1 3 0 02 1 0 00 0 6 50 0 8 1>> b=[4 0 0 30 3 2 00 6 5 04 0 0 2]b =4 0 0 30 3 2 00 6 5 04 0 0 2>> c=[4 2 0 3 00 1 0 4 00 0 0 5 07 11 6 5 7-1 8 1 8 4]c =4 2 0 3 00 1 0 4 00 0 0 5 07 11 6 5 7-1 8 1 8 4 >> x=det(a),y=det(b),z=det(c)x =170y =-12z =-3403-1.2.7>> a=[4 -1 -12 -1 01 0 2]a =4 -1 -12 -1 01 0 2>> x=inv(a)x =0.4000 -0.4000 0.20000.8000 -1.8000 0.4000-0.2000 0.2000 0.4000 > b=[4 1 -12 3 -22 1 0]b =4 1 -12 3 -22 1 0>> x=inv(b)x =0.2500 -0.1250 0.1250-0.5000 0.2500 0.7500-0.5000 -0.2500 1.2500>> c=[-4 -3 0 03 5 0 00 0 -2 30 0 -1 1]c =-4 -3 0 03 5 0 00 0 -2 30 0 -1 1>> x=inv(c)x =-0.4545 -0.2727 0 00.2727 0.3636 0 00 0 1.0000 -3.00000 0 1.0000 -2.0000 >> d=[4 1 1 11 1 -1 -11 -1 1 -11 -1 -1 4]d =4 1 1 11 1 -1 -11 -1 1 -11 -1 -1 4>> x=inv(d)0.1000 0.2500 0.2500 0.10000.2500 -0.1250 -0.6250 -0.25000.2500 -0.6250 -0.1250 -0.25000.1000 -0.2500 -0.2500 0.1000 >> a=[4 1-2 1]a =4 1-2 1>> b=[1 02 4]b =1 02 4>> c=[1 -10 -4]c =1 -10 -4>> d=a\b,e=inv(b)d =-0.1667 -0.66671.66672.6667e =1.0000 0-0.5000 0.2500>> x=e*dx =-0.1667 -0.66670.5000 1.0000>> a=[4 0 00 1 03 0 4]a =4 0 00 1 03 0 4>> b=[4 0 00 0 10 1 0]b =4 0 00 0 10 1 0>> c=[4 2 34 5 67 8 2]c =4 2 34 5 67 8 2>> d=a\b,e=inv(b)d =1.0000 0 00 0 1.0000-0.7500 0.2500 00.2500 0 00 0 1.00000 1.0000 0 >> x=e*dx =0.2500 0 0-0.7500 0.2500 00 0 1.0000 >> a=[4 -1 0-1 1 -1-1 0 1]a =4 -1 0-1 1 -1-1 0 1>> e=[1 0 00 1 00 0 1]e =1 0 00 1 00 0 1>> c=[1 11 23 0]c =1 11 23 0>> b=e-ab =-3 1 01 0 11 0 0>> x=b\cx =3 010 1-2 2>> a=[4 -1 -12 -1 -33 -2 -7]a =4 -1 -12 -1 -33 -2 -7>> b=[1 11 40 7]b =1 11 40 7>> x=a\bWarning: Matrix is singular to working precision.(Type "warning off MATLAB:singularMatrix" to suppress this warning.) x =Inf InfInf InfInf Inf4-2.1.>> a=[4 2 4 11 3 1 52 0 2 21 -1 -2 1]a =4 2 4 11 3 1 52 0 2 21 -1 -2 1>> b=[4 1 2 1 22 0 1 -1 50 2 3 3 -11 1 0 -1 1]b =4 1 2 1 22 0 1 -1 50 2 3 3 -11 1 0 -1 1 >> x1=rank(a),x2=rank(b)x1 =4x2 =4。