利用matlab由开环传递函数求闭环传递函数并求其单位冲击与阶跃响应

《利用matlab由开环传递函数求闭环传递函数并求其单位冲击与阶跃响应》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、利用matlab由开环传递函数求闭环传递函数并求其单位冲击和阶跃响应例：设有一个系统的开环传递函数如下函数, 并绘制输出阶跃响应曲线和脉冲响应曲线 解：编程（见：workCT_tchresp_2_20110522）clear all;close all;%%%%%%%%%%%%%%a0 = [0 0 0 0 0.8 ];            % a0: 开环传递函数分子多项式系数b1 = [1 0 ] ;                 %  sb2 = [0.3 1 ] ;          

2、      % b3 = [0.5 0.7 1] ;             %  (0.5 s2 + 0.7s + 1)b0 = conv(b1,conv(b2,b3));    % b0: 开环传递函数分母多项式系数%%%% 闭环传递函数aa = a0 ;                    %  aa：闭环传递函数分子多项式系数bb = b0 + a0;                %  bb：闭环传递函数分子多项式系数disp ('System Closed Loop Transfer Fu

3、nction is :')aabb%%%% 计算：阶跃响应t = 0 : 0.1 : 20y = step (aa, bb, t);           %  阶跃响应%%%% 绘制：阶跃响应figure(1)plot(t ,y) ;title('阶跃响应');            %  标题：阶跃响应xlabel('时间 /s');            %  横坐标ylabel('幅值');              %  纵坐标grid;                           % 

4、Íø¸ñ%%%% 计算：脉冲响应figure(2)yy = impulse (aa, bb, t);       %  标题：脉冲响应plot(t, yy);title('脉冲响应');xlabel('时间/s');ylabel('幅值');grid;                       %  网格%%%%  绘制：脉冲响应wt = logspace (-1,1);                % 对数空间(¨0.1, 10)[mag,phase] = bode (a0 ,b0 ,wt);  

5、     % 计算：Bode中的幅值和相位[ Gm,Pm,Wcg,Wcm] = margin(a0,b0);  % 计算：稳定裕度disp ('System Gain Margin and its associated frequency are :');Gm                               % 模值稳定裕度Wcg                              % 幅值穿越频率，剪切频率，1/sdisp ('System Phase Margin and its 

6、associated frequency are :');Pm                               % 相位稳定裕度Wcm                             % -π相位穿越频率，1/s%%%% 绘制：Bode图figure(3)Subplot (211);                     % 对数幅值-频率图amp = 20*log10(mag);              %  20*log(mag)，dBsemilogx(wt,amp);tit

7、le ('对数幅值-频率图');xlabel ('频率 / rad');ylabel ('幅值 / dB');grid;subplot (212);                    % 相位-频率图semilogx(wt,phase);title ('相位-频率图');xlabel ('频率/ rad');ylabel ('相位/ degree');grid;   运行该程序可得系统的单位阶跃和脉冲响应曲线如下， 图1 　系统的单位阶跃响应曲线 图2 　系统的脉冲响应曲线 系统的Bode图如下，图

8、3 　系统的Bode图