原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
數字信号试卷解答
7 Digital Filter Design
1 四种基本类型的滤波器（低通，高通，带通，带阻）
2 滤波器设计指标的描述
3 IIR数字滤波器设计方法
从模拟滤波器设计IIR数字滤波器：
脉冲响应不变法
双线性变换法
FIR数字滤波器设计方法
线性相位概念
具有线性相位的FIR滤波器的特点：（4种类型！）
设计方法有：窗函数法；频率采样法
1 Determine the period （58，2.21）
（1）
（2） N=18
（3） aperiodic
（4）
2 Determine whether or not the system is linear，casual，stable，and shift-invariable （20，2.26）
（1）
（2）（Page 108 T 2.26 （a））
（3）（Page 108 T 2.26 （e））
3 Determine the inverse DTFT of each the following DTFTs（page120，3.15）
（1）
（2）
Computations （Page 149，3.64）
（1）
（2）
（3）
% compute
x= [1, 3, 0, 0, 1]; y=[1,1,2];
L=length(x) + length(y)-1;
X=fft(x,L); Y=fft(y,L)?;
ZL=ifft(X,Y)；
% OR
ZL=cconv（x，y，L）；
% compute ZC
x= [1, 3, 0, 0, 1]; y= [1,1,2];
L=max（length(x)，length(y)）;
X=fft(x,L); Y=fft(y,L)?;
ZC=ifft(X,Y)；
% OR
ZC=cconv（x，y，L）；
Applications of MATLB functions
（1）page 129
Easily known，a=[1],b=[0.2, 0.06, 0.01, 0.3].
Matlab Programme
a=[1];b=[0.2, 0.06, 0.01, 0.3];
[h,w]=freqz(b,a,512);
subplot (2,2,1)
plot(w/pi,real(h));grid
subplot(2,2,2)
plot(w/pi,imag(h));grid
subplot(2,2,3)
plot(w/pi,abs(h));grid
subplot(2,2,4)
plot(w/pi,angle(h));grid
freqz(b,a,512)
page 85,173
Y = FILTER(B,A,X) filters the data in vector X with the
filter described by vectors A and B to create the filtered
data Y. The filter is a Direct Form II Transposed
implementation of the standard difference equation:
a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb)
a(2)*y(n-1) - ... - a(na+1)*y(n-na)
directly b=[0.8, -0.4, 0.36 ,0, 0.01]
a=[1, 0.26, -0.35 ,0,0,0,0,- 0.6]
2) Z-transform
b=[0.8, -0.4, 0.36 ,0, 0.01]
a=[1, 0.26, -0.35 ,0,0,0,0,- 0.6]
According to page 85：
b=[0.8, -0.4, 0.36 ,0, 0.01]
a=[1, 0.26, -0.35 ,0,0,0,0,- 0.6]
N=41;
x=[1 zeros(1,N-1)];
y=filter(b,a,x);
subplot(2,1,1)
k=0:1:N-1;
stem(k,y)
xlabel(‘Time index n’)；
ylabel(‘Amplitude’)；
x=[ones(1,N)];
y=filter(b,a,x);
subplot(2,1,2)
stem(k,y)
xlabel(‘Time index n’)；
ylabel(‘Amplitude’)；
文档评论（0）