matlab练习题和答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
matlab练习题与答案
控制系统仿真实验
Matlab部分实验结果
目录
实验一 MATLAB基本操作、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、 1 实验二 Matlab编程、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、 5 实验三Matlab底层图形控制、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、 6 实验四控制系统古典分析、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、12 实验五控制系统现代分析、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、15
实验六 PID控制器的设计、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、
、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、19 实验七系统状态空间设计、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、23 实验九直流双闭环调速系统仿真、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、25
实验一 MATLAB基本操作
1 用MATLAB可以识别的格式输入下面两个矩阵
1233,,1443678,i,,,,2357,,,,2335542,i,,,,
A,1357B,,,2675342,i,,3239,,,,189543,,,,1894,,
再求出它们的乘积矩阵C,并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后,调
用相应的命令查瞧MATLAB工作空间的占用情况。
A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];
B=[1+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;1,8,9,5,4,3];
C=A*B;
D=C(4:5,4:6);
whos
Name Size Bytes Class Attributes
A 5x4 160 double
B 4x6 384 double complex
C 5x6 480 double complex
D 2x3 96 double complex 2 选择合适的步距绘制出下面的图形
,其中 sin(/)1tt,,(,)11
t=[-1:0、1:1];
y=sin(1、/t);
plot(t,y)
1
3 对下面给出的各个矩阵求取矩阵的行列式、秩、特征多项式、范数、特征根、特征向量与
逆矩阵。
5765753500、、,,,,
,,,,71087833410、,,,,, B,A,,,,,681090910315,、
,,,,579100037193、、,,,,
12343324,,,,,,
,,,,56785518,,,,,, C,D,,,,,910111211857,
,,,,1,,,,,,,
A=[7、5,3、5,0,0;8,33,4、1,0;0,9,103,-1、5;0,0,3、7,19、3];
B=[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10]; C=[1:4;5:8;9:12;13:1rtf6];
D=[3,-3,-2,4;5,-5,1,8;11,8,5,-7;5,-1,-3,-1];
det(A);det(B);det(C);det(D);
rank(A);
rank(B);
rank(C);
rank(D);
a=poly(A);
b=poly(B);
c=poly(C);
d=poly(D);
norm(A);
2
norm(B);
norm(C);
norm(D);
[v,d]=eig(A,'nobalance');
[v,d]=eig(B,'nobalance');
[v,d]=eig(C,'nobalance');
[v,d]=eig(D,'nobalance');
m=inv(A);
n=inv(B);
p=inv(C);
q=inv(D);
4 求解下面的线性代数方程,并验证得出的解真正满足原方程。
72124,1321390,,,,,,,,
,,,,,,,,721264,915327,,,,,,,,,(a),(b) X,X,,,,,,,,,91532117,,,,221151 ,,,,,,,,,,,,2211521132130,,,,,,,,
(a)
A=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13];
B=[4;7;-1;0];
X=A\B;
C=A*X;
(b)
A=[1,3,2,13;7,2,1,-2;9,15,3,-2;-2,-2,11,5];
B=[9,0;6,4;11,7;-2,-1];
X=A\B;
C=A*X;
5、
(1) 初始化一10*10矩阵,其元素均为1 ones(10,10);
(2) 初始化一10*10矩阵,其元素均为0 zeros(10,10);
(3) 初始化一10*10对角矩阵
v=[1:10];
diag(v);