实验二 时域采样与频域采样打印.doc

HIGER海格实验二 时域采样与频域采样一、实验程序清单及曲线1、 时域采样定理的验证程序：clc;A=444.128;alph=50*2(0.5)*pi;w=50*2(0.5)*pi;tp=0.064;fs=1/1000; %采样周期是1000HZT=1/fs;m=tp*T;n=0:m-1;xn1=A*exp(-alph*n*fs).*sin(w*n*fs);subplot(3,2,1);stem(n,xn1,.);title(xn1)xlabel(n);ylabel(xn1);xk1=fft(xn1,m).*T;k=0:2*pi/m:2*pi/m*(m-1);subplot(3,2,2);plot(k,abs(xk1); %画出幅频特性title(xk1);xlabel(w);ylabel(xk1);fs=1/300; %采样周期是300HZT=1/fs;m=tp*T;n=0:m-1;xn2=A*exp(-alph*n*fs).*sin(w*n*fs);subplot(3,2,3);stem(n,xn2,.);title(xn2);xlabel(n);ylabel(xn2);xk2=fft(xn2,m).*T;k=0:2*pi/m:2*pi/m*(m-1); %归一化处理subplot(3,2,4);plot(k,abs(xk2); %画出幅频特性title(xk2);xlabel(w);ylabel(xk2);fs=1/200; %采样周期是200HZT=1/fs;m=tp*T;n=0:m-1;xn3=A*exp(-alph*n*fs).*sin(w*n*fs);subplot(3,2,5);stem(n,xn3,.);title(xn3)xlabel(n);ylabel(xn3);xk3=fft(xn3,m).*Tk=0:2*pi/m:2*pi/m*(m-1);subplot(3,2,6);plot(k,abs(xk3); %画出幅频特性title(xk3);xlabel(w);ylabel(xk3); 得到的图形：2、频域采样定理的验证程序：clc;xn1=1:14;n=0:13;subplot(4,2,1);stem(n,xn1,.);title(xn1);xlabel(n);ylabel(xn1);xn2=13:-1:1;n=14:26;subplot(4,2,2);stem(n,xn2,.);title(xn2);xlabel(n);ylabel(xn2);xn=xn1,xn2;n=0:length(xn)-1;subplot(4,2,3);stem(n,xn,.);title(a) xn三角波序列);xlabel(n);ylabel(xn);xk32=fft(xn,32); %32点FFTx(n)xn32=ifft(xk32); %32点IFFTX32(k)得到x32(n)xk16=xk32(1:2:32); %隔点抽取X32(k)得到X16(k)xn16=ifft(xk16); %16点IFFTX16(k)得到x16(n)subplot(4,2,4);k=0:15;wk=2*k/16;plot(wk,abs(xk16);title(b) FTx(n);xlabel(omega/pi);ylabel(|X(ejomega)|);subplot(4,2,5);k=0:15;stem(k,abs(xk16),.);title(c) 16点频率采样);xlabel(k);ylabel(|x_1_6(k)|);axis(0,32,0,200);subplot(4,2,6);n1=0:15;stem(n1,xn16,.);title(d) 16点IDFTx_1_6(k);xlabel(k);ylabel(|x_1_6(n)|);axis(0,32,0,20);subplot(4,2,7);k=0:31;stem(k,abs(xk32),.);title(e) 32点频率采样);axis(0,32,0,200);xlabel(k);ylabel(|x_3_2(k)|);subplot(4,2,8);n1=0:31;stem(n1,xn32,.);title(f) 32点IDFTx_3_2(k);xlabel(k);ylabel(|x_3_2(n)|);axis(0,32,0,20);图形： 结论：该图验证了频域采样理论和频域采样定理。对信号x(n)的频谱函数X(ej)在0，2上等间隔采样N=16时， N点IDFT得到的序列正是原序列x(n)以16为周期进行周期延拓后的主值区序列：由于NM，满足频域采样定理，所以不存在时域混叠失真，因此。与x(n)相同。二、思考