011100
9.已知A=[705;028]，B=[140;504]，则A&B=100，A|B=111。
001111
二、简答题
1.MATLAB软件系统主要有哪些部分组成？
答：
MATLAB系统主要由开发环境、数学函数库、MATLAB语言、MATLAB图形处理系统及外部应用程序接口组成。
2.启动MATLAB后，默认用户界面中包含哪些窗口元素？
答：
启动MATLAB程序后，进入MATLAB工作界面，包括用于管理文件的工具条（图形用户界面）、命令窗口（CommandWindow）、当前文件夹窗口（CurrentFolder）、工作空间（Workspace）及历史命令窗口（CommandHistory）等。
3.用哪些方法可以启动MATLAB的帮助系统？
使用help命令和lookfor命令有什么区别？
答：
点击桌面工具栏上的帮助按钮
，或在命令窗口中键入helpbrowser命令，或使用任何工具栏上的help菜单，都可打开帮助浏览器。
help：
在命令窗口中键入“help函数名”，将显示指定函数的M文件说明和语法，对于MDL文件，将显示模型的描述；
lookfor：
在命令窗口中键入“lookfor关键字”，将按关键字搜索相关的M文件。
4.如果文件保存在D:
\MyDocuments\MyWork中，而当前文件夹路径为C:
\MATLAB，如何操作才能执行该文件？
答：
选择菜单中的〖File〗→〖SetPath〗菜单项将D:
\MyDocuments\MyWork目录加入到MATLAB的搜索路径中，就可执行文件。
5.什么是M脚本文件？
什么是M函数文件？
两者之间有什么区别？
答：
脚本文件也称为命令文件，是由一系列MATLAB的命令行组成，既没有参数的输入也不返回输出参数。
MATLAB软件中提供了很多工具箱，工具箱中的每一个函数其实就是一个M函数文件。
M函数文件可以接受参数输入，并能返回输出参数。
脚本文件与函数文件的主要区别在于函数文件允许通过数值进行参数传递，函数文件使用局部变量而不是全局变量，函数文件的第一行必须是“函数声明行”，以关键词“function”进行声明，M函数文件名必须与function声明的名称相同。
三、编程题
1.创建三变量A、B、C，其值分别为eye（2,3），ones（1,2），magic（3），其中变量A、C显示，变量B不显示。
解：
clear
A=eye（2,3）,
B=ones（1,2）;
C=magic（3）
2.创建初始值为4，终止值为15，步长为4的向量D，元素个数为3的向量E，并对向量进行加、减、点积和叉积运算。
解：
clear
D=[4:
4:
15];
E=linspace（4,15,3）;
D+E
D-E
dot（D,E）
cross（D,E）
3.已知3×3维矩阵
，做如下操作：
（1）寻访F矩阵的第8个元素；
（2）将F矩阵的第二行重新赋值为[583]；
（3）分别生成F矩阵的对角矩阵和上三角矩阵；
（4）抽取F矩阵的第二行、第三行，第二列、第三列生成新的矩阵G。
解：
clear
F=[257;324;159];
F（8）
F（2,:
）=[583]
diag（F）
triu（F）
G=F（2:
3,2:
3）
4.已知矩阵
和
，查找出K大于L的元素，K不等于L的元素，及K|L、K&L的值。
解：
clear
K=[3012;1024;0-101];
L=[2013;0404;6001];
K>=L
K~=L
K|L
K&L
5.分别用for和while循环语句编写程序，求
的值。
解：
for语句：
clear
sum=0;
fori=1:
40
sum=sum+i^0.5;
end
y=sum
while语句：
clear
sum=0;i=1;n=40;
whilei<=n
sum=sum+i^0.5;
i=i+1;
end
y=sum
姓名：
朱剑
学号：
1111001261637
得分：
教师签名：
控制系统CAD作业2
第3章
一、填空题
1.在MATLAB的命令窗口中键入simulink命令可打开Simulink图形化设计环境。
2.引出信号支线的方法是按住鼠标右键，在需要增加分支的地方拉出引线；也可在按下Ctrl键的同时，在分支点按住鼠标左键拉出引线。
3.已知系统传递函数的增益为4.5，零点为-2±j，极点为0、-3和-1±j，则零极点模型模块的参数设置对话框内“Zeros”框应填入[-2+j-2-j]；“Poles”框应填入[0-3-1-j-1+j]；“Gain”框应填入4.5。
4.用零极点模型表示系统结构时若没有零点只有极点，需要在零极点模型模块的参数设置对话框将“Zeros”向量设为[]。
5.已知单位负反馈系统的开环传递函数为：
，若要完成以下的系统仿真分析：
（1）利用Simulink求系统的单位阶跃响应；
（2）用示波器显示仿真结果；（3）将结果输出到MATLAB的工作空间，则需要用到的模块有Sources库的Step模块、Continuous库的TransferFcn模块、MathOperations库的Sum模块、Sinks库的ToWorkspace模块和Scope模块。
6.利用Continuous库的TransferFcn模块表示传递函数
时，其参数设置对话框中“Numertorcoefficient（分子多项式）”框应填入[25]，“Denominatorcoefficient（分母多项式）”框应填入[120]。
7.利用已存在系统模型的部分结构建立子系统，则Simulink会自动添加输入模块和输出模块来表示子系统的输入端口和输出端口。
8.创建子系统、封装子系统、编辑已封装的子系统所用的菜单命令分别为〖CreateSubsystem〗、〖MaskSubsystem…〗、〖EditMask〗。
9.在仿真参数设置对话框中可选择求解器（Solver）类型，可供选择的类型有可变步长模式和固定步长模式。
10.选择可变步长模式时仿真步长随模型的动态特性变化，提供误差控制和过零检测两种功能，而固定步长模式则不提供。
二、选择题
1.已知单位负反馈系统的开环传递函数为：
，利用Simulink求系统的单位阶跃响应时，不需要下列模块中的哪一种？
（D）
A.
B.
C.
D.
2.能将Simulink中的结果数据导出到MATLAB工作空间的模块是：
（C）
A.
B.
C.
D.
三、判断题
1.Simulink将建模时所需的功能模块按其功能划分到不同的模块库中，因此任一模块都有唯一对应的模块库。
（错）
2.Simulink对子系统的封装功能不仅增强了模型的可读性，还大大简化了子系统的参数设置，在仿真前只需打开子系统编辑窗口，对其中的每个模块分别进行参数设置就可以了。
（错）
3.利用Simulink进行系统建模和仿真的一般步骤是调用各模块库中的模块建立仿真模型，设置各模块参数及全局仿真参数后运行仿真，并根据系统性能要求调整参数，输出结果。
（对）
第4章
一、填空题
1.在线性定常系统中，常用的数学模型有：
微分方程及传递函数的多项式模型、传递函数的零极点增益模型、状态空间模型和频率响应数据模型或者TF、ZPK、SS、FRD等。
（、、、）
2.说明以下函数的功能：
（1）tf2zp（）：
传递函数多项式模型转换为零极点模型
（2）zp2tf（）：
零极点模型转换为传递函数多项式模型
（3）tf2ss（）：
传递函数多项式模型转换为状态空间模型
（4）ss2tf（）：
状态空间模型转换为传递函数多项式模型
（5）zp2ss（）：
零极点模型转换为状态空间模型
（6）ss2zp（）：
状态空间模型转换为零极点模型
二、编程题
1.已知系统的传递函数为：
，在MATLAB中用两种方法建立系统的传递函数多项式模型。
解：
方法一
clear
num=[741];%分子多项式系数向量
den=[1050308];%分母多项式系数向量
G=tf（num,den）
方法二
clear
s=tf（'s'）;%将变量s定义为传递函数模型
G1=（7*s^2+4*s+1）/（s^6+5*s^4+3*s^2+8）
2.已知系统的传递函数为：
，在MATLAB中建立系统的传递函数多项式模型。
解：
clear
num=conv（conv（[32],[13]）,[241]）;
den=conv（conv（conv（[103],[20351]）,conv（[11],[11]））,[11]）;
G2=tf（num,den）
3.已知系统的传递函数为：
，在MATLAB中用两种方法建立系统的零极点增益模型。
解：
方法一
clear
z=[1];p=[-2,3,-1];k=2;%输入系统的零级、极点和增益
G3=zpk（z,p,k）
方法二
clear
s=zpk（'s'）;%将变量s定义为零级点增益模型
G4=2*（s-1）/（（s+2）*（s^2-2*s-3））
4.已知系统的传递函数为：
，在MATLAB中用两种方法将TF模型转换为ZPK模型。
解：
方法一
clear
num=[213];%分子多项式系数向量
den=[12431];%分母多项式系数向量
G=tf（num,den）;%创建传递函数模型
G5=zpk（G）
方法二
clear
num=[213];%分子多项式系数向量
den=[12431];%分母多项式系数向量
[z,p,k]=tf2zp（num,den）
G6=zpk（z,p,k）
5.已知系统的传递函数为：
，在MATLAB中用两种方法将ZPK模型转换为TF模型。
解：
方法一
clear
z=[-2,-2,+2];p=[1,1,-1,-5,-6^0.5,6^0.5];k=1;%输入系统的零点、极点和增益
G=zpk（z,p,k）;%生成ZPK模型
G1=tf（G）
方法二
clear
z=[-2,-2,+2]';p=[1,1,-1,-5,-6^0.5,6^0.5];k=1;%输入系统的零点、极点和增益
[num,den]=zp2tf（z,p,k）%模型转换，获得系统TF模型的分子分母系数
G=tf（num,den）%根据所求的参数生成TF模型
6.已知两个子系统的传递函数分别为：
，
，在MATLAB中用将两个子系统串联连接生成新系统。
解：
clear%清除工作空间中变量
z=[2];p=[-2i,2i,-2,-2,-1,3];k=2;%输入子系统的零点、极点和增益
G1=zpk（z,p,k）;%生成第一个子系统的ZPK模型
num=[10231];%子系统的分子多项式系数向量
den=[153241];%子系统分母多项式系数向量
G2=tf（num,den）;%生成第二个子系统的TF模型
G=series（G1,G2）
7.已知两个子系统的传递函数分别为：
，
。
在MATLAB中用将两个子系统并联连接生成新系统。
解：
clear%清除工作空间中变量
z=[1];p=[-1,-1,-1,2i,-2i];k=3;%输入子系统的零点、极点和增益
G1=zpk（z,p,k）;%生成第一个子系统的ZPK模型
num=[542];%子系统的分子多项式系数向量
den=[10421];%子系统分母多项式系数向量
G2=tf（num,den）;%生成第二个子系统的TF模型
G=parallel（G1,G2）
8.已知两个子系统的传递函数分别为:
，
。
在MATLAB中用将两个子系统正、负反馈连接生成新系统。
解：
clear%清除工作空间中变量
z=[-21];p=[314];k=1;%输入子系统的零点、极点和增益
G1=zpk（z,p,k）;%生成第一个子系统的ZPK模型
num=[11];%子系统的分子多项式系数向量
den=[112];%子系统分母多项式系数向量
G2=tf（num,den）;%生成第二个子系统的TF模型
Ga=feedback（G1,G2,1）
Gb=feedback（G1,G2）
姓名：
朱剑
学号：
1111001261637
得分：
教师签名：
控制系统CAD作业3
第5章
一、填空题
1.在控制系统时域分析时，常用的时域动态性能指标有：
上升时间tr、峰值时间tm、
调节时间ts和超调量σ%。
2.MATLAB中step（）函数是单位阶跃响应函数。
它可以用来绘制单位阶跃响应波形图或求取输出量的数值向量。
3.impulse函数不带输出变量时，可在当前窗口中直接绘制出系统的单位冲激响应曲线。
4.在MATLAB中，欲求取系统的时域性能，可以利用鼠标和菜单的操作，在单位阶跃响应曲线上求取，也可以编制简单程序求取。
5.要判别系统的稳定性，可使用的MATLAB函数有roots（）、zpkdata（）、tf2zp（）和
pzmap（）四个函数。
若要直接求出系统闭环特征方程的根，应使用roots（）函数。
6.若在Simulink仿真集成环境中分析系统在斜坡信号和阶跃信号作用下系统的稳态误差，在Simulink建模时，可分别选用Step和Ramp模块作为系统的输入信号模块。
二、写出下列程序（命令）的功能或执行结果
1.
>>G=zpk（[],[-1+3*I,-1-3*i],3）
>>step（G）
绘制出
的单位阶跃响应曲线
2.
num=7*[1,5];den=conv（[1,0,0],
conv（[1,10],[1,1]））;
[gm,pm,wg,wc]=margin（num,den）
求系统
的幅值裕度和相角裕度
3.
num=5*[1,5,6];den=[1,6,10,8];
impulse（num,den）;gridon;
求三阶系统的
单位阶跃响应和单位冲激响应。
三、编程题
1.典型二阶系统传递函数为：
，阻尼比ζ的值为0.4，试编制程序绘制出自然振荡角频率ωn从1变化到10（每次递增2）的单位阶跃响应曲线，并对系统的动态性能的变化作定性分析。
clc;clf;clear
wn=[1,3,5,7,9];%自然振荡角频率
dr=0.4;%阻尼比
t=linspace（0,20,200）;%将时间t在0-20间200等分
holdon;%打开图形叠加功能（同一个坐标上画多条曲线）
forj=1:
5
num=wn（j）^2;
den=[1,2*dr*wn（j）,0];%开环传递函数分母多项式系数
Gk=tf（num,den）;%求得开环传递函数
Gb=feedback（Gk,1）;%求得闭环传递函数
step（Gb）;grid;%画图
gtext（strcat（'wn=',num2str（j）））%添加鼠标文字注释，strcat（）为字符串连接函数
end
随着自然振荡角频率ωn的增大，系统的响应速度加快。
2.设单位反馈系统的开环传递函数：
，试求系统的单位阶跃响应的性能指标σ％、tS（2％误差带）。
clc;clear
num=49;den=[110.50];
Gk=tf（num,den）%求得开环传递函数
Gb=feedback（Gk,1,-1）%求得单位反馈的闭环传递函数
[y,t]=step（Gb）;%求取单位阶跃响应
%计算调节时间
Cend=dcgain（Gb）;%求取系统的终值
n=length（t）;%求得仿真时间t序列的长度
while（y（n）>0.98*Cend）&（y（n）<1.02*Cend）
%通过循环求取单位阶跃响应首次进入±2%h（∞）内的所需时间序列号
n=n-1;
end
Ts=t（n）%求得调节时间
%计算超调量
[Y,i]=max（y）;%求得系统响应的峰值
Os=100*（Y-Cend）/Cend%求得超调量
Ts=0.8074（s）
Os=2.8335（%）
3.单位负反馈系统的开环传递函数为：
，试求在单位斜坡输入信号作用下的稳态误差，试写出MATLAB命令序列（需首先判断系统的稳定性）。
>>num=[77];
>>den=[conv（conv（[10],[13]）,[145]）];
>>Gk=tf（num,den）;
>>Gb=feedback（Gk,1,-1）
Transferfunction:
7s+7
-------------------------------
s^4+7s^3+17s^2+22s+7
>>r=roots（[1717227]）
r=
-4.0287
-1.2632+1.5198i
-1.2632-1.5198i
-0.4449
闭环特征根的实部均为负，系统稳定。
>>rss=tf（[10],[01]）;
>>Kv=dcgain（rss*Gk）
>>essv=1/Kv
essv=2.1429
第6章
一、填空题
1.所谓根轨迹，是指控制系统开环传递函数的某一参数（如开环增益）从零变化到
无穷大时，闭环特征根在s平面上移动的轨迹。
2.根轨迹法的实质就是依据反馈系统中开、闭环传递函数间的关系，由开环传递函数的零、极点直接寻求闭环根轨迹的总体规律。
3.根据给定的根轨迹增益向量k，绘制系统sys的部分根轨迹的MATLAB函数是rlocus（sys,k）。
4.计算与根轨迹上任一闭环极点相对应的根轨迹增益的MATLAB函数是rlocfind（）。
二、编程分析题
1.已知单位反馈系统的开环传递函数，试用MALAB编程绘出系统根轨迹。
⑴
（2）
（1）num=1;den=[conv（conv（[10],[0.21]）,[0.51]）]
sys=tf（num,den）;
rlocus（sys）;
（2）num=[12];den=[conv（[10],[21]）]
sys=tf（num,den）;
rlocus（sys）;
2.已知控制系统的开环传递函数为
，试用MALAB分析该系统的根轨迹的起点和终点。
num=[1104];
den=[1370];
G=tf（num,den）;
rlocus（G）；
p=roots（den）
升级会员 