时域离散随机信号处理及实现

时域离散随机信号处理及实现（王雪梅老师）作业1按下表整理复习离散傅立叶变化相关内容FT积分变换FS级数DTFT时域连续，非周期连续，周期离散，非周期频域连续，非周期离散，非周期连续，周期拉普拉斯变换LTS平面Z变换ZTZ平面DFSDFT离散周期离散，非周期离散，周期离散，非周期作业2（1）用脉冲响应不变法设计数字低通滤波器，要求通带和阻带具有单调下降特性。具体指标参数：，，，.（2）用双线性变换法设计上述数字低通滤波器。（3）对上述两种方法所求结果进行分析、比较。(1)脉冲响应不变法相应的模拟滤波器的指标，取采样周期T=1s，题设要求通带和阻带具有单调下降特性，所以选择巴特沃斯滤波器。通带波纹幅度参数：阻带波纹幅度参数：过度比：偏离参数：巴特沃斯滤波器阶数：，可取N=43db截至频率：根据：查表可得p1=-0.3827+j0.9239,p2=-0.3827-j0.9239,p3=-0.9239+j0.3827,p4=-0.9239-j0.3827归一化系统函数：去归一化，得到模拟低通滤波器的系统函数：将Ha(s)化为：形式。其中A1=0.094+0.0939i,A2=0.094-0.0939i,A3=-0.227+0.2221i,A3=-0.227-0.2221i可以得到脉冲响应不变法所求数字滤波器脉冲响应法设计的数字滤波器频率无非线性畸变，耗损函数与模拟滤波器相似，在接近采样频率处衰减明显小于模拟滤波器，出现混叠失真。(2)双线性变换设计数字低通滤波器取T=2s非线性预畸变通带波纹幅度参数：阻带波纹幅度参数：过度比：偏离参数：巴特沃斯滤波器阶数：，可取N=33db截至频率：去归一化，得到模拟低通滤波器的系统函数：双线性变换：由图可见频率非线性畸变会造成数字滤波器的耗损函数与模拟滤波器耗损函数不同，双线性变换没有频谱混叠。作业3（1）对模拟信号进行低通滤波处理，要求通带内衰减小于，阻带上衰减大于。希望对模拟信号采样后用线性相位数字滤波器实现上述滤波要求，采样频率。用窗函数法设计满足要求的数字滤波器，求出，并画出损耗函数曲线和相频特性曲线。为了降低运算量，希望滤波器阶数尽量低。数字滤波器指标为了使滤波器的阶数尽量低，选择凯塞窗，由，则数字滤波器过度带宽度：滤波器阶数：，取M=23所以窗长度取N=24理想低通滤波器的通带截止频率希望逼近的频率响应函数其中傅里叶反变换得到对理想响应函数加窗可得满足要求的线性相位FIR数字滤波器作业4（1）3.17（用仿真）(1)可以看出fft快速算法计算计算自相关函数与直接公式计算自相关函数基本重合，但fft快速算法计算计算自相关函数用时0.000844秒，直接公式计算自相关函数用时0.035199秒。计算时间大幅减少。作业5用编程仿真分析：（1）3.17（3）参数谱（1）实际计算中，常对ω在2π内均匀采样，则可以用计算模型参数ak的FFT变换得到p阶AR其中，M为采样点数。作业6，6-15，用和算法实现的线性预测。并对两种方法进行分析、比较。(1)RLS算法单次实验结果如下图：(2)LMS算法单次实验结果如下图：作业71、考虑过程，其差分方程，其中是零均值、方差为的加性白噪声。参数，。（1）产生点的样本序列。（2）令为二阶线性预测器的输入，在、的情况下用滤波器来估计和。（3）在（2）的参数条件下，滤波器进行100次独立实验，通过平均预测误差的平均值，计算剩余均方误差和失