lx2.m

matlab与控制系统仿真的PPt多媒体教程和源码

%lx2.m
clc
clear
close all
%第一步：绘制未校正系统的根轨迹
num=[1 10];den=[1 2 20];%未校正系统传函
rlocus(num,den);

%第二步：根据时域调节时间、相角裕度可以得到：阻尼比zeta=0.42和自然振荡频率wn=1.67
%由此在未校正系统的根轨迹上确定校正后的预期主导极点；
hold on
zeta=0.42;wn=1.67;
x=[-10:0.1:0];y=-(sqrt(1-zeta^2)/zeta)*x;
plot(x,y,'r',x,-y,'y')

%第三步：计算与预期主导极点对应的系统增益K
%和未校正系统的位置误差系数Kpuncomp
[k,p]=rlocfind(num,den);
k
Kpuncomp=k/2;
%第四步：计算滞后校正网络参数a
Kpcomp=9;%Kpcomp为期望的系统位置误差系数＝9 (ess=1/(1+Kpcomp)=0.1)；
a=Kpcomp/Kpuncomp
%第五步：根据a的取值(a=3.6)，|z/p|=a，确定滞后校正网络的零极点，
%使得校正后的根轨迹经过预期的主导极点。
z=3.5;p=1;
numc=k*[1 z];denc=[1 p];
[nums,dens]=series(numc,denc,num,den);
%第六步：仿真计算系统闭环阶跃响应，检验设计结果
figure(2)
[numb,denb]=cloop(nums,dens);
step(numb,denb)
figure(3)
margin(nums,dens)