MATLAB仿真试题
《Matlab与系统仿真》期末考试试卷A卷(2013)
第1页 共2页 第2页 共2页
红河学院2013-2014学年秋季学期 《Matlab 与系统仿真》 课程期末考试试卷 卷别： A 卷 考试单位：工学院 考试日期： 年 月 日
2小题采用M 语言编程； 0分记； 3、提交的压缩包包括M 文件、mdl 模型、Word 文档（仿真输出曲线），压缩包以自己的姓名、学号、班级命名；发送至邮箱：791705181@
一、设计题(每小题50分，共100分) 1、精馏过程是在一定物料平衡和能量平衡基础上进行的。
在精馏塔的连续操作过程中应做到物料平衡、气－液平衡和热量平衡，这三个平衡互相影响，互相制约。
现以蒸汽流量作为副变量，精馏塔釜温度为主变量，构成精馏塔釜温度与蒸汽流量串级控制系统。
已知主被控对象的数学模型为G(s)=1.25e -44s /(25s+1)2，副被控对象的数学模型为G(s)=e -4s /(12s+1)，精馏塔釜设定温度为1300℃，副控制器选择Ｐ作用，主控制器选择PID 作用，试整定串级控制器的参数为最佳值，使系统的最大偏差≤±1.5℃。
（采用SIMLINK 仿真，参数自己整定） 2、设一被控对象17125.050)(2++=s s s G ，用增量式PID 控制算法编写仿真程序（输入分别为单位阶跃、正弦信号，采样时间为1ms ，控制器输出限幅:[-5,5]，仿真曲线包括系统输入、输出及误差曲线。
（采用M 语言编程，参数自己整定）。
MATLAB与控制系统仿真练习题及答案 (2)
MATLAB与控制系统仿真练习题1、在同一坐标系内，画出一条正弦曲线和一条余弦曲线，要求正弦曲线用红色实线、数据点用“+”号显示；余弦曲线用黑色点线、数据点用“*”显示。
x=0:1 /10:2*pi;y=sin(x);plot(x,y,'r--+')hold ont=cos(x);plot(x,t,'k-.*')>>2、绘制单位圆。
a=0:1/10:2*pi;x=cos(a);y=sin(a);plot(x,y)axis squaregrid on3、使用极坐标函数绘制)3/ρ=的曲线。
8θsin10(t=0:0.01:2*pi;>> y=8*sin((10*t)/3);>> polar(t,y)4、绘制方程⎪⎩⎪⎨⎧===)cos()sin(t z t y t x 在[]π20，=t 区间的三维曲线。
t=0:pi/50:2*pi;x=t;y=sin(t);z=cos(t);plot3(x,y,z)5、用不同的线型和颜色在同一坐标内绘制曲线x 5.0e 21y -=、)2sin(2y x π=的图形。
x=0:1/10:10;y1=2*exp(-0.5*x);y2=sin(2*pi*x);plot(x,y1,'r');>> hold onplot(x,y2,'b')6、绘制二维正态分布密度函数)(212221),(y x e y x f +-=π的三维图形。
syms x yz=(1/(2*pi))*exp(-0.5*(x.^2+y.^2));ezsurf(x,y,z)7、选择合适的步长绘制下列函数的图形。
（1）x y cos =，⎥⎦⎤⎢⎣⎡-∈2,2ππx （2）x xe y = x=-2/pi:1/20:2/pi; >> y=sqrt(cos(x)); >> plot(x,y)>> hold on>> y=x.*exp(x);>> plot(x,y)>>8、在同一坐标下绘制函数x 、2x 、3x -、)sin(x x 在()π,0的曲线。
MATLAB控制系统仿真试题B卷.doc
广东技术师范学院2010—2011学年度第 一学期期未考查试卷科 目：控制系统仿真 （B ）卷 考试形式：上机考试 考试时间: 120 分钟 所属学院：自动化学院 班级： 姓名： 学号：答题注意事项：以下题目请大家在计算机上完成，并将结果写成电子文档上交。
电子文档中应包含以下几项内容：仿真程序代码、仿真结果截图、仿真实验结论。
1．设单位负反馈系统的开环传递函数为)5)(2()(2++=s s s Ks G k ，试确定分别使闭环系统稳定的开环增益的取值范围。
（20分） 程序如下： num=1;den=conv([1 2 0 0],[1 5]); G=tf(num,den); k=0:0.05:100; rlocus(G) [K,P]=rlocfind(G)题 目 一 二 三 四 五 总分 标准分数 20 20 20 20 20 100 实得分数 评卷人装订线考 生 答 题 不 得 超 过 此 线Select a point in the graphics window selected_point =0.0024 + 0.3230iK =1.0607P =-4.9857-2.08380.0347 + 0.3176i0.0347 - 0.3176i运行结果：K=1.0607;t=0:0.05:100;G0=feedback(tf(K*num,den),1); step(G0,t)由上图分析可得，系统的临界稳定增益是1.0607，因此该闭环系统稳定的开环增益的取值范围是0~1.0607,随着增益的增加系统趋于发散统系统 2、试求以下单位负反馈系41)(2++=s s s G k 的单位阶跃响应曲线，并求其动态性能指标r p s t t t %σ。
num=1; den=[1 1 4]; G=tf(num,den); Gk=feedback(G ,1) figure(1);运行结果：2．已知单位负反馈系统的开环传递函数为)5(500)(+=s s Ks G k ，为使系统具有如下性能指标：速度误差系数1100-=s K v ，相角裕度 45≥γ。
MATLAB系统仿真试题复习资料,有答案
PART I 教材复习要点第一章 MATLAB 入门测试1.11.MATLAB 的命令窗口的作用是什么？编辑/调试窗口？图象窗口？ MATLAB 命令窗口是MATLAB 起动时第一个看到的窗口，用户可以在命令窗口提示符">>"后面输入命令，这些命令会被立即执行。
命令窗口也可以用来执行M 文件。
编辑/调试窗口是用来新建，修改或调试M 文件的。
图像窗口用来显示MATLAB 的图形输出。
2.列出几种不同的得到MATLAB 帮助的方法。
（1）在命令窗口中输入help <command_name>，本命令将会在命令窗口中显示关于些命令的有信息。
（2）在命令窗口中输入lookfor <keyword>，本命令将会在命令窗口中显示所有的在第一注释行中包含该关键字的命令和函数。
（3）通过在命令窗口输入helpwin 或helpdesk 启动帮助浏览器，或者是在启动板中选择"Help"。
帮助浏览器包含了基于超文本的MATLAB 所有特性的描述，HTML 或PDF 格式的在线手册，这是MATLAB 最全面的帮助资源。
3.什么是工作区？在同一工作区内，你如何决定它里面存储了什么？工作区是命令、M 文件或函数执行时被MATLAB 使用的变量或数组的收集器，所有命令都在命令窗口（所有的脚本文件也是从命令窗口执行）共享公共工作区，因此它们也共享所有变量，工作区的内容可以通过whos 命令来查看，或者通过工作区浏览器来图形化地查看。
4.你怎样清空MATLAB 工作区内的内容？要清除工作区的内容，只需在命令窗口中输入clear 或clear variables 即可。
第二章MATLAB 基础测试2.11.数组,矩阵,向量有什么区别?数组是在内存中被组织成行和列的数据集合，只有一个名称，数据要通过在数组名后面圆括号里加上表示数据所在行和列的数字来访问。
MATLAB与控制系统仿真大作业
>> G1=tf(num,den);
>> G1=zpk(G1)
G1 =
5 (s+0.6) ----------------(s+3) (s+2) (s+1)
>> num=[2 1]; den=[1 2.9 1]; G2=tf(num,den); >> G2=zpk(G2)
G2 =
2 (s+0.5)
>> xlabel('x') ylabel('exp(-2.*x)')
三、simulink 建模（每题 10 分，共 20 分）
1、已知单位负反馈的开环传递函数为 G(s) =
2 s2 + 4s
，试利用
simulink
建立系在单位阶跃输入作用下的模型。
要求答案包括：（1）simulink 建模结构图；
（2）在同一个坐标中的阶跃信号和响应曲线图。
2、已知系统的开环传递函数为
G1(s)
=
2s2 + 5s + 6 s2 + 2s + 3
、G2 (s)
=
s2
s+6 + 7s +1
，H
(s)
=
5(s + 2) s +10
求：建立 Simulink 仿真模型，并求出其系统在单位阶跃响应；
要求答案包括：（1）simulink 建模结构图；
x=
-2.9709
0.5491
3.6000
0.0509
2. 已知下列矩阵
2 3 1
−1 3 5
A
=
《MATLAB与系统仿真》期末考试试题
《MATLAB 与系统仿真》期末考试试题一、 填空题（1’ ×15=15’）1、 Matlab 预定义变量ans 表示： ；eps 表示： 。
2、 在Matlab 开发环境是一个集成的工作环境，包括 、 、 、 和 等。
3、 函数rand 表示： ；randn 函数表示： ；size 函数表示： 。
4、 Matlab 命令窗口中的“>>”符号表示 。
5、 实现Fourier 变换在Matlab 中的对应函数为： ；Laplace 变换的函数： 。
6、 建立结构数组或转换结构数组的函数为： 。
二、 简答题（3’×6=18’）1、 分析以下函数执行过程，回答问题：function c = testarg1(a, b)if (nargin == 1)c = a .^ 2;elseif (nargin == 2)c = a + b;end如果用户使用一个参数调用函数testarg1，函数执行什么功能？如果使用两个参数调用函数testarg1，函数又执行什么？2、 简述脚本M 文件和函数M 文件的区别？3、 简述绘制二维图形的一般步骤？4、 简述MATLAB 命令窗的主要作用？5、 打开和关闭文件的函数是什么？分别说明其作用。
6、 Matlab 提供了哪四类程序控制结构方式？三、程序设计及分析（1-9每题5’，10题22’，共67’）1、 编程实现函数3()25f x x x =--的根？2、 编程求解由方程2210,50,240x y z x y z x y z -+=+-=-+=构成的线性方程组？3、 将下列字符矩阵S 转换为数字矩阵（'1''2''3''4'S ⎡⎤=⎢⎥⎣⎦）？ 4、 编程实现323()5327()483b x x x x a x x x +-+=-++的部分分式展开？ 5、 把矩阵A ，B 分别在垂直方向，水平方向连接，A=[1 2 3]，B=[4 5 6]?6、 分析以下程序执行后的效果?x=-2:0.2:2;y=x;[x,y]=meshgrid(x,y);z=x.*exp(-x.^2-x.^2);mesh(x,y,z);7、 编写一个实现n 阶乘的函数。
Matlab与系统仿真试题及答案
Matlab 与系统仿真试题及答案一、填空题（每空2分，共30分）1、单位取样序列定义为________________。
2、单位阶跃序列定义为__________________。
3、取样定理中，采样频率和原始信号频率之间至少应该满足公式_______________。
4、实现IIR 滤波器可以采用三种不同的结构：_______________、__________________和_______________。
5、对于理想的滤波特性()jwH e ，相应的()d h n 一般为无限时宽，从而使非因果的。
为此，使显示一般通过如下公式__________实现加窗处理，相应的傅里叶表达形式为_________________。
6、设计IIR 滤波器的基础是设计模拟低通滤波器的原型，这些原型滤波器有：_________、________和______________低通滤波器。
7、实现FIR 滤波器可以采用以下4中结构：________________,_______________,______________和____________。
二、阅读下列Matlab 程序，并回答完成的功能(每题5分，共10分)1、f=[00.60.61];m=[1100];b=fir2(30,f,m);[h ,w]=freqz(b,1,128);Plot(f,m,w/pi,abs)本段程序的功能是：2、[b ,a]=butter(9,300/500,’high’);Freqz(b,a,128,1000)本段程序的功能是：三、用Matlab 语言编程实现下列各小题（每小题15分，共60分）1、产生序列0.3,16p s A dB πΩ==55n -≤≤，并绘出离散图，标注图形标题为“Sequence in Example 3.1a”。
实现程序如下：2、设11()23X z z z-=++，212()2435X z z z z -=+++，求312()()()X z X z X z =+。
MATLAB控制系统各种仿真例题(包括simulink解法)
一、 控制系统的模型与转换1． 请将下面的传递函数模型输入到matlab 环境。
]52)1)[(2(24)(32233++++++=s s s s s s s G )99.02.0)(1(568.0)(22+--+=z z z z z H ，T=0.1s >> s=tf('s');G=(s^3+4*s+2)/(s^3*(s^2+2)*((s^2+1)^3+2*s+5));GTransfer function:s^3 + 4 s + 2------------------------------------------------------s^11 + 5 s^9 + 9 s^7 + 2 s^6 + 12 s^5 + 4 s^4 + 12 s^3>> num=[1 0 0.56];den=conv([1 -1],[1 -0.2 0.99]);H=tf(num,den,'Ts',0.1)Transfer function:z^2 + 0.56-----------------------------z^3 - 1.2 z^2 + 1.19 z - 0.992． 请将下面的零极点模型输入到matlab 环境。
请求出上述模型的零极点，并绘制其位置。
)1)(6)(5()1)(1(8)(22+++-+++=s s s s j s j s s G )2.8()6.2)(2.3()(1511-++=----z z z z z H ，T=0.05s>>z=[-1-j -1+j];p=[0 0 -5 -6 -j j];G=zpk(z,p,8)Zero/pole/gain:8 (s^2 + 2s + 2)--------------------------s^2 (s+5) (s+6) (s^2 + 1)>>pzmap(G)>> z=[0 0 0 0 0 -1/3.2 -1/2.6];p=[1/8.2];H=zpk(z,p,1,'Ts',0.05)Zero/pole/gain:z^5 (z+0.3125) (z+0.3846)-------------------------(z-0.122)Sampling time: 0.05>>pzmap （H ）num=[0,7.1570,-6.4875 ];den=[1,-2.2326,1.7641,-0.4966];sysd=tf(num,den,0.05,'variable','z^-1')Transfer function:7.157 z^-1 - 6.487 z^-2-----------------------------------------1 - 2.233 z^-1 + 1.764 z^-2 - 0.4966 z^-3Sampling time: 0.05二、 线性系统分析1. 请分析下面传递函数模型的稳定性。
MATLAB与控制系统仿真练习题(含图)及答案
MATLAB 与控制系统仿真练习题（含图）1、已知函数x x e x f x sin cos )(-=，作出函数的大致图像。
>> syms x>> y=exp(x)*cos(x)-sin(x); >> ezplot(y)2、求下列极限：（1）30sin lim xx x x -→ >> syms x>> y=(x-sin(x))/(x^3);>> limit(y,x,0)ans =1/6(2) xx x ⎪⎭⎫ ⎝⎛+∞→11lim >> y=(1+1/x)^x;>> limit(y,x,inf)ans =exp(1)3、求下列函数的导数：(1)x e y x sin =>> syms x>> y=exp(x)*sin(x);>> diff(y,x)ans =exp(x)*sin(x)+exp(x)*cos(x)(2) x e x x y 22sin +=>> syms x>> y=sin(x)+x^2*exp(2*x);>> diff(y,x)ans =cos(x)+2*x*exp(2*x)+2*x^2*exp(2*x)4、求.)1(532⎰-dx x x 和.sin ⎰xdx e x（1）.)1(532⎰-dx x xsyms x>> int(x^2*(1-x^3)^5)ans =-1/18*x^18+1/3*x^15-5/6*x^12+10/9*x^9-5/6*x^6+1/3*x^3（2）.sin ⎰xdx e x>> int(exp(x)*sin(x))ans =-1/2*exp(x)*cos(x)+1/2*exp(x)*sin(x)5、求.)(102⎰-dx x x 和.1102⎰-dx x x (1) .)(102⎰-dx x x>> syms x>> int(x-x^2,0,1)ans =1/6(2) .1102⎰-dx x x>> syms x>> int(x*(1-x^2)^0.5,0,1)ans =1/36、已知二元函数),(cos )sin(2xy xy z +=试求y x z x z y z x z ∂∂∂∂∂∂∂∂∂222,,,。
MATlab 仿真部分习题答案
[4.1]控制系统结构如图4.1所示（1）利用MATLAB对以上单位负反馈控制系统建立传递函数；（2）将第一问中求得的传递函数模型转化为零极点增益形式和状态空间形式。
解:（1）num=[2 2];den=[1 2 1];[num1,den1]=cloop(num,den);sys=tf(num1,den1)程序运行结果如下：Transfer function:2 s + 2-------------s^2 + 4 s + 3（2）[z,p,k]=tf2zp(num1,den1);g_zpk=zpk(z,p,k);[A,B,C,D]=zp2ss(z,p,k);g_ss=ss(A,B,C,D)程序运行结果如下：z = -1 p = -3 -1 k = 2a = x1 x2x1 -4 -1.732x2 1.732 0b = u1x1 1x2 0c = x1 x2y1 2 1.155d = u1y1 0K[5.1]设单位负反馈控制系统的开环传递函数为G(S)=+s(s^2+7s17)（1）试绘制k=10、100时闭环系统的阶跃响应曲线，并计算稳态误差、上升时间、超调量和过渡过程时间；（2）绘制k=1000时闭环系统的阶跃响应曲线，与k=10、100时所得的结果相比较，分析增益系数与系统稳定性的关系；解：（1）k=10时，K=100时，K=10时，利用MA TLAB工作区输入程序：num=[10];den=[1,7,17,0];[z,p,k]=tf2zp( num,den);运行得z,p,k的值p= -3.5+2.1794*i -3.5-2.1794*i k=10G=zpk([ ],[-3.5+2.1794*i,-3.5-2.1794*i],10)；c=dcgain(G);[y,t]=step(G);plot(t,y)[Y,K]=max(y);timetopeak=t(k);percentovershoot=100*(y-c)/cn=1;while y(n)<c n=n+1;end risetime=t(n) i=length(t);while(y(i)>0.98*c&y(i)<1.02*c) i=i-1;End setllingtime=t(i) 运行程序结果为：稳定值c=0.5882，响应时间setllingtime=1.0096s,上升时间risetime=1.1989s,最大峰值时间timetopeak=1.4356,超调量percentovershoot=0.555% 同理得k=100时，稳定值c=5.882，响应时间setllingtime=1.0096s, 上升时间risetime=1.1989s,最大峰值时间timetopeak=1.4356,超调量percentovershoot=0.555%（3）k=1000时，由其响应曲线可知，增益系数越大，其稳定性越差。
MATLAB模拟试题
一、填空题1. MA TLAB命令窗口中可用clear命令清除工作区中的变量;用clc 命令清除命令窗口中的内容。
2. MA TLAB中的运算符包括数学运算符、关系运算符、逻辑运算符。
3 MA TLAB中的M文件有命令文件和函数文件。
4。
MATLAB中的函数类型有嵌套函数、内嵌函数、、私有函数和匿名函数。
5。
MATLAB中用系统预定义的特殊变量ans 表示函数输入参数个数。
6。
MATLAB中用系统预定义的特殊变量inf 表示无穷大。
7. MA TLAB中用myarray 符号定义单元数组。
8。
MATLAB中用hold 命令打开图形保持功能。
9 MA TLAB中如果一个变量名为123，该变量名是否合法否（填是或否）.10. MA TLAB中用inv 函数完成矩阵的求逆运算，用det 函数求矩阵的行列式。
二、阅读下述MATLAB脚本程序A=［3 5 7 9；0 2 4 1；6 4 4 8］;B=ones（3,3）；b=A([1 2],2：3）；B（4，4）=5；A（[1，2］，:）=［］；C=［A；B（:,1：3）'］；D=reshape（C,2，8)；写出b、A、B、C、D的结果。
〉〉 bb =5 72 4〉> AA =6 4 4 8〉〉BB =1 1 1 01 1 1 01 1 1 00 0 0 5〉〉 CC =6 4 4 81 1 1 01 1 1 01 1 1 0〉> DD =6 1 4 1 4 1 8 01 1 1 1 1 1 0 0三、已知矩阵A 为4×4的正态分布的随机矩阵，编写脚本程序实现矩阵B, B 矩阵为对角矩阵，主对角线元素为A矩阵的主对角线元素加1，并将矩阵B 存入文件result 。
mat 文件中。
>〉 A=randn(4，4)〉〉 v=diag(A ）〉> B=diag （v+1,0)save result.mat B四、已知多项式532)(3++=x x x a 和13)(24++=x x x b ,编写脚本程序计算 1） 多项式a(x ）+b （x ）2) a （x)＊b （x)2) 多项式a(x ）/b(x ）的部分分式展开。
matlab模拟试题及答案
matlab模拟试题及答案MATLAB模拟试题及答案1. 编写一个MATLAB函数，计算并返回一个向量中所有元素的和。
```matlabfunction sumVector = sumVectorElements(vector)sumVector = sum(vector);end```2. 给定一个矩阵A，编写一个MATLAB脚本，找出矩阵中的最大元素及其位置。
```matlabA = [1, 2, 3; 4, 5, 6; 7, 8, 9];[maxValue, maxIndex] = max(A(:));[maxRow, maxCol] = ind2sub(size(A), maxIndex);disp(['最大元素为: ', num2str(maxValue), '，位于第',num2str(maxRow), '行，第', num2str(maxCol), '列。
']);```3. 编写一个MATLAB函数，实现一个向量元素的排序。
```matlabfunction sortedVector = sortVector(vector)sortedVector = sort(vector);end4. 给定两个向量X和Y，编写一个MATLAB脚本，计算这两个向量的点积。
```matlabX = [1, 2, 3];Y = [4, 5, 6];dotProduct = dot(X, Y);disp(['X和Y的点积为: ', num2str(dotProduct)]);```5. 编写一个MATLAB函数，计算一个矩阵的行列式。
```matlabfunction determinant = calculateDeterminant(matrix)determinant = det(matrix);end```6. 给定一个向量，编写一个MATLAB脚本，找出向量中的最大值和最小值。
matlab仿真实例100题
matlab仿真实例100题Matlab是一种强大的数学软件，广泛应用于科学计算、数据分析和工程仿真等领域。
在学习和使用Matlab的过程中，通过实例的方式进行仿真练习是一种非常有效的学习方法。
下面将给出100个Matlab仿真实例题目，帮助读者更好地掌握Matlab的使用。
1. 编写一个程序，计算并输出1到100之间所有奇数的和。
2. 编写一个程序，计算并输出1到100之间所有偶数的乘积。
3. 编写一个程序，计算并输出1到100之间所有素数的个数。
4. 编写一个程序，计算并输出1到100之间所有整数的平方和。
5. 编写一个程序，计算并输出1到100之间所有整数的立方和。
6. 编写一个程序，计算并输出1到100之间所有整数的阶乘和。
7. 编写一个程序，计算并输出1到100之间所有整数的倒数和。
8. 编写一个程序，计算并输出1到100之间所有整数的平均值。
9. 编写一个程序，计算并输出1到100之间所有整数的中位数。
10. 编写一个程序，计算并输出1到100之间所有整数的标准差。
11. 编写一个程序，计算并输出1到100之间所有整数的方差。
12. 编写一个程序，计算并输出1到100之间所有整数的最大值。
13. 编写一个程序，计算并输出1到100之间所有整数的最小值。
15. 编写一个程序，计算并输出1到100之间所有整数的平方根和。
16. 编写一个程序，计算并输出1到100之间所有整数的立方根和。
17. 编写一个程序，计算并输出1到100之间所有整数的对数和。
18. 编写一个程序，计算并输出1到100之间所有整数的指数和。
19. 编写一个程序，计算并输出1到100之间所有整数的正弦和。
20. 编写一个程序，计算并输出1到100之间所有整数的余弦和。
21. 编写一个程序，计算并输出1到100之间所有整数的正切和。
22. 编写一个程序，计算并输出1到100之间所有整数的双曲正弦和。
23. 编写一个程序，计算并输出1到100之间所有整数的双曲余弦和。
2012Matlab 仿真考题A
2011-2012（二）MATLAB 仿真期末考试题（A 卷）考试说明：．．．．．1> 请在D 盘目录下建立一个以“学生姓名_班级_A ”为名字的文件夹，并在该文件夹下建立名为“学生姓名_班级_A.doc ”的WORD 文件，例如“张三_自动化083_A ” ，在本次考试中所有结果（包括源程序，图以及运行结果等）均须保存在该文本文档中，否则记零分；2> 答案中要求保留利用MATLAB 软件进行计算和仿真的过程；3> 要求在答卷上写清题号，并按顺序写答案。
一、已知矩阵A ＝⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡876654432， （18分） 求： 1）用MA TLAB 写出矩阵A ，2）获取矩阵的第三行、第三列的元素赋值为10，3）访问矩阵得到第1行所有元素，4）访问矩阵的最后一列元素，5）把矩阵按照列拉长为1行9列矩阵，6）把矩阵扩展为四行四列元素，扩展的元素全部为0。
二、用内部函数生成矩阵： （18分）1）对角阵B ，B 的主对角线元素为[-3 -2，-1，0，1，2 3]；2）4行4列对角阵C ，C 的主对角线元素均为4；3）2行3列的零矩阵D ；4）3行3列的全9阵E ；5）3行3列的伪随机阵F ；6）计算矩阵F 的行列式G 。
三、请在同一个坐标系内画出以下图形： （20分）1、y1=3*cos(x)；2、y2=sin(2x).*cos(x)；3、y3=x 2要求y1为蓝色、实线；y2为黑色、虚线；y3为红色、点线。
并对该图形定义坐标轴范围，[5,5]x ππ∈-，[7,7]y ∈-，x 的步长取0.02。
四、编写函数型M 文件，实现∑==N i i s 1，当和小于等于1000时，输出∑==Ni i s 1，当和>=1000时，s=1000，要求：1）编写程序，定义函数名称为mysum ，程序写在word 文档上，并存入D 盘\学生姓名_班级_A 文件夹下。
2）在工作空间中进行调用，当N=20，N=50分别调用两次，给出调用的语句和执行结果。
Matlab控制系统仿真试卷(横版)2013(精品文档)_共9页
6、已知某系统的实测输入数据 u 和输出数据 y，则用最小二乘法辨识出该系统模型
的命令为>>
7、已知某系统的传递函数为 G(s)
2s 1 ，则在 Matlab 空间中建立该系统
s3 5s2 3s 1
的数学模型 G=tf([2 1],[1 5 3 1])后，能够准确判断该系统稳定性的命令为>>
8、已知 z=5x+y2,在 matlab 中输入[x,y]=meshgrid(-2:0.1:2,-3:0.1:3)后，能够绘制该三
4
y(t)
2 0.3
0.5 0.3
0 0.2
0.8 1 x(t)
燕山大学试卷
4, x D 3、 用 Matlab 语言编写函数，实现如下分段函数： y 4 / Dx, x D （6 分）。
4, x D
密封线
共8 页 第 5 页
4、 考虑如下图所示的典型反馈控制系统的框图
（1）用 Matlab 语言绘制该系统的单位脉冲响应曲线以及在 u(t)=tsin(t)+t2 作用下的响应曲线；（8 分）
（2）绘制该系统的根轨迹、伯德图，并求取相角裕度和幅值裕度。（6 分）
5、
考虑简单的线性微分方程
y(3)
5&y&
4y&
2y
5e t
sin(2t
)
，且方程的
4
初值为 y(0) 1, y&(0) &y&(0) 0.4 。
（1） 试用 Simulink 搭建起系统的仿真模型，并将结果分别输出到示波器
考试时间 2012-
题号 一 二 三 四 分数
总分
燕山大学试卷
MATLAB考试试题及答案
MATLAB考试试题及答案一、选择题（每题5分，共25分）1. 在MATLAB中，下列哪个命令用于创建一个行向量？A. v = [1; 2; 3]B. v = [1 2 3]C. v = [1, 2, 3]D. v = (1, 2, 3)答案：B2. 在MATLAB中，下列哪个命令用于计算矩阵A的行列式？A. det(A)B. det(A')C. det(inv(A))D. det(A^2)答案：A3. 在MATLAB中，下列哪个命令用于计算矩阵A的逆？A. inv(A)B. A^(-1)C. pinv(A)D. A\B答案：A4. 在MATLAB中，下列哪个命令用于求解线性方程组Ax= b？A. A\bB. A/BC. B/AD. A^-1b答案：A5. 在MATLAB中，下列哪个命令用于绘制二维图形？A. plot(x, y)B. scatter(x, y)C. bar(x, y)D. pie(x, y)答案：A二、填空题（每题5分，共25分）6. 在MATLAB中，可以使用______命令创建一个等差数列。
答案：linspace7. 在MATLAB中，可以使用______命令创建一个等比数列。
答案：logspace8. 在MATLAB中，可以使用______命令计算矩阵A的特征值。
答案：eig(A)9. 在MATLAB中，可以使用______命令计算矩阵A的特征向量。
答案：eigenvector(A)10. 在MATLAB中，可以使用______命令计算矩阵A的奇异值。
答案：svd(A)三、解答题（每题25分，共75分）11. 编写MATLAB程序，求解以下线性方程组：2x + 3y - z = 1x - y + 2z = 03x + 2y - 4z = -3答案：```A = [2 3 -1; 1 -1 2; 3 2 -4];b = [1; 0; -3];x = A\b;disp('解为：');disp(x);```12. 编写MATLAB程序，绘制以下函数的图形：y = sin(x) + cos(x)，x ∈ [0, 2π]答案：```x = linspace(0, 2pi, 100);y = sin(x) + cos(x);plot(x, y);title('y = sin(x) + cos(x)');xlabel('x');ylabel('y');grid on;```13. 编写MATLAB程序，计算以下矩阵的特征值和特征向量：A = [1 2 3; 4 5 6; 7 8 9]答案：```A = [1 2 3; 4 5 6; 7 8 9];[V, D] = eig(A);disp('特征值：');disp(diag(D));disp('特征向量：');disp(V);```14. 编写MATLAB程序，使用牛顿迭代法求解方程f(x) = x^3 - 4x + 2 = 0在x = 1附近的根。
MATLAB仿真课后知识题
MATLAB仿真课后知识题第⼀章习题3.请指出以下的变量名（函数名、M⽂件名）中，哪些是合法的？Abc 2004x lil-1 wu_2004 a&b qst.u _xyz 解：合法的变量名有：Abc wu_20044．指令窗操作（1）求[12+2×（7-4）]÷32的运算结果解：>> [12+2*(7-4)]/3^2ans =2（2）输⼊矩阵A=[1，2，3；4，5，6；7，8，9]，观察输出。
解：>> A=[1,2,3;4,5,6;7,8,9]A =1 2 34 5 67 8 9（3）输⼊以下指令，观察运算结果；clear;x=-8:0.5:8;y=x';X=ones(size(y))*x;Y=y*ones(size(x));R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;mesh(X,Y,Z);colormap(hot)xlabel('x'),ylabel('y'),zlabel('z')解：7．指令⾏编辑（1）依次键⼊以下字符并运⾏：y1=2*sin(0.3*pi)/(1+sqrt(5)) 解：>>y1=2*sin(0.3*pi)/(1+sqrt(5))y1 =0.5000(2)通过反复按键盘的箭头键，实现指令回调和编辑，进⾏新的计算；y2=2*cos(0.3*pi)/(1+sqrt(5))解：>>y2=2*cos(0.3*pi)/(1+sqrt(5))y2 =0.363311.编写题4中（3）的M脚本⽂件，并运⾏之。
解：第⼆章习题1.在指令窗中键⼊x=1:0.2:2和y=2:0.2:1,观察所⽣成的数组。
解:>> x=1:0.2:2 x =1.0000 1.2000 1.4000 1.6000 1.80002.0000 >> y=2:0.2:1 y =Empty matrix: 1-by-02．要求在[0，2π]上产⽣50个等距采样数据的⼀维数组，试⽤两种不同的指令实现。
《MATLAB仿真技术》课程考试试题
《MATLAB仿真技术》课程考试试题年学期班级时量，皿分钟总分，皿分考试形式，团签一、填空题(每空2分，共10空，共20分)1、MATLAB常用操作界面包括、工作空间窗口...内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。
2、MATLAB 中Inf 或inf 表示、NaN 或nan 表示、nargout 表示。
3、MATLAB中逗号主要用作；用作输入量与输入量之间的分隔符；。
4、工作空间浏览器主要用于内存变量的、和。
5、MATLAB实现将全下标转换为单下标的指令为。
二、简答题(每题5分，共4题，共20分)1、简述MATLAB历史指令窗的主要作用。
2、简述空数组的功用。
3、简述MATLAB函数的基本结构。
4、简述绘制二维图形的一般步骤。
三、阅读程序并回答问题(每题4分，共28分)1、写出下列指令运行结果。
A=ZerOS(2,4);A(:)=l:8;s=[2 3 5];A(s)Sa=[10 20 30]'A(s)=Sa2、写出下列指令运行结果。
A=reshape(l: 16,2,8)reshape(A,4,4)s=[l 3 6 8 9 11 14 16];A(S)=O3、写出下列指令运行结果。
A=[l,2;3,4];B=[-l,-2;2;l];S=3;A.*BA*BS*AS*B4、下面的函数主要完成什么功能？function f=factor(n)if n<= 1f=l；elsef=factor(n-l)*n;end5、写出下列指令运行结果。
ch=t ABc 123d4e56Fg9,; subch=ch( 1:5)revch=ch(end:-l:l)k=find(ch>=t a,⅛ch<=t z');ch(k)=ch(k)-(*a,-4A*);char(ch)6、写出下列指令运行结果。
A( 1,1)= {'this is cell')；A{l,2}={[l 23;456]};A{2,l}=∣l+2*i];A{2,2}=A{ 1,2]{1 }+(A(l,2}{l }(1,1)+A{ 1,2}{1 }(2,2));celldisp(A)7、下面的程序完成功能是什么？t=0:pi/50:4*pi;yθ=exp(-t∕3);y=exp(-t∕3).*sin(3*t);plot(t,y,'-r',t,yOΛb∖t,-yOΛb')xlabel(ς∖bf∖it t,);ylabel(*∖bf∖it y,);grid on;四、编程题(32分)1、在同•图上分别用红色实线和绿色虚线绘制yl=sin(x)和y2=cos(x)在区间[0, 4*pi]的曲线，并用星号*标出两条曲线的交点以及建立图例。