数字信号处理DSP3-2.pptVIP

fft（x，N） 的点数设置： 设x序列长度为L; 若N值不设，则计算点数为L； 当N值小于L，取x序列前N点计算； 当N值大于L,则对x序列后续补零进行计算； MATLAB中的FFT函数 采用MATLAB函数fftshift可以实现频谱的对称化，该函数将频谱表现在 区间中，对应的频率坐标设置时应注意向下取整。 m=floor(-(N-1)/2):floor((N-1)/2); MATLAB中的FFT函数 [y,fs1,bit]=wav_read(audio1.wav); y1=y(198025:199048); n1=1:1024; figure(1), plot(n1,y1); axis([1,1024,-0.15,0.15]); Y=fftshift(fft(y1,1024)); w=-511:512;k=pi/512; figure(2),plot(k*w,abs(Y)); axis([-3.14,3.14,0,60]); 音频信号的分析仿真 音频信号的分析仿真 信号波形:1024点 信号频谱：幅频特性 第三章 数字信号的频谱分析 频谱分析概述 单频率信号的数字化 数字信号的付氏分析（DTFT） 频谱计算方法：DFT 快速计算方法：FFT 数字信号处理 DTFT的计算问题：连续频谱 如何由频谱表达式计算频谱？ DFT 光滑连接 DTFT DFT：频谱的离散计算 在 范围进行M点等距采样计算 M点DFT的采样点频率 频率点间隔（基本频率）为 采样点频率为 DFT：频谱离散化 DTFT： 从DTFT到DFT：解析表达 利用N个信号测量值，等距计算M个频率分量； 为避免信号发生混叠失真，M不得小于N; 在实际工程计算中，通常选取M=N。 从DTFT到DFT：解析表达 N点DFT的定义 1 对信号进行有限次测量，得到信号测量值表达的数字序列； 2 选择N等于或大于数字序列长度，将数字序列直接代入DFT定义式，得到频谱的解析表达； 3 利用解析表达计算各频率分量值。 N点DFT的计算 N点DFT的计算 16点 128点 N点DFT的计算 DTFT 16点 DFT 128点 DFT N点DFT的计算：MATLAB仿真程序 x=[0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 7/8 3/4 5/8 1/2 3/8 1/4 1/8]; N=4096;X=fftshift(fft(x,4096));T=1; D=2*pi/(N*T); m=floor(-(N-1)/2):floor((N-1)/2); subplot(3,1,1),plot(m*D,abs(X)); axis([-3.5,3.5,-0.5,10]); N=16;X=fftshift(fft(x,16));T=1; D=2*pi/(N*T); m=floor(-(N-1)/2):floor((N-1)/2); subplot(3,1,2),stem(m*D,abs(X)); axis([-3.5,3.5,-0.5,10]); N=128;X=fftshift(fft(x,128));T=1; D=2*pi/(N*T); m=floor(-(N-1)/2):floor((N-1)/2); subplot(3,1,3),stem(m*D,abs(X)); axis([-3.5,3.5,-0.5,10]); 第三章 数字信号的频谱分析 频谱分析概述 单频率信号的数字化 数字信号的付氏分析（DTFT） 频谱计算方法：DFT 快速计算方法：FFT 数字信号处理 N点DFT的计算量 数字信号的加法与乘法 N点DFT计算，需要N2次乘法和N(N-1)次加法； 对M位量化的数据，乘法运算量为加法的M倍。 N点DFT的旋转因子特点 单位圆上N等分点：只有N个不同的点！ N点DFT的旋转因子特点 并不是每个旋转因子都需要乘法： 8点DFT的计算 N点DFT的分解运算 N点DFT的分解运算 基于2点DFT分解的FFT 基于2点DFT的FFT 蝶形运算结构 考虑第k列的变换过程： N/2点DFT结果 N点DFT结果 2点DFT结果 4点DFT结果 8点DFT结果 基于2点DFT的FFT 需要安排K级运算； 每级需要N/2个蝶形运算模块； 累计模块使用量为： 基于2点DFT的FFT 1024点：模块使用量为5120 采用软件运算时，可以将带乘法的蝶形运算和不带乘法的蝶形运算区分对待。 基于2点DFT的FFT 乘法数量为： 运算量的对比： 16点 64点 1024点 整体计算（不考虑旋转因子差别） 256 4