基于MATLAB自动控制原理.pptVIP

基于MATLAB/Simulink的自动控制原理课程设计 强盛 2011-11-22 内容 用SIMULINK搭建系统模型 编辑M文件，将SIMULINK搭建的模型转换成传递函数模型 绘制伯德图(bode(),margin())、单位阶跃响应曲线(step()) 等图形 图形的提取及坐标轴的标注(title(),xlabel(),ylabel()) 举例：已知单位负反馈系统的固有部分的传递函数为设计指标： 2、MATLAB中M文件的建立 二、串联滞后校正 /*先画G=15/s(0.5s+1)渐近幅频特性 num=15;den=[0.5 1 0]; G=tf(num,den); w=logspace(-2,1,200); [wpos,ypos]=bd_asymp(G,w); /*求折点 semilogx(wpos,ypos) hold on /*再画实际幅频特性 [mag,phase,w]=bode(num,den,w); semilogx(w,20*log10(mag)) 四、闭环系统的单位阶跃响应 * 一、开环系统分析 1、Simulink中模型的建立 [a,b,c,d]=linmod2(sim1)g=tf(ss(a,b,c,d))margin(g) function [wpos,ypos]=bd_asymp(G,w) G1=zpk(G); Gtf=tf(G); if nargin==1, w=freqint2(Gtf.num{1},Gtf.den{1},100); end zer=G1.z{1}; pol=G1.p{1}; gain=G1.k; wpos=[ ]; pos1=[ ]; for i=1:length(zer); if isreal(zer(i)) wpos=[wpos, abs(zer(i))]; pos1=[pos1,20]; else if imag(zer(i))0 wpos=[wpos, abs(zer(i))]; pos1=[pos1,40]; end, end, end for i=1:length(pol); if isreal(pol(i)) wpos=[wpos, abs(pol(i))]; pos1=[pos1,-20]; else if imag(pol(i))0 wpos=[wpos, abs(pol(i))]; pos1=[pos1,-40]; end, end, end wpos=[wpos w(1) w(length(w))]; pos1=[pos1,0,0]; [wpos,ii]=sort(wpos); pos1=pos1(ii); ii=find(abs(wpos)eps); kslp=0; w_start=1000*eps; if length(ii)0, kslp=sum(pos1(ii)); ii=(ii(length(ii))+1):length(wpos); wpos=wpos(ii); pos1=pos1(ii); end while 1 [ypos1,pp]=bode(G,w_start); if isinf(ypos1), w_start=w_start*10; else, break;end end wpos=[w_start wpos]; ypos(1)=20*log10(ypos1); pos1=[kslp pos1]; for i=2:length(wpos) kslp=sum(pos1(1:i-1)); ypos(i)=ypos(i-1)+... kslp*log10(wpos(i)/wpos(i-1)); End ii=find(wpos=w(1)wpos=w(length(w))); wpos=wpos(ii); ypos=ypos(ii); 三、校正后的Bode图 *