matlab实验报告
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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