ex3_2.asv

我们上课的资料

% 关闭所有图形窗口并清除内存变量
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:3;
% 系统输出量的阶跃响应
    y=step(numc,denc,t);
% 系统输出量的脉冲响应
    yy=impulse(numc,denc,t);
% 输出量阶跃响应曲线绘制
plot(t, y);
title('System Step Response');
xlabel('Time—Sec');
ylabel('Response—value');
grid;
% 输出量脉冲响应曲线绘制
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]=step(A, B, C, D, 1, t);
% 系统状态变量的脉冲响应
[yys,xx]=impulse(A, B, C, D, 1, t);
% 状态变量阶跃响应曲线绘制
plot(t,x(:,1),t,x(:,2),'+',t,x(:,3)，，．—，，t，x(：，4)，，．r)
    title(，SystemState—VariablesStepResponse')；
    xlabel(，Time—secl)；    ?
    ylabel(，Response—valuet)；
    grid；
    ％状态变量脉冲响应曲线绘制
    吵t(t，xx(：，1)，t，xx(：，2)，，+，，t，xx(：，3)，，．—，，t，xx(：，4)，，．，)
    title(，SystemState—VariableslmpulseResponse')；    ‘
    xlabel(，Time—seef)；
    ylabel(，Response—value')；
    grid；