|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
在零状态时,MATLAB工具箱提供了一个filter函数,计算由差分方程描述的系统响应,其调用格式为* x: j5 K2 n6 s7 T, l: ^) a# p
- T+ v F3 s1 T* O. J
y=filter(b,a,f)
" W9 V, d8 s% Y" H8 r7 W# F0 e
+ z& A" ?2 Y1 a( x* M3 \* j0 e7 u7 @' ?, T5 i3 B, ^0 U8 x; G
其中b=[b0,b1,b2,…, bm],a=[a0,a1,a2,…, an]分别是差分方程左右的系数向量,f表示输入序列,y表示系统的零状态响应。注意输出和输入序列的长度相同。" j5 | a6 l3 y# p7 _. \! @2 W
) {4 b+ P0 U& o* [9 f/ _
/ c q$ c0 Z0 T6 C7 [# Z
/ s' H7 x4 [+ H: `+ hR=51;
+ I/ G; J/ x2 H: \% }' x" wd=rand(1,R)-0.5; S) w( ^9 Z$ }" w8 U- G+ z
k=0:R-1;( {1 v* ^8 M2 J! R6 f
s=2*k.*(0.9.^k);
/ K8 B: k% S; n3 n' I0 ]( b$ {f=s+d;
5 @* }0 W, K1 }+ o+ x% Zfigure(1);- S: G3 ~. B( T$ d* H/ y5 i9 b/ K
stem(k,f,'.');
5 I T, {0 A6 `( c% n( FM=5;
- u- H8 Q% D; p: S6 m) qb=ones(M,1)/M;
, r' S4 g( j8 p( sa=1;
4 Y' P/ D! S7 w0 u9 g( [9 ty=filter(b,a,f);
8 }* f! v% P+ ~: `% ifigure(2)+ y; p1 L' c8 o$ X- P. v
stem(k,y,'.'): M( X8 M5 A1 n( D# N
% p9 g& e: G, {3 h
& M! a& C9 y5 m {7 f: P$ ?. i' V1 L
+ u! T; } A; A$ L |
|