河北工业大学MATLAB实验二

1、1 / 21程序设计实验报告、实验目的：1、掌握矩阵和数组的表示与赋值方法2、了解字符串、结构和单元等数据类型在中的使用3、掌握中基本的数值运算，了解基本统计函数的使用4、掌握多项式基本运算以及线性方程组的求解方法、实验内容：1、利用“：”操作符和函数生成数组0,6,12,18,24,30,36,42,并将数组a分别转化为4*2和2*4的矩阵-3681 12、在中输入矩阵A =-572217，并使用回答以下的问题：6916-121513-210 一(1) 创建一个由A中第2列到第4列所有元素组成的4X3数组B;(2) 创建一个由A中第3行到第4行所有元素组成的2X4数组C;(3) 创建一个由A

2、中前两行和后三列所有元素组成的2X3数组D;-5 6 15和向量6 8(4) 根据A利用单下标方法和双下标方法分别创建向量-5 61,并利用向量a和b生成矩阵E= 6 8 ;15 1 一利用“”删除矩阵A的第二行和第三列3、利用()函数和()函数生成如下矩阵:0 0 0 0 00 1110A= 0 1 3 1 00 1110'0 0 0 0 0_|-4-1 24、生成矩阵A= -303:-214J(1) 取出所有大于0的元素构成一个向量a (可推广到从一个矩阵里筛选出 符合条件的元素组成一个向量)。(2) 将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩

3、阵中不符合条件的全用0来表示，符合条件的值不变)。5、已知a = 12|(4 5-111c =0，d = 8-2-331T2 3 J456712下列运算是否合法，为什么？如合法，写出计算结果?(1) 1'2* b34* d 5='* d 6.* b7(8) 8.*cb(10) 10.A2(11) 11A2(12) 1仁 2.Aa134.1-7846、已知 a=5 1244 ,b=122438,c=1,0,8,3,624,23,46,67 827 j.68-53j(1) 求*.*A2.A2的结果，并观察运算结果。(2) 求c中所有元素的平均值、最大值、中值，并将c中元素按从小到大

4、顺序排列。 解释b(2:3,1,3)的含义，并将生成的数组赋值给do关系运算与逻辑运算7、已知20201，计算下列表达式： r1>b r2»d r3*(-10)r4179215138、已知A =-2-211132-2-2513 一(1)求矩阵A的秩();求矩阵A的行列式()；求矩阵A的逆()；(4)求矩阵A的特征值及特征向量()。9、已知多项式fi(x)=2x42+32, f2(x)=32,禾用进行如下计算:(1)计算两多项式的乘法和除法；求多项式f1(x)的根；_131 分别计算多项式当3以及x = |时的值。2 4 一10、用两种方法求下列方程组的解，并比较两种方法执行的时

5、间。7% +14x2 9x3 2x4 +5xs =1003x1 - 15x2 - 13x3 - 6X4 - 4x5 - 200* -11x1 9x2 2x3 +5x4 +7卷=3005x1 +7x2 +14x3 +16x4 - 2x5 =400-2x1 5x2 12x3 -11x4 -4卷=50011、 建立一个元胞数组，要求第一个元素为4X4的魔术矩阵，第二个元素 为18.66,第三个元素为''，第四个元素为等差数列7,9,11,13-990计算第一个元 胞第4行第2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元胞的 第10个元素。12、建立一个结构体的数组，包括3个

6、人，字段有姓名，年龄，分数，其中 分数由随机函数产生一个3行10列的数据表示了有10门课程，每门课程有三个 阶段的分数。问题是：(1) 如何找到第 2 个人的分数并显示出来；(2) 如何找到第 2 个人的每门课程 3个阶段的平均分数并显示出来；(3) 所有同学的 10 门课程的每门课程的平均分如何计算出来？要求放到一 个数组里；(4) 找到这个班所有同学的姓名放到一个元胞数组里；13、完成以下操作：(1) 已有 20 10/20/10'，查找字符串中 20 出现的位置。(2) 有字符串s1= ，s2= 比较两个字符串的前8个字符是否相同，相 同输出 1 ，否则输出 0。(3) 统计字符

7、串s1=''中出现的次数，并替换为。三、源程序和实验结果1.>> 0:6:42a =06121824303642>> (0,42,8)06121824303642>> (a,4,2)02463012361842>> (a,2,4)012243661830422.>> 3 6 8 15 7 22 17;6 9 16 -12;15 13 -21 0 A =3681-5722176916-121513-210(1)>> (:,2:4)B =68172217916-1213-210(2)>> (3:4,

8、:)C =69 16 -121513 -210(3)>> (1:2,2:4)D =681722 17(4)>> A(2) A(5) A(4) a =-5 6 15>> A(5) A(9) A(13) b =>> A(2,1) A(1,2) A(4,1) a =-5 6 15 >> A(1,2) A(1,3) A(1,4) b =6 8 1>> a' b'E =-5 66 815 1 >> A(2,:)=A =36816916-121513-210>> A(:,3)=A =36169-

9、12151303.>> (5)A =0000000000000000000000000>> A(2:4,2:4)(1)7 / 210000001110011100111000000>> A(3,3)=3*(1)A =00000011100131001110000004.>> -4 -1 23 0 32 1 4 A =-4-12-303-214(1) >> (A>0) k =6789>> (k)a =4(2)>> (A<=0) m =12345>> A(m)=0A =0020030145.

10、>> 1 2 3;4 5 6a =1 23456>> 2 4 -1;1 3 5 b =2 4-1135>> 1;02 c =10-2 >> 1 4 7;8 5 2;3 6 0d =147852360(1)>> 1'1 =142536(2)不合法，两个矩阵的维数不一致>> 2*b ? >(3)>> 33 =3625811(4)>> 4*d4 =312222404913(5)>> 5='*d5 =312222404913-5-87(6)>> 6.*b6 =2

11、 8-341530(7)>> 77 =0.50000.5000-3.00004.00001.66671.2000(8)不合法，元素个数不一样 >> 8.*c? >(9)>> 9.b9 =2.00002.0000-0.33330.25000.60000.8333(10)>> 10.A210 =1 4 9162536(11)不合法，维数不一致>> 11A2? > (12)>> 12=2.Aa12 =2481632646.>> 1 3 4;5 12 44;7 8 27a =134512447827>

12、> -7 8 4;12 24 38;68 -5 3-712 / 21# / 21122438# / 21-568 >> 1,0,8,3,6,24,23,46,6# / 21# / 216 2 -4 2346 6(1)>>-611817368275330>> a*b301603101108>>a.*b188311313060841313 / 21-72416602881672476-4081>>0.09660.09450.0080-3.61251.5838-0.5778-1.99170.9414-0.2682>>-0.

13、14290.37501.00000.41670.50001.15790.1029-1.60009.0000>> aA2447124437351117362363331109>> a.A219162514419364964729(2)>> (c,2) y =9.1000>> (c)y =9.100014 / 21>> y1(c)y1 =9.1000>> y2(c)y2 =46>> y3(c)y3 =4.5000>> (c)Y =469-4012366823I =7 21645103 8(3)选出b矩阵

14、的2到3行，1列和3列的元素>> (2:3,1,3)d =12386837.>> 20201a =20b =-2c =0d =1(1)15 / 21>> r1>b r1 =1(2) >> r2>>d r2 =0(3) >> r3*(-10) r3 =1(4) >> r4 r4 =08.>> 7 2 1 -2;9 15 3 -22 -2 11 5;1 3 2 13 A =721-29153-2-2 -2 11 513213(1) >> (A)4(2) >> (A)1256

15、8(3)>> (A)0.1744-0.0303-0.01250.02700.10500.0789-0.01210.00060.00830.01730.0911-0.03110.0095-0.0185-0.01030.0795(4)>> (A) C =0.76290.0919 + 0.0640i0.0919 - 0.0640i-0.02990.62230.6087 + 0.0276i0.6087 - 0.0276i0.26370.0807-0.7474-0.74740.64340.15540.0342 - 0.2374i0.0342 + 0.2374i0.7180D =4

16、.8554000012.6460 + 1.8333i000012.6460 - 1.8333i000015.85269.>>f1=21,3,2;>>f2=3,2;(1)>>(f12)x =6 1 7 12 4 >> (f12)y =0.6667 -0.7778 1.518517 / 210 -0.00000 -1.037018 / 21# / 21(2)(f1)x =0.5000 + 1.3229i0.5000 - 1.3229i-0.5000(3)>> x1=3;>> X2=1 3;2 4X2 =1 32 4>&g

17、t; y1(f11)y1 =56>> y2(f12)y2 =6 5620 126>> y3(f21)y3 =11>> y4(f22)y4 =5 118 1410.>> 7 14 -9 -2 5;3 -15 -13 -6 -411 -9 -2 5 7;5 7 14 16 -22 5 12 -11 -4 A =# / 21714-9-253-15-13-6-411-9-257571416-2-2512-11-4>> 100;200;300;400;500B =100200300400500(1)>> B x =370.9455

18、-224.7276238.3333-138.1891503.6378(2)>> ()x =0001.000000 370.94440 -224.72730 238.33330 -138.18921.0000 503.63641.0000000 1.0000000 1.0000000000两个方法速度差不多，肉眼不能分辨11.19 / 21A =4x4 18.6600>> (4),18.66,'',7:2:99>> A1,11x47 20 / 21# / 2116231351110897612414151>> A1,1(4,2)1,

19、21,3(2)1,4(10)154.660012.>> ''>> 20;>> (3,10,60,100)H: 20: 3x10 >> (2)''>> (2)20;>> (2)(3,10,60,100)>> (2)H: 20: 3x10 >> (3)''>> (3)20;>> (3)(3,10,60,100);(1)>> (2)A =886188619180898764736163737792789086808371939875678671669969(2)>> (A)Y =1 873.3333 72.3333 86.3333 71.0000 83.3333 81.3333 83.333379.66679 1081.0000 75.0000( 3) >> (1)(2)(3