連續時間系統的時間分析

系統函數

clear all;
close all;
a=[1,7,16,12];
b=[1];
sys=tf(b,a);

直接表示法

H=tf([1],[1 2 1]);

p=[1 7 16 12];
a=roots(p)

衝激響應

impulse(sys);%其中sys可以是利用命令tf。zpk或者ss建立的系統函數
impulse(sys,t);%計算並畫出系統在向量t定義的時間內的衝激函數

階躍函數

step(sys);%其中sys可以是利用命令tf,zpk,或ss建立的系統函數
step(sys,t);%計算並畫出系統在向量定義的時間內的階躍函數

連續系統任意輸入下的時間響應
lsim函數
功能：計算並畫出系統在任意輸入下的零狀態響應
調用格式：
lism(sys,x,t)其中sys可以是利用命令tf,zpk或者ss建立的系統函數。x是系統的輸入，t定義的是時間範圍；
lism(sys,x,t,zi):計算系統在任意輸入的零狀態下的全響應，sys必須是狀態空間形式的系統函數，zi是系統的初始狀態。

b=[3 9];
a=[1 6 8];
sys=tf(b,a);
t=0:0.1:10;
y=impluse(sys,t);
plot(t,y);
xlabel('時間')；
ylabel('y(t)');
title('單位衝激響應');
a=[0 3 9];
b=[1 6 8];
sys=tf(a,b);
t=0:0.1:10;
y=step(sys,t);
plot(t,y);
xlabel('時間(t)');
ylabel('y(t)');
title('單位階躍響應');

a=[0 0 1];
b=[1 0 1];
sys=tf(a,b);
t=0:0.1:10;
x=cos(t);
y=lsim(sys,x,t);
plot(t,y);
xlabel('時間(t)');
ylabel('y(t)');
title('零狀態響應')；
b=[1];
a=[1 0 1];
[A B C D]=tf2ss(a,b);
sys=ss(A,B,C,D);
t=0:0.1:10;
x=cos(t);
zi=[-1 0];
y=lsim(sys,x,t,zi);
plot(t,y);
xlabel('時間(t)');
ylabel('y(t)');
title('全響應')；

離散時間系統