数字信号处理(王震宇张培珍编).pptVIP

（连续周期时间信号的傅里叶变换DFS） 时域连续函数对应频域是非周期的谱，而时域的周期性导致频域的离散。 （离散非周期时间信号的傅里叶变换DTFT） 时域的离散化造成频域的周期延拓，而时域的非周期对应于频域的连续。 离散周期时间信号的傅里叶变换DFS--DFT 离散傅里叶变换的时域和频域都是离散的和周期的。主值区间的变换为DFT变换 周期卷积与线性卷积的区别： (1) 线性卷积在无穷区间求和;周期卷积在一个主值周期内求和 (2) 两个不同长度的序列可以进行线性卷积;只有同周期的两个序列才能进行周期卷积，且周期不变 时域循环卷积过程： 1）补零 2）其中一个序列周期延拓 3）翻褶，截取计算区域 4）循环移位 5）被卷积两序列对应序号值相乘，再相加 6）取主值序列 线性卷积与循环卷积步骤比较 例 1.按抽取方法分： 时间抽取法(DIT Decimation-In-Time); 频率抽取法(DIF Decimation-In-Frequency) 2.按“基数”分：基-2FFT算法；基-4FFT算法；混合基FFT算法；分裂基FFT算法 3.其它方法：离散哈特莱变换；线性调频Z变换(Chrip-z法) 1）输入倒位序，输出自然序 2）蝶形运算 对N = 2L点FFT，共需L级蝶形运算，每级有N/2个蝶形运算组成，蝶形运算两节点的距离:2L-1（L表示级数） 每个蝶形运算有一次复乘和再次复加。如 相关和卷积的时域关系 x=[1 3 -1 1 2 3 3 1]; y=[2 1 -1 1 2 0 -1 3]; k=length(x); xk=fft(x,2*k); yk=fft(y,2*k); rm=real(ifft(conj(xk).*yk)); rm=[rm(k+2:2*k) rm(1:k)]; m=(-k+1):(k-1); stem(m,rm) xlabel(m); ylabel(幅度); 通过FFT来分析其信号频率成分 选择 对x1(n)(0≤n≤N1-1)和x2(n)(0 ≤ n≤N2-1)进行8点的循环卷积，其中( ) 的结果不等于线性卷积。 （A）N1=5，N2=5 （B）N1=3，N2=4 （C）N1=4，N2=4 （D）N1=5，N2=4 填空 1.已知h1(n)={1，2，3，4，4，3，2，1}，的傅里叶变换为H1(k)，h2(n)={4,3,2,1,1,2,3,4}，的傅里叶变换为H2(k)，则， H1(k)和H2(k)的幅频特性____(相等，不等写出关系)。 2.x(n)=cos[(π/6)n]R4(n),则x((9))4=__________。 3.计算N点的DFT，需要计算____次复数乘法，____复数加法。FFT算法需要_____次复数乘法，____复数加法。 4.基2DIF-FFT算法是在_____域进行的，N=2M点的FFT算法需要_____级蝶形运算，每一级要____个蝶，每个蝶形运算要_____次复数乘法，____复数加法。 判断 1.在DIT-FFT算法中，输入顺序为倒序排列，输出为自然顺序（ ） 2.时间抽取算法FFT对两个经时间抽取的N/2点离散序列X(2r)和X(2r+1)做DFT，并将结果相加，就得到一个N点的DFT（ ） 重点习题 5（4），12，13，14，15 用MATLAB FFT计算相关函数 4 -8 -6 -4 -2 0 2 4 6 8 -6 -4 -2 0 2 4 6 8 10 12 m 幅度 两个序列的自相关函数 用MATLAB FFT计算相关函数 4 实例 4 t=0:0.001:0.6; x=sin(2*pi*50*t)+sin(2*pi*120*t); %两个正弦信号合成一个输入 y=x+randn(1,length(t)); %加载噪声信号 Y=fft(y,512); %进行FFT变换 subplot(211);plot(y);title(受噪声污染的信号); n=0:511;f=1000*n/512; %将信号横向放大，由512放大为1000 subplot(212);plot(f,abs(Y));title(FFT); 实例 4 1点DFT x(0) 1点DFT x(4) X3(0) X3(1) 2个1点的DFT蝶形流图 -1 8点DFT为例 4 一个完整的按时间抽取的8点FFT 8点DFT为例 4 由以上按时间抽取法FFT的运算流图，可以看出，当N=2M时，共有M级蝶形，每级都由N/2个蝶形运算组成，每个蝶形有一次复乘、二次复加，因而每级运算都需N/2次复乘和N次复加，这样M级运算总共需要