|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为" N: s2 \1 j8 e, g
/ ^' w( B8 ~) }* E
y=filter(b,a,f)
W3 a* Q- x1 J
: P7 h$ ]3 ^; ~- B0 P& i) c& ?0 k
~& _9 Z( ~: Z8 ?其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。
& _5 I6 n9 L$ m$ Y
( k2 w& E! ?' z4 a5 b' i
3 K$ W, B& {# ]
1 c) q, u$ b5 |R=51;
" j# n- j) `. P- Q, Y0 B0 T; {d=rand(1,R)-0.5;
' D8 [( m: S5 {) Z) sk=0:R-1;% F4 O$ Z I2 L- P4 o ]
s=2*k.*(0.9.^k);( m6 L6 ~! ~7 `1 d
f=s+d;
5 _* U( r1 A' i# z+ ffigure(1);
; v& a( I% g8 ?$ s0 j* g' c) U/ M- d$ zstem(k,f,'.');
4 h A- F- o0 v. T# s& iM=5;
r! A* M) W5 Lb=ones(M,1)/M;
8 d+ z! W9 K! oa=1;
" l9 O8 [& F) N( D) zy=filter(b,a,f);
: |$ L, ^6 C8 q* e, @figure(2)9 F: W5 Y& z, ^3 w5 @& p7 `
stem(k,y,'.')* p. Z' `1 k6 r3 c1 D
2 |1 F U. C+ N, E! I: k+ G
) O+ M1 j) y5 A* e- D; s: q" H
1 }' F6 `4 X: t m4 X7 z) s3 ^8 K; |
|
|