《MATLAB程序设计》实验报告
maxC = 46
>> minC=min(c) minC = -4
>> mid=median(c) mid =
4.5000 >> sort(c) ans =
-4 0 1 2 3 6 6 8 (3) 解释b(2:3,[1,3])的含义，并将生成的数组赋值给d。
b(2:3,[1,3])表示取b矩阵的第2,3行、第1,3列的元素 >> d=b(2:3,[1,3]) d =
12 48 68 3 关系运算与逻辑运算
7、已知a=20,b=-2,c=0,d=1，计算下列表达式：
共64 页，第 21页
23 46 《MATLAB程序设计》实验报告
实验结果：
(1) r1=a>b
>> r1=a>b r1 = 1 (2) r2=a>b&c>d
>> r2=a>b&c>d r2 = 0 (3) r3=a==b*(-10) >> r3=a==b*(-10) r3 =
1
(4) r4=~b|c
>> r4=~b|c r4 = 0
??7218、已知A??9153???2?211?132(1) 求矩阵A的秩(rank)；
?2??2?5?? 13??共64 页，第 22页
《MATLAB程序设计》实验报告
(2) 求矩阵A的行列式(determinant)； (3) 求矩阵A的逆(inverse)；
(4) 求矩阵A的特征值及特征向量(eigenvalue and eigenvector)。
A的秩(rank) >> rankA=rank(A) rankA = 4 A的行列式 >> detA=det(A) detA =
12568 A的逆矩阵 >> invA=inv(A) invA =
0.1744 -0.0303 -0.1050 0.0789 0.0083 0.0173 0.0095 -0.0185 A的特征向量 >> eigA=eig(A) eigA =
4.8554
-0.0125 0.0270 -0.0121 0.0006 0.0911 -0.0311 -0.0103 0.0795 共64 页，第 23页
《MATLAB程序设计》实验报告
12.6460 + 1.8333i 12.6460 - 1.8333i 15.8526
9、已知多项式f1(x)=2x4-x2+3x+2，f2(x)=3x+2，利用MATLAB进行如下计算： 实验结果：
(1) 计算两多项式的乘法和除法；
>> a=[2,-1,3,2]; >> b=[3,2]; >> conv(a,b) ans =
6 1 7 12 4 (2) 求多项式f1(x)的根；
>> roots(a) ans =
0.5000 + 1.3229i 0.5000 - 1.3229i -0.5000
?13?(3) 分别计算多项式当x=3以及x???时的值。 24??X=3时两个多项式的值： >> x=3;
>> y1=polyval(a,x) y1 = 56
共64 页，第 24页
《MATLAB程序设计》实验报告
>> y2=polyval(b,x) y2 =
11
X为题中矩阵时，两个多项式的值： >> x2=[1,3;2,4]; >> y1=polyval(a,x2) y1 =
6 56 20 126
>> y2=polyval(b,x2) y2 =
5 11 8 14
10、用两种方法求下列方程组的解，并比较两种方法执行的时间。
?7x1 ?14x2?9x3?2x4?5x5?3x ?15x?13x?6x?4x12345????11x1 ?9x2?2x3?5x4?7x5?5x ?7x?14x?16x?2x2345?1???2x1 ?5x2?12x3?11x4?4x5?100?200?300 ?400?500>> A=[7,14,-9,-2,5;3,-15,-13,-6,-4;-11,-9,-2,5,7;5,7,14,16,-2;-2,5,12,-11,-4]; >> b=[100,200,300,400,500]'; >> tic,x=rref([A,b]),toc
共64 页，第 25页