通信 专业 ——dsp 实验五 频域.docVIP

实验五 频域抽样定理 牛聪 学号：200900121120 通信三班 实验内容 给定信号 1．计算信号的频谱，画出频谱图，标明坐标轴。 2．分别对信号的频谱等间隔抽样16点和32点，离散傅里叶反变换后得到时域信号。 3．画出相关信号的图形，计算均方误差。 4. 对比和分析，验证频域抽样定理。 编程原理、思路和公式 1）、首先表示出x(n)，画出其波形，并对其进行DTFT变换。 2）、将频谱分别等间隔的抽样16/32个点，画出其抽样后的频谱，然后分别进行IDTF变换. 3）、将两次抽样后的频谱和反变换得到的时域信号用图形表示出来，并计算反变换得到的时域信号与原信号的均方误差。 程序脚本，并注释 N=27; n=0:N-1; xn=(n+1).*(n=0n=13)+(27-n).*(n=14n=26); subplot(3,2,2);stem(n,xn,.);grid on %画出xn的图形 title((b) 序列x(n));xlabel(n);ylabel(x(n)); M=1024;i=0:M-1;w=2*pi*i/M; %确定数字角频率的范围 N=27;n=0:N-1; Xw=xn*exp(-j*n*w); Xw_abs=abs(Xw); %画出xn的频谱 subplot(3,2,1);plot(w,Xw_abs);grid on title((a)频域幅度谱); xlabel(\omega/\pi);ylabel(|X(e^j^\omega)|); N=16;k=0:N-1;m=M/N; %确定频谱抽样的间隔16及抽样点数 Xk1=Xw(m*k+1); %16点频谱抽样 Xk1_abs=abs(Xk1); Xk1_angle=angle(Xk1); subplot(3,2,3);stem(k,Xk1_abs,.);grid on title((c) 16点频域采样);xlabel(k);ylabel(|X_1_6(k)|); w1=w(m*k+1); %频谱抽样后各点对应的数字角频率 x16n=Xk1*exp(j*w1*n)./N; %将抽样后的频谱IDFT变换 subplot(3,2,4);stem(n,x16n,.);grid on title((d) 16点IDFT[X_1_6(k)]);xlabel(n);ylabel(x_1_6(n)); mes1=sum(abs(x16n-xn).^2)/27 %计算原信号与变换信号的均方误差 N=32;k=0:N-1;m=M/N; %确定频谱抽样的间隔16及抽样点数 Xk2=Xw(m*k+1); %16点频谱抽样 Xk2_abs=abs(Xk2); subplot(3,2,5);stem(k,Xk2_abs,.);grid on title((e) 32点频域采样);xlabel(k);ylabel(|X_3_2(k)|); w2=w(m*k+1); %频谱抽样后各点对应的数字角频率 x32n=Xk2*exp(j*w2*n)./N; %将抽样后的频谱IDFT变换 subplot(3,2,6);stem(n,x32n,.);grid on title((f) 32点IDFT[X_3_2(k)]);xlabel(n);ylabel(x_3_2(n)); mes2=sum(abs(x32n-xn).^2)/27 %计算原信号与变换信号的均方误差 仿真结果、图形 16点频域抽样后恢复的x16n与xn的均方误差 mes1 = 37.4815 32点频域抽样后恢复的x32n与xn的均方误差 mes2 = 4.7446e-028 结果分析和结论 当频谱抽样点数为16时，不能还原出原信号xn，16点抽样并反变换后的信号是原信号以16为周期延拓而来的，信号之间有交叠部分，所以不能还原成原信号。而当频谱抽样点数为32时，则可以还原出原信号。 频域抽样会造成时域的周期严拓。 对于M点的有限长序列，频域抽样不失真的条件是频域抽样点数N要大于或等于M，即满足 N=M 当NM时，x(n)以N为周期进行周期严拓，就会造成混叠，不能不失真地恢复出原信号x(n)来。 从均方误差也可以看出，32点抽样还原的信号与原信号的均方误差小于16点抽样还原的信号与原信号的均方误差，说明信号与原信号差别很小，信号复原的很好。 遇到的问题、解决方法及收获 基于频域采样定理的论证过程我们了解到了频域采样理论的要点，即频域采样点数N必须大于等于时域离散信号的长度M(