matlab上机操作.doc

Matlab上级操作及答案 解线性方程组 A=[5 1 –1 0;1 0 3 –1;-1 –1 0 5;0 0 2 4];b=[1;2;3;-1]; x=A\b ans= x1=1.4, x2= -5.9, x3=0.1, x4= -0.3. 设 f(x,y) = 4 sin (x 3y)，求 。 syms x y; f=diff(4*sin(x^3*y),x); f=diff(f,y); f=subs(f,x,2); f=subs(f,y,3) ans= 1063.6 求解下列微分方程 (只要求写出x=1.5时y的值) ex6fun=@(x,y)x+y; [x,y]=ode45(ex6fun,[1 1.5],2); [x,y] ans= x=1.5000时， y= 4.0949 4、已知：，分别计算a的数组平方和矩阵平方，并观察其结果。 a=[1 2 3;4 5 6;7 8 9]; a.^2 ans = 1 4 9 16 25 36 49 64 81 a^2 ans = 30 36 42 66 81 96 102 126 150 5、角度，求x的正弦、余弦、正切和余切。 x=[30 45 60]; x1=x/180*pi; sin(x1) ans = 0.5000 0.7071 0.8660 cos(x1) ans = 0.8660 0.7071 0.5000 tan(x1) ans = 0.5774 1.0000 1.7321 cot(x1) ans = 1.7321 1.0000 0.5774 6、矩阵，分别对a进行特征值分解、奇异值分解、LU分解、QR分解及Chollesky分解。 [v,d]=eig(a,b) v = -0.4330 -0.2543 -0.1744 -0.5657 0.9660 -0.6091 -0.7018 0.0472 0.7736 d = 13.5482 0 0 0 4.8303 0 0 0 3.6216 a=[9 1 2;5 6 3;8 2 7]; [u,s,v]=svd(a) u = -0.5601 0.5320 -0.6350 -0.4762 -0.8340 -0.2788 -0.6779 0.1462 0.7204 s = 15.5234 0 0 0 4.5648 0 0 0 3.3446 v = -0.8275 0.3917 -0.4023 -0.3075 -0.9156 -0.2592 -0.4699 -0.0907 0.8781 [l,u]=lu(a) l = 1.0000 0 0 0.5556 1.0000 0 0.8889 0.2041 1.0000 u = 9.0000 1.0000 2.0000 0 5.4444 1.8889 0 0 4.8367 [q,r]=qr(a) q = -0.6903 0.3969 -0.6050 -0.3835 -0.9097 -0.1592 -0.6136 0.1221 0.7801 r = -13.0384 -4.2183 -6.8260 0 -4.8172 -1.0807 0 0 3.7733 c=chol(a) c = 3.0000 0.3333 0.6667 0 2.4267 1.1447 0 0 2.2903 7、将矩阵、和组合成两个新矩阵： （1）组合成一个4(3的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素，即