Matlab前三次作业答案.pptVIP

作业答案 杨颖 中国农业大学 信息与电气工程学院 2009-10-31 题目 1、用函数ones和diag分别编写下列矩阵。 解法一 1、第一个矩阵的生成： a=ones(7) b=5*ones(5) a(2:6,2:6)=b(1:5,1:5) c=3*ones(3) a(3:5,3:5)=c a(4,4)=7 2、第二个矩阵的生成： a=[2 2 2], b=[3 3 3 3]; c=[4 4 4 4 4 ]; d=[5 5 5 5 ]; e=[6 6 6 ]; f=diag(a,-2)+diag(b,-1)+diag(c,0)+diag(d,1)+diag(e,2) 解法二 % 1、用函数ones和diag分别编写下列矩阵。 clear all;close all;clc; % 1.1---------------------------------------------------------------------- Result1 = ones(7); Result1(2:6,2:6) = Result1(2:6,2:6) + 2*ones(5); Result1(3:5,3:5) = Result1(3:5,3:5) + 2*ones(3); Result1(4,4) = Result1(4,4) + 2; % 1.2---------------------------------------------------------------------- Result2 = diag(ones(1,3)*2,-2) + diag(ones(1,4)*3,-1) + ... diag(ones(1,5)*4,0) + diag(ones(1,4)*5,1) + diag(ones(1,3)*6,2); % 输出结果------------------------------------------------------------------ Result1 Result2 题目 2、X为4阶随机矩阵，分别对其进行如下操作： （1）lu分解 （2）正交分解 （3）cholesky分解 （4）奇异值分解 解法 %a为4阶随机矩阵 a=rand(4) %（1）lu分解 [l,u,p]= lu(a) %（2）正交分解 [q,r]= qr(a) %（4）奇异值分解 s = svd(a) [u,s,v]=svd(a) %（3）cholesky分解，由于a不是正定矩阵，因此用pascal生成b正定矩阵，就可求得 b =pascal(4) T =chol(b) 题目 3、用两种方法求解下列矩阵的逆： 解法一 方法一： A=[1 1 1 1;1 2 3 4;1 3 6 10;1 4 10 20] inv(A) ans = 4.0000 -6.0000 4.0000 -1.0000 -6.0000 14.0000 -11.0000 3.0000 4.0000 -11.0000 10.0000 -3.0000 -1.0000 3.0000 -3.0000 1.0000 方法三： eye(4)/A ans = 4 -6 4 -1 -6 14 -11 3 4 -11 10 -3 -1 3 -3 1 解法二 B=[A,eye(4)] B = 1 1 1 1 1 0 0 0 1 2 3 4 0 1 0 0 1 3 6 10 0 0 1 0 1 4 10 20 0 0 0 1 rref(B) ans = 1 0 0 0 4 -6 4 -1 0 1 0 0 -6 14 -11 3 0 0 1 0 4 -11 10 -3 0 0 0 1 -1 3 -3 1 ans(:,5:8) ans = 0 4 -6 4 -1 0 -6 14 -