数字信号处理实验四离散傅里叶变换.docxVIP

一、实验目的： 加深对DFT性质的理解，拓展它们在DSP中的使用； 二、实验内容：实验原理： DFT的快速算法FFT利用了的三个固有特性：（1）对称性，，（2）周期性，，（3）可约性，和。FFT算法基本上可以分为两大类，即按时间抽选法(DIT,Decimation-In-Time)和按频率抽选法(DIF,Decimation-In-frequency)。 MATLAB中提供了进行快速傅里叶变换的fft函数：X=fft(x)，基2时间抽取FFT算法，x是表示离散信号的向量；X是系数向量；X=fft(x,N),补零或截断的N点DFT，当x得长度小于N时，对补零使其长度为N，当x的长度大于N时，对x截断使其长度为N。1 已知连续周期信号x(t)=cos(10πt)+2sin(18πt) (1)确定信号的基频Ω和基本周期Tp，以及分析时采用的采样点数N； (2)当分析长度取0.5Tp和1.5Tp时，对x(t)采样，利用FFT计算其幅度谱；对所得结果进行比较，总结应如何选取分析长度。2、设x（n）=R8(n)，分别计算X（ejw）在[0,2π]上的32点和64点等间隔采样，并绘制幅频和相频特性图。三、实验结果：1 已知连续周期信号x(t)=cos(10πt)+2sin(18πt) (1)确定信号的基频Ω和基本周期Tp，以及分析时采用的采样点数N； (2)当分析长度取0.5Tp和1.5Tp时，对x(t)采样，利用FFT计算其幅度谱；对所得结果进行比较，总结应如何选取分析长度。 （1）基频Ω=2πrad/s，基本周期Tp=1s，在一周期内有采样定理知采样频率应大于ω=36π，即fΩ18Hz。则定一周期内采样点数N=20。0.5T和1.5T内采样点数分别为10点和30点。 n=0:0.0025:2; x=cos(10*pi*n)+2*sin(18*pi*n); plot(n,x) grid n1=0:0.05:0.45; x1=cos(10*pi*n1)+2*sin(18*pi*n1); subplot(221),stem(n1,x1); subplot(222),plot(n1,x1); xk1=fft(x1); subplot(223),stem(n1,xk1); n2=0:0.05:1.5; x2=cos(10*pi*n2)+2*sin(18*pi*n2); subplot(221),stem(n2,x2); subplot(222),plot(n2,x2); xk2=fft(x2); subplot(223),stem(n2,xk2);2、设x（n）=R8(n)，分别计算X（ejw）在[0,2π]上的32点和64点等间隔采样，并绘制幅频和相频特性图。 x=[ones(1,8) zeros(1,22)]; n=0:499; L=length(n); n=n.*2/L; xk=fft(x,L); subplot(311), stem(x); grid,title(x(n)=R8(n)); xlabel(n); subplot(312), plot(n,abs(xk)); grid,title(MAGNITUDE RESPONSE) xlabel(NORMALIZED FREQUENCY),ylabel(|H(w)|) subplot(313),plot(n,180/pi*angle(xk)); grid ,xlabel(NORMALIZED FREQUENCY),ylabel(DEGREES), title(PHASE RESPONSE) x=[ones(1,8) zeros(1,22)]; n=0:31; L=length(n); n=n.*2/L; xk1=fft(x,L); subplot(211), stem(n,abs(xk1)); grid,title(MAGNITUDE RESPONSE) xlabel(NORMALIZED FREQUENCY),ylabel(|H(w)|) subplot(212),stem(n,180/pi*angle(xk1)); grid ,xlabel(NORMALIZED FREQUENCY),ylabel(DEGREES), title(PHASE RESPONSE) x=[ones(1,8) zeros(1,22)]; n=0:63; L=length(n); n=n.*2/L; xk2=fft(x,L); subplot(211), stem(n,abs(xk2)); grid,title(MAGNITUDE RESPONSE) xlabel(NORMALIZED FREQUENCY),ylabel(|H(w)|) subplot(212),stem(n,180/pi*angle(