Matlab实验报告终极版Word文档格式.docx

1、X = 0.4000 -1.4000 3.6000 0 1.5000 -2.5000 Y = BAY = -0.9500 1.3500 0.5500 -0.1500 0.9000 -0.70003、输入A=7 1 5；2 5 6；3 1 5，B=1 1 1; 2 2 2; 3 3 3，在命令窗口中执行下列表达式，掌握其含义：答：A（2, 3） %矩阵A第二行第三列的元素 A（:,2） %矩阵A第二列的数据A（3,:） %矩阵A第二行的数据,1:2:3） %矩阵A第一列和第三列的数据 ,3）.*B（:,2） %矩阵A的第三列和矩阵B的第二列点乘,3）*B（2,:） %矩阵A 的第三列和矩阵B的第

2、二行相乘 A*B %矩阵A和矩阵B相乘A.*B %矩阵A和矩阵B点乘A2 %矩阵A的乘方A.2 %矩阵A的点乘方B/A %矩阵B右除矩阵AB./A %矩阵B点除矩阵A A=7 1 5;2 5 6;3 1 5;B=1 1 1; 3 3 3; A（2, 3） ans = 6 A（:,2） ans = 1 5 1 A（3,:） ans = 3 1 53）ans = 7 5 2 6 3 5,2）ans = 5 12 15ans = 10 10 10 12 12 12 10 10 10 A*B ans = 24 24 24 30 30 30 20 20 20 A.*B ans = 7 1 5 4 10

3、12 9 3 15 A2 ans = 66 17 66 42 33 70 38 13 46 A.2 ans = 49 1 25 4 25 36 9 1 25 B/Aans = 0.1842 0.2105 -0.2368 0.3684 0.4211 -0.4737 0.5526 0.6316 -0.7105 B./Aans = 0.1429 1.0000 0.2000 1.0000 0.4000 0.3333 1.0000 3.0000 0.60004、输入C=1:20，则C（i）表示什么？其中i=1,2,3,10； C=1:20C = 1 3 5 7 9 11 13 15 17 19 C（i）

4、下标索引必须为正整数类型或逻辑类型。 C（1） C（2）ans = 3 C（3） C（4）ans = 7 C（5）ans = 9 C（6）ans = 11 C（7）ans = 13 C（8）ans = 15 C（9）ans = 17 C（10）ans = 195、先求下列表达式的值，然后显示Matlab工作空间的使用情况并保存全部变量。（1）（2），其中（3）提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。（1） z1=2*sin（85/360）/（1+（exp（1）2）z1 =0.0558 x=2 1+2*i;-0.45 5x = 2.0000 + 0.0000i 1.0000

5、+ 2.0000i -0.4500 + 0.0000i 5.0000 + 0.0000i（2） z2=（log（x+sqrt（1+x2）/2z2 = 0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i（3） a=-3:0.1:3a =Columns 1 through 9-3.0000 -2.9000 -2.8000 -2.7000 -2.6000 -2.5000 -2.4000 -2.3000 -2.2000Columns 10 through 18-2.1000 -2.0000 -1.9000 -1.800

6、0 -1.7000 -1.6000 -1.5000 -1.4000 -1.3000Columns 19 through 27 -1.2000 -1.1000 -1.0000 -0.9000 -0.8000 -0.7000 -0.6000 -0.5000 -0.4000 Columns 28 through 36-0.3000 -0.2000 -0.1000 0 0.1000 0.2000 0.3000 0.4000 0.5000Columns 37 through 450.6000 0.7000 0.8000 0.9000 1.0000 1.1000 1.2000 1.3000 1.4000C

7、olumns 46 through 541.5000 1.6000 1.7000 1.8000 1.9000 2.0000 2.1000 2.2000 2.3000Columns 55 through 612.4000 2.5000 2.6000 2.7000 2.8000 2.9000 3.0000 m1=linspace（0.3,0.3,61）m1 =Columns 1 through 90.3000 0.3000 0.3000 0.3000 0.3000 0.3000 0.3000 0.3000 0.3000Columns 28 through 360.3000 0.3000 0.300

8、0 0.3000 0.3000 0.3000 0.3000 m2=linspace（-0.3,-0.3,61）m2 =Columns 1 through 9-0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000-0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 -0.3000 z2=（power（exp（1）,m1.*a）

9、-power（exp（1）,m2.*a）/2.*sin（a+m1）+log（a+m1）/2）z2 =Columns 1 through 40.7388 + 3.1416i 0.7696 + 3.1416i 0.7871 + 3.1416i 0.7913 + 3.1416iColumns 5 through 80.7822 + 3.1416i 0.7602 + 3.1416i 0.7254 + 3.1416i 0.6784 + 3.1416iColumns 9 through 120.6196 + 3.1416i 0.5496 + 3.1416i 0.4688 + 3.1416i 0.3780

10、+ 3.1416iColumns 13 through 160.2775 + 3.1416i 0.1680 + 3.1416i 0.0497 + 3.1416i -0.0771 + 3.1416iColumns 17 through 20-0.2124 + 3.1416i -0.3566 + 3.1416i -0.5104 + 3.1416i -0.6752 + 3.1416iColumns 21 through 24-0.8536 + 3.1416i -1.0497 + 3.1416i -1.2701 + 3.1416i -1.5271 + 3.1416iColumns 25 through

11、 28-1.8436 + 3.1416i -2.2727 + 3.1416i -2.9837 + 3.1416i -37.0245 + 0.0000iColumns 29 through 32-3.0017 + 0.0000i -2.3085 + 0.0000i -1.8971 + 0.0000i -1.5978 + 0.0000iColumns 33 through 36-1.3575 + 0.0000i -1.1531 + 0.0000i -0.9723 + 0.0000i -0.8083 + 0.0000iColumns 37 through 40-0.6567 + 0.0000i -0

12、.5151 + 0.0000i -0.3819 + 0.0000i -0.2561 + 0.0000iColumns 41 through 44-0.1374 + 0.0000i -0.0255 + 0.0000i 0.0792 + 0.0000i 0.1766 + 0.0000iColumns 45 through 480.2663 + 0.0000i 0.3478 + 0.0000i 0.4206 + 0.0000i 0.4841 + 0.0000iColumns 49 through 520.5379 + 0.0000i 0.5815 + 0.0000i 0.6145 + 0.0000i

13、 0.6366 + 0.0000iColumns 53 through 560.6474 + 0.0000i 0.6470 + 0.0000i 0.6351 + 0.0000i 0.6119 + 0.0000iColumns 57 through 600.5777 + 0.0000i 0.5327 + 0.0000i 0.4774 + 0.0000i 0.4126 + 0.0000iColumn 610.3388 + 0.0000i6、已知：，求下列表达式的值。，其中I为单位矩阵；和（4）（5） A=12 34 -4;34 7 87;3 65 7;B=1 3 -1;2 0 3;3 -2 7;

14、A+6*Bans = 18 52 -10 46 7 105 21 53 49 I=eye（3）; A-B+Ians =12 31 -3 32 8 84 0 67 1 A*Bans = 68 44 62 309 -72 596 154 -5 241 A.*Bans = 12 102 4 68 0 261 9 -130 49 A3ans = 37226 233824 48604 247370 149188 600766 78688 454142 118820 A.3ans = 1728 39304 -64 39304 343 658503 27 274625 343 A/Bans = 16.400

15、0 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000 ABans = -0.0313 0.3029 -0.3324 0.0442 -0.0323 0.1063 0.0317 -0.1158 0.1558 A,Bans = 12 34 -4 1 3 -134 7 87 2 0 3 3 65 7 3 -2 7 A（1,3,:）;B2ans = 12 34 -4 3 65 7 4 5 1 11 0 19 20 -5 40实验二 Matlab的数值计算和符号运算一、上机验证实验A、开启PC机，进入MATLAB语言B、帮

16、助命令的使用，查找 sqrt函数的使用方法 help sqrtsqrt - Square rootThis MATLAB function returns the square root of each element of the array X.B = sqrt（X） C、矩阵运算（1）已知A=1 2;3 4;B=5 5;7 8;求A2*Bans =105 115229251（2）矩阵除法。已知A=1 2 3;4 5 6;7 8 9;B=1 0 0;0 2 0;0 0 3; 求A/B，ABA=1 2 3;A/Bans = 1.0000 1.0000 1.0000 4.0000 2.5000

17、 2.0000 7.0000 4.0000 3.0000AB警告: 矩阵接近奇异值，或者缩放错误。结果可能不准确。RCOND = 1.541976e-18。 ans =1.0e+16 * -0.4504 1.8014 -1.3511 0.9007 -3.6029 2.7022（3）矩阵的转置及共轭转置。已知A=15+i, 2-i ,1;6*i,4,9-i;求A.，A A=15+i, 2-i ,1; A.ans =15.0000 + 1.0000i 0.0000 + 6.0000i 2.0000 - 1.0000i 4.0000 + 0.0000i 1.0000 + 0.0000i 9.0000

18、 - 1.0000i Aans =15.0000 - 1.0000i 0.0000 - 6.0000i 2.0000 + 1.0000i 4.0000 + 0.0000i 1.0000 + 0.0000i 9.0000 + 1.0000i（4）使用冒号选出指定元素。 求A中第3列前2个元素；A中所有第2行的元素； A=1 2 3; A（1:2,3） 6 A（2,1:ans = 4 5 6（5）方括号。用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列。 A=magic（4）A = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 A（ : ,4）=A =16

19、2 3 5 11 10 9 7 6 4 14 15D、多项式（1）求多项式的根 P=1,0,-2,-4; roots（P）ans = 2.0000 + 0.0000i -1.0000 + 1.0000i -1.0000 - 1.0000i（2）已知A1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4, 求矩阵A的特征多项式； A=1.2 3 5 0.9;1 2 3 4; B=poly（A）B = 1.0000 -6.9000 -77.2600 -86.1300 604.5500（3）P（s）=1.0000s4-6.9000s3-77.2600s2-86.1300s+60

20、4.5500求矩阵多项式中未知数为20时的值； C=polyval（B,20）C = 7.2778e+04（4）把矩阵A作为未知数代入到多项式中； D=polyval（B,A）D = 1.0e+03 * 0.3801 -0.4545 -1.9951 0.4601 -1.9951 0.2093 -1.9951 -2.8880 -0.4545 -4.8978 0.6045 0.4353 0.4353 0.0841 -0.4545 -1.16172.上机练习实验（1）利用 help 检查 roots 指令的使用方法, 然后求 4x4 - 3x3 + 2x2 - x +1=0的所有根 。 help r

21、ootsroots - Polynomial rootsThis MATLAB function returns a column vector whose elements are the roots of the polynomial c. r = roots（c） P=4,-3,2,-1,1;ans = 0.5917 + 0.4864i 0.5917 - 0.4864i -0.2167 + 0.6158i -0.2167 - 0.6158i（2）令 A 是一个维度m x n的矩陣. 解释max（A） 和 min（max（A） 分別是什么意思？max（A）:表示由A矩阵每一列最大值组成的一

22、横向量;min（max（A）:表示得到的横向量之中的最小值。 A=rand（3,4）A = 0.8147 0.9134 0.2785 0.9649 0.9058 0.6324 0.5469 0.1576 0.1270 0.0975 0.9575 0.9706 max（A）ans = 0.9058 0.9134 0.9575 0.9706 min（max（A）ans = 0.9058（3）令 x 是一个维度 n 的向量. 解释find（x 0.8*max（x） 是什么意思？find（x 0.8*max（x）的意思是从向量中找出大于最大值0.8倍的所有数据得位置。 x=1;2;3;4;5;6; a

23、=find（x 0.8*max（x） a =（4）令 x=2 3 5 6 和 y=-1:2 而 z = x.y ，解释z 的值是什么？ x=2 3 5 6 ; y=-1:2y = -1 0 1 2 z = x.yz =0.5000 1.0000 5.0000 36.0000（5）试用解析解和数值解的方法求解微分方程x （t）= -2x（t）-3x（t）+exp（-5t）y（t）=2x（t）-3y（t）-4x（t）-4y（t）-sintx（0）=1,x（0）=2,y（0）=3,y（0）=4 解析解：x,y=dsolve（D2x=（-2）*x-3*Dx+exp（-5*t）,D2y=2*x-3*y-4*Dx-4*Dy-sint,x（0