|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为
6 m; a9 h7 s# ?* [# |- i2 q$ h* N% A- F( Z& X4 a
y=filter(b,a,f)) n, U: D8 c9 S4 I1 l) \
5 {1 Q9 L" E$ W! F0 L8 y6 [
% _6 K9 l, R" M3 L
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。' {# K" m7 f4 f3 n4 S" |
2 V! w$ `# s9 Q# e2 k5 t/ u% x
8 |- q0 M( Y$ o6 m$ `6 g, O5 n4 @ M% A- Y- f$ J
R=51;; k8 Z& A2 w$ c$ C0 ~
d=rand(1,R)-0.5;4 N; B3 f7 I, J7 s _3 l
k=0:R-1;* J' P5 o+ S5 e% W3 J
s=2*k.*(0.9.^k);
7 ?; }6 f% X7 cf=s+d;
2 y9 S6 H3 g$ _+ M; u( w6 sfigure(1);/ O. p- U& E( |& d
stem(k,f,'.');8 E5 \/ P, e w/ S9 f2 l! {+ Y
M=5;: Q- o; t) \+ d. }( L
b=ones(M,1)/M;
% r$ ?" m! l5 H; K0 F. Wa=1;: ^0 s# t) P7 ] \" a
y=filter(b,a,f);
$ v: B2 g" d. g% d+ @figure(2)! U1 s/ C5 V2 P
stem(k,y,'.')
+ |% Z( s8 s0 e- {
6 G6 ], i, ?1 F! F2 k" x
4 M& ^! R1 D. {, [% ?
7 s# \* n& J6 L6 R
% z1 a6 G1 `. c4 Q. D* [ |
|