数字信号处理三级项目.pptx

数字信号处理三级项目14-通信工程4班目录Part OnePart ThreePart FourPart Two01时域信号的采样及谱分析思考题设计FIR滤波器利用FIR滤波器滤波020304CONTENTS01Part One时域信号的采样及谱分析谱分析过程4.信号的频谱分布特性1.对模拟信号进行截取（选择合适的窗函数）3.对采样所得序列进行DFT(FFT)运算2.对截取的模拟信号进行采样（选择满足时域采样定理的频谱采样率）3241利用Tubor C绘制时域采样序列图及频谱图利用VC++6.0+easyx绘制时域采样序列图及频谱图利用MATLAB绘制时域采样序列图及频谱图信号的频谱分布特性取100点，补零至128点取100点，补零至256点，256点DFT高分辨率谱高密度谱取100点，补零至256点，256点DFT取256点，256点DFTDFT与FFT运算量 一般来说,FFT比DFT运算量小得多,N点的FFT需要做(N/2)log2N次乘法运算,而N点DFT需要做N2次乘法运算,由此看来N点DFT运算量大约是FFT的2N/log2N倍,例如对1 024点的变换,DFT大约是FFT的200倍.然而实际应用时存在下列情况: 点数或采样率的可选性 对DFT来讲,其变换点数可任意选定,.FFT的变换点数必须是有规律的 实时性 DFT运算可以用采一点后立即进行相乘、累加运算的方法FFT运算必须在全部点采集结束后才能开始进行计算,因此从某种角度讲DFT的实时性优于FFT.CONTENTS02Part Four设计FIR滤波器FIR滤波器的主要设计方法窗函数法切比雪夫等波纹逼近法频率取样法滤波利用窗函数设计FIR滤波器010203计算根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口长度N构造需要逼近的频率响应函数若已知通带边界频率和阻带边界频率，则取先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。然后根据过渡带宽度估计窗口长度N。04加窗得到设计结果：常用的窗函数的幅频特性汉宁窗海明窗矩形窗三角形窗布莱克曼窗汉宁窗布莱克曼窗三角形窗矩形窗海明窗窗函数设计FIR低通滤波器的频率特性不同窗函数对滤波性能的影响窗函数类型旁瓣峰值αn/dB过渡带宽度Bt阻带最小衰减αs/dB近似值精确值矩形窗-134π/N1.8π/N-21三角窗-258π/N6.1π/N-25汉宁窗-318π/N6.2π/N-44海明窗-418π/N6.6π/N-53布莱克曼窗-5712π/N11π/N-74窗函数的选择标准1、较低的旁瓣幅度2、旁瓣幅度要下降得快，以利于增加阻带衰减。3、主瓣宽度要窄，这样滤波器过渡带较窄。CONTENTS03Part five利用FIR滤波器滤波之3.1.滤波前后的频域波形（采样点数为200点）各滤波器截止频率选择低通滤波器高通滤波器带通滤波器带阻滤波器fc(Hz)204211fc1=204fc1=204fc2=211fc2=211ωc(rad)0.408π0.422π0.408π0.408π矩形窗低通FIR滤波器矩形窗高通FIR滤波器矩形窗——低通高通矩形窗带阻FIR滤波器矩形窗带通FIR滤波器矩形窗——带阻带通海明窗低通FIR滤波器汉宁窗低通FIR滤波器其他窗函数的低通滤波器汉宁窗高通FIR滤波器海明窗高通FIR滤波器其他窗函数的高通滤波器海明窗带阻FIR滤波器汉宁窗带阻FIR滤波器其他窗函数的带阻滤波器海明窗带通FIR滤波器汉宁窗带通FIR滤波器其他窗函数的带通滤波器之3.2.提高采样时间，对信号滤波前后时域频域分析取采样点为1000点滤波前的时域图滤波前的时域波形滤波后的时域图海明窗的设计，阶数N=475低通滤波器低通滤波器fc(Hz)204ωc(rad)0.408π滤波后的时域图海明窗的设计，阶数N=475带通滤波器带通滤波器fc(Hz)fc1=204fc2=211ωc(rad)0.408π滤波后的时域图海明窗的设计，阶数N=475带阻滤波器带阻滤波器fc(Hz)fc1=204fc2=211ωc(rad)0.408π滤波后的时域图海明窗的设计，阶数N=475高通滤波器高通滤波器fc(Hz)211ωc(rad)0.422πCONTENTS04Part Two思考题用DFT对信号进行频谱分析时，f1=201hz,f2=208hz,f3=214hz,采样频率fs=1000hz,观测时间为Tp=0.1s时的时域和频域波形：保持信号长度M不变，DFT变换的N=256时信号波形：频域波形 时域波形频域波形此时由f/fs=k/N可算得k1=20.1,k2=20.8,k3=21.4时域波形采样点N=100，补0至128采样点N=256此时采样分辨率为F=1/Tp=1/0.1s=10hz,而