matlab实验报告

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

MATLAB实验报告

一、实验过程

1.矩阵运算

(1)、矩阵的乘法

>> A=[1 2;3 4];

>> B=[5 5;7 8];

>> A^2*B

ans =

105 115

229 251

(2)、矩阵的除法

>> A=[1 2 3;4 5 6;7 8 9];

>> B=[1 0 0;0 2 0;0 0 3];

>> A\B

ans =

1.0e+16 *

-0.4504 1.8014 -1.3511

0.9007 -3.6029 2.7022

-0.4504 1.8014 -1.3511

>> A*inv(B)

ans =

1.0000 1.0000 1.0000

4.0000 2.5000 2.0000

7.0000 4.0000 3.0000

(3)、矩阵的转置及共轭转置

>> A=[5+i,2-i,1;6*i,4,9-i];

>> A.'

ans =

5.0000 + 1.0000i 0 +

6.0000i

2.0000 - 1.0000i 4.0000

1.0000 9.0000 - 1.0000i

>> A'

ans =

5.0000 - 1.0000i 0 -

6.0000i

2.0000 + 1.0000i 4.0000

1.0000 9.0000 + 1.0000i

(4)、使用冒号选出指定元素

>> A=[1 2 3;4 5 6;7 8 9];

>> A(1:2,3)

ans =

3

6

>> A(2:3,:)

ans =

4 5 6

7 8 9

(5)复数矩阵的生成

>> syms a1 a2 a3 a4 b1 b2 b3 b4 real;

>> Z1=[a1+j*b1 a2+j*b2;a3*exp(j*b3) a4*exp(j*b4)];

>> a1=3;a2=-2;a3=9;a4=23;b1=5;b2=3;b3=6;b4=33;

>> Z1=subs(Z1,'[a1 b1 a2 b2 a3 b3 a4 b4]',[a1 b1 a2 b2 a3 b3 a4 b4]) Z1 =

[ 3 + 5*i, - 2 + 3*i]

[ 9*exp(6*i), 23*exp(33*i)]

2 、多项式

(1)>> p=[1 0 2 4];

>> roots(p)

ans =

0.5898 + 1.7445i

0.5898 - 1.7445i

-1.1795

(2)

特征多项式

>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4]

A =

1.2000 3.0000 5.0000 0.9000

5.0000 1.7000 5.0000

6.0000

3.0000 9.0000 0 1.0000

1.0000

2.0000

3.0000

4.0000

>> B=sym(A);

>> poly(B)

ans =

x^4 - (69*x^3)/10 - (3863*x^2)/50 - (8613*x)/100 + 12091/20

未知数为20的值

>> x=20

x =

20

>> y=x^4 - (69*x^3)/10 - (3863*x^2)/50 - (8613*x)/100 + 12091/20 y =

7.2778e+04

矩阵A作为未知数代入到多项式中

>> p=[1 0 2 4];

>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4];

>> y=polyval(p,A)

y =

8.1280 37.0000 139.0000 6.5290

139.0000 12.3130 139.0000 232.0000

37.0000 751.0000 4.0000 7.0000

7.0000 16.0000 37.0000 76.0000

3、基本绘图命令

(1)

>> t=0:0.1:2*pi; >> y=cos(t);

>> plot(t,y);

>> title('y=cos(t)')

(2)

>> t=0:0.1:2*pi; >> y1=cos(t-0.25);

>> y2=sin(t-0.5);

>> plot(t,y1,t,y2)

>> legend('cos(t-0.25)','sin(t-0.5)')

4、基本绘图控制

>> t=0:0.5:4*pi;

>> y=10*sin(t);

>> plot(t,y,'r-.+');

>> xlabel('时间');

>> ylabel('正弦');

>> title('正弦曲线');

>> legend('sin(t)');

>> grid on

相关文档
最新文档