线代上机作业VIP专享VIP免费

上机作业（一） 随机生成5 阶方阵A,B,C 及5 维列向量b，求： ①A+B,A-B. 解：在 Matlab 中输入： A=rand(5,5); B=rand(5,5); C=rand(5,5); b=rand(5,1); 得： 输入 A+B；A-B. 得： ②A*B+B*A. 解：输入 A*B+B*A. 得： ③Ax=b 的解,并验证克莱姆法则. 解：输入：x=A\b. 得： 输入：D=A;D1=A;D2=A;D3=A;D4=A;D5=A;D1(:,1)=b;D2(:,2)=b;D3(:,3)=b;D4(:,4)=b;D5(:,5)=b; Y=[det(D1)/det(D);det(D2)/det(D);det(D3)/det(D);det(D4)/det(D);det(D5)/det(D)]. 得： 发现：x=y,故克莱姆法则成立. ④A,B 的行列式，逆，秩. 解：输入：det(A);det(B);inv(A);inv(B);rank(A);rank(B). 得： ⑤A*B 的行列式，逆，秩，并验证 det(A*B)=det(A)*det(B). 解：输入 det(A*B);det(A)*det(B). 得： 可见 det(A*B)=det(A)*det(B). ⑥验证(AB)ᵄ = ᵃᵄᵃᵄ, (AB)−1= ᵃ−1ᵃ−1, ᵃᵃ ≠ ᵃᵃ. 解：输入(A*B)’;B’*A’;inv(A*B);inv(B)*inv(A);A*B;B*A. 由此可见(AB)ᵄ = ᵃᵄᵃᵄ,(AB)−1= ᵃ−1ᵃ−1, ᵃᵃ ≠ ᵃᵃ. ⑦求矩阵 X 使得 AXB=C. 解：输入 X=(A\C)/B. 得： 上机作业（二） 验证：对于一般的方阵 A,B,C,D， 若 A,C 均为对角矩阵，且 A 可逆， 则： 解：输入：A=rand(3,3);B=rand(3,3);C=rand(3,3);D=rand(3,3);E(1:3,1:3)=A;E(1:3,4:6)=B;E(4:6,1:3)=C;E(4:6,4:6)=D;det(E);det(A)*det(D)-det(B)*det(C). ABA DB CCD ABADCBCD  由此可见 在Matlab 中输入：A=diag(diag(rand(3,3)));B=rand(3,3);C=diag(diag(rand(3,3)));D=rand(3,3);E(1:3,1:3)=A;E(1:3,4:6)=B;E(4:6,1:3)=C;E(4:6,4:6)=D;det(E);det(A*D-C*B). 由此可见 上机作业（三） N= 201465004 共9 位 a=最后两位 04. b=第4-5 位 46. c=第6-7 位 50. d=第4,8 位 40. e=第1,8 位 20. f=第5,9 位 64. g=第4,9 位 44. h=第5,7 位 60. 求A 列向量组的一个最大无关组，并把不属于 极大无关组的向量利用极大无关组表示. 解：由题意可得该矩阵为：A=[44650123403444312152220644417576080] . 在 MAtlab 中输入 A=[4,46,50,40,3,4;1,2,3,4,4,3;12,15,22,17,5,7;20,64,44,60,8,0];b=rref(A). 所以ᵯ1, ᵯ2, ᵯ3,ᵯ4是一个极大无关组，且有： ᵯ5 = −0.2945ᵯ1 − 1.4863ᵯ2 − 0.0062ᵯ3 + 1.8214ᵯ4 ᵯ6 = −0.0863ᵯ1 − 1.3001ᵯ2 + 0.3643ᵯ3 +...