FT实现频谱分析与线性卷积实验四.docVIP

实验报告 实验名称DFT实现频谱分析与线性卷积 课程名称 数字信号处理实验 院系部：电气学院 学生姓名：庞思颜 指导教师：范杰清 实验日期：2017-11-15 专业班级：电网1501 学 号：1151600115 成 绩： 华北电力大学 连续周期信号恥离散非周期信号％(/7)X(j/2) = J^xO) e 连续周期信号恥 离散非周期信号％(/7) X(j/2) = J^xO) e_j/2zdz w)= 火=一00 DSP4 : DFT实现频谱分析与线性卷积 一、 实验目的 应用离散傅里叶变换(DFT)，分析离散信号的频谱。 深刻理解DFT分析离散信号频谱的原理，掌握改善分析过程中产生的误差的方法。 应用离散傅里叶变换⑴H),实现离散信号与系统的线性卷积，并理解实现条件。 二、 实验原理一DFT实现频谱分析 根据信号傅里叶变换建立的时域与频域之间的对应关系，可以得到有限长序列的离散傅 里叶变换(DFT)与叫种确定信号傅里叶变换的之间的关系，实现由DFT分析其频谱。 四种信号的频谱函数 连续叩周期信号;vU) ? 2n ? 2n , -\—nk J N 离散周期信号%(/?)  /V-1 X(k) = ^x(n)^ 人=0 四种信号的频谱函数之间的相互关系 信号的傅里叶变换建立了信号的吋域与频域之间的一一对应关系，如果信号在吋域存在 某种联系，则在其频谱函数之间必然存在联系。 若离散非周期信号奴/;)是连续非周期信号Hf)的等间隔抽样序列，则信号奴/7)的频谱 函数是信号的频谱函数的周期化；若离散周期信号iOO是离散非周期信号奴/》的周 期化，则信号的频谱函数是信号的频谱函数的离散化。 信号在时域的离散化导致其频谱函数的周期化  信号在时域的周期化导致其频谱函数的离散化▲ TO通过离散周期序列的频谱可以近似分析连续非周期信号的频谱利用DH可以直接计算离散周期信号的频谱 DFT: X(k) = ^x(n)-e /Vz/=0々=0 信号在时域的周期化导致其频谱函数的离散化 ▲ TO 通过离散周期序列的频谱可以近似分析连续非周期信号的频谱 利用DH可以直接计算离散周期信号的频谱 DFT: X(k) = ^x(n)-e /V z/=0 々=0，1，…，TV —1 1DFT: 1 AM :⑻雌） k=Q .2jt kn Mat lab中提供了 fft函数，FFT是DFT的快速算法 X=fft (x):用于计算序列x的离散傅里叶变换(DFT) X=fft (x, n):对序列x补零或截短至n点的离散傅里叶变换。 当x的长度小于n时，在x的尾部补零使x的长度达到n点; 当x的长度大于n时，将x截短使x的长度成n点； x=ifft(X)和x=ifft(X, n)是相应的离散傅里叶反变换。 fftshift(x)将fft汁算输出的零频移到输出的中心位置。 利用DFT计算离散周期信号x(n)的频谱 分析步骤为： 确定离散周期序列x(zi)的基本周期光 利用fft函数求其一个周期的DFT，得到1(A); X(/d为％(幻的一个周期。 己知一个周期序列1(/7)二cos(—Z2 +—)+ 0.5cos(——n),利用FFT计算其频谱并与理论分 8 3 8 析相比较。 N=16; n=0:N-l; x=cos (n. *pi. /8+pi/3) +0. 5*cos (0. 875*pi. *n); X=fftshift(fft(x)); subplot (2,1，1); stem (n-N/2, abs (X)); title (’幅频特性’)； subplot (2,1, 2): stem(n-N/2, angle(X)); ^门以’相频特性’)； 幅频特性 III领特性 - A * ( ■ * o .T ? ? ■ 7 口 ■ a r 6 1 A ? 1 1 口 -8 -6 -4 -2 0 2 4 6 8 利用DFT计算离散非周期信号Yn)的频谱_ 分析步骤为： 确定序列的长度#当序列力无限长时，需要根据能量分布，利用窗函数进行截短。 确定作FFT的点数;V;根椐频域取样定理，为使时域波形不产生混叠，必须取;V 2 #。 使用fft函数作#点FFT计算J(A)。 利川DFT分析序列x(n) = Q.Sflu(n)的频谱。 k=0:30; x=0. 8. k; subplot (2, 1, 1): stem(k, x) ; %画出序列的时域波形 subplot (2, 1, 2); w=k~15; plot (w, abs (fftshi ft (fft (x)))) ; %画出序列频谱的幅度谱 二、实验原理一DFT实现线性卷积 x(n) (0nN}-l)有限长(7V1 点) /?() (0NN2-\)有限长(/V2 点) 当