ex3_2.m

我们上课的资料

% 关闭所有图形窗口并清除内存变量
close all;
clear;
% 系统开环传递函数初始化
numo=[0 0 0 0 200];
deno=[1 20 140 400 384];
% 求解系统的闭环传递函数
numc=numo;
n=length(deno);
denc=zeros(1,n);
denc=numo+deno;
% 结果显示
disp('System Closed Loop Transfer Function is：')
    numc
    denc
% 系统仿真数据初始化
    t=0:0.05:5;
% 系统输出量的阶跃响应
    y=step(numc,denc,t);
% 系统输出量的脉冲响应
    yy=impulse(numc,denc,t);
% 输出量阶跃响应曲线绘制
figure(1);
plot(t, y);
title('System Step Response');
xlabel('Time—Sec');
ylabel('Response—value');
grid;
% 输出量脉冲响应曲线绘制
figure(2);
plot(t, yy);
title('System impulse Response');
xlabel('Time—sec');
ylabel('Response—value');
grid;
% 求解系统的状态方程模型
[A, B, C, D]=tf2ss(numc, denc);
% 状态方程模型结果显示
disp('SystemState—Space Model is：');
    A
    B
    C
    D
% 系统状态变量的阶跃响应
[ys,x,t]=step(A, B, C, D);
% 系统状态变量的脉冲响应
[yys,xx,t]=impulse(A, B, C, D);
% 状态变量阶跃响应曲线绘制
figure(3);
plot(t,x(:,1),t,x(:,2),'+',t,x(:,3),'.-',t,x(:,4),'.')
title('System State—Variables Step Response');
xlabel('Time—sec');
ylabel('Response—value');
grid;
% 状态变量脉冲响应曲线绘制
figure(4);
plot(t,xx(:,1),t,xx(:,2),'+',t,xx(:,3),'.-',t,xx(:,4),'.')
title('System State—Variables impulse Response');
xlabel('Time—sec');
ylabel('Response—value');
grid;