matlab中怎么传函绘图,基于MATLAB的控制工程.ppt

mp1=(ymax-yss)/yss; % s=1001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;end ts1=(s-1)*0.001; % r=1;while y2(r)1-dta &y3(s)<1+dta;s=s-1;end ts2=(s-1)*0.001; % r=1;while y3(r)1-dta & y3(s)<1+dta;s=s-1;end ts3=(s-1)*0.001 % [tr1 tp1 mp1 ts1;tr2 tp2 mp2 ts2;tr3 tp3 mp3 ts3] 利用MATLAB进行频率特性分析 求取系统对数频率特性:bode( ) [mag,phase,w]= bode(num,den) [mag,phase,w]= bode(num,den,w) 格式： 幅值 相角 频率 求取系统对数频率特性: bode( ) bode(num,den)：可绘制出以连续时间多项式传递函数表示的系统的波特图。 bode(num,den,w)：利用指定的角频率矢量w绘制出系统的波特图。 当输出变量[mag,pha,w]或[mag,pha]引用函数时，可得到系统波特图相应的幅值mag、相角pha及角频率点矢量w或只是返回幅值与相角。相角以度为单位，幅值可转换为分贝单位： magdb=20×log10(mag) 二阶系统，已知 =2，阻尼比 =0.7。试求系统的传递函数，并绘制系统的Bode图。 [num,den]=ord2(2,0.7); G=tf(num,den); bode(G); grid on;%给Bode图加上网格 先利用ord2函数求出系统的传递函数模型，再利用bode函数绘制系统的Bode图。 例子 求取系统奈奎斯特图：nyquist( ) [re,im,w]= nyquist(num,den) [re,im,w]= nyquist(num,den,w) 格式： nyquist(num,den)可绘制出以连续时间多项式传递函数表示的系统的极坐标图。 nyquist(num,den,w)：可利用指定的角频率矢量w绘制出系统的极坐标图。 当不带返回参数时，直接在屏幕上绘制出系统的极坐标图(图上用箭头表示w的变化方向，负无穷到正无穷)。当带输出变量[re,im,w]引用函数时，可得到系统频率特性函数的实部re和虚部im及角频率点w矢量(为正的部分)。可以用plot(re,im)绘制出对应w从负无穷到零变化的部分。 nyquist(num,den) %注意不能直接用nyquist(G)，这样会造成输入向量不够， grid on 求系统的nyquist图. num=[2 0 3]; den=[5 2 3 1]; G=tf(num,den) 例子 s 作 的系统的nyquist和Bode图 %建模 k=24;nunG1=k*[0.25 1]; denG1=conv([5 2],[0.05 2]); %求实频和虚频函数 [re,im]=nyquist(nunG1,denG1); %画nyquist图 plot(re,im);grid on pause %画Bode图 w=logspace(-2,3 ,100); bode(nunG1,denG1,w);grid on 判断系统的稳定性 已知系统的开环传递函数为 绘制系统nyquist图并判断闭环系统稳定性 先用nyquist函数绘制出系统的nyquist图; 再利用nyquist判据判断闭环系统的稳定性. 思路： 例子 从上图可看出nyquist曲线逆时针包围(-1,j0)点1圈，而从开环系统的零极点图可知开环系统在右平面有1极 点，所以闭环系统稳定。 求相对稳定性 第一步: 先求出系统传递函数 G=tf(num,den) 第二步: 利用margin函数求出幅值裕量,相位裕度以及相应的频率. [Gm,Pm,Wcp,Wcg]=margin(G) 第三步: 将幅值裕量转化为幅值裕度 GmdB=20*log10(Gm) 利用MATLAB直接求相对稳定性 例子1: num=con

最后

以上就是沉默柜子为你收集整理的matlab中怎么传函绘图,基于MATLAB的控制工程.ppt的全部内容，希望文章能够帮你解决matlab中怎么传函绘图,基于MATLAB的控制工程.ppt所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。