基于MATLAB的控制工程资料.pptVIP

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)yss;r=r+1;end tr2=(r-1)*0.001;[ymax,tp]=max(y2); tp2=(tp-1)*0.001;mp2=(ymax-yss)/yss; s=1001;while y2(s)1-dta y3(s)1+dta;s=s-1;end ts2=(s-1)*0.001; % r=1;while y3(r)yss;r=r+1;end tr3=(r-1)*0.001;[ymax,tp]=max(y3); tp3=(tp-1)*0.001; mp3=(ymax-yss)/yss; s=1001;while y3(s)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( )