matlab题和答案,MATLAB练习题和答案

MATLAB练习题和答案

控制系统仿真实验 Matlab部分实验结果 目 录 实验一 MATLAB基本操作1 实验二 Matlab编程5 实验三 Matlab底层图形控制6 实验四 控制系统古典分析12 实验五 控制系统现代分析15 实验六 PID控制器的设计19 实验七 系统状态空间设计23 实验九 直流双闭环调速系统仿真25 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 再求出它们的乘积矩阵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 选择合适的步距绘制出下面的图形 ，其中 t=[-1:0.1:1]; y=sin(1./t); plot(t,y) 3 对下面给出的各个矩阵求取矩阵的行列式、秩、特征多项式、范数、特征根、特征向量和逆矩阵。 ， ， 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); 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 求解下面的线性代数方程，并验证得出的解真正满足原方程。 (a) ，(b) (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); (4) 输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，执行下列命令，理解其含义 A(2, 3) 表示取A矩阵第2行、第3列的元素； A(:,2) 表示取A矩阵的第２列全部元素； A(3,:) 表示取A矩阵第３行的全部元素； A(:,1:2:3)表示取A矩阵第1、3列的全部元素； A(:,3).*B(:,2)表示A矩阵第3列的元素点乘B矩阵第2列的元素 A(:,3)*B(2,:)表示A矩阵第3列的元素乘以B矩阵第2行 A*B 矩阵AB相乘 A.*B矩阵A点乘矩阵B A^2矩阵A的平方 A.^2矩阵表示求矩阵A的每一个元素的平方值 B/A 表示方程AX=B的解X B./A表示矩阵B的每一个元素点除矩阵A的元素 6在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]，用不同颜色，不同线的类型予以表示，注意坐标轴的比例控制。 t=[0:0.01:2*pi]; y1=cos(t-0.25); plot(t,y1, r-- ) hold on y2=sin(t-0.5); plot(t,y2, k ) 实验二 Matlab编程 1分别用for和while循环结构编写程序，求出 并考虑一种避免循环的简洁方法来进行求和。 (a)j=1;n=0;sum=1; for n=n+1:63 for i=1:n j=j*2; end sum=sum+j; j=1; end sum (b)j=1;n=1;sum=1; while n~=64 i=1; while i> t=[-pi:0.05:pi]; >>