Matlab函数选编.ppt

五、Routh稳定性分析 [rtab,msg]=routh(den) 超前滞后补偿器设计 Gc=leadlagc(sys1,wc,gamma,kv,key)---超前滞后补偿器设计 * 一、ＭＡＴＬＡＢ中连续系统模型表示方法 二、求连续系统的单位脉冲响应 三、求连续系统的单位阶跃响应 用ＭＡＴＬＡＢ求取瞬态响应 一、ＭＡＴＬＡＢ中连续系统模型表示方法 １、连续系统多项式模型 ＭＡＴＬＡＢ表示方法 分子多项式　num=[b0 b1 … bm-1 bm] 分母多项式　den=[a0 a1 … an-1 an] 系统表示方法　　　（num , den) sys1=tf(num,den) ２、连续系统零极点模型 ＭＡＴＬＡＢ表示方法 比例系数　Ｋ＝k 分子　Z=[-z1 , -z2 , … , -zm] 若无zi，则Z=[] 分母　Ｐ=[-p1 ,- p2 , … , -pm] 系统表示方法　sys1=zpk(Z,P,K) 模型转换　　　［num ,den]=zp2tf（Z , P , K ) [Z,P,K]=tf2zp(num,den) 二、连续系统的单位脉冲响应 例一：求如下系统的单位脉冲响应 %example 1 num=[1.9691 , 5.0395] den=[1,0.5572,0.6106] impulse(num,den) end 在ＭＡＴＬＡＢ的Ｅditor/Debugger输入程序 在ＴＯＯＬＳ菜单中选择ＲＵＮ得到结果 三、连续系统的单位阶跃响应 %example ２ num=[1.9691 , 5.0395] den=[1,0.5572,0.6106] step(num,den) end 例二：求如下系统的单位阶跃响应 在ＭＡＴＬＡＢ的Ｅditor/Debugger输入程序 在ＴＯＯＬＳ菜单中选择ＲＵＮ得到结果 Sys1=tf(num,den) y=step(Sys1,t) [y,t]=step(Sys1) n=100;d=[1,12,30,100]; g1=tf(n,d) t1=0:0.1:5; y1=step(n,d,t1); y1=step(g1,t1); n1=10;d1=[1,2,10]; g2=tf(n1,d1) t2=0:0.1:5; y2=step(n1,d1,t2); plot(t1,y1,‘-r’,t2,y2)； grid 主导极点： n=100;d=[conv ([1,10], [1,2,10])]; 4.5 用MATLAB绘制系统的根轨迹 一、求开环传递函数的零极点 例：已知系统的开环传递函数 求系统开环零、极点的位置 num=[2 5 1]; den=[1 2 3]; pzmap(num,den); title{‘pole-zero Map’} 分子多项式 分母多项式 求零极点函数 打印标题 二、绘制常规根轨迹 例：已知系统的开环传递函数 绘制该系统的根轨迹图 num=[2 5 1]; den=[1 2 3]; rlocus(num,den); Title（‘控制系统根轨图’） 分子多项式 分母多项式 绘制根轨迹函数 打印标题 三、绘制带阻尼比和自然振荡频率栅格 例：已知系统的开环传递函数 绘制系统的根轨迹图及带阻尼比和自然振荡频率栅格。 num=[2 5 1]; den=[1 2 3]; rlocus(num,den); sgrid Title（‘控制系统根轨图和栅格’） 分子多项式 分母多项式 绘制根轨迹函数 打印标题 绘制栅格 四、系统性能分析 例：已知系统的开环传递函数 绘制系统的根轨迹图并分析系统性能。 num=[2 5 1]; den=[1 2 3]; rlocus(num,den); Sgrid; [k,p]=rlocfind(num,den); Title（‘控制系统根轨图和栅格’）; [xx,yy]=rlasymp(num,den); 分子多项式 分母多项式 绘制根轨迹函数 打印标题 绘制栅格 求十字光标处的闭环极点 根轨迹渐近线 5.7 用MATLAB进行系统的频域分析 一、绘制BODE图 已知二阶系统的传递函数 绘制系统的Bode图。 num=1 den=[1,0.2,1] bode(num,den) grid 二、绘制Nyquist图 绘制系统 乃奎斯特图（Nyquist图）。 num= [2,5,1] den=[1,2,3] nyquist(num,den) 三、求