Matlab TransferFunction -- lism()

close all; clear; clc;

num=[1];
den=[1 4 13 0];
%sys=tf(num,den); %分子分母形式(传递函数形式)
[A,B,C,D] = tf2ss(num,den);
sys = ss(A,B,C,D); %状态空间形式

%[u,t] = gensig('square',4,10,0.1) %方波输入信号
%plot(t,u) %是一个方波
t=0:0.01:10; %采样时间
u=zeros(length(t),1); %对应t的输入信号(每行代表对应t时对应输入的值，每列代表每一个输入)
x0=[1 1 1]; %指定非零初始状态值
lsim(sys,u,t,x0); %线性时不变模型，给定任意输入，得到任意输出