数字信号处理-王雪梅老师

时域离散随机信号处理及实现（王雪梅老师） 作业1 （1）按下表整理复习离散傅立叶变化相关内容 FT积分变换 FS级数 DTFT

时域

连续，非周期 连续，周期 离散，非周期

12jjnxnXeed

2

1~0NjknNkkxnae

 1

2jnxnXne

频域

jjnnXexne



 

21~~0()NjnkNkXkxne



 

1~~0()NjnkXkxne



连续，非周期 离散，非周期 连续，周期 拉普拉斯变换 LT S平面 Z变换 ZT Z平面

DFS DFT 离散周期 离散，非周期

21~~01()NjnkNkxnXkeN





2101NjnkNkxnXkeN



21~~0()NjnkNkXkxne 210NjnkNkXkxne





离散，周期 离散，非周期 作业2 （1）用脉冲响应不变法设计数字低通滤波器，要求通带和阻带具有单调下降特性。具体指标参数：2.0p，dBp1，35.0s，dB

s10.

（2）用双线性变换法设计上述数字低通滤波器。 （3）对上述两种方法所求结果进行分析、比较。 (1)脉冲响应不变法 相应的模拟滤波器的指标，取采样周期T=1s，

0.2,1pppradsdbT

0.35,10sssradsdbT 题设要求通带和阻带具有单调下降特性，所以选择巴特沃斯滤波器。 通带波纹幅度参数：100.5081018p 阻带波纹幅度参数：20103.1623sA

过度比：47psk 偏离参数：2110.1696kA 巴特沃斯滤波器阶数：1lglg3.1704Nkk，可取N=4 3db截至频率：12210.2659NcsArads 根据：

121221,kNiNNckkkDssppe









查表 可得p1=-0.3827+j0.9239,p2=-0.3827-j0.9239,p3=-0.9239+j0.3827,p4=-0.9239-j0.3827 归一化系统函数：

1234

1Gppppppppp



去归一化，得到模拟低通滤波器的系统函数： 1234|ccapscccc

HsGpspspspsp



将Ha(s)化为：41kakkAHsss形式。 其中

kcksp A1=0.094+0.0939i, A2=0.094-0.0939i,A3=-0.227+0.2221i,A3=-0.227-0.2221i 可以得到脉冲响应不变法所求数字滤波器

324

143210.04560.10270.015411.91841.65460.68530.1127kksTkAzzzHzezzzzz







脉冲响应法设计的数字滤波器频率无非线性畸变，耗损函数与模拟滤波器相似，在接近采样频率处衰减明显小于模拟滤波器，出现混叠失真。 (2)双线性变换设计数字低通滤波器 取T=2s 非线性预畸变

tantan0.10.3249,1pppradsdbT





tantan0.1750.6128,10sssradsdbT





通带波纹幅度参数：100.5081018p 阻带波纹幅度参数：20103.1623sA 过度比：0.5302psk 偏离参数：2110.1696kA 巴特沃斯滤波器阶数：1lglg2.7966Nkk，可取N=3 3db截至频率：12210.4249NcsArads 去归一化，得到模拟低通滤波器的系统函数：



3

2220.07671

|0.42490.42490.1805ccapscccHsGpssssss



双线性变换： 11123123110.033530.10060.10060.03353|10.42490.88270.19azszzzzHzHszzz











由图可见频率非线性畸变会造成数字滤波器的耗损函数与模拟滤波器耗损函数不同，双线性变换没有频谱混叠。

作业3 （1）对模拟信号进行低通滤波处理，要求通带KHzf5.10内衰减小于dB1，阻带fKHz5.2上衰减大于dB40。希望对模拟信号采样后用线性相位FIR数字滤波器实现上述滤波要求，采样频率KHzF

s10

。

用窗函数法设计满足要求的FIR数字滤波器，求出)(nh，并画出损耗函数曲线和相频特性曲线。为了降低运算量，希望滤波器阶数尽量低。 数字滤波器指标 2/21.5/100.3ppsfF

2/22.5/100.5sssfF 为了使滤波器的阶数尽量低，选择凯塞窗，由40sdB，则 0.40.5842210.078863.395321ss

数字滤波器过度带宽度：0.2sp 滤波器阶数：82.28522.2887sMB，取M=23 所以窗长度取N=24

理想低通滤波器的通带截止频率20.4csp 希望逼近的频率响应函数

,0,jcjdceHe





其中1211.5N 傅里叶反变换得到

sin0.411.51.2151ccjindnhneedn







对理想响应函数加窗可得满足要求的线性相位FIR数字滤波器 2411.5sin0.1154.d

nhnhnnnn





 作业4 （1）3.17（用matlab仿真） (1)

可以看出fft快速算法计算计算自相关函数与直接公式计算自相关函数基本重合，但fft快速算法计算计算自相关函数用时0.000844 秒，直接公式计算自相关函数用时0.035199 秒。计算时间大幅减少。

作业5 用matlab编程仿真分析：

（1）3.17（3） AR参数谱 （1）实际计算中，常对ω在2π内均匀采样，则可以用计算模型参数𝑎 𝑘的FFT变换得到p阶AR模型的功率谱 22

2210ˆˆˆjlpMARMjlkMkkSeae







 其中，M为采样点数。 作业6 （1）242P，6-15，用RLS和LMS算法实现)(nu的线性预测。并对两种方法进行分析、比较。

(1)RLS算法 单次实验结果如下图：

(2)LMS算法 单次实验结果如下图： 作业7 1、考虑AR过程)(nu，其差分方程)()2()1()(

21nvnuanuanu

，其中

)(nv是零均值、方差为0731.02v的加性白噪声。AR参数975.01a，95.02a。 （1）产生512N点的),...,2,1)((Nnnu样本序列。 （2）令)(nu为二阶线性预测器)2(LP的输入，在05.0、005.0的情况下用LMS滤波器来估计1w和2w

。

（3）在（2）的参数条件下，滤波器进行100次独立实验，通过平均预测误差)(ˆ)()(nunune的平均值，计算剩余均方误差和失调参数，并画出学习曲线。 （4）改变005.0，其它参数不变，计算剩余均方误差和失调参数，并画出学习曲线，比较05.0和005.0是二者学习曲线的区别。