数字信号处理三级项目汇编.pptxVIP

数字信号处理 三级项目 14-通信工程4班 目录 CONTENTS 01 时域信号的采样及谱分析 Part One 谱分析过程 4 3 2 1 利用Tubor C绘制时域采样序列图及频谱图 利用VC++6.0+easyx绘制时域采样序列图及频谱图 利用MATLAB绘制时域采样序列图及频谱图 信号的频谱分布特性 取100点，补零至128点 取100点，补零至256点，256点DFT 高密度谱 取256点，256点DFT 高分辨率谱 取100点，补零至256点，256点DFT 运算量 一般来说,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. DFT与FFT CONTENTS 02 设计FIR滤波器 Part Four 窗函数法 频率取样法 切比雪夫等波纹逼近法 滤波 01 02 03 04 根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口长度N 先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。然后根据过渡带宽度估计窗口长度N。 构造需要逼近的频率响应函数 计算 若已知通带边界频率和阻带边界频率，则取 加窗得到设计结果： 矩形窗 三角形窗 汉宁窗 海明窗 布莱克曼窗 窗函数类型 旁瓣峰值αn/dB 过渡带宽度Bt 阻带最小衰减 αs/dB 近似值 精确值 矩形窗 -13 4π/N 1.8π/N -21 三角窗 -25 8π/N 6.1π/N -25 汉宁窗 -31 8π/N 6.2π/N -44 海明窗 -41 8π/N 6.6π/N -53 布莱克曼窗 -57 12π/N 11π/N -74 CONTENTS 03 利用FIR滤波器滤波 Part five 之 3.1.滤波前后的频域波形（采样点数为200点） 低通滤波器 高通滤波器 带通滤波器 带阻滤波器 fc(Hz) 204 211 fc1=204 fc1=204 fc2=211 fc2=211 ωc(rad) 0.408π 0.422π 0.408π 0.408π 之 3.2.提高采样时间，对信号滤波前后时域频域分析 低通滤波器 fc(Hz) 204 ωc(rad) 0.408π 海明窗的设计，阶数N=475 带通滤波器 fc(Hz) fc1=204 fc2=211 ωc(rad) 0.408π 海明窗的设计，阶数N=475 带阻滤波器 fc(Hz) fc1=204 fc2=211 ωc(rad) 0.408π 海明窗的设计，阶数N=475 高通滤波器 fc(Hz) 211 ωc(rad) 0.422π 海明窗的设计，阶数N=475 CONTENTS 04 思考题 Part Two 用DFT对信号进行频谱分析时，f1=201hz,f2=208hz, f3=214hz,采样频率fs=1000hz,观测时间为Tp=0.1s 时的时域和频域波形： 频域波形 时域波形 采样点N=100，补0至128 此时采样分辨率为F=1/Tp=1/0.1s=10hz,而f1、f2、f3之间互相差7hz、6hz，10hz不是它们的整数倍，所以不能将三个单倍频分辨出来。 保持信号长度M不变，DFT变换的N=256时信号波形： 频域波形 时域波形 采样点N=256 此时相对于N=128时频谱更加密集,其包络更加接近序列的FT。在分辨率F=1/Tp=1/0.1s=10hz,则在保持观测时间不变的条件下不能提高分辨率。 此时由f/fs=k/N可算得k1=20.1,k2=20.8, k3=21.4 采样点为20点（整周期） 采样点为25点（非整周期） 采样点为30点（非整周期） 采样点为80点（整周期） 当采样点N=32,并进行32点DFT，经过计算k=5,27的谱线代表f=40Hz； k=10,22的谱线代表f=80Hz；但对于f=60Hz，对应的k=7.5,是否该频率不能显示出来呢？ 60Hz以频谱泄露的形式出现 80Hz 当采样点N=64,并进行64点DFT，经过计算k=10,54的谱线代表f=40Hz； k=20,44的谱线代表f=80Hz； k=15,49的谱线代表f=60Hz； 40Hz 80Hz 60Hz 40Hz 结论：当输入频率f=k(Fs/N)