自动化_数字信号处理实验5答案
数字信号处理教程课后习题及答案
分析：已知边界条件，如果没有限定序列类型（例如因果序列、反因果序列等）， 则递推求解必须向两个方向进行（n ≥ 0 及 n < 0）。
解 : (1) y1 (0) = 0 时， (a) 设 x1 (n) = δ (n) ，
按 y1 (n) = ay1 (n − 1) + x1 (n) i) 向 n > 0 处递推，
10
T [ax1(n)+ bx2 (n)] =
n
∑
[ax1
(n
)
+
bx2
(n
)]
m = −∞
T[ax1(n) + bx2(n)] = ay1(n) + by2(n)
∴ 系统是线性系统
解：(2) y(n) =
[x(n )] 2
y1(n)
= T [x1(n)] = [x1(n)] 2
y2 (n) = T [x2 (n)] = [x2 (n)] 2
β α
n +1
β α β =
n +1− N −n0
N−
N
α −β
y(n) = Nα n−n0 ,
(α = β )
, (α ≠ β )
如此题所示，因而要分段求解。
2 .已知线性移不变系统的输入为 x( n ) ,系统的单位抽样响应
为 h( n ) ,试求系统的输出 y( n ) ,并画图。
(1)x(n) = δ (n)
当n ≤ −1时 当n > −1时
∑ y(n) = n a −m = a −n
m=−∞
1− a
∑ y(n) =
−1
a−m =
第五章 数字信号处理课后答案刘顺兰版
H α ( s ) = H αN (
ω c = 2πf c T = 2π × 400 HZ / 6000 HZ =
Ωc = 2 ωc 2 π 2 tg = tg ( ) = 0.2 × T T 2 T 15
2π 15
s=
2 1 − z −1 , T 1 + z −1
s = Ωc
1 − z −1 −1 π tg ( ) 1 + z 15 1 1
=
1 + 3z −1 + 3z −2 + z −3 0.005376(1 + 3z −1 + 3z −2 + z −3 ) = 186 − 412 z −1 + 318 z − 2 − 84 z −3 1 − 2.215 z −1 + 1.71z − 2 − 0.4516 z −3
5.24 用双线性变换设计一个三阶巴特沃思数字高通滤波器，采样频率为 f s = 6 KHZ ，截止 频率为 f c = 1.5 KHZ （不计 3KHZ 以上的频率分量） 。 解法 1：三阶巴特沃思低通模拟滤波器的原型函数：
按照冲激不变条件，可以写出
因此系统函数为
H ( z ) = ∑ h(n) z − n
n =0
∞
1 1 2 2 = + 1 − e − aT e − jbT z −1 1 − e −aT e jbT z −1 = 1 − (e − aT cos bT ) z −1 (1 − e − aT e − jbT z −1 )(1 − e −aT e jbT z −1 )
所以
ω1 + ω 2
H BP ( z ) = H αN ( s )
s=
1 1+ z − 2 3 1− z − 2
数字信号处理第5章答案史林赵树杰编著
数字信号处理第5章答案史林赵树杰编著第五章练习题答案%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%5.8 已知复序列()()()f n x n jy n =+的8点DFT 为()[()](07)F k DFTf n k =≤≤，其值为(0)13,(1)24,(2)37,(3)45,(4)25,(5)12,(6)48,(7)6,F j F j F j F j F j F j F j F j =-=-+=+=--=+=--=-=不计算()F k 的离散傅里叶逆变换（IFFT ），试求实序列()x n 和()y n 的8点DFT ()X k 和()Y k 。
解：利用DFT 的共轭对称性()()()f n x n jy n =+[]()()()()F k DFT f n X k jY k ==+[]Re ()()()f n Fep k x n ??[]Im ()()()j f n Fop k y n ??所以[][]*()()R e ()()1(())(())()2N N N X k D FT x n D FT f n Fep k F k F N k R k ??====+-?[][]*1()()Im ()()1(())(())()2N N N Y k D FT y n D FT f n Fep k j F k F N k R k j ??====--?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%5.9 设()x n 和()y n 是长度为N 的两个实序列。
已知()[()](01)X k DFT x n k N =≤≤-，()[()](01)Y k DFT y n k N =≤≤-。
现在希望根据()X k 和()Y k 求()x n 和()y n ，为了提高运算效率，试设计一种算法，用一次N 点IFFT 来完成。
数字信号处理,第5章课后习题答案
第五章习题与上机题5.1 已知序列12()(),0 1 , ()()()nx n a u n a x n u n u n N =<<=--，分别求它们的自相关函数，并证明二者都是偶对称的实序列。
解：111()()()()()nn mx n n r m x n x n m a u n au n m ∞∞-=-∞=-∞=-=-∑∑当0m ≥时，122()1mmnx n ma r m aaa∞-===-∑ 当0m <时，122()1m mnx n a r m aaa -∞-===-∑ 所以，12()1mx ar m a =-2 ()()()()N x n u n u n N R n =--=22210121()()()()()1,0 =1,00, =()(1)x NN n n N mn N n m N r m x n x n m Rn R n m N m N m N m m Nm N m R m N ∞∞=-∞=-∞--=-=-=-=-⎧=--<<⎪⎪⎪⎪=-≤<⎨⎪⎪⎪⎪⎩-+-∑∑∑∑其他从1()x r m 和2()x r m 的表达式可以看出二者都是偶对称的实序列。
5.2 设()e()nTx n u n -=，T 为采样间隔。
求()x n 的自相关函数()x r m 。
解：解：()()()()e()e ()nTn m T x n n r m x n x n m u n u n m ∞∞---=-∞=-∞=-=-∑∑用5.1题计算1()x r m 的相同方法可得2e()1e m Tx Tr m --=-5.3 已知12()sin(2)sin(2)s s x n A f nT B f nT ππ=+，其中12,,,A B f f 均为常数。
求()x n 的自相关函数()x r m 。
解：解：()x n 可表为)()()(n v n u n x +=的形式，其中)2sin()(11s nT f A n u π=，=)(n v 22sin(2)s A f nT π，)(),(n v n u 的周期分别为 s T f N 111=，sT f N 221=，()x n 的周期N 则是21,N N 的最小公倍数。
数字信号处理实验答案
实验一熟悉Matlab环境一、实验目的1.熟悉MATLAB的主要操作命令。
2.学会简单的矩阵输入和数据读写。
3.掌握简单的绘图命令。
4.用MATLAB编程并学会创建函数。
5.观察离散系统的频率响应。
二、实验内容认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。
在熟悉了MATLAB基本命令的基础上，完成以下实验。
上机实验内容：（1）数组的加、减、乘、除和乘方运算。
输入A=[1 2 3 4]，B=[3 4 5 6]，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B并用stem语句画出A、B、C、D、E、F、G。
clear all;a=[1 2 3 4];b=[3 4 5 6];c=a+b;d=a-b;e=a.*b;f=a./b;g=a.^b;n=1:4;subplot(4,2,1);stem(n,a);xlabel('n');xlim([0 5]);ylabel('A');subplot(4,2,2);stem(n,b);xlabel('n');xlim([0 5]);ylabel('B');subplot(4,2,3);stem(n,c);xlabel('n');xlim([0 5]);ylabel('C');subplot(4,2,4);stem(n,d);xlabel('n');xlim([0 5]);ylabel('D');subplot(4,2,5);stem(n,e);xlabel('n');xlim([0 5]);ylabel('E');subplot(4,2,6);stem(n,f);xlabel('n');xlim([0 5]);ylabel('F');subplot(4,2,7);stem(n,g);xlabel('n');xlim([0 5]);ylabel('G');（2）用MATLAB实现下列序列：a) x(n)=0.8n0≤n≤15b) x(n)=e(0.2+3j)n0≤n≤15c) x(n)=3cos(0.125πn+0.2π)+2sin(0.25πn+0.1π) 0≤n≤15d) 将c)中的x(n)扩展为以16为周期的函数x16(n)=x(n+16),绘出四个周期。
数字信号处理实验课后答案
数字信号处理实验课后答案数字信号处理实验课后答案【篇一：数字信号处理第三版课后实验程序(高西全)】txt>close all;clear all%======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== a=[1,-0.9];b=[0.05,0.05]; %系统差分方程系数向量b和ax1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=r8(n)x2n=ones(1,128); %产生信号x2(n)=u(n)hn=impz(b,a,58); %求系统单位脉冲响应h(n)subplot(2,2,1);y=h(n);tstem(hn,y); %调用函数tstem绘图title((a) 系统单位脉冲响应h(n));box ony1n=filter(b,a,x1n); %求系统对x1(n)的响应y1(n)subplot(2,2,2);y=y1(n);tstem(y1n,y);title((b) 系统对r8(n)的响应y1(n));box ony2n=filter(b,a,x2n); %求系统对x2(n)的响应y2(n)subplot(2,2,4);y=y2(n);tstem(y2n,y);title((c) 系统对u(n)的响应y2(n));box on%===内容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=r8(n)h1n=[ones(1,10) zeros(1,10)];h2n=[1 2.5 2.5 1 zeros(1,10)];y21n=conv(h1n,x1n);y22n=conv(h2n,x1n);figure(2)subplot(2,2,1);y=h1(n);tstem(h1n,y); %调用函数tstem绘图title((d) 系统单位脉冲响应h1(n));box onsubplot(2,2,2);y=y21(n);tstem(y21n,y);title((e) h1(n)与r8(n)的卷积y21(n));box onsubplot(2,2,3);y=h2(n);tstem(h2n,y); %调用函数tstem绘图title((f) 系统单位脉冲响应h2(n));box onsubplot(2,2,4);y=y22(n);tstem(y22n,y);title((g) h2(n)与r8(n)的卷积y22(n));box on%=========内容3：谐振器分析======================== un=ones(1,256); %产生信号u(n)n=0:255;xsin=sin(0.014*n)+sin(0.4*n); %产生正弦信号a=[1,-1.8237,0.9801];b=[1/100.49,0,-1/100.49]; %系统差分方程系数向量b和a y31n=filter(b,a,un); %谐振器对u(n)的响应y31(n) y32n=filter(b,a,xsin); %谐振器对u(n)的响应y31(n)figure(3)subplot(2,1,1);y=y31(n);tstem(y31n,y);title((h) 谐振器对u(n)的响应y31(n));box onsubplot(2,1,2);y=y32(n);tstem(y32n,y);title((i) 谐振器对正弦信号的响应y32(n));box on10.2.2 实验程序清单1 时域采样理论的验证程序清单% 时域采样理论验证程序exp2a.mtp=64/1000; %观察时间tp=64微秒%产生m长采样序列x(n)% fs=1000;t=1/fs;fs=1000;t=1/fs;m=tp*fs;n=0:m-1;a=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=a*exp(-alph*n*t).*sin(omega*n*t);xk=t*fft(xnt,m); %m点fft[xnt)]yn=xa(nt);subplot(3,2,1);tstem(xnt,yn); %调用自编绘图函数tstem绘制序列图box on;title((a) fs=1000hz);k=0:m-1;fk=k/tp;subplot(3,2,2);plot(fk,abs(xk));title((a) t*ft[xa(nt)],fs=1000hz);xlabel(f(hz));ylabel(幅度);axis([0,fs,0,1.2*max(abs(xk))])%===================================== ============ % fs=300hz和fs=200hz的程序与上面fs=1000hz完全相同。
数字信号处理习题及解答
只有在如上周期延拓序列中无混叠的点上， 才满足f(n)=fl(n),所以 f(n)=fl(n)=x(n)*y(n) 7≤n≤19
令
数字信号处理习题及解答
第二章Z变换及离散时间系统分析
3 解答
n≥0时， 因为c内无极点，x(n)=0; n≤－1时， c内有极点0 ， 但z=0是一个n阶极点， 改为求
圆外极点留数， 圆外极点有z1=0.5, z2=2， 那么
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 3 解答 (2) 收敛域0.5<|z|<2:
数字信号处理习题及解答
第三章信号的傅里叶变换 1 解答
(1) (2) (3)
数字信号处理习题及解答
第三章信号的傅里叶变换 2 试求如下序列的傅里叶变换：
(1) x1(n)=δ(n－3)
（2）
数字信号处理习题及解答
第三章信号的傅里叶变换 2 解答
(1) (2)
数字信号处理习题及解答
第三章信号的傅里叶变换
第一章离散时间信号与离散时间系统
4 解答
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 1
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 1 解答
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 1 解答
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 2
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 2 解答
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 2 解答
数字信号处理习题及解答
第二章Z变换及离散时间系统分析 3 已知
求出对应X(z)的各种可能的序列表达式。
数字信号处理课后习题答案(全)1-7章
1
4
(2m 5) (n m) 6 (n m)
m4
m0
第 1 章 时域离散信号和时域离散系统
（3） x1(n)的波形是x(n)的波形右移2位， 再乘以2， 画出图形如题2解图 （二）所示。
(4) x2(n)的波形是x(n)的波形左移2位， 再乘以2， 画出图形如题2解图（三） 所示。
因此系统是非线性系统。
第 1 章 时域离散信号和时域离散系统
（6） y(n)=x(n2)
令输入为
输出为
x(n－n0)
y′(n)=x((n－n0)2) y(n－n0)=x((n－n0)2)=y′(n) 故系统是非时变系统。 由于
T［ax1(n)+bx2(n)］=ax1(n2)+bx2(n2) =aT［x1(n)］+bT［x2(n)］
0≤m≤3
－4≤m≤n
非零区间如下:
第 1 章求解： ① n<0时， y(n)=0
② 0≤n≤3时， y(n)= ③ 4≤n≤7时， y(n)= ④ n>7时， y(n)=0
1=n+1
n
1=8－m n0
3
mn4
第 1 章 时域离散信号和时域离散系统
x(m)h(n－m)
m
第 1 章 时域离散信号和时域离散系统
题7图
第 1 章 时域离散信号和时域离散系统
y(n)={－2,－1,－0.5, 2, 1, 4.5, 2, 1; n=－2, －1, 0, 1, 2, 3, 4, 5}
第 1 章 时域离散信号和时域离散系统
解法（二） 采用解析法。 按照题7图写出x(n)和h(n)的表达式分别为
数字信号处理实验答案
第十章上机实验数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。
上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。
本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。
实验一系统响应及系统稳定性。
实验二时域采样与频域采样。
实验三用FFT对信号作频谱分析。
实验四IIR数字滤波器设计及软件实现。
实验五FIR数字滤波器设计与软件实现实验六应用实验——数字信号处理在双音多频拨号系统中的应用任课教师根据教学进度，安排学生上机进行实验。
建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。
学习完第六章进行；实验五在学习完第七章后进行。
实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。
10.1 实验一: 系统响应及系统稳定性1.实验目的（1）掌握求系统响应的方法。
（2）掌握时域离散系统的时域特性。
（3）分析、观察及检验系统的稳定性。
2.实验原理与方法在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。
已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。
在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB语言的工具箱函数filter函数。
也可以用MA TLAB语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。
系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。
重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。
系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。
或者系统的单位脉冲响应满足绝对可和的条件。
系统的稳定性由其差分方程的系数决定。
数字信号处理课后答案_史林版_科学出版社
第一章 作业题 答案############################################################################### 1.2一个采样周期为T 的采样器，开关导通时间为()0T ττ<<，若采样器的输入信号为()a x t ，求采样器的输出信号()()()a a x t x t p t ∧=的频谱结构。
式中()()01,()0,n p t r t n t r t ττ∞=-∞=-≤≤⎧=⎨⎩∑其他解：实际的采样脉冲信号为：()()n p t r t n τ∞=-∞=-∑其傅里叶级数表达式为：()000()jk tn p t Sa k T eTωωτω∞=-∞=∑采样后的信号可以表示为：()()()ˆa a xt x t p t δ= 因此，对采样后的信号频谱有如下推导：()()()()()()()()()()()()()0000000000000ˆˆsin 1j t a a jk t j t a n jk t j t a k j k ta k ak a k X j x t e dtx t Sa k T e e dtTSa k T x t e e dtTSa k T x t edtTSa k T X j jk Tk T X j jk T kωωωωωωωωτωωτωωτωωτωωωωωω∞--∞∞∞--∞=-∞∞∞--∞=-∞∞∞---∞=-∞∞=-∞∞=-∞Ω=====-=-⎰∑⎰∑⎰∑⎰∑∑%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 1.5有一个理想采样系统，对连续时间信号()a x t 进行等间隔T 采样，采样频率8s πΩ=rad/s ，采样后所得采样信号()a x t ∧经理想低通滤波器()G j Ω进行恢复，已知()41/4,,4G j ππ⎧Ω≤⎪Ω=⎨Ω>⎪⎩今有两个输入信号12()cos(2)()cos(5)a a x t t x t t ππ==和，对应的输出信号分别为12()()a a y t y t 和，如题1.5图所示，问12()()a a y t y t 、有没有失真，为什么？题1.5图 理想采样系统与恢复理想低通滤波器解：因为是理想采样系统，因此采样后的信号频谱可以表示为：()()1ˆa a s k X j X j jk T ∞=-∞Ω=Ω-Ω∑8s πΩ=，12πΩ=，25πΩ=，折叠频率为2s Ω，而滤波器对4πΩ≤的信号通过，因此有如下图：结论：1）1()a y t 不失真、2()a y t 失真。
数字信号处理课后习题答案(全)1-7章
x(n)=－δ(n+2)+δ(n－1)+2δ(n－3)
h(n)=2δ(n)+δ(n－1)+ δ(n－2)
由于
x(n)*δ(n)=x(n)
1
x(n)*Aδ(n－k)=Ax(n－k)
2
故
第 1 章 时域离散信号和时域离散系统
y(n)=x(n)*h(n)
=x(n)*［2δ(n)+δ(n－1)+ δ(n－2) 1 2
（5） 系统是因果系统， 因为系统的输出不取决于x(n)的未来值。 如果
|x(n)|≤M, 则|y(n)|=|ex(n)|≤e|x(n)|≤eM，
7． 设线性时不变系统的单位脉冲响应h(n)和输入序列x(n)如题7图所示，
要求画出y(n)输出的波形。
解： 解法（一）采用列表法。
y(n)=x(n)*h(n)=
0≤m≤3
－4≤m≤n
非零区间如下:
第 1 章 时域离散信号和时域离散系统
根据非零区间， 将n分成四种情况求解： ① n<0时， y(n)=0
② 0≤n≤3时， y(n)= ③ 4≤n≤7时， y(n)= ④ n>7时， y(n)=0
1=n+1
n
1=8－m n0
3
mn4
第 1 章 时域离散信号和时域离散系统
第 1 章 时域离散信号和时域离散系统
(3) 这是一个延时器， 延时器是线性非时变系统， 下面证明。 令输入为
输出为
x(n－n1)
y′(n)=x(n－n1－n0) y(n－n1)=x(n－n1－n0)=y′(n) 故延时器是非时变系统。 由于
T［ax1(n)+bx2(n)］=ax1(n－n0)+bx2(n－n0) =aT［x1(n)］+bT［x2(n)］
数字信号处理作业 第五章 参考答案
(1) 由极点构成 H a ( s ) 的分母多项式，分子为分母多项式的常数。 (2) H a ( s ) 展成部分分式。 (3) 据有理分式变换得到对应的 H ( z ) 各分式，整理得到最后的 H ( z ) 。 22、 取 T=1， 预畸， 由已知列出对模拟滤波器的衰减要求， 解出 N=6.04， 取 N=7， 得到
−0.5
Z −1
−1
0.9
−0.81
4、 H ( z ) = −4.9383 +
2.1572 4.7811 − 1.5959 z −1 + 1 + 0.5 z −1 1 − 0.9 z −1 + 0.81z −2
−4.9383
x ( n) y ( n)
2.1572 −0.5
Z −1
4.7811
Z
0.9 −0.81
= H 2 ( z)
α 02 + α12 z -1 -3.1986 + 0.2591z -1 = 1 +z 2 1 + 1.618 z - 4π 2 2 1 + r z 1 - 2rz -cos 5
频率取样型实现流程图：
−10.125
Z −1
18.3236
x ( n)
Z −1
x ( n)
Z −1
Z −1
+
Z −1
− 7 4
+
Z −1
− 69 8
+
y ( n) 4） 频率取样型：取 r=1，N=5，得到 DFT{h(n)}为：
{-10.1250 9.1618 + 6.6564i -1.5993 - 4.9221i -1.5993 + 4.9221i 9.1618 - 6.6564i}
数字信号处理_吴镇扬_第二版_第五章习题答案
5.7 （1）由于h2(n)是h1(n)圆周移位的序列，根据DFT的 2π 性质有： −j 4k − jπ k
H 2 (k ) = e
8
H 1 (k ) = e
H 1 (k )
~ ~ H1 ( k ) = H 2 ( k ) 成立 所以
(2)由于h1 (n ) 和h2 (n ) 均为偶对称序列，以其构成的低通滤波器
(3)若采用海明窗设计，则
⎡ ⎛ 2πn ⎞⎤ wHam ( n) = ⎢0.54 − 0.46 cos ⎜ ⎟ ⎥ RN ( n ) ⎝ N − 1 ⎠⎦ ⎣ 2 h( n) = sin[(n − α )ωc ]cos[(n − α )ω0 ]wHam (n) N 为奇数时， (n − α )π
h( n N 为偶数时， ) =
0 −ωc
e − jωα e jω nd ω
可见h(n)关于(N-1)/2偶对称，即 h( n) = h( N − 1 − n)
(1)当 N 为奇数时，为第一类滤波器。 (2)当N为偶数时，为第二类滤波器
⎧hd ( n) h( n) = hd ( n) ⋅ R(n ) = ⎨ ⎩0 0 ≤ n ≤ N −1
解：由经验公式可知若 不小于 At 40dB , 则
β = 0.5842 At - 21)0.4 + 0.07886(At - 21) ≈ 3.3953 （ At − 8 40 − 8 N= = ≈ 22.28 2.286∆ω 2.286× 0.2π ωc + ωr ωc′ = = 0.2π 2 ′ ⎧ωc ′ ⎪ π Sa[ωc (n − α )] n ≠ α ′ 1 ωc − jωα jωn ⎪ hd (n) = ∫ ′ e e dω = ⎨ ′ 2π −ωc ωc ⎪ n =α ⎪ ⎩ π
数字信号处理课后习题答案(全)1-7章
最后结果为 0
n<0或n>7
y(n)= n+1 0≤n≤3 8－n 4≤n≤7
y(n)的波形如题8解图（一）所示。 (2) y(n) =2R4(n)*［δ(n)－δ(n－2)］=2R4(n)－2R4(n－2)
=2［δ(n)+δ(n－1)－δ(n+4)－δ(n+5) y(n)的波形如题8解图（二）所示
因此系统是非线性系统。
第 1 章 时域离散信号和时域离散系统
（6） y(n)=x(n2)
令输入为
输出为
x(n－n0)
y′(n)=x((n－n0)2) y(n－n0)=x((n－n0)2)=y′(n) 故系统是非时变系统。 由于
T［ax1(n)+bx2(n)］=ax1(n2)+bx2(n2) =aT［x1(n)］+bT［x2(n)］
因此系统是非时变系统。
第 1 章 时域离散信号和时域离散系统
（5） y(n)=x2(n)
令输入为
输出为
x(n－n0)
y′(n)=x2(n－n0) y(n－n0)=x2(n－n0)=y′(n) 故系统是非时变系统。 由于
T［ax1(n)+bx2(n)］=［ax1(n)+bx2(n)］2 ≠aT［x1(n)］+bT［x2(n) =ax21(n)+bx22(n)
第 1 章 时域离散信号和时域离散系统
题4解图（一）
第 1 章 时域离散信号和时域离散系统
题4解图（二）
第 1 章 时域离散信号和时域离散系统
题4解图（三）
第 1 章 时域离散信号和时域离散系统
(4) 很容易证明： x(n)=x1(n)=xe(n)+xo(n)
数字信号处理第5章答案
第５章 无限脉冲响应(IIR)数字滤波器的设
计
特别是以理想滤波器特性作为Hd(ejω)时， 为了使ε2最小，
优化过程尽可能逼近Hd(ejω)的间断特性（即使过渡带最窄）， 而使通带出现较大过冲、 阻带最小衰减过小， 不能满足工
H(ejω)=|H(ejω)|ejθ(ω) 其中， |H(ejω)|称为幅频特性函数， θ(ω)称为相频特性函数。
常用的典型滤波器|H(ejω)|是归一化的， 即|H(ejω)|max=1， 下 的讨论一般就是针对归一化情况的。 对IIR数字滤波器， 通
常用幅频响应函数|H(ejω)|来描述设计指标， 而对线性相位特 性的滤波器， 一般用FIR数字滤波器设计实现。
计
图5.1.6
第５章 无限脉冲响应(IIR)数字滤波器的设
计
5.1.4 IIR-DF的直接设计法
所谓直接设计法， 就是直接在数字域设计IIR[CD*2]DF 的方法。 相对而言， 因为从AF入手设计DF是先设计相应的 AF， 然后再通过s-z平面映射， 将Ha(s)转换成H(z)， 所以 这属于间接设计法。 该设计法只能设计与几种典型AF相对 应的幅频特性的DF。 而需要设计任意形状幅频特性的DF时， 只能用直接设计法。 直接设计法一般都要借助于计算机进行 设计， 即计算机辅助设计（CAD）。 现在已有多种DF优化 设计程序。 优化准则不同， 所设计的滤波器特点亦不同。所 以最主要的是建立优化设计的概念， 了解各种优化准则的 特点， 并根据设计要求， 选择合适的优化程序设计DF。
≤≤
(5.1.1)
≤
(5.1.2)
数字信号处理课后答案第五章报告.ppt
1
1、用直接I型及典范结构实现以下系统函数：
H
z
=
3 2
4.2 0.6
z z
1 1
0.8 z 2 0.4 z 2
解：根据IIR滤波器的系统函数标准式
M
H
bm zm
z =
m0 N
1 an zn
Y z X z
n 1
将系统函数整理为：
H
z
=
1.5 2.1z1 0.4z2 1 0.3z1 0.2z2
H
z
=
4
z
z 1z 0.5 z2
2 1.4z 1 0.9z 0.8
试问一共能构成几种级联型网络。
解：H
z
A
k
1 1k z1 2k z2 1 1k z1 2k z2
4 1 z1 11.4z1 z2
1 0.5z1 1 0.9z1 0.8z2
4
考虑分子分母的组合及级联的次序，共有以下 四种级联型网络：
=
5
2z3 3z 1 z1
6
抽样点数 N 6，修正半径 r 0.9。
解：由N = 6，得频率抽样型结构：
H z= 1 6
1 r6z6
H0
z
H3
z
2 k 1
H
k
z
又
5 3z3 1 z3
H z=
1 z1
1 z1
1 z1 z2
5 3z3 1 z1 z2
5
h2 1
h1 h3 3 0.6
5
即 hn是偶对称，对称中心在n N 1 2处，
N为奇数 N 5 。
2
得线性相位结构：
17
1.5 2.1z1 0.4z2 1 0.3z1 0.2z2
数字信号处理课后习题答案(吴镇扬)(精编文档).doc
【最新整理，下载后即可编辑】习题一 (离散信号与系统)1.1周期序列，最小周期长度为5。
1.2 (1) 周期序列，最小周期长度为14。
(2) 周期序列，最小周期长度为56。
1.5()()()()()()()11s a s s s a n s s a s n X j x t p t X j ΩP j Ω2n τn τj sin j Ωjn e X 2n π2n n τj Sa X j jn e 2T 2πττ∞=-∞∞=-∞Ω==*⎡⎤⎣⎦ΩΩ⎛⎫-=-Ω ⎪⎝⎭ΩΩ⎛⎫-=Ω-Ω ⎪⎝⎭∑∑ 1.6 (1) )(ωj e kX (2) )(0ωωj n j e X e (3) )(21)(2122ωωj j e X e X -+ (4) )(2ωj e X1.7 (1) 0n z -(2) 5.0||,5.0111>--z z (3) 5.0||,5.0111<--z z (4)0||,5.01)5.0(11101>----z z z1.8 (1) 0,)11()(211>--=---z zz z z X N(2) a z az az z X >-=--,)1()(211 (3) a z az z a az z X >-+=---,)1()(311211.91.10 (1))1(2)(1----+n u n u n (2))1(24)()5.0(6--⋅--n u n u n n (3))()sin sin cos 1(cos 000n u n n ωωωω++(4) )()()(1n u a a a n a n ---+-δ 1.11(1))(1z c X - (2) )(2z X (3))()1(21z X z -+ (4)-+<<x x R z R z X /1/1),/1(1.12 (1) 1,11<-ab ab(2) 1 (3)00n a n1.13 (1) 该系统不是线性系统；该系统是时不变系统。
数字信号处理 Chapter05答案
2011/4/10
3
Notice: Skipped Sections
5.7 5.8 (all) 5.9 (all) 5.10 5.11
2011/4/10
4
5.1
Introduction
2011/4/10
5
5.1 Introduction
Digital processing of a real-world continuous-time signal involves the following basic steps:
The multiplication operation yields an impulse train
Note: Analog Filters are very important for digital systems. Since both the anti-aliasing filter and the reconstruction filter are analog lowpass filters, we review first the theory behind the design of such filters Also, the most widely used IIR digital filter design method is based on the conversion of an analog lowpass prototype
g[n ] = ga (nT ), -¥ <n < ¥
(5.1)
With T being the sampling period. The reciprocal of T is called the sampling frequency , i.e., 1⁄ . It is known that the frequency-domain of the analog signal is given by its FT:
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验五 IIR滤波器设计一、实验目的1. 掌握冲激响应不变法和双线性变换法设计IIR数字滤波器的原理和方法；2. 观察双线性变换法和冲激响应不变法设计的滤波器的频域特性，了解双线性变换法和冲激响应不变法的特点和区别。
二、实验仪器1. 计算机2. MATLAB软件三、预习要求理解冲激响应不变法，双线性变换法，掌握常用 matlab 程序。
四、实验原理1. IIR数字滤波器一般为线性移不变的因果离散系统，N阶IIR滤波器的传递函数为：，式中系数至少有一个不为0。
IIR数字滤波器的设计通常利用模拟滤波器作为原型滤波器，直接由模拟滤波器的频率响经过冲激响应不变法或双线性变换法转换成IIR数字滤波器。
2. 常用模拟原型滤波器（1）巴特沃斯滤波器巴特沃斯滤波器通带和阻带都单调衰减，所有结束的幅度函数-3dB点为同一点。
其幅度响应为：，式中N为滤波器阶数，通带截止频率。
（2）切比雪夫Ⅰ型滤波器切比雪夫Ⅰ型滤波器在通带呈现等波纹特性，阻带单调衰减，其幅度响应为：，式中N为滤波器阶数；，表示通带波纹大小，越大，波纹越大；为截止频率，不一定为3dB带宽；为N阶Chebyshev多项式。
（3）切比雪夫Ⅱ型滤波器切比雪夫Ⅱ型滤波器在阻带呈现等波纹特性，通带单调衰减，其幅度响应为：，式中是阻带衰减达到一定数值时的最低频率。
（4）椭圆滤波器椭圆滤波器在通带和阻带都呈现波纹特性，在带内均匀波动，具有最快的滚降。
，式中为椭圆函数。
3. 冲激响应不变法所谓冲激响应不变法就是使数字滤波器的单位冲激响应序列等于模拟滤波器的单位冲激响应和的采样值，即：，其中，T为采样周期。
令为模拟系统传递函数，且，则冲激响应不变法得到的数字滤波器传递函数。
冲激响应不变法的特点是：（1）时域逼近特性良好；（2）模拟频率和数字频率呈线性关系。
（3）存在频率混叠效应，故只适用于限带的模拟滤波器。
在 MATLAB 中，可用函数 impinvar 实现从模拟滤波器到数字滤波器的冲激响应不变映射，调用格式为：[bz,az]=impinvar(b,a,fs)[bz,az]=impinvar(b,a)其中，b、a 分别为模拟滤波器的分子和分母多项式系数向量；fs为采样频率（Hz），缺省值 fs=1Hz；bz、az分别为数字滤波器分子和分母多项式系数向量。
举例：若Ha(s)=，则b=[1,1];a=[1,5,6]。
4. 双线性变换法由于平面和平面的单值双线性映射关系为，其中 T 为采样周期。
因此，若已知模拟滤波器的传递函数，将上式代入即可得到数字滤波器的系统函数。
在双线性变换中，模拟角频率和数字角频率的变换关系为：，可见，模拟频率和数字频率之间的变换关系为非线性的。
双线性变换法的特点是：（1）由于映射是单值变换关系，避免了冲激响应不变法的频率混叠现象；（2）当频率增加时，变换关系存在着严重的非线性。
在 MATLAB 中，可用函数 bilinear 实现从模拟滤波器到数字滤波器的双线性变换映射，调用格式为：[bz,az]=bilinear(b,a,fs)5. 数字滤波器设计步骤（1）将技术指标转换为模拟低通原型设计性能指标。
（2）估计满足性能指标的模拟低通性能阶数和截止频率。
利用 MATLAB 中 buttord、cheb1ord、cheb2ord、ellipord 等函数，调用格式如：其中，和分别是滤波器的通带截止频率和阻止截止频率，单位为；和分别是通带最大衰减系数和阻带最小衰减系数，单位为dB；‘s’表示为模拟滤波器；函数返回值 n 为模拟滤波器的最小阶数；为模拟滤波器的截止频率(-3dB 频率)，rad/s。
函数适用低通、高通、带通、带阻滤波器。
（3）设计模拟低通原型滤波器。
利用 MATLAB 中 buttap、cheb1ap、cheb2ap、elliap等函数，调用格式为：[z,p,k]=buttap(n)。
采用上述函数所得到原型滤波器的传递函数为零点、极点、增益表达式，需要和函数[b,a]=zp2tf(z,p,k)配合使用，以转化为多项式形式。
（4）由模拟低通原型经频率变换获得模拟低通、高通、带通或带阻滤波器。
利用 MATLAB 中 lp2lp、lp2hp、lp2bp、lp2bs 等函数，调用格式如：[c1,d1]=lp2lp(c,d,Wn)。
（5）利用冲激响应不变法或双线性不变法，实现模拟滤波器到数字滤波器的映射。
五、实验内容1. 利用冲激响应不变法，设计一数字低通滤波器，并画出各滤波器的幅频和相频特性曲线。
滤波器指标、、、，采样频率fs=1kHz，要求：function [db,mag,pha,w]=freqz_m(b,a);%对FIR滤波器，b=h，a=[1][H,w]=freqz(b,a,1000,'whole');H=(H(1:1:501))';w=(w(1:1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);(1)用巴特沃斯滤波器设计；Wp=0.2*pi; Ws=0.3*pi;Rp=0.5; Rs=50;fs=1000;% T=1/fs;[n,Wn]=buttord(Wp*fs,Ws*fs,Rp,Rs,'s');[z,p,k]=buttap(n);[bb,aa]=zp2tf(z,p,k);[b,a]=lp2lp(bb,aa,Wn);[bz,az]=impinvar(b,a,fs);% [H,w]=freqz(b,a,fs);[db,mag,pha,w]=freqz_m(bz,az);subplot(211);plot(w/pi,mag);title(['巴特沃斯幅度响应 n=' int2str(n)])%标题里显示变量xlabel('');ylabel('|H|');axis([0,1,0,1.1]) ;%设置坐标变量及坐标轴范围% set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]); %在指定位置标记grid on; %显示格线% set(gca,'YTickMode','manual','YTick',[0,Attn,Ripple,1]);subplot(212);plot(w/pi,pha/pi);title('巴特沃斯相位响应');xlabel('');ylabel('单位：pi');axis([0,1,-1,1]);grid on; %显示格线% set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);% set(gca,'YTickMode','manual','YTick',[-1,0,1]);grid on;（2）用切比雪夫Ⅰ型设计。
Wp=0.2*pi; Ws=0.3*pi;Rp=0.5; Rs=50;fs=1000;% T=1/fs;[n,Wn]=cheb1ord(Wp*fs,Ws*fs,Rp,Rs,'s'); [z,p,k]=cheb1ap(n,Rp);[bb,aa]=zp2tf(z,p,k);[b,a]=lp2lp(bb,aa,Wn);[bz,az]=impinvar(b,a,fs);% [H,w]=freqz(b,a,fs);[db,mag,pha,w]=freqz_m(bz,az);subplot(211);plot(w/pi,mag);title(['切比雪夫Ⅰ型 n=' int2str(n)])%标题里显示变量xlabel('');ylabel('|H|');axis([0,1,0,1.1]) ;%设置坐标变量及坐标轴范围% set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]); %在指定位置标记grid on; %显示格线% set(gca,'YTickMode','manual','YTick',[0,Attn,Ripple,1]);subplot(212);plot(w/pi,pha/pi);title('切比雪夫Ⅰ型');xlabel('');ylabel('单位：pi');axis([0,1,-1,1]);grid on; %显示格线% set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);% set(gca,'YTickMode','manual','YTick',[-1,0,1]);grid on;提示：matlab中可用freqz( )函数得到数字系统的频率响应，调用格式为：[H,w]= freqz(bz,az,fs )，H为滤波器频率响应，w为H对应频率。
和分别是滤波器的通带截止频率和阻止截止频率，单位为；和分别是通带最大衰减系数和阻带最小衰减系数，单位为dB；1）利用双线性变换法，设计一数字低通滤波器，并画出各滤波器的幅频和相频特性曲线。
滤波器指标、、、，采样频率fs=1kHz，要求：（1）用巴特沃斯滤波器设计；wp=0.2*pi;ws=0.3*pi;fs=1000;T=1/fs;OmegaP=(2/T)*tan(wp/2);OmegaS=(2/T)*tan(ws/2);rp=0.5;rs=50;% ripple=10^(-rp/20);attn=10^(-rs/20);[n,wn]=buttord(OmegaP,OmegaS,rp,rs,'s');[z,p,k]=buttap(n);[b,a]=zp2tf(z,p,k);[bt,at]=lp2lp(b,a,wn);[bz,az]=bilinear(bt,at,fs);% [b,a]=impinvar(b,a,Fs);[db,mag,pha,w]=freqz_m(bz,az);%下面绘出各条曲线subplot(211);plot(w/pi,db);grid;title(['双线性变换法巴特沃斯幅度响应 N=' int2str(n)])%标题里显示变量。