资源描述：
《数字信号处理实验二.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、实用文档实验报告（本科）学号2015141443002 姓名柏冲 专业通信工程 日期2017/12/4 实验题目时域采样和频域采样一、实验目的实用文档时域采样理论与频域采样理论是数字信号处理中重要的理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使得采样后的信号不丢失信息；要求掌握频率采样会引起时域周期化的概念，以及频域采样定理及其对频域采样点数选择的指导作用。二、实验过程附：源程序(1)时域采样Tp=64/1000;%观察时间Tp=64毫秒%产生M长采样序列x(n)Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.12
2、8;a=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-a*n*T).*sin(omega*n*T);Xk=T*fft(xnt,M);%M点FFT[(xnt)]subplot(3,2,1);stem(xnt,'.');%调用编绘图函数stem绘制序列图boxon;title('(a)Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(3,2,2);stem(fk,abs(Xk),'.');title('(a)T*FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度
3、');axis([0,Fs,0,1.2*max(abs(Xk))]);%Fs=300Hz和Fs=200Hz的程序与上面Fs=1000Hz完全相同。Tp=64/1000;%观察时间Tp=64毫秒%产生M长采样序列x(n)Fs=300;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;a=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp(-a*n*T).*sin(omega*n*T);M1=fix(M);Xk=T*fft(xnt,M1);%M点FFT[(xnt)]subplot(3,2,3);stem(xnt,'.');
4、%调用自编绘图函数stem绘制序列图boxon;title('(b)Fs=300Hz');k=0:M-1;fk=k/Tp;subplot(3,2,4);stem(fk,abs(Xk),'.');title('(b)T*FT[xa(nT)],Fs=300Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]);Tp=64/1000;%观察时间Tp=64毫秒%产生M长采样序列x(n)Fs=200;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;a=pi*50*2^0.5;om
5、ega=pi*50*2^0.5;xnt=A*exp(-a*n*T).*sin(omega*n*T);M2=fix(M);实用文档Xk=T*fft(xnt,M2);%M点FFT[(xnt)]yn='xa(nT)';subplot(3,2,5);stem(xnt,'.');%调用自编绘图函数stem绘制序列图boxon;title('(c)Fs=200Hz');k=0:M-1;fk=k/Tp;subplot(3,2,6);stem(fk,abs(Xk),'.');title('(c)T*FT[xa(nT)],Fs=200Hz');xlabel('f(Hz)');yl
6、abel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))]);程序运行结果（2）频域采样M=27;N=32;n=0:M-1;%产生M长三角波序列x(n)xa=1:ceil(M/2);xb=fix(M/2):-1:1;xn=[xa,xb];Xk=fft(xn,1024);%1024点FFT[x(n)],用于近似序列x(n)的FtX32k=fft(xn,32);%32点FFT[x(n)]x32n=ifft(X32k);%32点IFFT[X32(k)]得到x32(n)X16k=X32k(1:2:N);%隔点抽取X32k得到X16(K)x16n=
7、ifft(X16k,N/2);%16点IFFT[X16(k)]得到x16(n)subplot(3,2,2);stem(n,xn,'.');boxontitle('(b)三角波序列x(n)');xlabel('n');ylabel('x(n)');axis([0,32,0,20])实用文档k=0:1023;wk=2*k/1024;%subplot(3,2,1);plot(wk,abs(Xk));title('(a)FT[x(n)]');xlabel('omega/pi');ylabel('
8、X(e^j^omega)
9、');axis([0,1,0,200])k
10、=0:N/2-1;sub
