北京邮电大学_dsp_matlab实验报告卢跃凯.docVIP

数字信号处理实验报告 PAGE \* MERGEFORMAT8 数字信号处理实验报告 实验名称： 数字信号处理实验二 学生姓名： 班 级： 班内序号： 1.实验要求 数字信号处理的核心算法是离散傅立叶变换(DFT)，是DFT使信号在数字域和频域都实现了离散化，从而可以用 HYPERLINK /view/1472448.htm \t _blank 通用计算机处理离散信号。而使数字信号处理从理论走向实用的是 HYPERLINK /view/1006229.htm \t _blank 快速傅立叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。 本次试验有两个内容: 1、用Matlab编程上机练习。已知： N=25。这里Q=0.9+j0.3。可以推导出 ， 首先根据这个式子计算X(k)的理论值，然后计算输入序列x(n)的32个值，再利用基2时间抽选的FFT算法，计算x(n)的DFT X(k)，与X(k)的理论值比较（要求计算结果最少6位有效数字）。 2、假设信号 x(n) 由下述信号组成： 请选择合适的长度 N 和窗函数，用 DFT 分析其频谱，得到清楚的三根谱线。 2.实验代码和实验结果 实验一： 代码 format long; %显示15位小数 Q=0.9+0.3i;n=0:24;f=Q.^n; %构建25点序列 y1=(1-Q^25)./(1-Q.*exp(-j*2*pi*n/25)); %用公式计算25点DFT f1=[f,zeros(1,7)]; %补0构建32点序列 y2=fft(f1,32); %使用基2FFT算法计算 subplot(2,1,1);stem(n,abs(y1));title(N=25,理论值);xlabel(n); n2=0:31;subplot(2,1,2);stem(n2,abs(y2));title(N=32, FFT);xlabel(n); %显示两组结果，进行对比 f %25点序列 ans = 1.000000000000000 + 0.000000000000000i 0.900000000000000 - 0.300000000000000i 0.720000000000000 - 0.540000000000000i 0.486000000000000 - 0.702000000000000i 0.226800000000000 - 0.777600000000000i -0.029160000000000 - 0.767880000000000i -0.256608000000000 - 0.682344000000000i -0.435650400000000 - 0.537127200000000i -0.553223520000000 - 0.352719360000000i -0.603716976000000 - 0.151480368000000i -0.588789388800000 + 0.044782761600000i -0.516475621440000 + 0.216941302080000i -0.399745668672000 + 0.350189858304000i -0.254714144313600 + 0.435094573075200i -0.098714357959680 + 0.467999359061760i 0.051556885554816 + 0.450813730543488i 0.181645316162381 + 0.390265291822695i 0.280560372092951 + 0.296745167791711i 0.341527885221169 + 0.182902539384655i 0.362245858514449 + 0.062153919879838i 0.344667448626955 - 0.052735229662480i 0.294380134865516 - 0.150861941284319i 0.219683538993669 - 0.224089787615542i 0.130488248809639 - 0.267585870552088i 0.0