15年13级matlab实验报告(5)

《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页