MATLAB控制系统仿真作业1
一、 控制系统的模型与转换
1． 请将下面的传递函数模型输入到matlab 环境。
]52)1)[(2(24)(322
33++++++=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)); G
Transfer 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.99
2． 请将下面的零极点模型输入到matlab 环境。请求出上述模型的零极点，并绘制其位置。
)1)(6)(5()1)(1(8)(22
+++-+++=s s s s j s j s s G )
2.8()
6.2)(2.3()(1
511-++=----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)
Zero/pole/gain:
z^5 (z+0.3125) (z+0.3846) ------------------------- (z-0.122)
Sampling time: 0.05
>>pzmap （H ）
二、 线性系统分析
1. 请分析下面传递函数模型的稳定性。
2
21)(23+++=
s s s s G 1
3)50600300(1
3)(22+++++=
s s s s s s G
>> num=[1];
den=[1 2 1 2]; G=tf(num,den); eig(G)' ans =
-2.0000
0.0000 - 1.0000i 0.0000 + 1.0000i
可见，系统有两个特征根在虚轴上，一个特征根在虚轴左侧，所以系统是临界稳定的。
G=tf(num,den); eig(G)' ans =
-1.9152 -0.1414
0.0283 - 0.1073i 0.0283 + 0.1073i
可见，有两个特征根在虚轴右侧，所以系统是不稳定的。
2. 请判定下面离散系统的稳定性。
)
05.025.02.0(2
3)(23+--+-=
z z z z z H
)340
39.804.10215.20368.791
.1576.1112.2)(1
2345
12-++--++=-------z z z z z z z z H
>> num=[-3 2];
den=[1 -0.2 -0.25 0.05]; H=tf(num,den,'Ts',0.1); [eig(H) abs(eig(H))] ans =
-0.5000 0.5000 0.5000 0.5000
0.2000 0.2000
可以看出，由于各个特征根的模均小于1，所以可以判定闭环系统是稳定的。
[eig(H) abs(eig(H))] ans =
0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 4.1724 4.1724 0.3755 + 0.1814i 0.4170 0.3755 - 0.1814i 0.4170 -0.5292 0.5292 -0.2716 0.2716
0.1193 0.1193
可以看出，由于4.1724这个特征根的模大于1，所以可以判定闭环系统是不稳定的。
3. 设描述系统的传递函数为
40320
10958411812467284224494536546364032018576022208812266436380598251418)(2345678
234567+++++++++++++++=s s s s s s s s s s s s s s s s G ，假定系统
具有零初始状态，请求出单位阶跃响应曲线和单位脉冲响应曲线。
>> num=[18 514 5982 36380 122664 22088 185760 40320]; den=[1 36 546 4536 22449 67284 118124 109584 40320]; G=tf(num,den)
Transfer function:
18 s^7 + 514 s^6 + 5982 s^5 + 36380 s^4 + 122664 s^3 + 22088 s^2 + 185760 s + 40320 -----------------------------------------------------------------------------------------
s^8 + 36 s^7 + 546 s^6 + 4536 s^5 + 22449 s^4 + 67284 s^3 + 118124 s^2 + 109584 s + 40320
>> step(G ,10) >> impulse(G ,10)
单位脉冲响应：
1. 请分析下面传递函数模型阶跃响应。
2
21)(2
3
+++=
s s s s G
利用Simulink 建模，建立系统仿真模型如下：
单击启动仿真按钮，双击示波器得到系统的阶跃响应如下：
2. 请分析下面离散系统的脉冲响应。
13)50600300(1
43)(2
223++++++=s s s s s s s G
利用Simulink 建模，建立系统仿真模型如下：
单击启动仿真按钮，双击示波器得到系统的脉冲响应如下：
3. 对离散采样系统进行分析，并求出其阶跃响应。
其中：2
373
)(2
3++++=
s s s s s G
利用Simulink 建模，建立系统仿真模型如下：
单击启动仿真按钮，双击示波器得到系统的阶跃响应如下：
4. 设计控制器，使得下列系统稳定。
)2.1)(2)(3()3.2)(1()(+-+++=
s s s s s s G
利用Simulink 建模，未连入控制器时，仿真模型和响应如下：
利用Simulink 建模，设计控制器：
从响应输出图形可以看出，连入控制器后系统稳定，性能明显提高。
四、 基于MATLAB 的PID 控制器设计
设计题目：
1. 应用Ziegler — Nichols 算法设计PID 控制器，实现系统的闭环稳定，并比较对各个系统的控制效果。
3
)1s (1G(s)+=
未连入PID 控制器时的系统仿真及其性能指标如下：
可见，未调节时的系统性能有待提高，需设计PID控制器连入。
输入：
>> num=1;
den=conv( [1,1],conv( [1,1],[1,1] ));
Step(num,den);
K=dcgain (num,den)
得出：
K =1
根据图形，得出：
L=1.86 T=4.4
利用自定义的Ziegler_std函数求出Kp、Ti、Td
输入：
>> K=1;
L=1.86;
T=4.4;
[num,den,Kp,Ti,Td]=Ziegler_std (3,[K,L,T])
得出：
num =
2.6400 2.8387 1.5262
den =
根据得出的Kp 、Ti 、Td 值，设计PID 控制器，并利用利用Simulink 仿真建模。 仿真模型及其响应如下：
可见，加入PID 控制器调节后，系统性能明显改善。
5
)1s (1G(s)+=
未连入PID 控制器时的系统仿真及其性能指标如下：
输入：
>> num=1;
den=conv([1,1],conv([1,1],…;
conv( [1,1],conv( [1,1],[1,1] ))));
Step(num,den);
K=dcgain (num,den)
得出：
K = 1
根据图形，得出：
L=3.4 T=6.8
利用自定义的Ziegler_std函数求出Kp、Ti、Td
输入：
>> K=1;
L=3.4;
T=6.8;
[num,den,Kp,Ti,Td]=Ziegler_std (3,[K,L,T])
得出:
num =
4.0800 2.4000 0.7059
den =
1 0
Kp =
2.4000
Ti =
6.8000
Td =
1.7000
根据得出的Kp、Ti、Td值，设计PID控制器，并利用利用Simulink仿真建模。仿真模型及其响应如下：
可见，加入PID控制器调节后，系统性能明显改善。
3)1 s(1
-1.5s G(s)
++
=
利用Simulink建模，未连入控制器时，仿真模型和响应如下：
输入:
>>num=[-1.5 1];
den=conv( [1,1],conv( [1,1],[1,1] ));
Step(num,den);
K=dcgain (num,den)
得出:
K =1
根据图形，得出：
L=1.8 T=5.7
利用自定义的Ziegler_std函数求出Kp、
Ti、Td
输入:
>> K=1;
L=1.8
T=5.7;
[num,den,Kp,Ti,Td]=Ziegler_std (3,[K,L,T])
得出:
num =
3.4200 3.8000 2.1111
den =
10
Kp =
3.8000
Ti =
3.6000
Td =
0.9000
根据得出的Kp、Ti、Td值，设计PID控制器，并利用利用Simulink仿真建模。仿真模型及其响应如下：
可见，加入PID控制器调节后，系统性能明显改善。
五、模糊控制器设计
设计任务：试设计一个模糊控制器，实现对室内温度的控制的模拟。
参考输入：
（1）温度18-40℃范围内分为七个论域，NB NM NS ZE PS PM PB；隶属度函数满足高斯分布；
（2）温度变化率-2 ~2℃范围内分为七个论域，NB NM NS ZE PS PM PB；隶属度函数满足高斯分布；
参考输出：
变频空调输出的控制信号。在一定范围内分为七个论域，NB NM NS ZE PS PM PB，隶属度函数为常数1。
模糊推理过程，output=输入隶属度函数值*输出论域的中心值。
注：本模糊程序采用PAM控制方式的压缩机，则其输出的转速范围为：0~10500转/分。
控制规则：
% % if input is NB and errorinput is PS, then output is NM; % % if input is NB and errorinput is PM, then output is NM; % % if input is NB and errorinput is PB, then output is NS; % % if input is NM and errorinput is NB, then output is NB; % % if input is NM and errorinput is NM, then output is NM; % % if input is NM and errorinput is NS, then output is NM; % % if input is NM and errorinput is ZE, then output is NM; % % if input is NM and errorinput is PS, then output is NM; % % if input is NM and errorinput is PM, then output is NS; % % if input is NM and errorinput is PB, then output is NS; % % if input is NS and errorinput is NB, then output is NM; % % if input is NS and errorinput is NM, then output is NS; % % if input is NS and errorinput is NS, then output is NS; % % if input is NS and errorinput is ZE, then output is NS; % % if input is NS and errorinput is PS, then output is NS; % % if input is NS and errorinput is PM, then output is ZE; % % if input is NS and errorinput is PB, then output is ZE; % % if input is ZE and errorinput is NB, then output is NS; % % if input is ZE and errorinput is NM, then output is ZE; % % if input is ZE and errorinput is NS, then output is ZE; % % if input is ZE and errorinput is ZE, then output is ZE; % % if input is ZE and errorinput is PS, then output is ZE;
% % if input is ZE and errorinput is PM, then output is PS; % % if input is ZE and errorinput is PB, then output is PS;
% % if input is PS and errorinput is NB, then output is ZE; % % if input is PS and errorinput is NM, then output is PS; % % if input is PS and errorinput is NS, then output is PS;
% % if input is PS and errorinput is ZE, then output is PS;
% % if input is PS and errorinput is PS, then output is PS;
% % if input is PS and errorinput is PM, then output is PM; % % if input is PS and errorinput is PB, then output is PM; % % if input is PM and errorinput is NB, then output is PS; % % if input is PM and errorinput is NM, then output is PS; % % if input is PM and errorinput is NS, then output is PM; % % if input is PM and errorinput is ZE, then output is PM; % % if input is PM and errorinput is PS, then output is PM; % % if input is PM and errorinput is PM, then output is PM; % % if input is PM and errorinput is PB, then output is PB; % % if input is PB and errorinput is NB, then output is PS; % % if input is PB and errorinput is NM, then output is PM; % % if input is PB and errorinput is NS, then output is PM; % % if input is PB and errorinput is ZE, then output is PM; % % if input is PB and errorinput is PS, then output is PB;
>> x1 = (18:0.1:40)';
y0 = gaussmf(x1, [1 18]);
y1 = gaussmf(x1, [1 21]);
y2 = gaussmf(x1, [1 25]);
y3 = gaussmf(x1, [1 29]);
y4 = gaussmf(x1, [1 33]);
y5 = gaussmf(x1, [1 37]);
y6 = gaussmf(x1, [1 40]);
plot(x1,[y0 y1 y2 y3 y4 y5 y6])
2.误差图：
程序为：
>> x1 = (-2:0.1:2)';
y0 = gaussmf(x1, [0.5 -2]);
y1 = gaussmf(x1, [0.5 -1.3]); y2 = gaussmf(x1, [0.5 -0.7]); y3 = gaussmf(x1, [0.5 0]);
y4 = gaussmf(x1, [0.5 0.7]);
y5 = gaussmf(x1, [0.5 1.3]);
y6 = gaussmf(x1, [0.5 2]);
plot(x1,[y0 y1 y2 y3 y4 y5 y6])
3.程序为；
x=35;
ex=-0.8;
% define input type in fuzzy zone
y0 = gaussmf(x, [1 18]);
y1 = gaussmf(x, [1 21]);
y2 = gaussmf(x, [1 25]);
y3 = gaussmf(x, [1 29]);
y4 = gaussmf(x, [1 33]);
y5 = gaussmf(x, [1 37]);
y6 = gaussmf(x, [1 40]);
a=[y0 y1 y2 y3 y4 y5 y6];
b=max(a);
% caculate input in fuzzy zone,get input_type and input_authorityvalue if x<=40 & x>=18
if b==a(1)
type='NB';
authorityvalue=y0;
elseif b==a(2)
type='NM';
authorityvalue=y1;
elseif b==a(3)
type='NS';
authorityvalue=y2;
elseif b==a(4)
type='ZE';
authorityvalue=y3;
elseif b==a(5)
type='PS';
authorityvalue=y4;
elseif b==a(6)
type='PM';
authorityvalue=y5;
elseif b==a(7)
type='PB';
authorityvalue=y6;
end
else
if x>40
type='PB';
authorityvalue=1;
elseif x<18
type='NB';
authorityvalue
%error calculate.
ey0 = gaussmf(x, [0.5 -2]);
ey1 = gaussmf(x, [0.5 -1.3]);
ey2 = gaussmf(x, [0.5 -0.7]);
ey3 = gaussmf(x, [0.5 0]);
ey4 = gaussmf(x, [0.5 0.7]);
ey5 = gaussmf(x, [0.5 1.3]);
ey6 = gaussmf(x, [0.5 2]);
a=[ey0 ey1 ey2 ey3 ey4 ey5 ey6];
b=max(a);
% caculate input in fuzzy zone,get input_type and input_authorityvalue if x<=2 & x>=-2
if b==a(1)
etype='NB';
eauthorityvalue=y0;
elseif b==a(2)
etype='NM';
eauthorityvalue=y1;
elseif b==a(3)
etype='NS';
eauthorityvalue=y2;
elseif b==a(4)
etype='ZE';
eauthorityvalue=y3;
elseif b==a(5)
etype='PS';
eauthorityvalue=y4;
elseif b==a(6)
etype='PM';
eauthorityvalue=y5;
elseif b==a(7)
etype='PB';
eauthorityvalue=y6;
end
else
if x>2
etype='PB';
eauthorityvalue=1;
elseif x<-2
etype='NB';
eauthorityvalue=1;
end
if type=='NB'& etype=='NB'
out=authorityvalue*0.001;
% if input is NB and errorinput is NB, then output is NB;
elseif type=='NB'& etype=='NM'
out=authorityvalue*0.001;
% if input is NB and errorinput is NM, then output is NB;
elseif type=='NB'& etype=='NS'
out=authorityvalue*0.001;
% if input is NB and errorinput is NS, then output is NB;
elseif type=='NB'& etype=='ZE'
out=authorityvalue*1750;
% if input is NB and errorinput is ZE, then output is NM;
elseif type=='NB'& etype=='PS'
out=authorityvalue*1750;
% if input is NB and errorinput is PS, then output is NM;
elseif type=='NB'& etype=='PM'
out=authorityvalue*1750;
% if input is NB and errorinput is PM, then output is NM;
elseif type=='NB'& etype=='PB'
out=authorityvalue*3500;
% if input is NB and errorinput is PB, then output is NS;
elseif type=='NM'& etype=='NB'
out=authorityvalue*0.001;
% if input is NM and errorinput is NB, then output is NB;
elseif type=='NM'& etype=='NM'
out=authorityvalue*1750;
% if input is NM and errorinput is NM, then output is NM;
elseif type=='NM'& etype=='NS'
out=authorityvalue*1750;
% if input is NM and errorinput is NS, then output is NM;
elseif type=='NM'& etype=='ZE'
out=authorityvalue*1750;
% if input is NM and errorinput is ZE, then output is NM;
elseif type=='NM'& etype=='PS'
out=authorityvalue*1750;
% if input is NM and errorinput is PS, then output is NM;
elseif type=='NM'& etype=='PM'
out=authorityvalue*3500;
% if input is NM and errorinput is PM, then output is NS;
elseif type=='NM'& etype=='PB'
out=authorityvalue*3500;
% if input is NM and errorinput is PB, then output is NS;
elseif type=='NS'& etype=='NB'
Matlab结构图控制系统仿真
图5. 利用 SIMULINK仿
4. 建立如图11-54所示的仿真模型，其中PID控 制器采用Simulink子系统封装形式，其内部 结构如图11-31（a)所示。试设置正弦波信号 幅值为5、偏差为0、频率为10πHz\始终相位 为0，PID控制器的参数为Kp=10.75、 Ki=1.2、Kd=5，采用变步长的ode23t算法、 仿真时间为2s，对模型进行仿真。 （6)观察仿真结果。系统放着结束后，双击仿真模型中的示波器模块，得到仿真结果。单击示波器窗口工具栏上的Autoscale按钮，可以自动调整坐标来 使波形刚好完整显示，这时的波形如图所示。 图3 2. 题操作步骤如下： (1） 打开一个模型编辑窗口。 (2） 将所需模块添加到模型中。在模块库浏览器中单击Sources,将 Clock（时钟）拖到模型编辑窗口。同样，在User-Defined Functions（用户定义模块库）中把Fcn(函数模块）拖到模型编辑窗口，在Continuous（连续系统模块库）中把 Integrator（积分模块）拖到模型编辑窗口，在Sinks中把Display模块编辑窗口。 (3） 设置模块参数并连接各个模块组成仿真模型。双击Fcn模块，打开Function Block operations中把Add模块拖到模型编辑窗口，在Sinks中把Scope模块拖到模型编辑窗口。 （3） 设置模块参数并连接各个模块组成仿真模型。先双击各个正弦源，打开其Block Parameters对话框，分别设置Frequency（频率）为2*pi、 6*pi、10*pi、 14*pi、18*pi，设置Amplitude（幅值）为1、1/3、1/5、1/7和1/9，其余参数不改变。对于求和模块，將符号列表List of signs设置为 +++++。 （4） 设置系统仿真参数。单击模型
控制系统MATLAB仿真基础
系统仿真 § 4.1控制系统的数学模型 1、传递函数模型（tranfer function） 2、零极点增益模型（zero-pole-gain） 3、状态空间模型（state-space） 4、动态结构图（Simulink结构图） 一、传递函数模型(transfer fcn-----tf) 1、传递函数模型的形式 传函定义：在零初始条件下，系统输出量的拉氏变换C（S）与输入量的拉氏变换R（S）之比。 C(S) b1S m+b2S m-1+…+b m G（S）=----------- =- -------------------------------- R(S) a1S n + a2S n-1 +…+ a n num(S) = ------------ den(S) 2、在MATLAB命令中的输入形式 在MATLAB环境中，可直接用分子分母多项式系数构成的两个向量num、den表示系统： num = [b1, b2, ..., b m]; den = [a1, a2, ..., a n]; 注：1）将系统的分子分母多项式的系数按降幂的方式以向量的形式输入两个变量，中间缺项的用0补齐，不能遗漏。 2）num、den是任意两个变量名，用户可以用其他任意的变量名来输入系数向量。 3）当系统种含有几个传函时，输入MATLAB命令状态下可用n1,d1;n2,d2…….。 4）给变量num，den赋值时用的是方括号；方括号内每个系数分隔开用空格或逗号；num,den方括号间用的是分号。 3、函数命令tf( ) 在MATLAB中，用函数命令tf( )来建立控制系统的传函模型，或者将零极点增益模型、状态空间模型转换为传函模型。 tf( )函数命令的调用格式为： 圆括号中的逗号不能用空格来代替 sys = tf ( num, den ) [G= tf ( num, den )]
Matlab通信系统仿真实验报告
Matlab通信原理仿真 学号： 2142402 姓名：圣斌
实验一Matlab 基本语法与信号系统分析 一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。 二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);
春MATLAB仿真期末大作业
MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：
2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。
仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。
控制系统的MATLAB仿真与设计课后答案
控制系统的MATLAB仿真与设计课后答案
>>z=-4*sqrt(2)*sin(t); >>plot3(x,y,z,'p'); >>title('Line in 3-D Space'); >>text(0,0,0,'origin'); >>xlabel('X'),ylable('Y'),zlable('Z');grid; 4>>theta=0:0.01:2*pi; >>rho=sin(2*theta).*cos(2*theta); >>polar(theta,rho,'k'); 5>>[x,y,z]=sphere(20); >>z1=z; >>z1(:,1:4)=NaN; >>c1=ones(size(z1)); >>surf(3*x,3*y,3*z1,c1); >>hold on >>z2=z; >>c2=2*ones(size(z2)); >>c2(:,1:4)=3*ones(size(c2(:,1:4))); >>surf(1.5*x,1.5*y,1.5*z2,c2); >>colormap([0,1,0;0.5,0,0;1,0,0]); >>grid on >>hold off 第四章 1>>for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 disp(m) end end 2M文件：function[s,p]=fcircle(r) s=pi*r*r; p=2*pi*r; 主程序： [s,p]=fcircle(10) 3>>y=0;n=100; for i=1:n y=y+1/i/i; end >>y
电机学matlab仿真大作业报告
. 基于MATLAB的电机学计算机辅助分析与仿真 实验报告
一、实验内容及目的 1.1 单相变压器的效率和外特性曲线 1.1.1 实验内容 一台单相变压器，N S =2000kVA, kV kV U U N N 11/127/21=，50Hz ，变压器的参数 和损耗为008.0* ) 75(=C k o R ，0725.0*=k X ，kW P 470=，kW P C KN o 160)75(=。 (1)求此变压器带上额定负载、)(8.0cos 2滞后=?时的额定电压调整率和额定效率。 (2)分别求出当0.1,8.0,6.0,4.0,2.0cos 2=?时变压器的效率曲线，并确定最大效率和达到负载效率时的负载电流。 (3)分析不同性质的负载（),(8.0cos 0.1cos ),(8.0cos 222超前，滞后===???）对变压器输出特性的影响。 1.1.2 实验目的 （1）计算此变压器在已知负载下的额定电压调整率和额定效率 （2）了解变压器效率曲线的变化规律 （3）了解负载功率因数对效率曲线的影响 （4）了解变压器电压变化率的变化规律 （5）了解负载性质对电压变化率特性的影响 1.1.3 实验用到的基本知识和理论 （1）标幺值、效率区间、空载损耗、短路损耗等概念 （2）效率和效率特性的知识 （3）电压调整率的相关知识 1.2串励直流电动机的运行特性 1.2.1实验内容 一台16kw 、220V 的串励直流电动机，串励绕组电阻为0.12Ω，电枢总电阻为0.2Ω。电动势常数为.电机的磁化曲线近似的为直线。其中为比例常数。假设电枢电流85A 时，磁路饱和(为比较不同饱和电流对应的效果，饱和电流可以自己改变）。
matlab控制系统仿真.
课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称龙蟠学院 专业自动化 班级M10自动化 学生姓名 学号 课程设计地点 C208 课程设计学时一周 指导教师应明峰 金陵科技学院教务处制成绩
一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。
（d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统PID控制器取参数分别为：50 2 5 有积分作用单回路控制系统PID控制器取参数分别为：50 0 5
大比例作用单回路控制系统PID控制器取参数分别为：50 0 0 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响应，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏；
MATLAB仿真实验全部
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间围0---Tn 。 3、),(T sys step ；表示时间围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：)()()()(1 )(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
matlab 大作业
上海电力学院 通信原理Matlab仿真 实验报告 实验名称： 8QAM误码率仿真 试验日期： 2014年 6月3日 专业：通信工程 姓名：罗侃鸣 班级： 2011112班 学号： 20112272
一、实验要求 写MATLAB程序，对图示的信号星座图完成M=8的QAM通信系统Monte Carlo仿真，在不同SNRindB=0：15时，对N=10000（3比特）个符号进行仿真。画出该QAM系统的符号误码率。 二、实验原理 1 QAM调制原理 QAM(Quadrature Amplitude Modulation)正交幅度调制技术，是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽内的正交性，实现两路并行的数字信息的传输。该调制方式通常有8QAM，16QAM，64QAM。 QAM调制实际上就是幅度调制和相位调制的组合，相位+ 幅度状态定义了一个数字或数字的组合。QAM的优点是具有更大的符号率，从而可获得更高的系统效率。通常由符号率确定占用带宽。因此每个符号的比特（基本信息单位）越多，频带效率就越高。 调制时，将输入信息分成两部分：一部分进行幅度调制；另一部分进行相位调制。对于星型8QAM信号，每个码元由3个比特组成，可将它分成第一个比特和后两个个比特两部分。前者用于改变信号矢量的振幅，后者用于差分相位调制，通过格雷编码来改变当前码元信号矢量相位与前一码元信号矢量相位之间的相位差。 QAM是一种高效的线性调制方式，常用的是8QAM，16QAM，64QAM等。当随着M 的增大，相应的误码率增高，抗干扰性能下降。 2 QAM星座图 QAM调制技术对应的空间信号矢量端点分布图称为星座图。QAM的星座图呈现星状分
MATLAB Simulink系统建模与仿真 实验报告
MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************
实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：
1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：
matlab机电系统仿真大作业
一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环
3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置
M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置
实验一 典型环节的MATLAB仿真汇总
实验一 典型环节的MATLAB 仿真 一、实验目的 1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。 2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。 3．定性了解各参数变化对典型环节动态特性的影响。 二、SIMULINK 的使用 MATLAB 中SIMULINK 是一个用来对动态系统进行建模、仿真和分析的软件包。利用SIMULINK 功能模块可以快速的建立控制系统的模型，进行仿真和调试。 1．运行MATLAB 软件，在命令窗口栏“>>”提示符下键入simulink 命令，按Enter 键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK 仿真 环境下。 2．选择File 菜单下New 下的Model 命令，新建一个simulink 仿真环境常规模板。 3．在simulink 仿真环境下，创建所需要的系统 三、实验内容 按下列各典型环节的传递函数，建立相应的SIMULINK 仿真模型，观察并记录其单位阶跃响应波形。 ① 比例环节1)(1=s G 和2)(1=s G 实验处理：1)(1=s G SIMULINK 仿真模型
波形图为： 实验处理：2)(1=s G SIMULINK 仿真模型 波形图为： 实验结果分析：增加比例函数环节以后，系统的输出型号将输入信号成倍数放大. ② 惯性环节11)(1+= s s G 和15.01)(2+=s s G 实验处理：1 1 )(1+=s s G SIMULINK 仿真模型
波形图为： 实验处理：1 5.01 )(2+= s s G SIMULINK 仿真模型 波形图为： 实验结果分析：当1 1 )(1+= s s G 时，系统达到稳定需要时间接近5s,当
电机大作业(MATLAB仿真-电机特性曲线)
电机大作业 专业班级：电气XXXX 姓名：XXX 学号：XXX 指导老师：张威
一、研究课题（来源：教材习题 4-18 ） 1. 74 、R 2 0.416 、X 2 3.03 、R m 6. 2 X m 75 。电动机的机械损耗p 139W,额定负载时杂散损耗p 320W, 试求额定负载时的转差率、定子电流、定子功率因数、电磁转矩、输出转矩和效 率。 二、编程仿真 根据T 形等效电路: 3D - R Q 运用MATLAB 进行绘图。MATLAB 文本中，P N PN ，U N UN ，尺 R 1， X 1 X1 , R 2 R 2，X 2 X 2，R m Rm, X m Xm ，p pjixiesunh ao ， p pzasansunhao 。定子电流I11，定子功率因数 Cosangle1，电磁转矩Te ， 效率 Xiaolv 。 1.工作特性曲线绘制 MATLA 文本： R1=0.715;X 仁1.74;Rm=6.2;Xm=75;R2=0.416;X2=3.03;pjixiesu nhao=139; pzasa nsu nhao=320;p=2;m 仁 3; ns=1500;PN=17000;UN=380;fN=50; Z1=R1+j*X1; Zm=Rm+j*Xm; for i=1:2500 s=i/2500; nO=n s*(1-s); Z2=R2/s+j*X2; Z=Z1+Zm*Z2/(Zm+Z2); 有一台三相四极的笼形感应电动机， 参数为P N 17kW 、U N 380V (△联 Rm 结）、尺 0. 715 、X j lcr S
U1=UN; I1=U1/Z; l110=abs(l1); An gle 仁an gle(ll); Cosa ngle10=cos(A ngle1); P仁3*U1*l110*Cosa ngle10; l2=l1*Zm/(Zm+Z2); Pjixie=m1*(abs(I2))A2*(1-s)/s*R2; V=(1-s)*pi*fN; Te0=Pjixie/V; P20=Pjixie-pjixies un hao-pzasa nsun hao; Xiaolv0=P20/P1; P2(i)=P20; n (i)=n0; l11(i)=l110; Cosa ngle1(i)=Cosa ngle10; Te(i)=Te0; Xiaolv(i)=Xiaolv0; hold on; end figure(1) plot(P2, n); xlabel('P2[W]');ylabel(' n[rpm]'); figure(2) plot(P2,l11); xlabel('P2[W]');ylabel('l1[A]'); figure(3) plot(P2,Cosa nglel); xlabel('P2[W]');ylabel('go nglvyi nshu'); figure(4) plot(P2,Te); xlabel('P2[W]');ylabel('Te[Nm]'); figure(5) plot(P2,Xiaolv); xlabel('P2[W]');ylabel('xiaolv');
MATLAB控制系统与仿真设计
MATLAB控制系统与仿真 课 程 设 计 报 告 院（系）：电气与控制工程学院 专业班级：测控技术与仪器1301班 姓名：吴凯 学号：1306070127
指导教师：杨洁昝宏洋 基于MATLAB的PID恒温控制器 本论文以温度控制系统为研究对象设计一个PID控制器。PID控制是迄今为止最通用的控制方法，大多数反馈回路用该方法或其较小的变形来控制。PID控制器（亦称调节器）及其改进型因此成为工业过程控制中最常见的控制器(至今在全世界过程控制中用的84%仍是纯PID调节器，若改进型包含在内则超过90%)。在PID控制器的设计中，参数整定是最为重要的,随着计算机技术的迅速发展，对PID参数的整定大多借助于一些先进的软件，例如目前得到广泛应用的MATLAB仿真系统。本设计就是借助此软件主要运用Relay-feedback法，线上综合法和系统辨识法来研究PID控制器的设计方法，设计一个温控系统的PID控制器，并通过MATLAB中的虚拟示波器观察系统完善后在阶跃信号下的输出波形。 关键词：PID参数整定；PID控制器；MATLAB仿真。 Design of PID Controller based on MATLAB Abstract This paper regards temperature control system as the research object to design a pid controller. Pid control is the most common control method up until now; the great majority feedback loop is controlled by this method or its small deformation. Pid controller (claim regulator also) and its second generation so become the most common controllers in the industry process control (so far, about 84% of the controller being used is the pure pid controller, it’ll exceed 90% if the second generation included). Pid parameter setting is most important in pid controller designing, and with the rapid development of the computer technology, it mostly recurs to some advanced software, for example, mat lab simulation software widely used now. this design is to apply that soft mainly use Relay feedback law and synthetic method on the line to study pid
MATLAB与系统仿真实验指导书
实验指导书（标准格式） 《MATLAB与系统仿真》 实验指导书 课程编号：ME2121025课程名称：MATLAB与系统仿真适应专业：电气工程及其自动化课程类别：限选 实验教学种类：上机课程总学时：30＋6 实验学时：6 执笔人：勾燕洁 西安电子科技大学机电工程学院 2006 年9 月
目录 实验一MA TLAB的基本使用及矩阵操作------------------- 3 实验二MA TLAB编程与图形处理---------------------------- 9 实验三MA TLAB中的数值运算与系统仿真---------------- 15
实验一MA TLAB的基本使用及矩阵操作 一、实验目的和要求 熟悉MA TLAB的界面和基本操作，掌握矩阵的建立方法及各种运算。 1．熟悉MA TLAB软件的界面和帮助系统。 2．掌握MA TLAB软件中关于矩阵建立、矩阵初等变换以及矩阵算术、关系、逻辑运算的各种命令。 3．掌握MA TLAB软件中M函数和M文件的编写，以及程序结构与控制，学会编写一般程序。 二、实验内容 1．启动与退出 2．数、数组、矩阵的输入 3．MA TLAB的基本命令 4．矩阵大小的测试 5．矩阵元素的操作 6．特殊矩阵的产生 7．矩阵的算术运算 8．矩阵的关系运算 9．矩阵的逻辑运算 10．常用函数 三、实验仪器、设备（软、硬件）及仪器使用说明 PC机一台，注意正确开、关机及打开软件。 四、实验原理 无 五、实验方法与步骤 1．启动与退出 双击MA TLAB图标，进入MA TLAB命令窗口，即可输入命令，开始运算； 观察各个窗口以及菜单； 单击File菜单中的Exit，或使用MATLAB的Exit命令退出。
运动控制MATLAB仿真
大作业: 直流双闭环调速MATLAB仿真 运动控制技术课程名称： 名：姓电气学院院：学 自动化业：专 号：学 孟濬指导教师： 2012年6月2日
------------------------------------- -------------学浙大江 李超 一、Matlab仿真截图及模块功能描述 Matlab仿真截图如下，使用Matlab自带的直流电机模型： 模块功能描述： ⑴电机模块（Discrete DC_Machine）：模拟直流电机 ⑵负载转矩给定（Load Torque）：为直流电机添加负载转矩 ⑶Demux：将向量信号分离出输出信号 ⑷转速给定（Speed Reference）：给定转速 ⑸转速PI调节（Speed Controller）：转速PI调节器，对输入给定信号与实际信号
的差值进行比例和积分运算，得到的输出值作为电流给定信号。改变比例和积分运算系数可以得到不同的PI控制效果。 ⑹电流采样环节（1/z）：对电流进行采样，并保持一个采样周期 ⑺电流滞环调节（Current Controller）：规定一个滞环宽度，将电流采样值与给定值进行对比，若：采样值＞给定值+0.5*滞环宽度，则输出0； 若：采样值＜给定值—0.5*滞环宽度，则输出1； 若：给定值—0.5*滞环宽度＜采样值＜给定值+0.5*滞环宽度，则输出不变 输出值作为移相电压输入晶闸管斩波器控制晶闸管触发角 ：根据输入电压改变晶闸管触发角，从而改变电机端电压。GTO⑻晶闸管斩波．⑼续流二极管D1：在晶闸管关断时为电机续流。 ⑽电压传感器Vd：测量电机端电压 ⑾示波器scope：观察电压、电流、转速波形 系统功能概括如下：直流电源通过带GTO的斩波器对直流电机进行供电，输出量电枢电流ia和转速wm通过电流环和转速环对GTO的通断进行控制，从而达到对整个电机较为精确的控制。 下面对各个部分的功能加以详细说明： （1）直流电机 双击电动机模块，察看其参数：
基于MATLAB的系统仿真实验
基于MATLAB的系统仿真实验 实验指导书 新乡学院 二○○八年三月
说明 MATLAB是MATrix LABoratory的缩写，是一种基于矩阵的数学与工程计算系统，可以用作动态系统的建模与仿真。研究系统的结构和参数的变化对系统性能的影响可利用MATLAB强大的计算和作图功能，因此本实验采用MATLAB仿真研究连续控制系统和离散控制系统的性能分析过程。通过该实验提高学生对控制系统的分析与设计能力，加深对《自动控制原理》课程内容的理解。 由于学时有限，该实验由同学们参考有关MATLAB应用的书籍利用课后时间完成。
第一部分 基于MATLAB 连续系统的仿真 一、 实验目的 在研究系统的结构和参数的变化对系统性能的影响时，采用解析和作图的方法比较麻烦，而且误差也大，用MATLAB 仿真实现则简单方便，精度高。本实验采用MATLAB 实现控制系统的数学描述、控制系统的时域分析及根轨迹和频率特性分析。通过该实验，加深学生对系统阶次，型号，参数与系统性能的关系的理解。 二、实验环境 在计算机Windows 环境下安装好MATLAB6.3以上版本后，双击MATLAB 图标或成“开始”菜单打开MATLAB ，即可进入MATLAB 集成环境。 三、MATLAB 应用实例 1．拉氏变换和反变换 例 求22)(2++=t t t f 的拉氏变换 解 键入 syms s t; ft=t^2+2*t+2; st=laplace(ft,t,s) 运行结果为 st= 2/s^3+2/s^2+2/s 例 求) 2)(34(6 )(2 ++++= s s s s s F 的拉氏反变换 解 键入 syms s t; Fs=(s+6)/(s^2+4*s+3)/(s+2); ft=ilaplace(Fs,s,t) 运行结果为 ft= 3/2*exp(-3*t)+5/2*exp(-t)-4*exp(-2*t) 2。求根运算 例 求多项式 43)(2 3 ++=s s s p 的根，再由根建多项式。 解 键入 p=[1 3 0 4]; r=root(p) 运行结果为 r= -3.3553 0.1777+1.0773i 1.7777-1.0773i
华科matlab大作业
MATLAB语言、控制系统分析与设计 大作业 升降压斩波电路MATLAB仿真 专业：电气工程及其自动化 班级： 设计者： 学号： 评分： 华中科技大学电气与电子工程学院 2016 年1月
评分栏
目录 一、引言 (4) 二、电路设计与仿真 (4) 三、仿真结果及分析 (7) 四、深入讨论 (10) 五、总结 (10) 六、参考资料 (11)
升降压斩波电路MATLAB 仿真 一．引言 Buck/Boost 变换器是输出电压可低于或高于输入电压的一种单管直流变换器，其主电路与Buck 或Boost 变换器所用元器件相同，也有开关管、二极管、电感和电容构成。与Buck 和Boost 电路不同的是，电感L 在中间，不在输出端也不在输入端，且输出电压极性与输入电压相反。开关管也采用PWM 控制方式。Buck/Boost 变换器也有电感电流连续喝断续两种工作方式，在此只讨论电感电流在连续状态下的工作模式。 二．电路设计与仿真 1、电路原理 当可控开关V 处于通态时，电源E 经V 向电感L 供电使其储存能量，此时电流为i1。同时，电容C 维持输出电压基本恒定并向负载R 供电。此后，使V 关断，电感L 中储存的能量向负载L 释放，电流为i2。负载电压极性为上正下负，与电源电压极性相反。 稳态时，一个周期T 内电感L 两端电压UL 对时间积分为零，即 当V 处于通态期间，UL=E ；而当V 处于断态期间，UL=-Uo 。于是 所以，输出电压为 =?dt T L U off o on t U Et =E E t T t E t t U on on off on o α α -=-== 1
《MATLAB与控制系统。。仿真》实验报告
《MATLAB与控制系统仿真》 实验报告 班级： 学号： 姓名： 时间：2013 年 6 月
目录实验一MATLAB环境的熟悉与基本运算（一）实验二MATLAB环境的熟悉与基本运算（二）实验三MATLAB语言的程序设计 实验四MATLAB的图形绘制 实验五基于SIMULINK的系统仿真 实验六控制系统的频域与时域分析 实验七控制系统PID校正器设计法 实验八线性方程组求解及函数求极值
实验一MATLAB环境的熟悉与基本运算（一） 一、实验目的 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验基本原理 1.熟悉MATLAB环境: MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MATLAB常用命令 表1 MATLAB常用命令 变量与运算符 3．1变量命名规则 3．2 MATLAB的各种常用运算符 表3 MATLAB关系运算符 表4 MATLAB逻辑运算符
| Or 逻辑或 ~ Not 逻辑非 Xor逻辑异或 符号功能说明示例符号功能说明示例 ：1:1:4;1:2:11 . ；分隔行.. ，分隔列… （）% 注释 [] 构成向量、矩阵！调用操作系统命令 {} 构成单元数组= 用于赋值 的一维、二维数组的寻访 表6 子数组访问与赋值常用的相关指令格式 三、主要仪器设备及耗材 计算机 四．实验程序及结果 1、新建一个文件夹（自己的名字命名，在机器的最后一个盘符） 2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。 3、学习使用help命令。
《控制系统MATLAB仿真》实验讲义88
《自动控制原理实验》 目录 第一部分实验箱的使用 第二部分经典控制实验 第一章基本实验 实验一典型环节及其阶跃响应 实验二二阶系统阶跃响应 实验三控制系统的稳定性分析 实验四控制系统的频率特性 实验五连续控制系统的串联校正 实验六数字PID控制实验 第二章综合实验 第三部现代控制理论实验 第一章基本实验 第二章综合实验
实验一 典型环节及其阶跃响应 预习要求： 1、复习运算放大器的工作原理；了解采用A μ741运算放大器构成各种运算电路的方法； 2、了解比例控制、微分控制、积分控制的物理意义。 一、实验目的 1、学习自动控制系统典型环节的电模拟方法，了解电路参数对环节特性的影响。 2、学习典型环节阶跃响应的测量方法； 3、学会根据阶跃响应曲线计算确定典型环节的传递函数。 二、实验内容 1、比例环节 电路模拟： 图1-1 传递函数： 2211 ()()()U s R G s U s R ==- 2、惯性环节 电路模拟： 图1-2 传递函数： 22112()/()()11 U s R R K G s U s Ts R Cs = =-=- ++ 3、积分环节 电路模拟： A/D1 D/A1 A/D1
图1-3 传递函数： 21()11 ()()U s G s U s Ts RCs = =-=- 4、微分环节 电路模拟： 图1-4 传递函数： 211() ()() U s G s s RC s U s τ= =-=- 5、比例微分 电路模拟： 图1-5 传递函数： 222111 ()()(1)(1)()U s R G s K s R C s U s R τ= =-+=-+ 6、比例积分 电路模拟： 图1-6 A/D1 2 R D/A1 A/D1 A/D1 A/D1 C
MATLAB与系统仿真实验手册
MATLAB与控制系统仿真实验指导书 河北大学电子信息工程学院 2006年9月
目录 实验一MATLAB基本操作与基本运算 (1) 实验二M文件及数值积分仿真方法设计 (3) 实验三MATLAB 的图形绘制 (4) 实验四函数文件设计和控制系统模型的描述 (6) 实验五控制系统的分析与设计 (7) 实验六连续系统离散化仿真方法设计 (8) 实验七SIMULINK 仿真 (9) 实验八SIMULINK 应用进阶 (10) 附录MATLAB常用函数 (12)
实验一MATLAB基本操作与基本运算 一、实验目的及要求： 1．熟悉MATLAB6.5的开发环境； 2．掌握MATLAB6.5的一些常用命令； 3．掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容： 1.熟悉MATLAB6的开发环境： ①MATLAB的各种窗口： 命令窗口、命令历史窗口、启动平台窗口、工作空间窗口、当前路径窗口。 图1 MA TLAB界面窗口 ②路径的设置： 建立自己的文件夹，加入到MA TLAB路径中，并保存。 设置当前路径，以方便文件管理。 ③改变命令窗口数据的显示格式 >> format short >> format long
然后键入特殊变量：pi （圆周率），比较显示结果。 2.掌握MATLAB 常用命令 >> who %列出工作空间中变量 >> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量 >>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称 如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。 3. 在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ 4. 设计M 文件计算: 已知 求出： 1) a+b a-b a+b*5 a-b+I （单位阵） 2) a*b a.*b a/b 3) a^2 a.^2 注意：点运算 . 的功能，比较结果。 5. 设计M 文件计算: x=0:0.1:10 当sum>1000时停止运算，并显示求和结果及计算次数。 三、实验报告要求： 1）体会1、2的用法； 2）对3、4、5写出程序及上机的结果。 ?? ??? ?????=987654321a ?? ?? ? ?????=300120101b ) 2(100 2 i i i x x sum -=∑=