IIR及FIR数字滤波器
一题干
对模拟信号进行低通滤波处理，要求通带0≤f≤4kHz，通带衰减小于0.5dB，阻带4.5k Hz≤f<∞，阻带衰减大于50dB，设采样频率Fs=20kHz。
(1)设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。
(2)分别用脉冲响应不变法和双线性变换法设计IIR低通数字滤波器，求出Ha(z) 的分子、分母多项式系数Bz和Az，并画出幅频响应损耗函数曲线
(3)采用窗函数法（分别用汉宁窗、哈明窗、布莱克曼窗函数）设计满足要求的FIR 低通滤波器，求出h(n)，并画出幅频响应损耗函数曲线.
(4)用频率采样法设计满足要求的FIR低通滤波器，求出h(n),并画出幅频响应损耗函数曲线。
二求解过程
具体内容如下:
（1）设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。
程序：
wp=2*pi*4000;
ws=2*pi*5800;
Rp=0.5;
As=50;
[N,wc]=buttord(wp,ws,Rp,As,'s');
[B,A]=butter(N,wc,'s');
k=0:511;
fk=0:20000/512:20000;
wk=2*pi*fk;
Hk=freqs(B,A,wk);
plot(fk/1000,20*log10(abs(Hk)));
grid on
xlabel('频率/kHz');
ylabel('幅度/dB');
axis([0,6,-65,5]);
波形图：
A = 1.0e+207 *
0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0020 2.1576
B = 1.0e+207 *
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.1576
N = 46
（2）分别用脉冲响应不变法和双线性变换法设计IIR低通数字滤波器，求出Ha(z) 的分子、分母多项式系数Bz和Az，并画出幅频响应损耗函数曲线
脉冲响应不变法
程序：
Fs=20000;
wp=10000*pi;
ws=11600*pi;
Rp=0.5;
As=50;
[N,wc]=buttord(wp,ws,Rp,As,'s');
[B,A]=butter(N,wc,'s');
[Bz,Az]=impinvar(B,A);
k=0:511;
fk=0:20000/512:20000;
wk=2*pi*fk;
Hk=freqs(B,A,wk);
plot(fk/1000,20*log10(abs(Hk)));
grid on;
xlabel('频率/kHz');
ylabel('幅值/dB');
axis([0,6,-65,5]);
波形图:
Bz = 1.0e-007 *
0 -0.0000 0.0000 -0.0000 0.0001 -0.0006 0.0035 -0.0116 0.0279 -0.0745 0.1490 -0.2235 0.3353 -0.3725 0.4470 -0.4098 0.3353 -0.2235 0.1304 -0.0698 0.0291 -0.0093 0.0026 -0.0006 0.0001 -0.0000 0.0000 -0.0000
Az = 1.0e+007 *
0.0000 -0.0000 0.0000 -0.0003 0.0018 -0.0081 0.0296 -0.0888 0.2219 -0.4684 0.8431 -1.3028 1.7370 -2.0041 2.0040 -1.7367 1.3024 -0.8427 0.4681 -0.2217 0.0887 -0.0296 0.0081 -0.0018 0.0003 -0.0000 0.0000 -0.0000
N = 46
双线性变换法
程序：
Fs=20000;wpz=10000 /Fs;
wsz=11600/Fs;
Rp=0.5;
As=50;
wp=2*tan(wpz*pi /2);
ws=2*tan(wsz*pi
/2);[N,wc]=buttord(wp,ws,Rp,As,'s');[B,A]=butter(N,wc,'s');[Bz,Az]=bilinear(B,A,Fs); [Nd,wdc]=buttord(wpz,wsz,Rp,As);
[Bdz,Adz]=butter(Nd,wdc);
k=0:511;
fk=0:20000/512:20000;
wk=2*pi*fk;Hk=freqs(B,A,wk);plot(fk/1000,20*log10(abs(Hk)));
grid on;xlabel('频率/kHz');
ylabel('幅值/dB');axis([0,16,-2800,5])
波形图:
Bz = 1.0e-007 *
0 -0.0000 0.0000 -0.0000 0.0001 -0.0006 0.0035 -0.0116 0.0279 -0.0745 0.1490 -0.2235 0.3353 -0.3725 0.4470 -0.4098 0.3353 -0.2235 0.1304 -0.0698 0.0291 -0.0093 0.0026 -0.0006 0.0001 -0.0000 0.0000 -0.0000
Az = 1.0e+007 *
0.0000 -0.0000 0.0000 -0.0003 0.0018 -0.0081 0.0296 -0.0888 0.2219 -0.4684 0.8431 -1.3028 1.7370 -2.0041 2.0040 -1.7367 1.3024 -0.8427 0.4681 -0.2217 0.0887 -0.0296 0.0081 -0.0018 0.0003 -0.0000 0.0000 -0.0000
N = 27
(3)采用窗函数法（分别用汉宁窗、哈明窗、布莱克曼窗函数）设计满足要求的FIR
低通滤波器，求出h(n)，并画出幅频响应损耗函数曲线.
a ）汉宁窗
程序：
Fs=20000;
fp=4000;
fs=5800;
m=[1 1 0 0];
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
Rp=0.5;
As=50;
Bt=ws-wp;
N0 = ceil (6.6* pi /Bt);
N=N0+mod(N0+1,2);
wc=(wp+ws)/2/pi;
hn=fir1(N-1,wc,hanning(N));
freqz(hn,1,512);
plot (w,20*log(abs(hn)));
grid;
axis ([0 ,1, -1000 , 100]);
xlabel ('频率/kHz');
ylabel ('幅值/dB' );
波形图：
hn =
0.0000 -0.0000 -0.0000 0.0002 -0.0000 -0.0004 0.0002 0.0007 -0.0006 -0.0010 0.0012 0.0012 -0.0021 -0.0010 0.0032 0.0005 -0.0045 0.0006 0.0057 -0.0025 -0.0066 0.0050 0.0070 -0.0083 -0.0065 0.0123 0.0047 -0.0168 -0.0012 0.0215 -0.0046 -0.0262 0.0134 0.0307 -0.0270 -0.0345 0.0497 0.0375 -0.0974 -0.0394 0.3154 0.5400 0.3154 -0.0394 -0.0974 0.0375 0.0497 -0.0345 -0.0270 0.0307 0.0134 -0.0262 -0.0046 0.0215 -0.0012 -0.0168 0.0047 0.0123 -0.0065 -0.0083 0.0070 0.0050 -0.0066 -0.0025 0.0057 0.0006 -0.0045 0.0005 0.0032 -0.0010 -0.0021 0.0012 0.0012 -0.0010 -0.0006 0.0007 0.0002 -0.0004 -0.0000 0.0002 -
0.0000 -0.0000 0.0000
b )哈明窗
程序：
Fs=20000;
fp=4000;
fs=5800;
m=[1 1 0 0];
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
Rp=0.5;
Rs=50;
Bt=ws-wp;
N0 = ceil (6.6* pi /Bt);
N=N0+mod(N0+1,2);
wc=(wp+ws)/2/pi;
hn=fir1(N-1,wc,hamming(N)); freqz(hn,1,512);
plot (w,20*log(abs(hn)));
grid on;
axis ([0 ,1, -1000 , 100]); xlabel ('频率/kHz');
ylabel ('幅值/dB' );
波形图：
hn =
0.0003 -0.0006 -0.0001 0.0008 -0.0001 -0.0010 0.0004 0.0012 -0.0008 -0.0014 0.0016 0.0015 -0.0025 -0.0012 0.0037 0.0005 -0.0049 0.0007 0.0061 -0.0026 -0.0069 0.0052 0.0072 -0.0086 -0.0066 0.0125 0.0048 -0.0170 -0.0012 0.0217 -0.0046 -0.0264 0.0135 0.0308 -0.0271 -0.0346 0.0498 0.0375 -0.0975 -0.0394 0.3154 0.5401 0.3154 -0.0394 -0.0975 0.0375 0.0498 -0.0346 -0.0271 0.0308 0.0135 -0.0264 -0.0046 0.0217 -0.0012 -0.0170 0.0048 0.0125 -0.0066 -0.0086 0.0072 0.0052 -0.0069 -0.0026 0.0061 0.0007 -0.0049 0.0005 0.0037 -0.0012 -0.0025 0.0015 0.0016 -0.0014 -0.0008 0.0012 0.0004 -0.0010 -0.0001 0.0008 -
0.0001 -0.0006 0.0003
c )布莱克曼窗
程序：
Fs=20000;
fp=4000;
fs=5800;
m=[1 1 0 0];
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
Rp=0.5;
Rs=50;
Bt=ws-wp;
N0 = ceil (6.6* pi /Bt);
N=N0+mod(N0+1,2);
wc=(wp+ws)/2/pi;
hn=fir1(N-1,wc,blackman(N)); freqz(hn,1,512);
plot (w,20*log(abs(hn))); grid;
axis ([0 ,1, -1000 , 100]); xlabel ('频率/kHz');
ylabel ('幅值/dB' );
波形图：
hn =
-0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0001 0.0001 0.0002 -0.0002 -0.0004 0.0005 0.0005 -0.0009 -0.0005 0.0016 0.0002 -0.0024 0.0004 0.0034 -0.0015 -0.0043 0.0034 0.0049 -0.0060 -0.0049 0.0095 0.0037 -0.0137 -0.0010 0.0186 -0.0040 -0.0237 0.0124 0.0288 -0.0257 -0.0333 0.0485 0.0369 -0.0965 -0.0392 0.3150 0.5400 0.3150 -0.0392 -0.0965 0.0369 0.0485 -0.0333 -0.0257 0.0288 0.0124 -0.0237 -0.0040 0.0186 -0.0010 -0.0137 0.0037 0.0095 -0.0049 -0.0060 0.0049 0.0034 -0.0043 -0.0015 0.0034 0.0004 -0.0024 0.0002 0.0016 -0.0005 -0.0009 0.0005 0.0005 -0.0004 -0.0002 0.0002 0.0001 -0.0001 -0.0000 0.0000 -
0.0000 -0.0000 -0.0000
（4）用频率采样法设计满足要求的FIR低通滤波器，求出h(n),并画出幅频响应损耗函数曲线。
程序：
T=input('T=')
fp=4000;
fs=5800;
Fs=20000;
wp=2*pi*fp/Fs;
ws=2*pi*fs/Fs;
Rp=0.5;
As=50;
Bt=ws-wp;
m=1;
N=ceil((m+1)*2*pi/Bt);
N=N+mod(N+1,2);
Np=fix(wp/(2*pi/N));
Ns=N-2*Np-1;
Hk=[ones(1,Np+1),zeros(1,Ns),ones(1,Np)];
Hk(Np+2)=T;
Ak(N-Np)=T;
thetak=-pi*(N-1)*(0:N-1)/N;
hdk=Hk.*exp(j*thetak);
hn=real(ifft(hdk));
hw=fft(hn,1024);
wk=2*pi*[0:1023]/1024;
hgw=hw.*exp(j*wk*(N-1)/2);
Rp=max(20*log10(abs(hgw)));
hgmin=min(real(hgw));
As=20*log10(abs(hgmin));
[N,wc]=buttord(wp,ws,Rp,As,'s');
[B,A]=butter(N,wc,'s');
k=0:511;fk=0:20000/512:20000;wk=2*pi*fk; Hk=freqs(B,A,wk);
plot(fk/1000,20*log10(abs(Hk)));
grid on;
xlabel('频率/kHz');
ylabel('幅度/dB');
T=1
T =
1
Rp =
0.2801
As =
-29.6745
波形图：
hn =
0.0004 0.0012 -0.0023 -0.0027 0.0045 0.0040 -0.0069 -0.0051 0.0098 0.0061 -0.0133 -0.0070 0.0175 0.0078 -0.0230 -0.0084 0.0305 0.0089 -0.0416 -0.0093 0.0609 0.0096 -0.1044 -0.0097 0.3178 0.5098 0.3178 -0.0097 -0.1044 0.0096 0.0609 -0.0093 -0.0416 0.0089 0.0305 -0.0084 -0.0230 0.0078 0.0175 -0.0070 -0.0133 0.0061 0.0098 -0.0051 -0.0069 0.0040 0.0045 -0.0027 -0.0023 0.0012 0.0004
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形
第一章数字信号处理概述 简答题： 1．在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器，它们分别起什么作用？ 答：在A/D变化之前为了限制信号的最高频率，使其满足当采样频率一定时，采样频率应大于等于信号最高频率2倍的条件。此滤波器亦称为“抗混叠”滤波器。 在D/A变换之后为了滤除高频延拓谱，以便把抽样保持的阶梯形输出波平滑化，故又称之为“平滑”滤波器。 判断说明题： 2．模拟信号也可以与数字信号一样在计算机上进行数字信号处理，自己要增加一道采样的工序就可以了。 （） 答：错。需要增加采样和量化两道工序。 3．一个模拟信号处理系统总可以转换成功能相同的数字系统，然后基于数字信号处理理论，对信号进行等效的数字处理。（） 答：受采样频率、有限字长效应的约束，与模拟信号处理系统完全等效的数字系统未必一定能找到。因此数字信号处理系统的分析方法是先对抽样信号及系统进行分析，再考虑幅度量化及实现过程中有限字长所造成的影响。故离散时间信号和系统理论是数字信号处
理的理论基础。 第二章 离散时间信号与系统分析基础 一、连续时间信号取样与取样定理 计算题： 1．过滤限带的模拟数据时，常采用数字滤波器，如图所示，图中T 表示采样周期（假设T 足够小，足以防止混叠效应），把从)()(t y t x 到的整个系统等效为一个模拟滤波器。 （a ） 如果kHz T rad n h 101,8)(=π截止于，求整个系统的截止频 率。 （b ） 对于kHz T 201=，重复（a ）的计算。 采样（T） () n h () n x () t x () n y D/A 理想低通T c πω=() t y 解 （a ）因为当0)(8=≥ω πωj e H rad 时，在数 — 模变换中 )(1)(1)(T j X T j X T e Y a a j ωω=Ω= 所以)(n h 得截止频率8πω=c 对应于模拟信号的角频率c Ω为 8 π = ΩT c 因此 Hz T f c c 625161 2==Ω= π
数字信号处理作业
DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~1k X 是周期性的，周期为N ，而)(~2k X 也是周期性的，周期为N 2。试利用)(~1k X 确定)(~2k X 。（76-4）
2. 研究两个周期序列)(~n x 和)(~n y 。)(~n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000) ()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
一、单项选择题 1．数字信号的特征是( ) A.时间离散、幅值连续 B.时间离散、幅值量化 C.时间连续、幅值量化 D.时间连续、幅值连续 2．若一线性移不变系统当输入为x(n)=δ(n)时，输出为y(n)=R 2(n)，则当输入为u(n)-u(n-2)时，输出为( ) A.R 2(n)-R 2(n-2) B.R 2(n)+R 2(n-2) C.R 2(n)-R 2(n-1) D.R 2(n)+R 2(n-1) 3．下列序列中z 变换收敛域包括|z|=∞的是( ) A.u(n+1)-u(n) B.u(n)-u(n-1) C.u(n)-u(n+1) D.u(n)+u(n+1) 4．下列对离散傅里叶变换（DFT ）的性质论述中错误的是( ) A.DFT 是一种线性变换 B.DFT 具有隐含周期性 C.DFT 可以看作是序列z 变换在单位圆上的抽样 D.利用DFT 可以对连续信号频谱进行精确分析 5．若序列的长度为M ，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混叠现象，则频域抽样点数N 需满足的条件是( ) A.N ≥M B.N ≤M C.N ≥M/2 D.N ≤M/2 6．基-2 FFT 算法的基本运算单元为( ) A.蝶形运算 B.卷积运算 C.相关运算 D.延时运算 7．以下对有限长单位冲激响应（FIR ）滤波器特点的论述中错误的是( ) A.FIR 滤波器容易设计成线性相位特性 B.FIR 滤波器的单位冲激抽样响应h(n)在有限个n 值处不为零 C.系统函数H(z)的极点都在z=0处 D.实现结构只能是非递归结构 8．下列结构中不属于IIR 滤波器基本结构的是( ) A.直接型 B.级联型 C.并联型 D.频率抽样型 9．下列关于用冲激响应不变法设计IIR 滤波器的说法中错误的是( ) A.数字频率与模拟频率之间呈线性关系 B.能将稳定的模拟滤波器映射为一个稳定的数字滤波器 C.使用的变换是s 平面到z 平面的多值映射 D.可以用于设计低通、高通和带阻等各类滤波器 10．离散时间序列x (n )=cos(n 73π-8 π)的周期是( ) A.7 B.14/3 C.14 D.非周期 11．下列系统(其中y(n)是输出序列，x(n)是输入序列)中______属于线性系统。( ) A.y (n )=x 2(n ) B.y (n )=4x (n )+6 C.y (n )=x (n -n 0) D.y (n )=e x (n )
2019~2020年度《数字信号处理》大作业题目与要求 大作业要求： 本学期大作业总分40分，学生可选择任意数量的题目完成，只要所选题目总分达到40分即可，所选题目总分如果超过40分，超过的部分不计入大作业总分。大作业以电子版的形式提交，内容应包括详细的程序设计思路与题目分析（题目分析指的是对该题目中所用到的知识点的说明，不要照搬书上或网上的内容，写出你自己对该知识点的理解。），程序截图，程序源码，其中设计思路和程序截图可写在同一个文档中，程序源码可以是.txt或.m 文件，并在源码中标注代码注释。另：题目中有GUI设计要求的部分占该题目分值的20%，功能实现部分占该题目分值的80%。 注：以下题目均用MATLAB完成。 大作业题目： 1、实现有限长序列的基本运算（包括：加法、乘法、累加、移位、翻褶、抽取、插值、卷积和），并以GUI的形式将这些运算整合起来，使用者可通过向GUI输入任意有限长序列得到对应的运算结果。（5分） 2、设计一个GUI，实现奈奎斯特采样定理，要求：1、在GUI中输入任意一个模拟信号，显示该模拟信号的时域和频域谱图；2、在GUI中设置任意采样频率，对输入的模拟信号进行采样处理，显示采样信号的时域和频域谱图； 3、在GUI中实现采样信号向模拟信号的恢复功能，要求显示恢复后的模拟信号的时域和频域谱图。（10分） 3、通过GUI动态展示z变换与s变换之间的所有关系。（5分） 4、设计一个GUI，通过向GUI输入任意系统函数，得到其对应系统的相关信息（包括：系统频率响应中的幅度响应和相位响应、系统零极点的分布、系统的稳定性判定）。（10分） 5、设计一个GUI，实现利用DFT（或FFT）完成任意时域信号的频谱分析，要求：1、可在GUI中输入时域数字或模拟信号；2、可设置DFT点数；3、在GUI中显示输入信号经DFT（或FFT）处理后的频谱图；3、若输入信号为模拟信号，需完成对该模拟信号的采样，采样频率可在GUI中设置。（10分） 6、在GUI中，实现IIR滤波器的直接型、级联型和并联型三种结构之间的任意转换，要求：在GUI中输入任意一型的系统函数后可在该GUI中显示出对应的另外两型的系统函数。（10分） 7、实现巴特沃斯样本模拟低通滤波器及其对应的数字低通滤波器的设计，以GUI的形式给出。要求：输入所需的模拟低通滤波器参数指标后，程序能将该指标转化为数字低通滤波器指标（在GUI中应能选择转化方式：冲激响应不变法、双线性变换法），并在GUI中显示出所给参数下巴特沃斯样本模拟低通滤波器及其对应的数字低通滤波器的频率响应中幅度响应的频谱图。（15分） 8、已知某组数字信号（见大作业数据压缩包中HWDATA.mat文件），该信号中除了目标信号之外还掺杂有强噪声，但噪声与目标信号的频率不重叠，要求采用本学期已学的知识对该信
武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02
一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分
实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列');
实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列');
正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列');
随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列');
第二章 2.25 已知线性时不变系统的差分方程为 若系统的输入序列x(x)={1,2,3,4,2,1}编写利用递推法计算系统零状态响应的MATLAB程序，并计算出结果。 代码及运行结果： >> A=[1,-0.5]; >> B=[1,0,2]; >> n=0:5; >> xn=[1,2,3,4,2,1]; >> zx=[0,0,0];zy=0; >> zi=filtic(B,A,zy,zx); >> yn=filter(B,A,xn,zi); >> figure(1) >> stem(n,yn,'.'); >> grid on;
2.28图所示系统是由四个子系统T1、T2、T3和T4组成的，分别用单位脉冲响应或差分方程描述为 T1： 其他 T2: 其他 T3： T4: 编写计算整个系统的单位脉冲响应h(n)，0≤n≤99的MATLAB程序，并计算结果。 代码及结果如下： >> a=0.25;b=0.5;c=0.25; >> ys=0; >> xn=[1,zeros(1,99)]; >> B=[a,b,c]; >> A=1; >> xi=filtic(B,A,ys); >> yn1=filter(B,A,xn,xi); >> h1=[1,1/2,1/4,1/8,1/16,1/32]; >> h2=[1,1,1,1,1,1]; >> h3=conv(h1,h2); >> h31=[h3,zeros(1,89)]; >> yn2=yn1+h31; >> D=[1,1];C=[1,-0.9,0.81]; >> xi2=filtic(D,C,yn2,xi); >> xi2=filtic(D,C,ys); >> yn=filter(D,C,yn2,xi); >> n=0:99; >> figure(1) >> stem(n,yn,'.'); >> title('单位脉冲响应'); >> xlabel('n');ylabel('yn');
数字信号处理上机作业 学院：电子工程学院 班级：021215 组员：
实验一：信号、系统及系统响应 1、实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 (2) 熟悉时域离散系统的时域特性。 (3) 利用卷积方法观察分析系统的时域特性。 (4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 2、实验原理与方法 (1) 时域采样。 (2) LTI系统的输入输出关系。 3、实验内容及步骤 (1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。 (2) 编制实验用主程序及相应子程序。 ①信号产生子程序，用于产生实验中要用到的下列信号序列： a. xa(t)=A*e^-at *sin(Ω0t)u(t) b. 单位脉冲序列：xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③有限长序列线性卷积子程序 用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数conv。 conv 用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。调用格式如下： y=conv (x, h) 4、实验结果分析 ①分析采样序列的特性。 a. 取采样频率fs=1 kHz,，即T=1 ms。 b. 改变采样频率，fs=300 Hz，观察|X(e^jω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的|X(e^j ω)|曲线。 程序代码如下： close all;clear all;clc; A=50; a=50*sqrt(2)*pi; m=50*sqrt(2)*pi; fs1=1000; fs2=300; fs3=200; T1=1/fs1; T2=1/fs2; T3=1/fs3; N=100;
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
DSP的 数字信号处理（Digital Signal Processing，简称DSP）是一门涉及许多学科而又广泛应用于许多领域的新兴学科。DSP有两种含义：Digital Signal Processing（数字信号处理）、Digital Sign al Processor（数字信号处理器）。我们常说的DSP指的是数字信号处理器。数字信号处理器是一种适合完成数字信号处理运算的处理器。20世纪60年代以来，随着计算机和信在过去的二十多年时间里，数字信号处理已经在通信等领域得到息技术的飞速发展，数字信号处理技术应运而生并得到迅速的发展。极为广泛的应用。数字信号处理是利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以得到符合人们需要的信号形式。 数字信号处理是围绕着数字信号处理的理论、实现和应用等几个方面发展起来的。数字信号处理在理论上的发展推动了数字信号处理应用的发展。反过来，数字信号处理的应用又促进了数字信号处理理论的提高。而数字信号处理的实现则是理论和应用之间的桥梁。数字信号处理是以众多学科为理论基础的，它所涉及的范围极其广泛。例如，在数学领域，微积分、概率统计、随机过程、数值分析等都是数字信号处理的基本工具，与网络理论、信号与系统、控制论、通信理论、故障诊断等也密切相关。近来新兴的一些学科，如人工智能、模式识别、神经网络等，都与数字信号处理密不可分。可以说，数字信
号处理是把许多经典的理论体系作为自己的理论基础，同时又使自己成为一系列新兴学科的理论基础。顾名思义,DSP主要应用在数字信号处理中，目的是为了能够满足实时信号处理的要求，因此需要将数字信号处理中的常用运算执行的尽可能快，这就决定了DSP的特点和关键技术。适合数字信号处理的关键技术：DSP包含乘法器、累加器、特殊地址产生器、领开销循环等；提高处理速度的关键技术：流水线技术、并行处理技术、超常指令（VLIW）、超标量技术、DMA等。从广义上讲，DSP、微处理器和微控制器（单片机）等都属于处理器，可以说DSP是一种CPU。DSP和一般的CPU又不同，最大的区别在于：CPU是冯.诺伊曼结构的；DSP是数据和地址空间分开的哈佛结构。 世界上第一个单片 DSP 芯片应当是1978年 AMI公司发布的 S2 811，1979年美国Intel公司发布的商用可编程器件2920是DSP芯片的一个主要里程碑。这两种芯片内部都没有现代DSP芯片所必须有的单周期乘法器。1980 年，日本 NEC 公司推出的μP D7720是第一个具有乘法器的商用 DSP 芯片。在这之后，最成功的DSP 芯片当数美国德州仪器公司（Texas Instruments，简称TI）的一系列产品。TI 公司在1982年成功推出其第一代 DSP 芯片 TMS32010及其系列产品TMS32011、TMS320C10/C14/C15/C16/C17等，之后相继推出了第二代DSP芯片TMS32020、TMS320C25/C26/C28，第三代DSP芯片TMS 320C30/C31/C32，第四代DSP芯片TMS320C40/C44，第五代 DSP 芯片TMS320C5X/C54X，第二代DSP芯片的改进型TMS320C2XX，集多片DSP芯片于一体的高性能DSP芯片TMS320C8X以及目前速度最快的第
数字信号处理作业 哈尔滨工业大学 2006.10
DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~ 2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~ 1k X 是周期性的，周期为N ，而)(~ 2k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2k X 。（76-4）
2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列 )(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~ k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~ n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~ k X 和)(~ k Y 求)(~ k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
数字信号处理习题集 第一章习题 1、已知一个5点有限长序列，如图所示，h (n )=R 5(n )。（1）用写出的 ()n δ()x n 函数表达式；（2）求线性卷积*。 ()y n =()x n ()h n 2、已知x (n )=(2n +1)[u (n +2)-u (n -4)],画出x (n )的波形，并画出x (-n )和x (2n )的波形。 3、判断信号是否为周期信号，若是求它的周期。3()sin 7 3x n n π π??=+ ???4、判断下列系统是否为线性的，时不变的，因果的，稳定的？ （1），（2）2()(3)y n x n =-0()()cos() y n x n n ω=5、已知连续信号。()2sin(2),3002 a x t ft f Hz π π=+=（1）求信号的周期。 ()a x t （2）用采样间隔T=0.001s 对进行采样，写出采样信号的表达式。()a x t ?()a x t （3）写出对应于的时域离散信号的表达式，并求周期。?()a x t ()x n 6、画出模拟信号数字处理的框图，并说明其中滤波器的作用。
第二章习题 1、求下列序列的傅立叶变换。 （1）， （2）11()333n x n n ?? =-≤ ? ?? [] 2()()()n x n a u n u n N =--2、已知理想低通滤波器的频率响应函数为：为整数，000(),0j n j e H e n ωωωωωωπ-?≤≤?=? <≤?? c c 求所对应的单位脉冲响应h (n )。 3、已知理想高通滤波器的频率响应函数为：，求所对应 0()1j H e ω ωωωωπ ?≤≤?=? <≤?? c c 的单位脉冲响应h (n )。 4、已知周期信号的周期为5,主值区间的函数值=,求该周期信号的 ()(1)n n δδ+-离散傅里叶级数和傅里叶变换. 5、已知信号的傅立叶变换为，求下列信号的傅立叶变换。 ()x n ()j X e ω（1） （2）(3)x n -*() x n -6、已知实因果信号如图所示，求和。 ()x n ()e x n ()o x n 7、已知实因果信号的偶分量为{-2,-3,3,4,1,4,3,-3,-2},求信号。 ()x n ()x n 8、已知信号,对信号采样,得到时域采样信号和时()cos(2100),300a s x t t f Hz π==?()a x t 域离散信号x(n),求: (1)写出信号的傅里叶变换. ()a x t
数字信号处理作业-答案
数字信号处理作业
DFT 习题 1. 如果)(~ n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1 k X 表示)(~ n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2 k X 表示)(~ n x 的离散傅里叶级数之系数。当然，)(~ 1 k X 是周期性的，周期为N ，而)(~ 2 k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2 k X 。（76-4）
2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~ n w 定义为)()()(~~ ~ n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~ n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地，由于)(~ n y 的周期为M ，其离散傅里叶级数之系数)(~ k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~ k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=0 0)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
第二章上机作业 1、ljdt(A,B)函数定义 function ljdt(A,B) p=roots(A); q=roots(B); p=p'; q=q'; x=max(abs([p q 1])); x=x+0.1; y=x; clf hold on axis([-x x -y y]) w=0:pi/300:2*pi; t=exp(i*w); plot(t) axis('square') plot([-x x],[0 0]) plot([0 0],[-y y]) text(0.1,x,'jIm[z]') text(y,1/10,'Re[z]') plot(real(p),imag(p),'x') plot(ral(q),imag(q),'o') title('pole-zero diagram for discrete system') hold off 例2.26 a=[3 -1 0 0 0 1]; b=[1 1]; ljdt(a,b) p=roots(a) q=roots(b) pa=abs(p) 程序运行结果如下： P= 0.7255+0.4633i 0.7255+0.4633i -0.1861+0.7541i -0.1861-0.7541i -0.7455 q=
-1 pa= 0.8608 0.8608 0.7768 0.7768 0.7455 例2.27 b=[0 1 2 1];a=[1 -0.5 -0.005 0.3]; subplot 311 zplane(b,a);xlabel('实部');ylabel('虚部'); num=[0 1 2 1];den=[1 -0.5 -0.005 0.3]; h=impz(num,den); subplot 312
《数字信号处理Ⅰ》作业 姓名： 学号： 学院： 2012 年春季学期
第一章 时域离散信号和时域离散系统 月 日 一 、判断： 1、数字信号处理和模拟信号处理在方法上是一样的。（ ） 2、如果信号的取值和自变量都离散，则称其为模拟信号。（ ） 3、如果信号的取值和自变量都离散，则称其为数字信号。（ ） 4、时域离散信号就是数字信号。（ ） 5、正弦序列都是周期的。（ ） 6、序列)n (h )n (x 和的长度分别为N 和M 时，则)n (h )n (x *的长度为N+M 。（ ） 7、如果离散系统的单位取样响应绝对可和，则该系统稳定。（ ） 8、若满足采样定理，则理想采样信号的频谱是原模拟信号频谱以s Ω（采样频率）为周期进行周期延拓的结果。（ ） 9、序列)n (h )n (x 和的元素个数分别为21n n 和，则)n (h )n (x *有（1n n 21-+）个元素。（ ） 二、选择 1、R N (n)和u(n)的关系为（ ）： A. R N (n)=u(n)-u(n-N) B. R N (n)=u(n)+u(n-N) C. R N (n)=u(n)-u(n-N-1) D. R N (n)=u(n)-u(n-N+1) 2、若f(n)和h(n)的长度为别为N 、M ，则f(n)*h(n)的长度为 （ ）： A.N+M B.N+M-1 C.N-M D.N-M+1 3、若模拟信号的频率范围为[0，1kHz],对其采样，则奈奎斯特速率为（ ）： A.4kHz B. 3kHz C.2kHz D.1kHz 4、LTIS 的零状态响应等于激励信号和单位序列响应的（ ）： A.相乘 B. 相加 C.相减 D.卷积 5、线性系统需满足的条件是（ ）： A.因果性 B.稳定性 C.齐次性和叠加性 D.时不变性 6、系统y(n)=f(n)+2f(n-1)(初始状态为0)是（ ）： A. 线性时不变系统 B. 非线性时不变系统 C. 线性时变系统 D. 非线性时变系统
IIR及FIR数字滤波器 一题干 对模拟信号进行低通滤波处理，要求通带0≤f≤4kHz，通带衰减小于0.5dB，阻带4.5k Hz≤f<∞，阻带衰减大于50dB，设采样频率Fs=20kHz。 (1)设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。 (2)分别用脉冲响应不变法和双线性变换法设计IIR低通数字滤波器，求出Ha(z) 的分子、分母多项式系数Bz和Az，并画出幅频响应损耗函数曲线 (3)采用窗函数法（分别用汉宁窗、哈明窗、布莱克曼窗函数）设计满足要求的FIR 低通滤波器，求出h(n)，并画出幅频响应损耗函数曲线. (4)用频率采样法设计满足要求的FIR低通滤波器，求出h(n),并画出幅频响应损耗函数曲线。
二求解过程 具体内容如下: （1）设计巴特沃斯模拟低通滤波器，求出Ha(s)的分子、分母多项式系数B和A，并画出幅频响应损耗函数曲线。 程序： wp=2*pi*4000; ws=2*pi*5800; Rp=0.5; As=50; [N,wc]=buttord(wp,ws,Rp,As,'s'); [B,A]=butter(N,wc,'s'); k=0:511; fk=0:20000/512:20000; wk=2*pi*fk; Hk=freqs(B,A,wk); plot(fk/1000,20*log10(abs(Hk))); grid on xlabel('频率/kHz'); ylabel('幅度/dB'); axis([0,6,-65,5]); 波形图：
A = 1.0e+207 * 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0020 2.1576 B = 1.0e+207 * 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.1576 N = 46
数字信号处理第三章作业 1．(第三章习题3)在图P3-2中表示了两个周期都为6的周期性序列，确定这个两个序列的周期卷积的结果3()x n ，并画出草图。 2．(第三章习题5)如果()x n 是一个具有周期为N 的周期性序列，它也是具有周期为2N 的周期性序列。令~1()X k 表示当()x n 看做是具有周期为N 的周期性序列的DFS 系数。而~2()X k 表示当()x n 看作是具有周期为2N 的周期性序列的DFS 系数。当然~1()X k 是具有周期为N 的周期性序列，而~2()X k 是具有周期为2N 的周期性序列，试根据~1()X k 确定~2()X k 。 3．(第三章习题6) （a ）试证明下面列出的周期性序列离散傅里叶级数的对称特性。在证明中，可以利用离散傅里叶级数的定义及任何前面的性质，例如在证明性质③时可以利用性质①和②。 序列 离散傅里叶级数 ① *()x n ~*()X k - ②*()x n - ~*()X k ③Re ()x n ???? ~ e ()X k ④Im ()j x n ???? ~()o X k
（b ）根据已在（a ）部分证明的性质，证明对于实数周期序列()x n ，离散傅里叶级数的下列对称性质成立。 ①~~Re ()Re ()X k X k ????=-???????? ②~~Im ()Im ()X k X k ????=--???????? ③~~()()X k X k =- ④~~arg ()arg ()X k X k ????=--???????? 4．(第三章习题7)求下列序列的DFT (a) {}11 1-，，，-1 (b) {}1 j 1j -，，，- (c) ()cn 0n 1x n N =≤≤-， (d) 2n ()sin 0n 1x n N N π??=≤≤- ??? ， 5．(第三章习题8)计算下列各有限长序列的离散傅立叶变换（假设长度为N ） 1 0)()(0) ()()() ()()(00-≤≤=<<-==N n a n x c N n n n n x b n n x a n δδ 6．(第三章习题9)在图P3-4中表示了一有限长序列)(n x ，画出序列)(1n x 和)(2n x 的草图。（注意：)(1n x 是)(n x 圆周移位两个点） )())(()() ())2(()(442441n R n x n x n R n x n x -=-=