matlab实验一 数值运算实验Word格式.docx

1、,1） 1 4 7键入t=0:10;u=0:0.1:观察矩阵变量t、u的值。t = 0 1 2 3 4 5 6 7 8 9 10u = 1 至 15 列 0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000 1.1000 1.2000 1.3000 1.4000 16 至 30 列 1.5000 1.6000 1.7000 1.8000 1.9000 2.0000 2.1000 2.2000 2.3000 2.4000 2.5000 2.6000 2.7000 2.8000 2.9000 31 至 45

2、列 3.0000 3.1000 3.2000 3.3000 3.4000 3.5000 3.6000 3.7000 3.8000 3.9000 4.0000 4.1000 4.2000 4.3000 4.4000 46 至 60 列 4.5000 4.6000 4.7000 4.8000 4.9000 5.0000 5.1000 5.2000 5.3000 5.4000 5.5000 5.6000 5.7000 5.8000 5.9000 61 至 75 列 6.0000 6.1000 6.2000 6.3000 6.4000 6.5000 6.6000 6.7000 6.8000 6.900

3、0 7.0000 7.1000 7.2000 7.3000 7.4000 76 至 90 列 7.5000 7.6000 7.7000 7.8000 7.9000 8.0000 8.1000 8.2000 8.3000 8.4000 8.5000 8.6000 8.7000 8.8000 8.9000 91 至 101 列 9.0000 9.1000 9.2000 9.3000 9.4000 9.5000 9.6000 9.7000 9.8000 9.9000 10.0000,3）=2;3;4;a观察矩阵a的变化。 0 0 2 0 0 3 0 0 4b=1 1+2i;3+4i 3;观察复数矩阵

4、。b = 1.0000 + 0.0000i 1.0000 + 2.0000i 3.0000 + 4.0000i 3.0000 + 0.0000i2）创建特殊矩阵a=ones（3,3）;b=zeros（2,2）;c=eye（4）;magic（4）;观察特殊矩阵。 1 1 1 0 0c = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 13）练习矩阵运算a=0 1 0;0 0 1;-6 -11 -6;b=1 2;3 4;5 6;c=1 1 0;0 1 1;作矩阵乘运算v1=c*av1 = 0 1 1 -6 -11 -5v2=a*bv2 = 3 4 5 6 -69 -92v3=c*a*

5、bv3 = 8 10 -64 -86v4=b*cv4 = 1 3 2 3 7 4 5 11 6v5=c*bv5 = 4 6矩阵乘方运算a2 0 0 1 -6 -11 -6 36 60 25a（1/2） 0.0000 + 0.4894i 0.0000 - 0.5588i 0.0000 - 0.0482i 0.0000 + 0.2891i 0.0000 + 1.0195i 0.0000 - 0.2696i 0.0000 + 1.6179i 0.0000 + 3.2553i 0.0000 + 2.6374i矩阵加减运算a1=a+b*ca1 = 1 4 2 3 7 5 -1 0 0a2=c*b-a（1

6、:2,1:2）a2 = 4 5a3=a（1:2,2:3）+c*ba3 = 8 11矩阵右除ar=c/aar = -0.8333 -1.0000 -0.1667 1.0000 1.0000 0矩阵左除al=abal = -5.6667 -8.6667 1.0000 2.0000 3.0000 4.00004）练习矩阵特征运算完成以下特征运算inv（a）,tril（a）,rank（a） -1.8333 -1.0000 -0.1667 1.0000 0 0 0 1.0000 0 0 0 0 3（2）MATLAB语言的点运算1）练习点乘与点除a1=1 2;3 4;a2=0.2*a1;观察a1 a2 1

7、.0000 2.0000 0.2000 0.4000 3.0000 4.0000 0.6000 0.8000a1.*a2 a1./a2 0.2000 0.8000 5.0000 5.0000 1.8000 3.2000 5.0000 5.00002）由点运算完成标量函数运算与作图正余弦函数的点运算2*pi/180:2*pi;y1=sin（t）;y2=cos（t）;y=y1.*y2;plot（t,y y1 y2）;复变函数的点运算。w=0.1:2;g1=（1+0.5*w*i）/（1-0.5*w*i）;g1g1 = 0.4719 + 0.7728ig2=（1+0.5*w*i）./（1-0.5*w*

8、i）;g2g2 = 1 至 7 列 0.9950 + 0.0998i 0.9802 + 0.1980i 0.9560 + 0.2934i 0.9231 + 0.3846i 0.8824 + 0.4706i 0.8349 + 0.5505i 0.7817 + 0.6236i 8 至 14 列 0.7241 + 0.6897i 0.6632 + 0.7484i 0.6000 + 0.8000i 0.5355 + 0.8445i 0.4706 + 0.8824i 0.4060 + 0.9139i 0.3423 + 0.9396i 15 至 20 列 0.2800 + 0.9600i 0.2195 +

9、 0.9756i 0.1611 + 0.9869i 0.1050 + 0.9945i 0.0512 + 0.9987i 0.0000 + 1.0000iplot（g2）;xlabel（real g2（w）;ylabel（inag g2（w）axis（square）;（3）多项式运算1）建立多项式向量ap=1 2 2 1;b=-1 -2 -3;bp=poly（b）bp = 1 6 11 62）练习多项式乘与求根p=conv（ap,bp）p = 1 8 25 41 40 23 6roots（p） -3.0000 + 0.0000i -2.0000 + 0.0000i -0.5000 + 0.8660i -0.5000 - 0.8660i -1.0000 + 0.0000i3）练习多项式运算a=1 2 3 4;b=1 -1;c=a+zeros（1,length（a）-length（b） b;poly2str（c,x） x3 + 2 x2 + 4 x + 3polyvalm（a,3） 58（4）代数方程组1）恰定方程组A=1 2;2 3;b=8;13;方法1：逆矩阵求解X=inv（a）*bX = 2.0000 3.0000方法2：矩阵左除求解X=ab3 实验报告要求按上述步骤完成数值运算实验，并按练习实例做好记录。