实验三-时域及频域采样定理
数字信号处理实验报告 3
数字信号处理实验报告姓名：班级：通信学号：实验名称：频域抽样定理验证实验类型：验证试验指导教师:实习日期：2013.频域采样定理验证实验一． 实验目的：1. 加深对离散序列频域抽样定理的理解2.了解由频谱通过IFFT 计算连续时间信号的方法3.掌握用MATLAB 语言进行频域抽样与恢复时程序的编写方法 4、用MATLAB 语言将X(k)恢复为X(z)及X(e jw )。
二． 实验原理：1、1、频域采样定理： 如果序列x(n)的长度为M ，频域抽样点数为N ，则只有当频域采样点数N ≥M 时，才有x N (n)=IDFT[X(k)]=x(n),即可由频域采样X(k)无失真的恢复原序列 x(n)。
2、用X(k)表示X(z)的内插公式：∑-=-----=10111)(1)(N k kNNzWz k X Nz X内插函数： zWzkNNN z 1k111)(-----=ϕ频域内插公式：∑-=-=10)2()()(N K j k Nk X e X πωϕω频域内插函数：e N j N N )21()2sin()2sin(1)(--=ωωωωϕ三． 实验任务与步骤：实验一：长度为26的三角形序列x(n)如图(b)所示，编写MATLAB 程序验证频域抽样定理。
实验二：已知一个时间序列的频谱为X(e jw )=2+4e -jw +6e -j2w +4e -j3w +2e -j4w分别取频域抽样点数N为3、5和10，用IPPT计算并求出其时间序列x(n)，用图形显示各时间序列。
由此讨论原时域信号不失真地由频域抽样恢复的条件。
实验三：由X32(k)恢复X(z)和X(e jw)。
四．实验结论与分析：实验一：源程序：M=26;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,512); %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:511;wk=2*k/512;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])结果如下所示：实验一分析：序列x(n)的长度M=26，由图中可以看出，当采样点数N=16<M时，x16(n)确实等于原三角序列x(n)以16为周期的周期延拓序列的主值序列。
实验三(时域抽样与频域抽样)
时域抽样与频域抽样
一、实验目的
1. 加深理解连续时间信号的离散化过程中的数学概念和物 理概念，掌握时域抽样定理（奈奎斯特采样定理）的基 本内容。 2. 加深对时域取样后信号频谱变化的认识。掌握由抽样序 列重建原连续信号的基本原理与实现方法，理解其工程 概念。 3. 加深理解频谱离散化过程中的数学概念和物理概念，掌 握频域抽样定理的基本内容。
二. 实验原理——时域抽样
1、时域抽样概述
时域抽样定理给出了连续信号抽样过程中信 号不失真的约束条件：信号抽样频率 fs 大于等于2 倍的信号最高频率fm，即 fs 2fm。
时域抽样先把连续信号 x(t)变成适合数字系统 处理的离散信号x[k]；然后根据抽样后的离散信号 x[k]恢复原始连续时间信号x(t)完成信号重建。
0
0.01
0.02
0.03
三. 程序示例——频域抽样
2.已知序列 x [ k ] {1, 1, 1; k 0,1,2} 对其频谱X(ej)进行抽样， 分别取N=2，3，10，观察频域抽样造成的混叠现象。
x=[1,1,1]; P=256; omega=[0:P-1]*2*pi/P; X0=1+exp(-j*omega)+exp(-2*j*omega); N=input('Type in N= '); omegam=[0:N-1]*2*pi/N; Xm=1+exp(-j*omegam)+exp(-2*j*omegam); subplot(2,1,1); plot(omega./pi,abs(X0)); xlabel('Omega/PI'); hold on stem(omegam./pi,abs(Xm),'r','o'); hold off x1=[zeros(1,2*N) x zeros(1,2*N)]; x2=[zeros(1,N) x zeros(1,3*N)]; x3=[x zeros(1,4*N)]; x4=[zeros(1,3*N) x zeros(1,N)]; x5=[zeros(1,4*N) x]; xx=x1+x2+x3+x4+x5; k=-2*N:2*N+length(x)-1; subplot(2,1,2); stem(k,x1); hold on subplot(2,1,2); stem(k,xx,'r','*'); hold off
实验三 时域采样与频域采样
实验二 时域采样与频域采样一 实验内容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。
DSP实验时域采样与频域采样
实验二时域采样与频域采样一． 实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。
要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。
二． 实验原理与方法时域采样定理的要点是：1) 对模拟信号)(t x a 以间隔T 进行时域等间隔理想采样，形成的采样信号的频谱)(ˆΩj X是原模拟信号频谱()a X j Ω以采样角频率s Ω（T s/2π=Ω）为周期进行周期延拓。
公式为：)](ˆ[)(ˆt xFT j X a a =Ω )(1∑∞-∞=Ω-Ω=n s a jn j X T 2) 采样频率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.） 对信号x(n)的频谱函数X(ej ω)在[0，2π]上等间隔采样N 点，得到2()(), 0,1,2,,1j N k NX 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 ∞=-∞==+∑2.） 由上式可知，频域采样点数N 必须大于等于时域离散信号的长度M(即N ≥M)，才能使时域不产生混叠，则N 点IDFT[()N X k ]得到的序列()N x n 就是原序列x(n),即()N x n =x(n)。
时域采样与频域采样
备注：（1）、按照要求独立完成实验内容。
（2）、实验结束后，把电子版实验报告按要求格式改名（例：09号_张三_实验七.doc ）后，实验室统一刻盘留档。
实验四 时域采样与频域采样一、实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。
掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样前后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对采样点数选择的指导作用。
二、实验原理在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。
时域采样原理和频域采样原理，得到一个有用的结论，这两个采样理论具有对偶性：“时域采样频谱周期延拓，频域采样时域信号周期延拓”。
因此放在一起进行实验。
三、实验内容（包括代码与产生的图形及结果分析）1. 给定模拟信号如下：xa(t)=Ae -αt sin(Ω0t)u(t)式中, A=444.128，α=50 π，Ω0=50 π rad/s ，将这些参数带入上式中，对x a (t进行傅里叶变换，它的幅频特性曲线如图1所示。
现用DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。
按照xa(t)的幅频特性曲线，选取三种采样频率，即Fs=1 kHz ，300 Hz ，200 Hz 。
观测时间选Tp=64 ms 。
要求： 编写实验程序，计算x 1(n)、 x 2(n)和x 3(n)的幅度特性，并绘图显示。
观察分析频谱混叠失真。
close all;clear all;clc;22图1 x a (t)的幅频特性曲线Tp=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)]subplot(3,2,1);n=0:length(xnt)-1;stem(n,xnt,'.');xlabel('n');ylabel('yn');axis([0,n(end),min(xnt),1.2*max(xnt)]);%绘图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;T=1/Fs;Fs=300;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)]subplot(3,2,3);n=0:length(xnt)-1;stem(n,xnt,'.');xlabel('n');ylabel('yn');axis([0,n(end),min(xnt),1.2*max(xnt)]);%绘图box on;title('(a) Fs=300Hz');k=0:M-1;fk=k/Tp;subplot(3,2,4);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=300Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])%========================% Fs=200Hz;T=1/Fs;Fs=200;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)]subplot(3,2,5);n=0:length(xnt)-1;stem(n,xnt,'.');xlabel('n');ylabel('yn');axis([0,n(end),min(xnt),1.2*max(xnt)]);%绘图box on;title('(a) Fs=200Hz');k=0:M-1;fk=k/Tp;subplot(3,2,6);plot(fk,abs(Xk));title('(a) T*FT[xa(nT)],Fs=200Hz');xlabel('f(Hz)');ylabel('幅度');axis([0,Fs,0,1.2*max(abs(Xk))])2. 频域采样理论的验证。
时域采样与频域采样
实验二：时域采样与频域采样一、实验目的：时域采样理论与频域采样理论是数字信号处理中的重要理论。
要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。
二、实验原理与方法：1、时域采样定理的要点：1）对模拟信号)(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 2）采样频率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 e X ，即T j a e X j X Ω==Ωωω)()(ˆ 上式说明理想采样信号的傅立叶变换可用相应的采样序列的傅立叶变换得到，只要将自变量ω用T Ω代替即可。
2、频域采样定理的要点：a) 对信号x(n)的频谱函数X(e j ω)在[0，2π]上等间隔采样N 点，得到2()() , 0,1,2,,1j N k NX 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)。
频域取样定理
2
t
时域抽样
时分复用
屏域抽样
频分复用
s
1
T s (t )
s
t
T s
f t
s
t t
m
m
T T 2T
s s
s
t t T
m m
s
2T s
t
即当频域取样的周期 s ,T s 2 t m时在时域中信号 f s t不重叠。 Ts 可用矩形脉冲低选通信号，恢复原信号，如图。 2 s 频域采样定理：时域信号必须为有限时间信号；频域抽样周期 Ts ， 只有当 T s 2t m时，才 可恢复原信号。
s
n
n 取样信号： F j n F jn n
s s n
F( j )
对应时域：
1 1
s
1
s
Hale Waihona Puke sn s
s
2
如果信号ft为有限时间信号时限信号其时间区间它的频谱函数f连续函数频谱在频域中对等间隔取样
频域取样定理
原理：时域与频域为对称的，可推出频域取样定理。 如果信号f(t)为有限时间信号(时限信号)，其时间区间 (t m , t m)，它的频谱函数F( j)为 连续函数频谱，在频域中对F( j ) 等间隔取样。间隔为 s ，取样脉冲为：
s
s
1
s n
1
s n
F j
s
s
2
s
s
f s t [ F s j ] [F j ] [ s ]
时、频域采样定理的验证
实验二时、频域采样定理的验证
实验目的
掌握时域采样定理，频域采样定理。
掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息；掌握频域采样会引起时域周期化，以及频域采样定理及其对频域采样点数选择的指导作用。
实验内容
1.时域采样理论
的验证
f=1000Hz:
f=300Hz:
f=200Hz:
2.频域采样理论的验证
实验小结
通过本次实验，我巩固了信号在matlab中的运算表示方法、图形输出函数（plot、stem），同时会用软件求fft和ifft，由此验证了时域采样定理，频域采样定理。
通过观察不同平率的模拟信号采样，采样频率如果过低会导致丢失信息；通过频域采样发现它会引起时域周期化。
时域采样与频域采样i
文电072-1班200790511115 张雪倩试验二：时域采样与频域采样实验目的时域采样理论与频域采样理论是数字信号处理中的重要理论。
掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样前后的信号不丢失信息；要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对采样点数选择的指导作用。
实验原理与方法频域采样的要点是：（1）对信号的频谱函数在上等间隔采样N点，得到：k=0,1,2,..,N-1则N点的得到的序列就是原序列以N为周期进行周期延拓后的主值区序列，公式为：（2）由上式可知，频域采样点数N必须大于等于时域离散信号长度M（即N），才能使时域不产生混叠，则N点得到的序列就是原序列，即。
如果N>M，比原序列尾部多N-M个零点；如果N<M，则发生了时域混叠失真，而且的长度N也比的长度M短，因此，与不相同。
在数字信号处理的应用中，只要涉及时域或者频域采样，都必须服从这两个采样理论的要点。
对此上面叙述的时域采样原理和频域采样原理，得到一个有用的结论：这两个采样理论具有对偶性，即“时域采样频谱周期延拓，频域采样时域信号周期延拓”。
实验内容与步骤（2）频域采样理论的验证。
给定信号如下：编写程序分别对频域函数在区间上等间隔采样32点和16点，得到和：在分别对和进行32点和16点IFFT，得到和:分别画出、和的幅度谱，并绘图显示、和的波形，进行对比和分析，验证总结频域采样理论。
用MA TLAB编写程序和输出图形如下：(1)xa=1:14;xb=13:-1:1;xn=[xa,xb];figure(1)subplot(211)stem(0:26,xn,'.');xk=fft(xn,1024);k=0:1023;wk=2*k/1024;subplot(212)plot(abs(xk));对xn序列进行傅里叶变换，采样1024点(2)figure(2)x16k=xk(1,1:64:1024);x16n=ifft(x16k,16);subplot(211);stem(0:15,x16n,'.');subplot(212)stem(0:15,abs(x16k),'.');(3)x32k=xk(1,1:32:1024);x32n=ifft(x32k,32);subplot(211);stem(0:31,x32n,'.');subplot(212);stem(0:31,abs(x32k),'.');xn=1:10;a=1:14;xb=13:-1:1;subplot(231)stem(0:26,xn,'.');xk=fft(xn,1024);subplot(232)plot(abs(xk));x16k=xk(1,1:64:1024);x16n=ifft(x16k,16);subplot(233);stem(0:15,x16n,'.');subplot(234)stem(0:15,abs(x16k),'.');x32k=xk(1,1:32:1024);x32n=ifft(x32k,32);subplot(235);stem(0:31,x32n,'.');subplot(236);stem(0:31,abs(x32k),'.');32点采样时,采样点数为N=32，时域信号的长度M=27，，时域不产生混叠；16点采样时，采样点数为N=16，N<M，发生了时域混叠失真。
实验三-时域及频域采样定理
实验三-时域及频域采样定理⼴州⼤学学⽣实验报告开课学院及实验室:电⼦楼３１7 ?２０13年⽉⽇ｎ）以Ｎ为周期进⾏周期延拓后的主值区序列，--三、实验⽤ＭＡTL ＡＢ函数介绍Xk=fft(x ｎ,N)：采⽤FFT 算法计算序列向量x ｎ的Ｎ点DFT ，缺省N 时，ｆft 函数⾃动按xn 的长度计算xn 的ＤFT 。
当N 为2的整数次幂时，ff ｔ按基2算法计算,否则⽤混合基算法。
2. i ｆf ｔ功能:⼀维快速逆傅⽴叶变换（IFFT)。
调⽤格式：与f ｆt 相同。
四、实验内容和步骤（⼀）时域采样定理实验１. 给定模拟信号如下:0()sin()()at a x t Ae t u t -=Ω假设式中A=4４4.１28,250π=a , 2500π=Ωｒad ／s ，将这些参数代⼊上式中,对()a x t 进⾏傅⽴叶变换,得到()a X j Ω,画出它的幅频特性()~a X jf f,如图３．1所⽰。
根据该曲线可以选择采样频率。
图3.１()a x t 的幅频特性曲线实验过程及原始数据：ｃlf ；A=4４4.128；a=50*ｐi ＊sqrt(２);ｗ０＝5０＊pi ＊sqrt （2）; ｆs=１000; %采样频率10０0HZ T ＝1/fs;n=０：０.05＊fs －1; ％产⽣的长度区间nx ｔ=A*ｅxp(－ａ*n*T).*sin(w0＊ｎ*Ｔ）; %产⽣采样序列xt(n) f ＝f ｆt(xt,l ｅng ｔh(n)); %采样序列xt(n)的FFT 变换 k1=0:ｌe ｎg ｔh （f)-1；ｆk1=k1/0.0５; %设置xt(ｎ）的频谱的横坐标的取值ｓubplot （1，２,1)st ｅｍ(n,xt,＇.＇) %xt 离散图 tit ｌe('(a)fs=1000Hz');xlab ｅl （'n ＇);ylab ｅl('x(n ）'）; subplot （１,２，2) plot(fk1,ab ｓ（ｆ))title(＇(a ） FT[x(ｎＴ）]，Fs ＝100０Hz ＇)； xlabel('ｆ(H ｚ)');ylabel('幅度'）2. 按照选定的采样频率对模拟信号进⾏采样,得到时域离散信号()x n :0()()sin()()anT a x n x nT Ae nT u nT ==Ω这⾥给定采样频率如下：1s f kHz =，300Ｈz ，20０Hz 。
时域采样定理实验报告
一、实验目的1. 理解时域采样定理的基本原理。
2. 掌握信号的采样过程，并分析采样频率对信号的影响。
3. 通过实验验证时域采样定理的正确性。
二、实验原理时域采样定理（Nyquist-Shannon采样定理）指出：一个频带限制在（0，fM）内的信号，如果以不低于2fM的采样频率进行采样，则采样信号能够无失真地恢复原信号。
三、实验设备1. 信号发生器2. 采样器3. 数据采集器4. 计算机5. 信号处理软件（如MATLAB）四、实验步骤1. 设置信号发生器，产生一个频带限制在（0，fM）内的信号，例如正弦波信号，频率为fM。
2. 设置采样器，选择合适的采样频率fS。
根据时域采样定理，fS应满足fS≥2fM。
3. 采集信号，记录采样数据。
4. 利用信号处理软件对采集到的数据进行处理，分析采样频率对信号的影响。
5. 对比不同采样频率下的信号，验证时域采样定理的正确性。
五、实验结果与分析1. 采样频率为fS=2fM时，采样信号能够无失真地恢复原信号。
此时，信号处理软件分析结果显示，信号频谱在fM以下没有出现混叠现象。
2. 采样频率为fS=fM时，采样信号出现失真。
此时，信号处理软件分析结果显示，信号频谱在fM以下出现混叠现象，导致信号失真。
3. 采样频率为fS=1.5fM时，采样信号失真较大。
此时，信号处理软件分析结果显示，信号频谱在fM以下出现较严重的混叠现象，信号失真明显。
六、实验结论通过本次实验，我们验证了时域采样定理的正确性。
实验结果表明，在满足时域采样定理的条件下，采样信号能够无失真地恢复原信号。
同时，实验也表明，采样频率对信号的影响较大，应选择合适的采样频率以保证信号质量。
七、实验总结本次实验使我们深入理解了时域采样定理的基本原理，掌握了信号的采样过程，并分析了采样频率对信号的影响。
通过实验验证了时域采样定理的正确性，提高了我们的信号处理能力。
在今后的学习和工作中，我们将继续关注信号处理技术，不断提高自己的专业知识水平。
数字信号处理实验三时域及频域采样定理
Xk1=fft(x1,length(n1)); %采样序列x1(n)的FFT变换
Xk2=fft(x2,length(n2)); %采样序列x2(n)的FFT变换
Xk3=fft(x3,length(n3)); %采样序列x3(n)的FFT变换
k1=0:length(Xk1)-1;
fk1=k1/Tp; %x1(n)的频谱的横坐标的取值
这里给定采样频率如下： ，300Hz，200Hz。分别用这些采样频率形成时域离散信号，按顺序分别用 、 、 表示。选择观测时间 。
3.计算 的傅立叶变换 ：
（3.6）
式中， ，分别对应三种采样频率的情况 。采样点数用下式计算：
（3.7）
（3.6）式中， 是连续变量。为用计算机进行数值计算，改用下式计算：
下面分析频域采样定理。对信号x(n)的频谱函数 ，在[0，2π]上等间隔采样N点，得到
（3.4）
则N点IDFT[ ]得到的序列就是原序列x(n)以N为周期进行周期延拓后的主值区序列，公式为：
(3.5)
由上式可知，频域采样点数N必须大于等于时域离散信号的长度M(即N≥M)，才能使时域不产生混叠，则N点IDFT[ ]得到的序列 就是原序列x(n),即 =x(n)。如果N>M， 比原序列尾部多N-M个零点；如果N<M，z则 =IDFT[ ]发生了时域混叠失真，而且 的长度N也比x(n)的长度M短，因此。 与x(n)不相同。
第四章3-时域采样与频域采样
抽样间隔(周期) 抽样角频率 抽样频率
T
(s)
wsam=2/T (rad/s)
fsam=1/T (Hz)
例 已知实信号x(t)的最高频率为fm (Hz)，试计 算对各信号x(2t)， x(t)*x(2t)， x(t)x(2t),
x(t)x(2t)抽样不混叠的最小抽样频率。
解：根据信号时域与频域的对应关系得：
F ( jw)
1
f (t) F( jw)
最高频率： fm
wm
w
0 wm
最小抽样频率： 2fm
17
根据信号时域与频域的对应关系及抽样定理得：
f (2t) 1 F( jw )
22
最高频率： 2fm 最小抽样频率： 4fm
根据信号时域卷积与频域的对应关系及抽样定理得：
X (jw)
1 X (j ) T
X (e j )
W sam
ω
0 wm s
W sam
离散序列x[k]频谱与抽样间隔T之间的关系
X ( jw)
wsam 2wm
1
w
wm 0 wm
X (e jwT )
X [ j(w wsam )]
1 X ( jw)
X [ j(w wsam )]
...
f (t)* f (2t) F( jw) 1 F( jw )
22
最高频率： fm
最小抽样频率： 2fm
根据信号时域与频域的对应关系及抽样定理得：
f (t) f (2t) 1 F( jw)* 1 F( j w )
2
22
1
1 F ( jw)* 1 F( jw )
2
实验三 时域采样与频域采样
实验三程序代码及实验结果图：（1）时域采样理论的验证。
给定模拟信号，x a t=Ae−αt sinΩ0t u(t),现用DFT(FFT)求该模拟信号的幅频特性，以验证时域采样理论。
按照x a t的幅频特性曲线，选取三种采样频率，即F s=1kHz，300Hz，200Hz。
观测时间选T p=50ms。
要求：编写实验程序，计算x1n、x2n和x3n的幅度特性，并绘图显示，观察分析频谱混叠失真。
实验程序代码及结果如下：Tp=64/1000; %观察时间Tp=64msFs=1000; %采样率1khzT=1/Fs; %采样间隔M=Tp*Fs; %fft给定点数n=0:M-1; %序列从0开始，至少要达到fft的点数%产生模拟信号对应的离散序列x1(n)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);%无需乘u(t)，因为序列从0开始Xk=T*fft(xnt,M); %频谱函数n1 =0: (length(xnt)-1); %求出序列长度n2 =0: (length(Xk)-1);subplot(3,2,1); %位置为左上stem(n1,xnt); %时域波形title('时域1000hz采样波形'); %标题fk=n2/Tp; %求出频率subplot(3,2,2); %位置为右上stem(fk,abs(Xk)); %幅频特性曲线title('频域1000hz采样'); %标题%-----------300hz-------------Tp=64/300; %观察时间Tp=64msFs=300; %采样率300hzT=1/Fs; %采样间隔M=Tp*Fs; %fft给定点数n=0:M-1; %序列从0开始，至少要达到fft的点数%产生模拟信号对应的离散序列x1(n)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);%无需乘u(t)，因为序列从0开始Xk=T*fft(xnt,M); %频谱函数n1 =0: (length(xnt)-1); %求出序列长度n2 =0: (length(Xk)-1);subplot(3,2,3); %位置为左中stem(n1,xnt); %时域波形title('时域300hz采样'); 标题fk=n2/Tp; %求出频率subplot(3,2,4); %位置为右中stem(fk,abs(Xk)); %幅频特性曲线title('频域300hz采样'); 标题%-----------200hz-------------Tp=64/200; %观察时间Tp=64msFs=200; %采样率200hzT=1/Fs; %采样间隔M=Tp*Fs; %fft给定点数n=0:M-1; %序列从0开始，至少要达到fft的点数%产生模拟信号对应的离散序列x1(n)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);%无需乘u(t)，因为序列从0开始Xk=T*fft(xnt,M); %频谱函数n1 =0: (length(xnt)-1); %求出序列长度n2 =0: (length(Xk)-1);subplot(3,2,5);stem(n1,xnt); %时域波形title('时域200hz采样');fk=n2/Tp;subplot(3,2,6); %位置在右下stem(fk,abs(Xk)); %幅频特性曲线title('频域200hz采样'); %标题（2）频域采样理论的验证。
验证时域采样定理和频域采样定理__数字信号处理.doc
恍恍惚惚课程设计报告Fra bibliotek课程名称
系 别：
专业班级：
学 号：
姓 名：
数字信号课程设计
工程技术系
电子信息工程
09XXXXXX7
课程题目： 验证时域采样定理和频域采样定理
完成日期：
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力根保通据护过生高管产中线工资敷艺料设高试技中卷术资配0料不置试仅技卷可术要以是求解指，决机对吊组电顶在气层进设配行备置继进不电行规保空范护载高高与中中带资资负料料荷试试下卷卷高问总中题体资，配料而置试且时卷可，调保需控障要试各在验类最；管大对路限设习度备题内进到来行位确调。保整在机使管组其路高在敷中正设资常过料工程试况中卷下，安与要全过加，度强并工看且作护尽下关可都于能可管地以路缩正高小常中故工资障作料高；试中对卷资于连料继接试电管卷保口破护处坏进理范行高围整中，核资或对料者定试对值卷某，弯些审扁异核度常与固高校定中对盒资图位料纸置试，.卷保编工护写况层复进防杂行腐设自跨备动接与处地装理线置，弯高尤曲中其半资要径料避标试免高卷错等调误，试高要方中求案资技，料术编试交写5、卷底重电保。要气护管设设装线备备置敷4高、调动设中电试作技资气高，术料课中并中3试、件资且包卷管中料拒含试路调试绝线验敷试卷动槽方设技作、案技术，管以术来架及避等系免多统不项启必方动要式方高，案中为；资解对料决整试高套卷中启突语动然文过停电程机气中。课高因件中此中资，管料电壁试力薄卷高、电中接气资口设料不备试严进卷等行保问调护题试装，工置合作调理并试利且技用进术管行，线过要敷关求设运电技行力术高保。中护线资装缆料置敷试做设卷到原技准则术确：指灵在导活分。。线对对盒于于处调差，试动当过保不程护同中装电高置压中高回资中路料资交试料叉卷试时技卷，术调应问试采题技用，术金作是属为指隔调发板试电进人机行员一隔，变开需压处要器理在组；事在同前发一掌生线握内槽图部内纸故，资障强料时电、，回设需路备要须制进同造行时厂外切家部断出电习具源题高高电中中源资资，料料线试试缆卷卷敷试切设验除完报从毕告而，与采要相用进关高行技中检术资查资料和料试检，卷测并主处且要理了保。解护现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
时域采样和频域采样实验报告
时域采样和频域采样实验报告一、实验目的本次实验旨在掌握时域采样和频域采样的原理、方法和技巧，研究它们在信号处理中的应用。
二、实验原理1. 时域采样时域采样是指将连续时间信号转换为离散时间信号的过程。
其原理是在一定时间间隔内对连续时间信号进行采样，得到离散时间信号。
采样定理规定：如果一个连续时间信号没有高于Nyquist频率两倍以上的频率分量，那么它可以通过等间隔采样来完全恢复。
2. 频域采样频域采样是指将连续频率信号转换为离散频率信号的过程。
其原理是对连续频率信号进行傅里叶变换，得到其频谱，并按照一定间隔取出其中若干个点，得到离散频率信号。
三、实验步骤1. 时域采样实验步骤：（1）使用函数发生器产生正弦波信号；（2）将正弦波信号输入示波器，并设置合适的水平和垂直尺度；（3）调整示波器触发方式为单次触发，同时设置触发电平和触发边沿；（4）按下示波器的单次触发按钮，记录采样到的离散时间信号；（5）将离散时间信号输入计算机，并进行处理和分析。
2. 频域采样实验步骤：（1）使用函数发生器产生正弦波信号；（2）将正弦波信号输入示波器，并设置合适的水平和垂直尺度；（3）通过示波器自带的FFT功能，对正弦波信号进行傅里叶变换，并得到其频谱图；（4）选取频谱图中若干个点，记录其幅值和相位信息；（5）将记录的幅值和相位信息输入计算机，并进行处理和分析。
四、实验结果与分析1. 时域采样实验结果与分析：在本次实验中，我们使用函数发生器产生了一个频率为1kHz、幅度为5V的正弦波信号，并将其输入示波器。
通过调整示波器触发方式为单次触发，同时设置触发电平和触发边沿，我们成功地对正弦波信号进行了时域采样，并得到了一组离散时间信号。
将这些离散时间信号输入计算机，并进行处理和分析，我们得到了正弦波信号的时域图像。
2. 频域采样实验结果与分析：在本次实验中，我们使用函数发生器产生了一个频率为1kHz、幅度为5V的正弦波信号，并将其输入示波器。
验证时域采样定理和频域采样定理—
课程设计报告课程名称数字信号课程设计系别：专业班级：学号：姓名：课程题目：验证时域采样定理和频域采样定理完成日期：2013年5月23日指导老师：2013年5 月23 日验证时域采样定理和频域采样定理摘要数字信号处理是将信号以数字方式表示并处理的理论和技术。
数字信号处理与模拟信号处理是信号处理的子集。
数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。
因此在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数转换器实现。
而数字信号处理的输出经常也要变换到模拟域，这是通过数模转换器实现的。
编制Matlab程序，完成以下功能，对给定模拟信号进行时域采样，观察不同采样频率对采样信号频谱的影响，验证时域采样定理；对给定序列进行傅里叶变换，并在频域进行采样，观察不同采样点数对恢复序列的影响，验证频域采样定理；绘制相关信号的波形。
关键字：时域采样，频域采样，数字信号处理，matlab目录一、摘要 (4)二、绪论 (6)三、方案 (6)1.验证时域采样定理 (6)2.详细程序及仿真波形分析 (7)3.频域采样理论的验证 (15)4.频域采样定理程序 (16)5.频域采样定理信号波形 (17)四、结论 (17)致谢 (18)参考文献 (18)一、绪论数字信号处理是将信号以数字方式表示并处理的理论和技术。
数字信号处理与模拟信号处理是信号处理的子集。
数字信号处理的目的是对真实世界的连续模拟信号进行测量或滤波。
因此在进行数字信号处理之前需要将信号从模拟域转换到数字域，这通常通过模数转换器实现。
而数字信号处理的输出经常也要变换到模拟域，这是通过数模转换器实现的。
编制Matlab程序，完成以下功能，对给定模拟信号进行时域采样，观察不同采样频率对采样信号频谱的影响，验证时域采样定理；对给定序列进行傅里叶变换，并在频域进行采样，观察不同采样点数对恢复序列的影响，验证频域采样定理；绘制相关信号的波形。
二、方案1.验证时域采样定理基本要求：①掌握数字信号处理的基本概念、基本理论和基本方法；②学会 MATLAB 的使用，掌握 MATLAB 的程序设计方法；③学会用 MATLAB 对信号进行分析和处理；④信号的各参数需由键盘输入，输入不同参数即可得不同的x(t) 和x(n)；⑤撰写课程设计论文，用数字信号处理基本理论分析结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
广州大学学生实验报告开课学院及实验室:电子楼３１7 ﻩ２０13年月日ｎ）以Ｎ为周期进行周期延拓后的主值区序列，--三、实验用ＭＡTL ＡＢ函数介绍１． fft功能:一维快速傅立叶变换(FFT)。
Xk=fft(x ｎ,N)：采用FFT 算法计算序列向量x ｎ的Ｎ点DFT ，缺省N 时，ｆft 函数自动按xn 的长度计算xn 的ＤFT 。
当N 为2的整数次幂时，ff ｔ按基2算法计算,否则用混合基算法。
2. i ｆf ｔ功能:一维快速逆傅立叶变换（IFFT)。
调用格式：与f ｆt 相同。
四、实验内容和步骤 （一） 时域采样定理实验１. 给定模拟信号如下:0()sin()()at a x t Ae t u t -=Ω假设式中A=4４4.１28,250π=a , 2500π=Ωｒad ／s ，将这些参数代入上式中,对()a x t 进行傅立叶变换,得到()a X j Ω,画出它的幅频特性()~a X jf f,如图３．1所示。
根据该曲线可以选择采样频率。
图3.１()a x t 的幅频特性曲线实验过程及原始数据： ｃlf ；A=4４4.128；a=50*ｐi ＊sqrt(２);ｗ０＝5０＊pi ＊sqrt （2）; ｆs=１000; %采样频率10０0HZ T ＝1/fs;n=０：０.05＊fs －1; ％产生的长度区间nx ｔ=A*ｅxp(－ａ*n*T).*sin(w0＊ｎ*Ｔ）; %产生采样序列xt(n) f ＝f ｆt(xt,l ｅng ｔh(n)); %采样序列xt(n)的FFT 变换 k1=0:ｌe ｎg ｔh （f)-1；ｆk1=k1/0.0５; %设置xt(ｎ）的频谱的横坐标的取值 ｓubplot （1，２,1)st ｅｍ(n,xt,＇.＇) %xt 离散图 tit ｌe('(a)fs=1000Hz');xlab ｅl （'n ＇);ylab ｅl('x(n ）'）; subplot （１,２，2) plot(fk1,ab ｓ（ｆ))title(＇(a ） FT[x(ｎＴ）]，Fs ＝100０Hz ＇)； xlabel('ｆ(H ｚ)');ylabel('幅度'）2. 按照选定的采样频率对模拟信号进行采样,得到时域离散信号()x n :0()()sin()()anT a x n x nT Ae nT u nT ==Ω这里给定采样频率如下：1s f kHz =，300Ｈz ，20０Hz 。
分别用这些采样频率形成时域离散信号,按顺序分别用1()x n 、2()x n 、3()x n 表示。
选择观测时间50p T ms=。
实验过程及原始数据： c ｌf;A=444.12８;a=50＊ｓqrt(2)*pi;ｗ０=５０*sqrt(2)*p ｉ; Ｔｐ=50/１０00;F1=1000;F2=3０0;Ｆ3=200; Ｔ１=１/F1;T2=１／F ２;T3=1/F ３；n1＝0:0．5： T ｐ*F1-1;n2＝0:0.5:T ｐ*F2-1;n3=0:0．５:T ｐ＊Ｆ3－1; x1=A*exp （-ａ＊n1*T1).＊s ｉn （w0*n １*T １)；ｘ2=A*ｅxp （－ａ*n2＊T ２).*ｓin(w0*n ２*T2）； x3=A*exp （-a*n3*T3).*ｓｉn(w ０*ｎ３*T3）； f1=fft （x １，length （n1)); f2＝fft(x2，len ｇth(n2））； ｆ3＝fft(x3，l ｅng ｔｈ(n3)); k １=0:ｌｅｎgth （f1)-1;f ｋ1=ｋ1/Ｔp; ｋ２=0：length(f2）－1；fk2=ｋ2/Tp; k3＝0：le ｎg ｔh(ｆ3)-１;--ｆk3=ｋ3/Tp; sub ｐlot(3,2，1) stem(ｎ1,x1,＇.') title('(ａ)Fs ＝1000Hz');ｘlabe ｌ(＇ｎ＇）;ylabe ｌ('x1(n)'); sub ｐlot(３,2,3) ｓtem(n2,x2，'．') tit ｌｅ('(b)Ｆｓ=3０0Hz');ｘlab ｅｌ('ｎ'）;ｙl ａbel('x2(n)＇); ｓｕｂpl ｏt(3,2,5） ｓtem(n3,x ３,＇．') ti ｔｌe('(c)Ｆs=200Ｈz');ｘｌａbel （＇n');yla ｂel （'x ３(n ）'); su ｂｐlot （３,2,2) pl ｏｔ(fk １,ａbs(ｆ１))tit ｌe('(ａ) FT[x ａ(ｎＴ)]的频谱特性,Fs=1000Hz'）; xlabe ｌ（'f(Ｈz)'）;ylabel(＇幅度') ｓubpl ｏt （3,2,4) pl ｏt(ｆｋ2，a ｂs(f2))ｔitl ｅ（'(ｂ） F Ｔ[xa(ｎT)] 的频谱特性,Ｆｓ=300Hz'); xlab ｅl('f(Ｈz)'）;yl ａbe ｌ('幅度') subp ｌot(3,2,6） plot(fk3,ab ｓ(ｆ３))ti ｔle('（c ） F Ｔ[xa(ｎT ）］ 的频谱特性,F ｓ=2０0Hz ＇）; xlabel('f(H ｚ）＇);ｙlab ｅl('幅度＇)3． 计算()x n 的傅立叶变换()jwX e ： 100()[()]sin()i in anT jwj ni n X e FT x n AenT eω--===Ω∑ （3.6)式中，1,2,3i =，分别对应三种采样频率的情况123111(,,)1000300200T s T s T s ===。
采样点数用下式计算：pi i T n T =（3.７）（3.6）式中,ω是连续变量。
为用计算机进行数值计算，改用下式计算：100()[()]sin()i ki k n jw anT jw n M i n X eDFT x n Ae nT e --===Ω∑ （3．8）式中，2k kM πω=，0,1,2,3...k =,1M -；64M =。
可以调用ＭATLA Ｂ函数fft 计算3.8式。
4． 打印三种采样频率的幅度曲线()~k jw kX e w ，0,1,2,3...k =，1M -;64M =。
f=inp ｕｔ(＇ｆ= '); %输入取样频率f=１000 T ＝1/f;n=0:50/（１00０*T);Ｍ=64;%设置信号有关参数 Ａ=444.１28; a=50*sq ｒt(2.０);ｗ0=５0*sqrt （２.0)*pi; xn=A*ex ｐ（-ａ*n*T).*sin(w0*n*T); % 产生ｘ（n)Ｘk=f ｆt(xn,1024)； %1０24 点FF Ｔ[x （n)],用于近似序列 ｘ(n)的FT Ｘ6４ｋ=fft(ｘｎ,M)； %ＦF Ｔ［ｘ（n)]x ６４ｎ＝ｉfft(X6４k ）； ％６４ 点IFFT[X64（ｋ)] 得到64(ｎ） ｓubpl ｏｔ(2,２,1）；ste ｍ(n,x ｎ,'.');ｇrid ；box on title （＇(a) x1（n)＇)；xlab ｅl('n');y ｌabel （'x １(n)'); k=０:M －1;s ｕｂplo ｔ（2,2,2)；stem （k,abs(X ６4k ）,'.＇）；ｇｒid; t ｉｔle('(ｂ) 64 点频域采样 输入取样频率ｆ=10０0'); xlab ｅl('k');ylabel('|X ６4(ｋ)｜＇)； k ＝0:1023;ｗk ＝2*k/1024;sub ｐｌｏt （2,2，3);plot(wk,abs(Xk ）);gr ｉｄ;ｔi ｔle('(c)FT[x １(ｎ)]＇)；xlabel('\o ｍega/\pi ＇);ylabel('|X （e^j ＾\ome ｇa)｜'); ｎ1＝0:M-1;su ｂplo ｔ(２,2,4);ste ｍ（n １，x6４n,'.');grid ；box o ｎ titl ｅ('(d ） ６4 点'）ifft(X64（k))；ｘｌa ｂe ｌ('n'）;y ｌａb ｅｌ('ｘ64（ｎ)')； ｆ＝input('f= '）; % 输入取样频率ｆ＝300 Ｔ=1/f;n=0:50/(1000*T);Ｍ＝64;%设置信号有关参数Ａ=4４4.１28； a=50*ｓqrt(2.0）;ｗ0=5０*ｓqrt （２．０)*pi; x ｎ=Ａ＊exp(－a*ｎ*T).*s ｉn(w0*ｎ*Ｔ); % 产生x(ｎ) X ｋ＝f ｆt(xn,１０24); %1024 点F ＦT[x(n)], 用于近似序列x(n)的F ＴＸ６4ｋ=fft （ｘn,Ｍ); %F ＦＴ［x(n)]ｘ６4n ＝ｉf ｆt(X64ｋ); %64点 IFFT[X64(k)］得到 x6４(n) subplot （2,2,１);ｓtem （n ，xn,＇.');g ｒｉd ；ti ｔｌe('(ａ) x2（ｎ）＇);xlabe ｌ（'n ＇)；ylab ｅl('ｘ２(ｎ)'); k=0:M －1;subplot(2,2，2）;st ｅｍ(k,abs （X ６4ｋ),＇．'）;gr ｉｄ;box o ｎtit ｌｅ('(b) 64 点频域采样输入取样频率f=30０'); ｘlabel('ｋ＇);yla ｂｅｌ(＇|Ｘ6４（k)|')； k=0:10２3;wk=２＊k/102４;ｓubplot(２,２,3)；plot(wk ，abs(X ｋ));ｇrid;title （＇（c ）F Ｔ[x2(n)］＇);xlabel （'\ome ｇa/\pi');ylabel('|X （e^j^＼ome ｇa)|');n1=0:M-1;s ｕbplot(2，2，4）;s ｔe ｍ(n1,x64n,'．＇);g ｒｉd;ｂox on titl ｅ('(d) 64 点IFF Ｔ[X64(k)]')；-- xlaｂeｌ(＇n')；ylａbｅl('x6４(n）'）;f＝inｐｕt('f＝ ');% 输入取样频率f=200T=1/f;n=0:5０/(10０0*T)；M=64; %设置信号有关参数A=44４．128; a=50*sqrt(２.０);w0=5０*ｓqｒt(2.0)*pi；xn=A*eｘp(-a*ｎ＊T)．*sin(w0*n*T）; %产生 x（n)Xk=fft（xn,102４)；%10２4 点FFT[ｘ(ｎ)]，用于近似序列x(ｎ)的FTX６４k＝ffｔ（xｎ,Ｍ);％FＦT［ｘ(ｎ）]x６4n＝iｆft(X64k）; %64 点ＩＦFT[X64(k)］得到x64(n）subplｏt(2,2，１）;ｓtｅm(n,ｘn,'.'）;gｒid;boｘontitｌe('(a）ｘ3（n)＇);xlabｅl(＇n＇)；yｌabeｌ（＇x3(ｎ)＇);k=０:Ｍ-1；suｂploｔ(2，２,2);stｅｍ（ｋ,ａbs（X64k),'．')；grid;tｉtle('(b） 6４点频域采样输入取样频率f=20０＇);xlabel('k')；ylａｂel（'｜X64(k）|');k=0：１02３;wk＝2*k／10２４;subplｏｔ(2,２,３）;plｏｔ(wｋ,abs(Xk)）;grid;ｔitｌｅ(＇（c)FT[x3(n)]')；ｘｌabｅl（'＼omegａ/\ｐi')；yｌabel('｜Ｘ(e^ｊ＾＼ｏmega)｜');ｎ1=0:M－１;subpｌot(２,2,4);ｓtem(ｎ1，x6４ｎ,'．')；grｉd;bｏx oｎtitlｅ('(d） 64 µã')ifft（X64（ｋ）);xlabｅl(＇ｎ'）;yｌａbeｌ('x６4(n)');采样频率ｆs=１０00采样频率ｆs=３0０采样频率fs=200 (二)频域采样定理实验给定信号如下:⎪⎩⎪⎨⎧≤≤-≤≤+=其它261427131)(nnnnnx编写程序分别对频谱函数()FT[()]jX e x nω=在区间]2,0[π上等间隔采样３２和16点,得到)()(1632kXkX和:32232()() , 0,1,2,31jkX k X e kωπω===16216()() , 0,1,2,15jkX k X e kωπω===再分别对)()(1632kXkX和进行３2点和16点IFFT,得到)()(1632nxnx和：323232()IFFT[()] ,0,1,2,,31x n X k n==161616()IFFT[()] ,0,1,2,,15x n X k n==分别画出()jX eω、)()(1632kXkX和的幅度谱,并绘图显示x(n）、)()(1632nxnx和的波形，进行对比和分析,验证总结频域采样理论。