shiyan3.m

系统校正,对MATLAB学习者帮助很大

clear
num=100
den1=[1 0]
den2=[0.1 1]
den3=[0.01 1]
den=conv(den1,conv(den2,den3))
%系统闭环传递函数
[numc,denc]=cloop(num,den,-1)
wt=0:0.5:60
t=0:0.01:1.5
%系统频率响应数据
g=freqs(numc,denc,wt)
mag=abs(g)
figure(1)
subplot(121)
bode(num,den)
pause
subplot(222)
step(numc,denc,t)
grid
pause

subplot(224)
plot(wt,mag)
title('C-loop Freguency Response-Amplitude')
xlabel('Frequency-rad')
ylabel('Amplitude')
grid
%Gc(s)=(a*T*s+1)/(T*s+1);a>1
gama=55;wc=50
%超前校正装置传递函数
a=(1+sin(gama*pi/180))/(1-sin(gama*pi/180))
T=1/(wc*sqrt(a))
num=[0 0 0 100]
numa=[a*T 1]
dena=[T 1]
numao=conv(num,numa)
denao=conv(den,dena)
figure(2)
bode(numao,denao)
[numac,denac]=cloop(numao,denao,-1)
figure(3)
step(numac,denac)
%显示超前校正装置的参数
disp('a T a*T')
a 
T
a*T
%显示超前校正装置传递函数多项式
numa
dena
%求取校正后系统的相角裕量、截止角频率
[Gm1,Pm1,Wcg1,Wcm1]=margin(numao,denao)
Pm1
Wcm1
wc=5;g=25
beta=10^(g/20)
T=4/wc
wt=logspace(-2,2)
numb=[T 1]
denb=[beta*T 1]
numbo=conv(num,numb)
denbo=conv(den,denb)
figure(4)
bode(numbo,denbo,wt)
[numbc,denbc]=cloop(numbo,denbo,-1)
figure(5)
step(numbc,denbc)
%
beta
T
numb
denb
[Gm2,Pm2,Wcg2,Wcm2]=margin(numbo,denbo)
Pm2
Wcm2
%
%相角裕量、截止角频率
gama=50;wc=20
a=(1+sin(gama*pi/180))/(1-sin(gama*pi/180))
T1=1/(wc*sqrt(a))
numa=[a*T1 1]
dena=[T1 1]
numao=conv(numa,num)
denao=conv(dena,den)
figure(6)
bode(numao,denao,wt)
%相角裕量
g=14
beta=10^(g/20)
T=5/wc
wt=logspace(-1,3)
numb=[T 1]
denb=[beta*T 1]
numo=conv(conv(num,numa),numb)
deno=conv(conv(den,dena),denb)
figure(7)
bode(numo,deno,wt)
a
T1
numa
dena
beta
T
numb
denb
[numc,denc]=cloop(numo,deno,-1)
figure(8)
step(numc,denc)