jz_4.m

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

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

%第二步：根据阻尼比zeta=0.59和自然振荡频率wn=2.26
%确定预期主导极点的允许区域，并进一步在未校正系统的根轨迹上确定校正后的预期主导极点；
hold on
zeta=0.59;wn=2.26;
x=[-10:0.1:0];y=-(sqrt(1-zeta^2)/zeta)*x;
plot(x,y,'r',x,-y,'y')

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