lx1_3.m

来自「matlab与控制系统仿真的PPt多媒体教程和源码」· M 代码 · 共 27 行

%lx1_3.m
%由此未校正系统Bode图可以得到wm=4.5rad/s
%第五步：计算超前校正环节的零极点
clc
close all
clear
wm=4.4;
alpha=4;
p=wm*sqrt(alpha);
z=p/alpha;

%第六步：增大系统增益，补偿由超前校正网络带来的增益衰减（1/a）
%取K=10*alpha
%所以有校正网络传函为：Gc(S)=(1+s/z)/(a(1+s/p))

K=10*alpha;
%第七步：绘制校正后的 Bode图，校验系统的相角裕度是否满足设计要求
num=[K];den=[1 1 0];%未校正系统
numc=[1/z 1];denc=[alpha/p alpha];%超前校正网络
[nums dens]=series(numc,denc,num,den);
figure(2)
margin(nums,dens);
%若不满足相角要求，重复前面设计步骤
figure(3)
[numcc dencc]=cloop(nums,dens);
step(numcc,dencc)