Matlab第5章线性代数相关运算程序.doc

第五章 线性代数相关运算 习题1、 A=[1 0 0;0 2 -1;0 -1 3] B=[3 5 7;0 1 0] C=[1 1;1 1;1 1] D=[0 0 0;0 0 0] 运行结果： A = 1 0 0 0 2 -1 0 -1 3 B = 3 5 7 0 1 0 C = 1 1 1 1 1 1 D = 0 0 0 0 0 0 习题2、 （1）、q=rand(5,1) （2）、m=0+(100-0)*rand(3,4) 运行结果： q =0.3922 0.6555 0.1712 0.7060 0.0318 m =27.6923 82.3458 95.0222 38.1558 4.6171 69.4829 3.4446 76.5517 9.7132 31.7099 43.8744 79.5200 q =0.1869 0.4898 0.4456 0.6463 0.7094 m =75.4687 65.5098 49.8364 58.5268 27.6025 16.2612 95.9744 22.3812 67.9703 11.8998 34.0386 75.1267 习题3、 a=magic(5) 运行结果： a =17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 习题4、 a1=[-2 -2 -2 -2 -2],x=diag(a1,0) %%%这一行也可改为：x=-2*eye(5) a2=[1 3 5 7],y=diag(a2,1) a3=[2 4 6 8],z=diag(a3,-1) A=x+y+z 运行结果： A =-2 1 0 0 0 2 -2 3 0 0 0 4 -2 5 0 0 0 6 -2 7 0 0 0 8 -2 习题5、 q=[1 2 3 4 5 6] w=[2 4 6 8 10 12] t=[1 1 1 1 1 1] e=-1*q,r=-1*w,y=-1*t A=[q;w;e;r;t;y] 将程序保存为xiti5，在命令窗口输入xiti5，会直接出现运行结果 运行结果： A =1 2 3 4 5 6 2 4 6 8 10 12 -1 -2 -3 -4 -5 -6 -2 -4 -6 -8 -10 -12 1 1 1 1 1 1 -1 -1 -1 -1 -1 -1 书130页——7 x=0:0.1:2*pi; y=cos(x); p2=polyfit(x,y,2); p3=polyfit(x,y,3); p4=polyfit(x,y,4); p6=polyfit(x,y,6); disp(二次拟合函数),f2=poly2str(p2,x) disp(三次拟合函数),f3=poly2str(p3,x) disp(四次拟合函数),f4=poly2str(p4,x) disp(六次拟合函数),f6=poly2str(p6,x) x1=0:0.1:2*pi; y2=polyval(p2,x1); y3=polyval(p3,x1); y4=polyval(p3,x1); y6=polyval(p6,x1); plot(x,y,rp,x1,y2,*,x1,y3,k-,x1,y4,y.,x1,y6,y-); legend(拟合点,二次拟合,三次拟合,四次拟合,六次拟合) 二阶拟合函数 f2 = 0.22992 x^2 - 1.4381 x + 1.4909 三阶拟合函数 f3 = 0.0038744 x^3 + 0.19389 x^2 - 1.3494 x + 1.4469 四阶