基于Matlab的计算机控制技术仿真实验Word格式文档下载.docx

基于Matlab的计算机控制技术仿真实验Word格式文档下载.docx

《基于Matlab的计算机控制技术仿真实验Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《基于Matlab的计算机控制技术仿真实验Word格式文档下载.docx（32页珍藏版）》请在冰豆网上搜索。

f=5*exp（-10*t）subplot（2,1,1）plot（t,f）

gridsubplot（2,1,2）stem（t,f）

grid

请改变采样周期T，观察不同的采样周期下的采样效果。

幅频曲线：

w=-50:

1:

50

F=5./sqrt（100+w.^2）plot（w,F）

若|F（jmax）|=0.1|F（0）|，选择合理的采样周期

T并验加以证

w=-400:

20:

400

ws=200Ts=2*pi/ws

F0=5/Ts*（1./sqrt（100+（w）.^2））F1=5/Ts*（1./sqrt（100+（w-ws）.^2））F2=5/Ts*（1./sqrt（100+（w+ws）.^2））plot（w,F0,w,F1,w,F2）

请改变采样频率ws，观察何时出现频谱混叠？

2.拉式变换和Z变换

使用Matlab求函数的拉氏变换和Z变换拉式变换：

symsawtf1=exp（-a*t）laplace（f1）f2=tlaplace（f2）

f3=t*exp（-a*t）laplace（f3）f4=sin（w*t）

Z变换：

symsakTf1=exp（-a*k*T）ztrans（f1）f2=k*Tztrans（f2）

f3=k*T*exp（-a*k*T）ztrans（f3）f4=sin（a*k*T）

laplace（f4）

f5=exp（-a*t）*cos（w*t）laplace（f5）

反拉式变换symssaf1=1/silaplace（f1）f2=1/（s+a）ilaplace（f2）f3=1/s^2ilaplace（f3）f4=w/（s^2+w^2）ilaplace（f4）

f5=1/（s*（s+2）^2*（s+3））ilaplace（f5）

ztrans（f4）f5=a^kztrans（f5）

反Z变换symszaTf1=z/（z-1）iztrans（f1）

f2=z/（z-exp（-a*T））iztrans（f2）f3=T*z/（z-1）^2iztrans（f3）

f4=z/（z-a）iztrans（f4）

f5=z/（（z+2）^2*（z+3））iztrans（f5）

3.控制系统模型的建立与转化

num

bsm+bsm-1++b

n

传递函数模型：

num=[b1,b2,…bm]，den=[a1,a2,…an]，G（s）==12m

1

2

den

asn+a

sn-1++b

零极点增益模型：

z=[z1,z2,……zm]，p=[p1,p2……pn]，k=[k]，G（s）=k

（s-z1）（s-z2

）（s-zm）

（s-p1）（s-p2）（s-pn）

建立系统模型

s（s+1）s2+s

G（s）==和

（s+2）（s+3）s2+5s+6

z（z+1）z2+z

G（z）==

（z+2）（z+3）z2+5z+6

num=[1,1,0]

den=[1,5,6]T=0.1

Gs1=tf（num,den）Gz1=tf（num,den,T）

z=[0,-1]

p=[-2,-3]k=[1]T=0.1

Gs2=zpk（z,p,k）Gz2=zpk（z,p,k,T）

传递函数模型和零极点增益模型相互转化传递函数模型转化零极点增益模型：

num=[1,1,0]

Gs1=tf（num,den）Gz1=tf（num,den,T）[z,p,k]=tf2zp（num,den）Gs2=zpk（z,p,k）Gz2=zpk（z,p,k,T）

零极点增益模型转化传递函数模型：

p=[-2,-3]k=[2]T=0.1

Gs1=zpk（z,p,k）Gz1=zpk（z,p,k,T）[num,den]=zp2tf（z'

p'

k）Gs2=tf（num,den）Gz2=tf（num,den,T）

（s+1）（s2+2s+2）（z+1）（z2+2z+2）

建立系统模型G（s）=和G（z）=

（s2+2）（s2+4s+8）（z2+2）（z2+4z+8）

num1=[1,1]num2=[1,2,2]

den1=[1,0,2]

den2=[1,4,8]

num=conv（num1,num2）den=conv（den1,den2）T=0.1

四、实验步骤

1.根据参考程序，验证采样定理、拉氏变换和Z变换、控制系统模型建立的方法

2.观察记录输出的结果，与理论计算结果相比较

3.自行选则相应的参数，熟悉上述的各指令的运用方法

五、实验数据及结果分析

记录输出的数据和图表并分析

六、总结

实验二基于Matlab的离散控制系统仿真

1.学习使用Matlab的命令对控制系统进行仿真的方法

2.学习使用Matlab中的Simulink工具箱进行系统仿真的方法

1.控制系统命令行仿真

建立如图所示一阶系统控制模型并进行系统仿真。

一阶系统闭环传递函数为

G（s）=3/s=3，转换为离散系统脉冲传

1+3/ss+3

递函数并仿真。

%模型建立

num=[3]%传递函数分子den=[1,3]%传递函数分母T=0.1%采样周期

gs=tf（num,den）%传递函数模型建立gz=c2d（gs,T,'

zoh'

）%转化为离散系统脉冲传递函数模型

%'

零阶保持器变换

foh'

三角变换（一阶保持器）

tustin'

双线性变换

prewarp'

带频率预畸的双线性变换

matched'

零极点匹配变换

%模型特性

[z,p,k]=tf2zp（num,den）%求零极点pzmap（gs）%零极点图

pzmap（gz）%零极点图grid

rlocus（gs）%根轨迹图

rlocus（gz）%根轨迹图grid

%时间响应

impulse（gs）%单位脉冲响应impulse（gz）%离散单位脉冲响应step（gs）%单位阶跃响应

step（gz）%离散单位阶跃响应

%频率响应

freqs（num,den）%频率响应freqz（num,den）%频率响应

close

bode（gs）%波特图bode（gz）%波特图

nyquist（gs）%奈奎斯特曲线nyquist（gz）%奈奎斯特曲线

nichols（gs）%尼科尔斯曲线nichols（gz）%尼科尔斯曲线

25

R（s）

R（z）

+

-

E（s）

C（s）

C（z）

二阶系统闭环传递函数为G（s）=

52

s（s+2⨯0.4⨯5）

=

s2+2⨯0.4⨯5s+52

，请转换为离散系

1+

统脉冲传递函数并仿真，改变参数，观察不同的系统的仿真结果。

2.控制系统的Simulink仿真

按图建立系统的Simulink模型，对不同的输入信号进行仿真，改变参数，观察不同的仿真结果。

将上述系统离散化并仿真，观察仿真结果。

1.根据实验原理对控制系统进行软件仿真

2.观察记录输出的结果，与理论计算值相比较

3.自行选择参数，练习仿真方法，观察不同的仿真结果

记录输出的数据和图表并分析。

实验三连续系统PID控制器设计及其参数整定

一、实验目的

（1）掌握PID控制规律及控制器实现。

（2）对给定系统合理地设计PID控制器。

（3）掌握对给定控制系统进行PID控制器参数在线实验工程整定的方法。

二、实验原理

在串联校正中，比例控制可提高系统开环增益，减小系统稳态误差，提高系统的控制精度，但会降低系统的相对稳定性，甚至可能造成系统闭环系统不稳定；

积分控制可以提高系统的型别（无差度），有利于提高系统稳态性能，但积分控制增加了一个位于原点的开环极点。

使信号产生90°

的相位滞后，于系统的稳定不利，故不宜采用单一的积分控制器；

微分控制规律能反映输入信号的变化趋势，产生有效的早期修正信号，以增加系统的阻尼程度，从而改善系统的稳定性，但微分控制增加了一个-1/τ的开环零点，使系统的相角裕度提高，因此有助于系统稳态性能的改善。

在串联校正中，PI控制器增加了一个位于原点的开环极点，同时也增加了一个位于s左半平面的开环零点。

位于原点的开环极点可以提高系统的型别（无差度），减小稳态误差，有利于提高系统稳态性能；

负的开环零点可以减小系统的阻尼，缓和PI极点对系统产生的不利影响。

只要积分时间常数Ti足够大，PI控制器对系统的不利影响可大为减小。

PI控制器主要用来改善控制系统的稳态性能。

在串联校正中，PID控制器增加了一个位于原点的开环极点，和两个位于s左半平面的开环零点。

除了具有PI控制器的优点外，还多了一个负实零点，动态性能比PI更具有优越性。

通常应使积分发生在低频段，以提高系统的稳态性能，而使微分发生在中频段，以改善系统的动态性能。

PID控制器传递函数为Ge（s）=Kp（1+1/Tis+Tds），注意工程PID控制器仪表中比例参数整定常用比例度δ%，δ%=1/Kp*100%.

三、实验内容

（1）Ziegler-Nichols——反应曲线法

反应曲线法适用于对象传递函数可以近似为e-Ls的场合。

先测出系统处于开环状态下的对象动态特性（即先输入阶跃信号，测得控制对象输出的阶跃响应曲线），如图6-25所示，然后根据动态特性估算出对象特性参数，控制对象的增益K、等效滞后时间L和等效时间常数T，然后根据表6-4中的经验值选取控制器参