时域采样与频域采样实验报告

实验二 时域采样与频域采样一 实验内容1 时域采样定理的验证给定模拟信号0()sin()()t a x t Ae t u t α-=Ω，式中，A=444.128，α=，0/rad s Ω=选取三种采样频率，即1s F kH z =，300Hz ，200Hz ，对()a x t 进行理想采样，得到采样序列：0()()sin()()nT a x n x nT Ae nT u nT α-==Ω。

观测时间长度为64p T m s =。

分别绘出三种采样频率得到的序列的幅频特性曲线图，并进行比较。

注：为与课本中幅频特性曲线比较，将纵坐标进行了归一化。

实验结果：由实验结果发现，采样频率为1000HZ 时，时域采样后的频谱函数可以较好的表现出原模拟信号的幅频特性，且是原幅频特性的周期延拓。

当采样频率为300HZ和200HZ时，其频谱函数与原幅频特性相比，有较大的误差，且在fs/2的位置误差最大。

实验分析：理想采样信号的频谱是原模拟信号的频谱沿频率轴，每间隔采样角频率2*pi*fs重复出现一次，并叠加形成的周期函数，所以只有当采样角频率2*pi*fs大于等于原模拟信号的角频率时才不会发生混叠。

2 频域采样定理的验证给定信号：1013()271426n nx n n nothers+≤≤⎧⎪=-≤≤⎨⎪⎩，对()x n的频谱函数()jX eω在[0，2π]上分别等间隔采样16点和32点，得到16()X k和32()X k，再分别对16()X k和32()X k进行IDFT，得到16()x n和32()x n。

分别画出()jX eω、16()X k和32()X k的幅度谱，并绘图显示()x n、16()x n和32()x n的波形，进行对比和分析。

实验结论：由上图分析知，频域采样32点时，其逆变换得到的xn32能较好的还原xn，只是尾部多了几个0而已，而对于频域采样16点时，逆变换之后已经产生较大的误差，不能等效为xn。

实验二 时域采样与频域采样学校:西南大学 班级:通信工程班一、实验目的时域采样理论与频域采样理论就是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化,以及如何选择采样频率才能使采样后的信号不丢失信息;要求掌握频率域采样会引起时域周期化的概念,以及频率域采样定理及其对频域采样点数选择的指导作用。

二、实验原理时域采样定理的要点就是采样频率s Ω必须大于等于模拟信号最高频率的两倍以上,才能使采样信号的频谱不产生频谱混叠。

频域采样定理的要点就是:a) 对信号x(n)的频谱函数X(e j ω)在[0,2π]上等间隔采样N 点,得到2()() , 0,1,2,,1j N k N X k X e k N ωπω===-则N 点IDFT[()N X k ]得到的序列就就是原序列x(n)以N 为周期进行周期延拓后的主值区序列,公式为 ()IDFT[()][()]()N N N N i x n X k x n iN R n ∞=-∞==+∑b) 由上式可知,频域采样点数N 必须大于等于时域离散信号的长度M(即N ≥M),才能使时域不产生混叠,则N 点IDFT[()N X k ]得到的序列()N x n 就就是原序列x(n),即()N x n =x(n)。

如果N>M,()N x n 比原序列尾部多N-M 个零点;如果N<M,z 则()N x n =IDFT[()N X k ]发生了时域混叠失真,而且()N x n 的长度N 也比x(n)的长度M 短,因此。

()N x n 与x(n)不相同。

三、实验程序(1)时域采样理论的验证。

Tp=64/1000;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444、128;alph=pi*50*2^0、5;omega=pi*50*2^0、5;xnt=A*exp(-alph*n*T)、*sin(omega*n*T);Xk=T*fft(xnt,M);yn='xa(nT)';subplot(3,2,1);tstem(xnt,yn);box on;title('(a) Fs=1000Hz');k=0:M-1;fk=k/Tp;subplot(3,2,2);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1、2*max(abs(Xk))])(Fs=300Hz与Fs=200Hz的程序与上面Fs=1000Hz完全相同。

数字信号处理及实验实验报告实验题目时域抽样与频域抽样姓名MYT 组别班级学号【实验目的】加深理解连续时间信号离散化过程中的数学概念和物理概念，掌握时域抽样定理的基本内容。

掌握由抽样序列重建原连续信号的基本原理与实现方法，理解其工程概念。

加深理解频谱离散化过程中的数学概念和物理概念，掌握频域抽样定理的基本内容。

【实验原理】离散系统在处理信号时，信号必须是离散的序列。

因此，在利用计算机等离散系统分析处理连续时间信号时必须对信号进行离散化处理。

时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件；对于基带信号，信号抽样频率f大于等于2倍的信号最高频率fm，即 f ≥ fm 。

信号的重建使信号抽样的逆过程。

非周期离散信号的频谱是连续的。

计算机在分析离散信号的频谱时，必须将其连续频谱离散化。

频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。

【实验结果与数据处理】1、为了观察连续信号时域抽样时抽样频率对抽样过程的影响，在[0，0.1]区间上以50HZ的抽样频率对下列3个信号分别进行抽样，试画出抽样后序列的波形，并分析产生不同波形的原因，提出改进措施。

（1）x1(t)=cos(2π*10t)（2）x2(t)=cos(2π*50t)（3）x3(t)=cos(2π*100t)程序代码如下:clc,clear,close allt0=0:0.001:0.1;Fs=50;t=0:1/Fs:0.1;figure(1)x1=cos(2*pi*10*t0);plot(t0,x1,'r')hold onx=cos(2*pi*10*t);stem(t,x);hold offfigure(2)x2=cos(2*pi*50*t0);plot(t0,x2,'r')hold onx=cos(2*pi*50*t);stem(t,x);hold offfigure(3)x3=cos(2*pi*100*t0);plot(t0,x3,'r')hold onx=cos(2*pi*100*t);stem(t,x);hold off图 1 x1(t)=cos(2π*10t)图 2 x2(t)=cos(2π*50t)图 3 x3(t)=cos(2π*100t)2、产生幅度调制信号X(t)=cos(2πt)cos（200πt），推导其频率特性，确定抽样频率，并绘出波形。

时域采样与频域采样实验报告一、实验目的：1.理解采样定理的原理和应用；2.掌握时域采样和频域采样的方法和步骤；3.学习使用MATLAB软件进行采样信号的分析和处理。

二、实验原理：采样是指将连续时间信号转换为离散时间信号的过程。

采样过程中，时间轴被分成若干个时间间隔，每个时间间隔内只有一个采样值，即取样点，采样信号的幅度就是该时间间隔内对应连续时间信号的幅度，称为采样值。

时域采样：利用采样定理进行抽样，采样时将模拟信号保持在一个固定状态下，以等间隔时间取样，实现模拟信号的离散化。

时域采样的反变换为恢复成为原连续时间信号，称为重构。

在数字信号中，通过离散时间信号构建模拟信号。

频域采样：首先通过傅里叶变换将时域信号转换到频域，然后在频域对其进行采样，将频域采样结果再进行反傅里叶变换恢复成时域信号。

三、实验内容及步骤：1.时域采样实验①模拟信号的采样：在MATLAB软件中设计一个三角波信号和正弦波信号，并画出其时域图像。

分别设定采样频率为1.5kHz和3kHz，进行采样。

重构时域信号，并画出重构信号的时域图像。

比较原信号和重构信号，在时域和频域上进行对比和分析。

②数字信号的量化：对采集的信号进行量化处理，设量化步长分别为1、2、3。

计算量化误差和信噪比，并作图进行比较分析。

2.频域采样实验设计一个具有3kHz频率的信号，并绘制其频域图像。

设定采样率为10kHz，进行采样，同时对采样信号进行降采样处理。

恢复实验所得到的采样信号，绘制重构后的时域图像，并分析其质量。

四、实验结果与分析：1.时域采样实验：①模拟信号的采样：通过MATLAB软件设计得到的三角波和正弦波信号及其时域图像如下所示：其中，Fs1 = 1.5kHz，Fs2 = 3kHz，信号的采样频率与信号频率的比值应大于2，以保证采样后的信号不失真。

通过采样得到的信号及其重构图像如下所示：可以看出，采样和重构得到的信号与原信号的时域图像是相似的，重构后的信号和原信号之间的误差可以忽略不计。

数字信号处理实验二时域采样和频域采样 Document serial number【KKGB-LBS98YT-BS8CB-BSUT-BST108】实验二-时域采样和频域采样一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

二、实验原理及方法1、时域采样定理的要点：a)对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X是原模拟信号频谱()a X j Ω以采样角频率s Ω（T s /2π=Ω）为周期进行周期延拓 b)采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。

2、频域采样定理的要点：a)对信号x(n)的频谱函数X(ej ω)在[0，2π]上等间隔采样N 点则N 点IDFT[()N X k ]得到的序列就是原序列x(n)以N 为周期进行周期延拓后的主值区序列。

三、实验内容及步骤1、时域采样理论的验证程序：clear;clcA=;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;Tp=50/1000;F1=1000;F2=300;F3=200;T1=1/F1;T2=1/F2;T3=1/F3;n1=0:Tp*F1-1;n2=0:Tp*F2-1;n3=0:Tp*F3-1;x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);f1=fft(x1,length(n1)); f2=fft(x2,length(n2)); % f3=fft(x3,length(n3)); %k1=0:length(f1)-1;fk1=k1/Tp; %k2=0:length(f2)-1;fk2=k2/Tp; %k3=0:length(f3)-1;fk3=k3/Tp; %subplot(3,2,1)stem(n1,x1,'.')title('(a)Fs=1000Hz');xlabel('n');ylabel('x1(n)');subplot(3,2,3)stem(n2,x2,'.')title('(b)Fs=300Hz');xlabel('n');ylabel('x2(n)');subplot(3,2,5)stem(n3,x3,'.')title('(c)Fs=200Hz');xlabel('n');ylabel('x3(n)');subplot(3,2,2)plot(fk1,abs(f1))title('(a) FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('·ùè')subplot(3,2,4)plot(fk2,abs(f2))title('(b) FT[xa(nT)],Fs=300Hz');xlabel('f(Hz)');ylabel('·ùè')subplot(3,2,6)plot(fk3,abs(f3))title('(c) FT[xa(nT)],Fs=200Hz');xlabel('f(Hz)');ylabel('·ùè')结果分析：由图可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。

实验二：时域采样与频域采样1. 实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

2. 实验原理与方法对模拟信号X a(t)以间隔T进行时域等间隔理想采样，形成的采样信号的频谱)?(j-)是原模拟信号频谱X a(j fl )以采样角频率门sC1 s二2 /T)为周期进行周期延拓。

公式为：fa(j「)二FT[X a(t)] X a(j「- jn 丄)T采样频率「I s必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

3. 实验内容及步骤%物联一班胡洪2%2015年10 月24 日澈验二：程序1Tp=64/1000;Fs=1000;T=1/Fs;M=ceil(Tp*Fs); n=0:M-1;A=444.128;a=pi*50*2A0.5;w=pi*50*2A0.5;xnt=A*exp(-a*n*T).*sin(w*n*T);Xk=fft(xnt,M);subplot(3,2,1);stem(n,xnt, '.' );axis([1,65,-5,150]);title( ' 图1 Fs=1000Hz' );subplot(3,2,2);plot(n/Tp,abs(Xk));title( '图2 Fs=1000Hz 幅度'); Fs=300;T=1/Fs;M=ceil(Tp*Fs);n=0:M-1;A=444.128;a=pi*50*2A0.5;w=pi*50*2A0.5;xnt=A*exp(-a*n*T).*sin(w*n*T);Xk=fft(xnt,M);subplot(3,2,3);stem(n,xnt, '.' );axis([0,M,-10,150])title( ' 图3 Fs=300Hz' );subplot(3,2,4);plot(n/Tp,abs(Xk));title( '图4 Fs=300Hz 幅度'); Fs=200;T=1/Fs;M=ceil(Tp*Fs);n=0:M-1;A=444.128;a=pi*50*2A0.5;w=pi*50*2A0.5;xnt=A*exp(-a*n*T).*sin(w*n*T);Xk=fft(xnt,M);subplot(3,2,5);图 6 Fs=200Hz 幅度');stem(n,xnt, '.' );axis([0,M,-10,150])title('图 5 Fs=200Hz');subplot(3,2,6);plot (n /Tp,abs(Xk));title(图1%物联一班胡洪2%2015年 10 月 24 日澈验二：程序2n=0:13;xa=n+1;n=14:26;xb=27-n;xn=[xa,xb]; n=0:26;subplot(3,2,2);stem(n,xn, '.');title('三角波序列 x(n)' );axis([0,32,0,15])Xk=fft(x n,1024);k=0:1023;wk=2*k/1024;subplot(3,2,1);plot(wk,abs(Xk));ylabel('|X(e A (j*w))|' );xlabel( 'w/ pi') Si P 沪lO&k W 3M SB 凹血 M Fffl) M TO lOBEliiF^3nQH £tt*axis([0,1,0,200]);Xk32=fft(xn,32);subplot(3,2,3);k=0:31;stem(k,abs(Xk32), '.' );axis([0,16,0,200]);xlabel( 'k' );ylabel( '| Xk3 2|' )xn32=ifft(Xk32);subplot(3,2,4);k=0:31;stem(k,xn32, '.' );axis([0,32,0,15]);Xk16=Xk32(1:2:32);subplot(3,2,5);k=0:15;stem(k,abs(Xk16), '.' );axis([0,8,0,200]);xlabel( 'k' );ylabel( '|X k16 ' )xn16=ifft(Xk16);subplot(3,2,6);k=0:15;stem(k,xn16, '.' );axis([0,32,0,15]);(j0.5w/pi200 r10010 15200 r1004•思考题:1.如果序列x (n)的长度为M希望得到其平铺X(e A(jw))在［0,2pi］上的N点等间隔采样，当N<M如何用一次最小的点数DFT得到该频谱采样？答：先对原序列x(n)以N为周期进行周期延拓后取主值区序列，再计算N 点DFT则得到N点频域采样。

一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论.要求掌握模拟信号采样前后频谱的变化,以与如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念,以与频率域采样定理与其对频域采样点数选择的指导作用.二、实验原理与方法时域采样定理的要点是：a.对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样,形成的采样信号的频谱)(ˆΩj X是原模拟信号频谱()a X j Ω以采样角频率s Ω〔T s /2π=Ω〕为周期进行周期延拓.公式为： )](ˆ[)(ˆt x FT j X a a =Ω )(1∑∞-∞=Ω-Ω=n s a jn j X T b.采样频率s Ω必须大于等于模拟信号最高频率的两倍以上,才能使采样信号的频谱不产生频谱混叠.利用计算机计算上式并不方便,下面我们导出另外一个公式,以便用计算机上进行实验.理想采样信号)(ˆt xa 和模拟信号)(t x a 之间的关系为： 对上式进行傅立叶变换,得到：在上式的积分号内只有当nT t =时,才有非零值,因此：上式中,在数值上)(nT x a ＝)(n x ,再将T Ω=ω代入,得到：上式的右边就是序列的傅立叶变换)(ωj e X ,即 T j a e X j X Ω==Ωωω)()(ˆ上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到,只要将自变量ω用T Ω代替即可.频域采样定理的要点是：a) 对信号x<n>的频谱函数X<e j ω>在[0,2π]上等间隔采样N 点,得到则N 点IDFT[()N X k ]得到的序列就是原序列x<n>以N 为周期进行周期延拓后的主值区序列,公式为：b) 由上式可知,频域采样点数N 必须大于等于时域离散信号的长度M<即N ≥M>,才能使时域不产生混叠,则N 点IDFT[()N X k ]得到的序列()N x n 就是原序列x<n>,即()N x n =x<n>.如果N>M,()N x n 比原序列尾部多N-M 个零点；如果N<M,z 则()N x n =IDFT[()N X k ]发生了时域混叠失真,而且()N x n 的长度N 也比x<n>的长度M 短,因此.()N x n 与x<n>不相同.在数字信号处理的应用中,只要涉与时域或者频域采样,都必须服从这两个采样理论的要点. 对比上面叙述的时域采样原理和频域采样原理,得到一个有用的结论,这两个采样理论具有对偶性："时域采样频谱周期延拓,频域采样时域信号周期延拓〞.因此放在一起进行实验.三、 实验内容与步骤〔1〕时域采样理论的验证给定模拟信号,)()sin()(0t u t Ae t x t a Ω=-α式中A =444.128,α=502π,0Ω=502πrad/s,它的幅频特性曲线如图图)(t x a 的幅频特性曲线现用DFT<FFT>求该模拟信号的幅频特性,以验证时域采样理论.安照)(t x a 的幅频特性曲线,选取三种采样频率,即s F =1k Hz ,300Hz ,200Hz .观测时间选ms T p 50=.为使用DFT,首先用下面公式产生时域离散信号,对三种采样频率,采样序列按顺序用)(1n x ,)(2n x ,)(3n x 表示.因为采样频率不同,得到的)(1n x ,)(2n x ,)(3n x 的长度不同, 长度〔点数〕用公式s p F T N ⨯=计算.选FFT 的变换点数为M=64,序列长度不够64的尾部加零.X <k >=FFT[x <n >] , k =0,1,2,3,-----,M -1式中k 代表的频率为 k Mk πω2=. 要求： 编写实验程序,计算)(1n x 、)(2n x 和)(3n x 的幅度特性,并绘图显示.观察分析频谱混叠失真.〔2〕频域采样理论的验证给定信号如下：编写程序分别对频谱函数()FT[()]j X e x n ω=在区间]2,0[π上等间隔采样32和16点,得到)()(1632k X k X 和：再分别对)()(1632k X k X 和进行32点和16点IFFT,得到)()(1632n x n x 和：分别画出()j X e ω、)()(1632k X k X 和的幅度谱,并绘图显示x <n>、)()(1632n x n x 和的波形,进行对比和分析,验证总结频域采样理论.提示：频域采样用以下方法容易变程序实现.1〕直接调用MATLAB函数fft计算就得到在的32点频率域采样2〕抽取的偶数点即可得到在的16点频率域采样,即.3〕当然也可以按照频域采样理论,先将信号x<n>以16为周期进行周期延拓,取其主值区〔16点〕,再对其进行16点DFT<FFT>,得到的就是在的16点频率域采样 .四、实验程序1 时域采样理论的验证程序清单% 时域采样理论验证程序exp2a.mTp=64/1000; %观察时间Tp=64微秒%产生M长采样序列x<n>% Fs=1000;T=1/Fs;Fs=1000;T=1/Fs;M=Tp*Fs;n=0:M-1;A=444.128;alph=pi*50*2^0.5;omega=pi*50*2^0.5;xnt=A*exp<-alph*n*T>.*sin<omega*n*T>;Xk=T*fft<xnt,M>; %M点FFT[xnt>]yn='xa<nT>';subplot<3,2,1>;tstem<xnt,yn>; %调用自编绘图函数tstem绘制序列图box on;title<'<a> Fs=1000Hz'>;k=0:M-1;fk=k/Tp;subplot<3,2,2>;plot<fk,abs<Xk>>;title<'<a> T*FT[xa<nT>],Fs=1000Hz'>;xlabel<'f<Hz>'>;ylabel<'幅度'>;axis<[0,Fs,0,1.2*max<abs<Xk>>]>%=================================================% Fs=300Hz和Fs=200Hz的程序与上面Fs=1000Hz完全相同.2 频域采样理论的验证程序清单%频域采样理论验证程序exp2b.mM=27;N=32;n=0:M;%产生M长三角波序列x<n>xa=0:floor<M/2>; xb= ceil<M/2>-1:-1:0; xn=[xa,xb];Xk=fft<xn,1024>; %1024点FFT[x<n>], 用于近似序列x<n>的TFX32k=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=ifft<X16k,N/2>; %16点IFFT[X16<k>]得到x16<n>subplot<3,2,2>;stem<n,xn,'.'>;box ontitle<'<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<'|X<e^j^\omega>|'>;axis<[0,1,0,200]>k=0:N/2-1;subplot<3,2,3>;stem<k,abs<X16k>,'.'>;box ontitle<'<c> 16点频域采样'>;xlabel<'k'>;ylabel<'|X_1_6<k>|'>;axis<[0,8,0,200]>n1=0:N/2-1;subplot<3,2,4>;stem<n1,x16n,'.'>;box ontitle<'<d> 16点IDFT[X_1_6<k>]'>;xlabel<'n'>;ylabel<'x_1_6<n>'>;axis<[0,32,0,20]>k=0:N-1;subplot<3,2,5>;stem<k,abs<X32k>,'.'>;box ontitle<'<e> 32点频域采样'>;xlabel<'k'>;ylabel<'|X_3_2<k>|'>;axis<[0,16,0,200]>n1=0:N-1;subplot<3,2,6>;stem<n1,x32n,'.'>;box ontitle<'<f> 32点IDFT[X_3_2<k>]'>;xlabel<'n'>;ylabel<'x_3_2<n>'>;axis<[0,32,0,20]>五、实验程序运行结果与分析1、时域采样理论的验证程序运行结果exp2a.m如图所示.由图可见,采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓.当采样频率为1000Hz时频谱混叠很小；当采样频率为300Hz时,在折叠频率150Hz附近频谱混叠很严重；当采样频率为200Hz时,在折叠频率110Hz附近频谱混叠更很严重.图2 时域采样理论的验证程序exp2b.m运行结果如图所示.图该图验证了频域采样理论和频域采样定理.对信号x<n>的频谱函数X<e jω>在[0,2π]上等间X k]得到的序列正是原序列x<n>以16为周期进行周期延拓隔采样N=16时,N点IDFT[()N后的主值区序列：x n与x<n>不相同,如图图<c>和<d>所示.当由于N<M,所以发生了时域混叠失真,因此.()NN=32时,如图图<c>和<d>所示,由于N>M,频域采样定理,所以不存在时域混叠失真,因x n与x<n>相同.此.()N。

10.2 实验二 时域采样与频域采样 10.2.1 实验指导1. 实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

2. 实验原理与方法时域采样定理的要点是：a) 对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X 是原模拟信号频谱()aX j Ω以采样角频率s Ω（T s /2π=Ω）为周期进行周期延拓。

公式为：)](ˆ[)(ˆt xFT j X a a =Ω )(1∑∞-∞=Ω-Ω=n s a jn j X T b) 采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。

理想采样信号)(ˆt xa 和模拟信号)(t x a 之间的关系为： ∑∞-∞=-=n a a nT t t x t x)()()(ˆδ对上式进行傅立叶变换，得到：dt e nT t t x j X t j n a a Ω-∞∞-∞-∞=⎰∑-=Ω])()([)(ˆδ dt e nT t t x t j n a Ω-∞-∞=∞∞-∑⎰-)()( δ＝在上式的积分号内只有当nT t =时，才有非零值，因此：∑∞-∞=Ω-=Ωn nT j aae nT xj X )()(ˆ上式中，在数值上)(nT x a ＝)(n x ，再将T Ω=ω代入，得到：∑∞-∞=-=Ωn nj aen x j X ω)()(ˆ上式的右边就是序列的傅立叶变换)(ωj eX ，即T j a e X j X Ω==Ωωω)()(ˆ上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用T Ω代替即可。

实验二-时域采样和频域采样一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。

要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。

二、实验原理及方法1、时域采样定理的要点：a)对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X是原模拟信号频谱()a X j Ω以采样角频率s Ω（T s /2π=Ω）为周期进行周期延拓 b)采样频率s Ω必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。

利用计算机计算上式并不方便，下面我们导出另外一个公式，以便用计算机上进行实验。

2、频域采样定理的要点：a)对信号x(n)的频谱函数X(ej ω)在[0，2π]上等间隔采样N 点则N 点IDFT[()N X k ]得到的序列就是原序列x(n)以N 为周期进行周期延拓后的主值区序列。

三、实验内容及步骤1、时域采样理论的验证程序：clear;clcA=;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;Tp=50/1000;F1=1000;F2=300;F3=200;T1=1/F1;T2=1/F2;T3=1/F3;n1=0:Tp*F1-1;n2=0:Tp*F2-1;n3=0:Tp*F3-1;x1=A*exp(-a*n1*T1).*sin(w0*n1*T1);x2=A*exp(-a*n2*T2).*sin(w0*n2*T2);x3=A*exp(-a*n3*T3).*sin(w0*n3*T3);f1=fft(x1,length(n1));f2=fft(x2,length(n2));%f3=fft(x3,length(n3));%k1=0:length(f1)-1;fk1=k1/Tp;%k2=0:length(f2)-1;fk2=k2/Tp;%k3=0:length(f3)-1;fk3=k3/Tp;%subplot(3,2,1)stem(n1,x1,'.')title('(a)Fs=1000Hz');xlabel('n');ylabel('x1(n)');subplot(3,2,3)stem(n2,x2,'.')title('(b)Fs=300Hz');xlabel('n');ylabel('x2(n)');subplot(3,2,5)stem(n3,x3,'.')title('(c)Fs=200Hz');xlabel('n');ylabel('x3(n)');subplot(3,2,2)plot(fk1,abs(f1))title('(a)FT[xa(nT)],Fs=1000Hz');xlabel('f(Hz)');ylabel('·ù?è')subplot(3,2,4)plot(fk2,abs(f2))title('(b)FT[xa(nT)],Fs=300Hz');xlabel('f(Hz)');ylabel('·ù?è')subplot(3,2,6)plot(fk3,abs(f3))title('(c)FT[xa(nT)],Fs=200Hz');xlabel('f(Hz)');ylabel('·ù?è')结果分析：由图可见，采样序列的频谱的确是以采样频率为周期对模拟信号频谱的周期延拓。