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

《北京邮电大学_dsp_matlab实验报告卢跃凯》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、数字信号处理实验报告数字信号处理实验报告实验名称：数字信号处理实验二学生姓名：班级：班内序号：1.实验要求数字信号处理的核心算法是离散傅立叶变换(DFT)，是DFT使信号在数字域和频域都实现了离散化，从而可以用通用计算机处理离散信号。而使数字信号处理从理论走向实用的是快速傅立叶变换(FFT)，FFT的出现大大减少了DFT的运算量，使实时的数字信号处理成为可能、极大促进了该学科的发展。本次试验有两个内容:1、用Matlab编程上机练习。已知：N=25。这里Q=0.9+j0.3。可以推导出，首先根据这个式子计算X(k)的理论值，然后计算输入序列x(n)的32个值，再利用基2时间抽选的FFT算法，计

2、算x(n)的DFTX(k)，与X(k)的理论值比较（要求计算结果最少6位有效数字）。2、假设信号x(n)由下述信号组成：请选择合适的长度N和窗函数，用DFT分析其频谱，得到清楚的三根谱线。2.实验代码和实验结果实验一：代码>>formatlong;%显示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点DFTf1=[f,zeros(1,7)];%补0构建32点序列y2=fft(f1,32);%使用基2FFT算法计算subplot(2,1,1);stem(n,abs(y1));

3、title('N=25,理论值');xlabel('n');n2=0:31;subplot(2,1,2);stem(n2,abs(y2));title('N=32,FFT');xlabel('n');8数字信号处理实验报告%显示两组结果，进行对比>>f'%25点序列ans=1.000000000000000+0.000000000000000i0.900000000000000-0.300000000000000i0.720000000000000-0.540000000000000i0.486000000000000-0.702000000000000i0.226800000000000-0

4、.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.350

5、189858304000i-0.254714144313600+0.435094573075200i-0.098714357959680+0.467999359061760i0.051556885554816+0.450813730543488i0.181645316162381+0.390265291822695i0.280560372092951+0.296745167791711i0.341527885221169+0.182902539384655i0.362245858514449+0.062153919879838i0.344667448626955-0.0527352296624

6、80i0.294380134865516-0.150861941284319i0.219683538993669-0.224089787615542i0.130488248809639-0.267585870552088i0.037163662763049-0.279973758139771i>>y1'%公式计算结果ans=1.839921274419313-2.888509011710853i10.079297058328029-7.634101810228918i0.751947550413978+5.860554178694223i0.231477875261307+2.56282253

7、4924723i0.280180444813830+1.629731693733070i0.333453299816167+1.185686071198745i0.374735699173569+0.920169346971339i0.407012287168865+0.738970117586464i0.433274806079868+0.603824818398893i8数字信号处理实验报告0