《MATLAB 程序设计与应用.doc

西安建筑科技大学 陈1.(1)eye(3) ans = 1 0 0 0 1 0 0 0 1 (2)A=round(rand(5,6)*100+100) A = 142 146 130 155 138 192 186 145 105 140 101 184 149 141 169 120 142 137 182 190 165 163 175 162 146 101 198 173 179 173 randn(1,500)*0.2+1 ones(size(A) ans = 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 C=diag(30 30 30 30 30 30)C(1:5,:)+A a=diag(A)B=diag(a) a = 142 145 169 163 179B = 142 0 0 0 0 0 145 0 0 0 0 0 169 0 0 0 0 0 163 0 0 0 0 0 179 2.使用函数，实现方阵左旋90或右旋90的功能A=1 4 7 10;2 5 8 11;3 6 9 12;B=rot90(A,1) B = 10 11 12 7 8 9 4 5 6 1 2 3 C=rot90(A,3) C = 3 2 1 6 5 4 9 8 7 12 11 10 3.A=rand(3);B=inv(A) B = -0.3872 1.2413 -0.0844 -5.1444 -5.1549 9.9467 4.6429 3.2377 -6.7282 det(A) ans = 0.0782 A*B ans = 1.0000 -0.0000 0.0000 -0.0000 1.0000 0.0000 -0.0000 -0.0000 1.0000 B*A ans = 1.0000 0.0000 0.0000 0.0000 1.0000 -0.0000 -0.0000 0.0000 1.0000 4.A=4 2 -1;3 -1 2;12 3 0;b=2 10 8;x=Ab x = -6.0000 26.6667 27.3333 或者用x=inv(A)*b x = -6.0000 26.6667 27.3333 5.A=1 -1 2 3;5 1 -4 2;3 0 5 2;11 15 0 9diag(A) A = 1 -1 2 3 5 1 -4 2 3 0 5 2 11 15 0 9ans = 1 1 5 9 triu(A) ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ans = 4 norm(A) ans = 21.3005 norm(A,1) ans = 20 norm(A,inf) ans = 35 cond(A) ans = 11.1739 cond(A,1) ans = 14.4531 cond(A,inf) ans = 22.0938 trace(A) ans = 16 B=0.43 43 2;-8.9 4 21 B = 0.4300 43.0000 2.0000 -8.9000 4.0000 21.0000 diag(B) ans = 0.4300 4.0000 triu(B) ans = 0.4300 43.0000 2.0000 0 4.0000 21.0000 tril(B) ans = 0.4300 0 0 -8.9000 4.0000 0 rank(B) ans = 2 norm(B) ans = 43.4271 norm(B,1) ans = 47 norm(B,inf) ans = 45.4300 cond(B) ans = 1.9354 trace(B) ans =4.4300 6.A=1 1 0.5;1 1 0.25;0.5 0.25 2 A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.25000.5000 0.2500 2.0000 V,D=eig(A) V = 0.7212 0.444