实验一：低通采样定理和内插与抽取实现a
实验一：低通采样定理和内插与抽取实现一．实验目的1. 连续信号和系统的表示方法，以及坊真方法。
2.用MATLAB实现连续信号采用与重构的方法，3. 采样信号的插值和抽取等重采样实现方法。
4. 用时域采样信号重构连续时域信号的原理和方法。
5. 用MATLAB绘图函数表示信号的基本方法，实验数据的可视化表示。
二．原理1 、时域抽样定理令连续信号xa(t)的傅里叶变换为Xa(jΩ),抽样脉冲序列p(t)傅里叶变换为P(jΩ),抽样后的信号x^(t)的傅里叶变换为X^(jΩ)若采用均匀抽样,抽样周期Ts,抽样频率为Ωs=2πfs,由前面分析可知:抽样的过程可以通过抽样脉冲序列p(t)与连续信号xa(t)相乘来完成,即满足:x^(t)=xa(t) p(t)，又周期信号f(t)傅里叶变换为：故可以推得p(t)的傅里叶变换为:其中:根据卷积定理可知:得到抽样信号x(t)的傅里叶变换为:其表明:信号在时域被抽样后,他的频谱X(jΩ)是连续信号频谱X(jΩ)的形状以抽样频率Ω为间隔周期重复而得到,在重复过程中幅度被p(t)的傅里叶级数Pn 加权。
因为Pn只是n的函数,所以X(jΩ)在重复的过程中不会使其形状发生变化。
假定信号x(t)的频谱限制在-Ωm~+Ωm的范围内, 若以间隔Ts对xa(t)进行抽样,可知抽样信号X^(t)的频谱X^(jΩ)是以Ωs为周期重复。
显然,若在抽样的过程中Ωs<2Ωm,则X^(jΩ)将发生频谱混叠现象,只有在抽样的过程中满足Ωs>=2Ωm条件,X^(jΩ)才不会产生频谱的混叠,接收端完全可以由x^(t)恢复原连续信号xa(t),这就是低通信号抽样定理的核心内容。
2、信号的重建从频域看,设信号最高频率不超过折叠频率:Xa(jΩ)=Xa(jΩ) |Ω|<Ωs/2Xa(jΩ)=0 |Ω|>Ωs/2则理想取样后的频谱就不会产生混叠,故有:让取样信号x^(t)通过一带宽等于折叠频率的理想低通滤波器:H(jΩ)=T |Ω|<Ωs/2H(jΩ)=0 |Ω|>Ωs/2滤波器只允许通过基带频谱,即原信号频谱,故:Y(jΩ)=X^(jΩ)H(jΩ)=Xa(jΩ)因此在滤波器的输出得到了恢复的原模拟信号:y(t)=xa(t)从时域上看,上述理想的低通滤波器的脉冲响应为:根据卷积公式可求得理想低通滤波器的输出为:由上式显然可得:则:上式表明只要满足取样频率高于两倍信号最高频率,连续时间函数xa(t)就可用他的取样值xa(nT)来表达而不损失任何信息，这时只要把每一个取样瞬时值与内插函数式相乘求和即可得出xa(t),在每一取样点上,由于只有该取样值所对应的内插函数式不为零,所以各个取样点上的信号值不变。
三． 内容1.设计连续时间信号线性滤波器分离信号组份已知信号()()990(1)cos 2(10050)m x t m m t π==++∑，试设计滤波器，分离出如下信号：（1）m 1,2,3......50= （2）m 51,52,53......100= （3）m 40,41,42......60= （4）m 1,2,3......40,61,62......100=据以下采样频率：(a) 20000s f Hz = (b) 10000s f Hz = (c)30000s f Hz = 求信号频谱及相应的滤波器。
参考程序如下：设计一个Butterworth 模拟带通滤波器，设计指标为：通带频率：1000- 2000Hz ，两侧过渡带宽500Hz ，通带波纹1dB ，阻带衰减100dB 。
假设一个信号，其中f1=100Hz,f2=1500Hz,f3=2900Hz 。
信号的采样频率为10000Hz 。
试将原信号与通过该滤波器的模拟信号进行比较。
参考程序如下：wp=[1000 2000]*2*pi;ws=[500 2500]*2*pi;Rp=1;Rs=100; %滤波器设计参数，对于给定Hz 应乘以2[N,Wn]=buttord(wp,ws,Rp,Rs,'s'); %求得滤波器的最小阶数和截止频率 w=linspace(1,3000,1000)*2*pi; %设置绘制频率响应的频率点 [b,a]=butter(N,Wn,'s'); %设计模拟Butterworth 滤波器H=freqs(b,a,w); %计算给定频率点的复数频率响应magH=abs(H);phaH=unwrap(angle(H)); %计算幅频响应和相频响应 plot(w/(2*pi),20*log10(magH)); %以频率为横坐标绘制幅频响应 xlabel('频率/Hz');ylabel('振幅/dB'); title('Butterworth 模拟带通滤波器');hold on;plot([1000 1000],ylim,'r');plot([2000 2000],ylim,'r');%绘带边界grid on figure(2)dt=1/10000; %模拟信号采样间隔f1=100;f2=1500;f3=2900;%输入信号的三个频率成分 t=0:dt:0.04; %给定模拟时间段x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t); %输入信号 H=[tf(b,a)]; %滤波器在MATLAB 系统中的表示 [y,t1]=lsim(H,x,t); %模拟输出subplot(2,1,1),plot(t,x),title('输入信号') %绘出输入信号 subplot(2,1,2),plot(t1,y) %绘制输出信号 title('输出信号'),xlabel('时间/s')2. 连续时间信号的采样和重建已知信号()()990(1)cos 2(10050)m x t m m t π==++∑，试以20000s f Hz =采样频率对该信号采样，并用插值公式重建该信号。
参考程序：①、分别用150HZ 及300HZ 对信号采样源信号为：fa=5*sin(2*pi*40*t1)+1.8*sin(4*pi*40*t1)+0.8*sin(5*pi*40*t1)，用150Hz 的频率对f(t)进行采样，其采样图如图1所示；用300Hz 的频率对f(t)进行采样，其采样图如图2所示。
程序如下： fs1=150;t1=-0.1:1/fs1:0.1;fa=5*sin(2*pi*40*t1)+1.8*sin(4*pi*40*t1)+0.8*sin(5*pi*40*t1); figure(1);plot(t1,fa),xlabel('fs1=150Hz 时,fa 采样时域图'); hold off; fs2=300;t2=-0.1:1/fs2:0.1;fb=5*sin(2*pi*40*t2)+1.8*sin(4*pi*40*t2)+0.8*sin(5*pi*40*t2); figure(2);plot(t2,fb),xlabel('fs2=300Hz 时,fb 采样时域图');fs1=150Hz时,fa采样时域图图1 150HZ采样频率对信号采样图fs2=300Hz时,fb采样时域图图2 300HZ采样频率对信号采样图②、对信号进行快速离散傅里叶变换将两个采样信号进行快速离散傅里叶变换(FFT),用150Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图3所示；用300Hz的频率对f(t)进行采样，其采样后快速傅立叶变换频谱图图4所示。
程序如下：f=40;fs=150;N=300;k=0:N-1;t=-0.1:1/fs:0.1;w1=150*k/N;fa=5*sin(2*pi*f*t)+1.8*sin(4*pi*f*t)+0.8*sin(5*pi*f*t);xfa=fft(fa,N);xf1=abs(xfa);figure(1);plot(w1,xf1),xlabel('fs=150Hz时,fa经fft后频谱图.单位：Hz');f=40;fs=300;N=300;k=0:N-1;t=-0.1:1/fs:0.1;w2=300*k/Nfb=5*sin(2*pi*f*t)+1.8*sin(4*pi*f*t)+0.8*sin(5*pi*f*t);xfb=fft(fb,N);xf2=abs(xfb);figure(2);plot(w2,xf2),xlabel('fs=300Hz时,fb经fft后频谱图.单位：Hz ');fs=150Hz时,fa经fft后频谱图.单位：Hz图3 150HZ采样后经FFT后频谱图fs=300Hz时,fb经fft后频谱图.单位：Hz图4 300HZ采样后经FFT后频谱图③、信号的重建我们可以通过利用内插法把原信号从采样信号中恢复出来，观察信号在满足怎样的采样条件下能够恢复为原信号，图5和图6分别为恢复后的原信号。
程序如下：Wm=180*pi;Wc=Wm;fs1=150;Ws=2*pi*fs1;n=-800:800;nTs1=n/fs1;fa=5.1*sin(2*pi*40*nTs1)+1.8*sin(4*pi*40*nTs1)+0.8*sin(5*pi*40*nTs1);Dt=1/fs1;t1=-0.1:Dt:0.1;fa1=fa/fs1*Wc/pi*sinc((Wc/pi)*(ones(length(nTs1),1)*t1-nTs1'*ones(1,length( t1)));figure(1);plot(t1,fa1);axis([-0.1 0.1 -8 8]);xlabel('fs=150Hz时,fa利用内插由样本重建原信号图.');Wm=180*pi;Wc=Wm;fs2=300;Ws=2*pi*fs2;n=-800:800;nTs2=n/fs2;fb=5.1*sin(2*pi*40*nTs2)+1.8*sin(4*pi*40*nTs2)+0.8*sin(5*pi*40*nTs2);Dt=1/fs2;t1=-0.1:Dt:0.1;fb1=fb/fs2*Wc/pi*sinc((Wc/pi)*(ones(length(nTs2),1)*t1-nTs2'*ones(1,length( t1)));figure(2);plot(t1,fb1);axis([-0.1 0.1 -8 8]);xlabel('fs=300Hz时,fb利用内插由样本重建原信号图.');grid;fs=150Hz时,fa利用内插由样本重建原信号图.图5 150HZ采样后的信号的重建信号fs=300Hz 时,fb 利用内插由样本重建原信号图.图6 300HZ 采样后的信号的重建信号3.采样信号的抽取和插值⑴.已知信号()()990(1)cos 2(10050)m x t m m t π==++∑，以20000s f Hz =采样频率采样后，设计抽取因子D 和滤波器，分离出如下信号m 1,2,3......50= 信号组分。
南理工通信原理实验报告
目录实验一抽样定理实验 (3)实验七HDB3码型变换实验 (14)实验十一BPSK调制与解调实验 (21)实验十九滤波法及数字锁相环法位同步提取实验 (29)实验一抽样定理实验一、实验目的1.了解抽样定理在通信系统中的重要性。
2.掌握自然抽样与平顶抽样的实现方法。
3.理解低通采样定理的原理。
4.理解实际的采样系统。
5.理解低通滤波器的幅频特性和对抽样信号恢复的影响。
6.理解带通采样定理的原理。
二、实验器材1.主控&信号源、3号模块。
各一块2.双踪示波器一台3.连接线若干三、实验原理1.实验原理框图2.实验框图说明抽样信号由抽样电路产生。
将输入的被抽样信号与抽样脉冲相乘就可以得到自然抽样信号，自然抽样信号经过保持电路得到平顶抽样信号。
平定抽样和自然抽样信号是通过S1切换输出的。
抽样信号的恢复是将抽样信号经过低通滤波器，即可得到恢复的信号。
这里滤波器可以选用抗混叠滤波器（8阶的巴特沃斯低通滤波器）或fpga数字滤波器（有FIR、IIR两种）。
反sinc滤波器不是用来恢复抽样信号的，而是用来应对孔径失真现象。
要注意，这里的数字滤波器是借用的信源编译码部分的端口。
在做本实验室与信源编译码的内容没有联系。
四、实验结果与波形观测实验项目一抽样信号观测及抽样定理验证概述：通过不同频率的抽样时钟，从时域与频域两方面观测自然抽样和平顶抽样的输出波形，以及信号恢复的混叠情况，从而了解不同抽样方式的输出差异和联系，验证抽样定理。
注：通过观测频谱可以看到当抽样脉冲小于2倍被抽样信号频率时，信号会产生混叠。
源端口目标端口连线说明信号源：MUSIC模块3：TH1(被抽样信号) 将被抽样信号送入抽样单元信号源：A-OUT 模块3：TH2(抽样脉冲)提供抽样时钟模块3：TH3(抽样输出)模块3：TH5(LPF-IN)送入模拟低通滤波器2. 开电，设置主控菜单，选择【主菜单】→【通信原理】→【抽样定理】。
调节主控模块的W1使A-out输出峰峰值为3V。
图解低通抽样定理
即采样频率至少 是基带信号最高 频率的2倍，这就 是低通抽样定理
f
低通信号抽样定理图解：
m(t)
M( f )
t
fH
fH
f
T (t )
T
t
卷积
T( f )
fs
f
t
f
第9章模拟信号的数字传输
我们假设信号m(t)的最高频率小于fH，所以若频率 间隔fs 2fH，则Ms(f)中包含的每个原信号频谱M(f) 之间互不重叠。这样就能够从Ms(f)中用一个低通滤 波器分离出信号m(t)的频谱M(f)，也就是能从抽样 信号中恢复原信号。
这里，恢复原信号的条件是：
fs 2 fH
即抽样频率fs应不小于fH的两倍。这一最低抽样速 率2fH称为奈奎斯特速率。与此相应的最小抽样时 间间隔称为奈奎斯特间隔。
M( f )
如果f s 2 f H T ???
fH
T( f )
fs
fH 混叠，无法提 取出纯净的 M(f)信号了
如果即采样频率至少是基带信号最高频率的2倍这就是低通抽样定理卷积信号频谱发生混叠无法提取出纯净的mf信号了
第9章模拟信号的数字传输
9.2.1 低通模拟信号的抽样定理
抽样定理：设一个连续模拟信号m(t)中的最 高频率 ≤ fH，则以间隔时间为T 1/2fH的周 期性冲激脉冲对它抽样时，m(t)将被这些抽 样值所完全确定。
实验一低通采样定理和内插与抽取实现
实验一低通采样定理和内插与抽取实现一、实验目的用Matlab 编程实现自然采样与平顶采样过程，根据实验结果给出二者的结论；掌握利用MATLAB 实现连续信号采样、频谱分析和采样信号恢复的方法。
二、实验原理 1．抽样定理若)(t f 是带限信号，带宽为m ω,)(t f 经采样后的频谱)(ωs F 就是将)(t f 的频谱)(ωF 在频率轴上以采样频率s ω为间隔进行周期延拓。
因此，当s ω≥m ω时，不会发生频率混叠；而当s ω<m ω时将发生频率混叠。
2．信号重建经采样后得到信号)(t f s 经理想低通)(t h 则可得到重建信号)(t f ，即：)(t f =)(t f s *)(t h 其中：)(t f s =)(t f ∑∞∞--)(s nT t δ=∑∞∞--)()(s s nT t nT f δ，)()(t Sa T t h c csωπω= 所以：)(t f =)(t f s *)(t h =∑∞∞--)()(s s nT t nT f δ*)(t Sa T c csωπω =πωcs T ∑∞∞--)]([)(s csnT t Sa nT f ω上式表明，连续信号可以展开成抽样函数的无穷级数。
利用MATLAB 中的t t t c ππ)sin()(sin =来表示)(t Sa ，有)(s i n )(πt c t Sa =，所以可以得到在MATLAB 中信号由)(s nT f 重建)(t f 的表达式如下：)(t f =πωcs T ∑∞∞--)]([sin )(s cs nT t c nT f πω 我们选取信号)(t f =)(t Sa 作为被采样信号，当采样频率s ω=2m ω时，称为临界采样。
我们取理想低通的截止频率c ω=m ω。
下面程序实现对信号)(t f =)(t Sa 的采样及由该采样信号恢复重建)(t Sa ： 三、实验内容已知信号()()99(1)cos 2(10050)m x t m m t π==++∑，试以以下采样频率对信号采样：(a)20000s f Hz =;(b)10000s f Hz =;(c)30000s f Hz =，求x(t)信号原信号和采样信号频谱，及用采样信号重建原信号x’(t)时序图。
采样定理
sin[ xa (nT )
(t nT )
S 2C
T (t nT ) T
]
采样内插公式
采样内插公式说明：只要满足采样频率高于两倍
信号最高截止频率，则整个连续时间信号就可以 用它的采样值来完全代表，而不会丢失任何信息。
sin[
(t nT )
T (t nT ) T
即
1 ˆ X a ( j) X a ( j) P ( j) 2
ˆa (t ) xa (t ) P (t ) , x
P ( j)
由于
[P (t )] （傅里叶变换）
是周期函数,即
P (t )
n
(t nT )
jk S t a e k
n
理想采样
T
t
(t nT )
0
xa(t)
ˆa (t ) xa (t ) P (t ) x
t
ˆa (t ) x
n
x (nT ) (t nT )
a
二、频谱的周期延拓
xa (t ) X a ( j)
即
-1
ˆ ( j) ˆa (t ) X x a
t
2 s T
称为内插函数
ˆa (t ) g (t ) ya (t ) x
n
[ xa (nT ) ( nT )]g (t )dt
ˆ ( ) g (t )dt x
a
n
由上面的分析有，频谱发生混叠的原因有两个： 1.采样频率低 2.连续信号的频谱没有被限带
通信原理抽样定理实验报告
通信原理实验（五）实验一抽样定理实验项目一、抽样信号观测及抽样定理实验1、观测并记录抽样前后的信号波形，分别观测music和抽样输出。
由分析知，自然抽样后的结果如图，很明显抽样间隔相同，且抽样后的波形在其包络严格被原音乐信号所限制加权，与被抽样信号完全一致。
2、观测并记录平顶抽样前后信号的波形。
此结果为平顶抽样结果，仔细观察可发现与上一实验中的自然抽样有很大差距，即相同之处，其包络也由原信号所限制加权，但是在抽样信号的每个频率分量呈矩形，顶端是平的。
3、观测并对比抽样恢复后信号与被抽样信号的波形，并以100HZ为步进，减小A-OUT的频率，比较观测并思考在抽样脉冲频率为多少的情况下恢复信号有失真。
(1)9.0KHZ(2)7.7KHZ(3)7.0KHZ实验二 PCM 编译码实验实验项目一 测试W681512的幅频特性1、将信号源频率从50HZ 到4000HZ ，用示波器接模块21的音频输出，观测信号的幅频特性。
在频率为9HZ 时的波形如上图，低通滤波器恢复出的信号与原信号基本一致，只是相位有了延时，约1/4个Ts ； 逐渐减小抽样频率可知在7.7KHZ 左右，恢复信号出现了幅度的失真，且随着fs 的减小，失真越大。
上述现象验证了抽样定理，即，在信号的频率一定时，采样频率不能低于被采样信号的2倍，否则将会出现频谱的混(1)、4000HZ (2)、3500HZ(3)120HZ (4)50HZ在实验中仔细观察结果，可知，当信号源的频率由4000HZ不断下降到3000HZ 的过程中，信号的频谱幅度在不断地增加；在3000HZ~1500HZ的过程中，信号的幅度在一定范围内变化，但是没有特别大的差距；在1500HZ~50HZ的过程中，信号的幅度有极为明显的下降。
实验项目二 PCM编码规则实验1、以FS为触发，观测编码输入波形。
示波器的DIV档调节为100微秒。
图中分别为输入被抽样信号和抽样脉冲，观察可发现正弦波与编码对应。
低通抽样定理验证实验
实验二低通抽样定理验证实验一、实验目的1、熟悉使用System View软件，了解各部分功能模块的操作和使用方法。
2、通过实验进一步掌握低通抽样定理的原理。
二、实验内容}用System View建立一个低通抽样定理仿真电路，通过观察各个模块输出波形变化，理解低通抽样定理原理。
三、电路构成图1 低通抽样定理验证实验原理图参数设置：Token3：产生模拟信号（参数设置：Source――Periodic――Sinusoid，幅度1V，频率50HZ，相位0度）《Token4：MultiplierToken5：产生抽样信号（参数设置：Source——Periodic——Pulse Train，幅度1V，频率100Hz，脉冲宽度，偏移0V，相位0度，抽样速率可调）Token6：产生一个模拟低通滤波器，滤除高频信号，保留低频信号（参数设置：Operator——Filters/Systems——Linear Sys Filters，选择：Analog——Lowpass——Butterworth，Lowcuttoff＝50Hz，No of Poles＝3，截止频率＝模拟信号最高频率）四、实验结果(1)原始的输入信号波形图）图2 原始的输入信号波形图(2)原始的输入信号的频谱图图3 原始的输入信号频谱图。
(3)被抽样以后的图形图4 被抽样以后的图形>(4)被抽样以后的频谱图图5 被抽样以后的频谱图分析：由于原始输入波形的离散化，使得输出频谱周期化。
输出频谱如图5所示。
\(5)经过低通滤波器后，还原出波形如图6】图6 还原出的波形(6)经过低通滤波器后，还原后的频谱图!图7 还原后的频谱图可以发现频谱图基本和图3所示相同，但是由于滤波器不是理想低通，使得使得输出频谱周期化的现象仍然存在。
但是基本上已被滤波器滤除，不影响输出波形。
五、思考题#1、观察仿真电路中各个模块输出波形变化，理解低通抽样定理原理。
答：输出波形如上图2至7所示。
抽样定理——精选推荐
抽样定理实验⼀抽样定理实验⼀、实验⽬的1、了解抽样定理在通信系统中的重要性2、掌握⾃然抽样及平顶抽样的实现⽅法3、理解低通采样定理的原理4、理解实际的抽样系统5、理解低通滤波器的幅频特性对抽样信号恢复的影响6、理解低通滤波器的相频特性对抽样信号恢复的影响7、理解平顶抽样产⽣孔径失真的原理8、理解带通采样定理的原理⼆、实验内容1、验证低通采样定理原理2、验证低通滤波器幅频特性对抽样信号恢复的影响3、验证低通滤波器相频特性对抽样信号恢复的影响4、验证带通抽样定理原理5、验证孔径失真的原理三、实验原理抽样定理原理：⼀个频带限制在（0，H f ）内的时间连续信号()m t ，如果以T ≤H f 21秒的间隔对它进⾏等间隔抽样，则()m t 将被所得到的抽样值完全确定。
（具体可参考《信号与系统》）我们这样开展抽样定理实验：信号源产⽣的被抽样信号和抽样脉冲经抽样/保持电路输出抽样信号，抽样信号经过滤波器之后恢复出被抽样信号。
抽样定理实验的原理框图如下：抽样/保持被抽样信号抽样脉冲低通滤波器抽样恢复信号图1抽样定理实验原理框图抽样/保持被抽样信号抽样脉冲低通滤波器抽样恢复信号低通滤波器图2实际抽样系统为了让学⽣能全⾯观察并理解抽样定理的实质，我们应该对被抽样信号进⾏精⼼的安排和考虑。
在传统的抽样定理的实验中，我们⽤正弦波来作为被抽样信号是有局限性的，特别是相频特性对抽样信号恢复的影响的实验现象不能很好的展现出来，因此，这种⽅案放弃了。
另⼀种⽅案是采⽤较复杂的信号，但这种信号不便于观察，如错误!未找到引⽤源。
所⽰：被抽样信号抽样恢复后的信号图3复杂信号抽样恢复前后对⽐你能分辨错误!未找到引⽤源。
中抽样恢复后信号的失真吗因此，我们选择了⼀种不是很复杂，但⼜包含多种频谱分量的信号：“3KHz正弦波”+“1KHz正弦波”，波形及频谱如所⽰：图1被抽样信号波形及频谱⽰意图对抽样脉冲信号的考虑⼤家都知道，理想的抽样脉冲是⼀个⽆线窄的冲激信号，这样的信号在现实系统中是不存在的，实际的抽样脉冲信号总是有⼀定宽度的，很显然，这个脉冲宽度（简称脉宽）对抽样的结果是有影响的，这就是课本上讲的“孔径失真”，⽤不同的宽度的脉冲信号来抽样所带来的失真程度是不⼀样的，为了让⼤家能很好地理解和观察孔径失真现象，我们将抽样脉冲信号设计为脉宽可调的信号，在实验中⼤家可以⼀边调节脉冲宽度，⼀边从频域和时域两个⽅⾯来观察孔径失真现象。
低通信号的抽样定理
实验一抽样定理一．概述抽样的分类：(1) 根据信号是低通型的还是带通型的，抽样定理分低通抽样定理和带通抽样定理。
(2) 根据用来抽样的脉冲序列是等间隔的还是非等同隔的，又分均匀抽样定理和非均匀抽样。
(3) 根据抽样的脉冲序列是冲击序列还是非冲击序列，又可分理想抽样和实际抽样。
二．实验原理及其框图抽样定理是通信原理中十分重要的定理之一，是模拟信号数字化的理论基础。
低通型连续信号的抽样定理一个频带限制在内的时间连续信号，若以的间隔对它进行等间隔抽样，则将被所得到的抽样值完全确定。
原理框图图1 抽样说明：抽样过程中满足抽样定理时，PCM系统应无失真。
这一点与量化过程有本质区别。
量化是有失真的,只不过失真的大小可以控制。
三．实验步骤1、根据抽样原理，用Systemview 软件建立一个仿真电路，如下图所示：图2 仿真电路元件参数配制Token 0: 被采样的模拟信号—正弦波（频率=100Hz，电平=1V，相位=0）Token 2: 乘法器Token 5 抽样脉冲——窄脉宽矩形脉冲（脉宽=1us ）Token1，3: 模拟低通滤波器（截止频率=100 Hz ）Token 4,6,7: 观察点—分析窗（6频率=100Hz 电压=-1V）2、运行时间设置运行时间=0.3 秒采样频率=10,00 赫兹3、运行系统在Systemview 系统窗内运行该系统后，转到分析窗观察Token 5,6,8三个点的波形。
4、功率谱在分析窗绘出该系统调制后的功率谱。
四、实验报告1）观察实验波形：Token 0－被采样的模拟信号波形；Token 2－采样后波形；Token 3－恢复信号的波形。
2）整理波形，存入文档。
3）观察采样前后各信号的功率谱，结果存入文档，进行比较。
4）分析说明实验结果与理论值之间的差别。
5）改变参数配置，将所得不同结果存档后，与实验结果进行比较，说明参数改变对结果的影响。
通信原理MATLAB验证低通抽样定理实验报告
通信原理实验报告一、实验名称MATLAB验证低通抽样定理二、实验目的1、掌握抽样定理的工作原理。
2、通过MATLAB编程实现对抽样定理的验证，加深抽样定理的理解。
同时训练应用计算机分析问题的能力。
3、了解MATLAB软件，学习应用MATLAB软件的仿真技术。
它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。
4、计算在临界采样、过采样、欠采样三种不同条件下恢复信号的误差，并由此总结采样频率对信号恢复产生误差的影响，从而验证时域采样定理。
三、实验步骤及原理1、对连续信号进行等间隔采样形成采样信号，采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。
2、设连续信号的的最高频率为Fmax，如果采样频率Fs>2Fmax，那么采样信号可以唯一的恢复出原连续信号，否则Fs<=2Fmax会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。
四、实验内容1、画出连续时间信号的时域波形及其幅频特性曲线，信号为x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t)2、对信号进行采样，得到采样序列，画出采样频率分别为10Hz，20 Hz，50 Hz时的采样序列波形；3、对不同采样频率下的采样序列进行频谱分析，绘制其幅频曲线，对比各频率下采样序列和的幅频曲线有无差别。
4、对信号进行谱分析，观察与3中结果有无差别。
5、由采样序列恢复出连续时间信号，画出其时域波形，对比与原连续时间信号的时域波形。
五、实验仿真图(1) x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t)的时域波形及幅频特性曲线。
clear;close all;dt=0.05;t=-2:dt:2x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t);N=length(t);Y=fft(x)/N*2;fs=1/dt;df=fs/(N-1);f=(0:N-1)*df;plot(t,x)title('抽样时域波形')xlabel('t')grid;subplot(2,1,2)plot(f,abs(Y));title('抽样频域信号 |Y|');xlabel('f');grid;(2)采样频率分别为10Hz时的采样序列波形, 幅频特性曲线,以及由采样序列恢复出连续时间信号时域、频域波形；clear;close all;dt=0.1;t0=-2:0.01:2t=-2:dt:2ts1=0.01x0=cos(4*pi*t0)+1.5*sin(6*pi*t0)+0.5*cos(20*pi*t0);x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t);B=length(t0);Y2=fft(x0)/B*2;fs2=1/0.01;df2=fs2/(B-1);f2=(0:B-1)*df2;N=length(t);Y=fft(x)/N*2;fs=1/dt;df=fs/(N-1);f=(0:N-1)*df;tm=-50:ts1:50gt=sinc(fs*tm)st=sigexpand(x,dt/ts1)x3=conv(st,gt)A=length(tm(5001:5401));Y1=fft(x3(5001:5401))/A*2;fs1=1/ts1;df1=fs1/(A-1);f1=(0:A-1)*df1;subplot(3,2,1)plot(t0,x0)title('原始时域波形')xlabel('t')subplot(3,2,2)title('原始频域波形')xlabel('t')subplot(3,2,3)plot(t,x)title('抽样时域波形')xlabel('t')grid;subplot(3,2,4)plot(f,abs(Y));title('抽样频域信号 |Y|');xlabel('f');subplot(3,2,5)plot(t0,x3(5001:5401))title('恢复后的信号');xlabel('tm')subplot(3,2,6)plot(f1,abs(Y1));title('恢复频域信号 |Y1|');xlabel('f1');grid;（3）采样频率分别为20 Hz时的采样序列波形，幅频特性曲线,以及由采样序列恢复出连续时间信号时域、频域波形；clear;close all;dt=0.05;t0=-2:0.01:2t=-2:dt:2ts1=0.01x0=cos(4*pi*t0)+1.5*sin(6*pi*t0)+0.5*cos(20*pi*t0); x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t);B=length(t0);Y2=fft(x0)/B*2;fs2=1/0.01;df2=fs2/(B-1);f2=(0:B-1)*df2;N=length(t);Y=fft(x)/N*2;fs=1/dt;df=fs/(N-1);f=(0:N-1)*df;tm=-50:ts1:50gt=sinc(fs*tm)st=sigexpand(x,dt/ts1)x3=conv(st,gt)A=length(tm(5001:5401));Y1=fft(x3(5001:5401))/A*2;fs1=1/ts1;df1=fs1/(A-1);f1=(0:A-1)*df1;subplot(3,2,1)plot(t0,x0)title('原始时域波形')xlabel('t')subplot(3,2,2)plot(f2,abs(Y2))title('原始频域波形')xlabel('t')subplot(3,2,3)plot(t,x)title('抽样时域波形')xlabel('t')grid;subplot(3,2,4)plot(f,abs(Y));title('抽样频域信号 |Y|');xlabel('f');subplot(3,2,5)plot(t0,x3(5001:5401))title('恢复后的信号');xlabel('tm')subplot(3,2,6)plot(f1,abs(Y1));title('恢复频域信号 |Y1|');xlabel('f1');grid;(4)采样频率分别为50 Hz时的采样序列波形，幅频特性曲线,以及由采样序列恢复出连续时间信号时域、频域波形；；clear;close all;dt=0.02;t0=-2:0.01:2t=-2:dt:2ts1=0.01x0=cos(4*pi*t0)+1.5*sin(6*pi*t0)+0.5*cos(20*pi*t0);x=cos(4*pi*t)+1.5*sin(6*pi*t)+0.5*cos(20*pi*t);B=length(t0);Y2=fft(x0)/B*2;fs2=1/0.01;df2=fs2/(B-1);f2=(0:B-1)*df2;N=length(t);Y=fft(x)/N*2;fs=1/dt;df=fs/(N-1);f=(0:N-1)*df;tm=-50:ts1:50gt=sinc(fs*tm)st=sigexpand(x,dt/ts1)x3=conv(st,gt)A=length(tm(5001:5401));Y1=fft(x3(5001:5401))/A*2; fs1=1/ts1;df1=fs1/(A-1);f1=(0:A-1)*df1;subplot(3,2,1)plot(t0,x0)title('原始时域波形') xlabel('t')subplot(3,2,2)plot(f2,abs(Y2))title('原始频域波形') xlabel('t')subplot(3,2,3)plot(t,x)title('抽样时域波形') xlabel('t')grid;subplot(3,2,4)plot(f,abs(Y));title('抽样频域信号 |Y|'); xlabel('f');subplot(3,2,5)plot(t0,x3(5001:5401)) title('恢复后的信号'); xlabel('tm')subplot(3,2,6)plot(f1,abs(Y1));title('恢复频域信号 |Y1|'); xlabel('f1');grid;六、实验结论实验中对模拟信号进行采样，需要根据最高截止频率Fmax，按照采样定理的要求选择采样频率的两倍，即 Fs>2Fmax。
低通采样定理
成绩西安邮电大学《通信原理》软件仿真实验报告实验名称：低通型采样定理院系：通信与信息工程学院专业班级：通工1102学生姓名：董坤坤学号：03111057（班内序号）13指导教师：张明远报告日期：2013年9月26日●实验目的：1、掌握低通型采样定理；2、掌握理想采样、自然采样和瞬时采样的特点；3*、掌握混叠失真和孔径失真。
●仿真设计电路及系统参数设置：1、理想采样：2、平顶采样3、自然采样建议时间参数：No. of Samples =4096；Sample Rate = 20000Hz1、记录理想采样时信源、样值序列和恢复信号的波形和频谱；信源为截止频率200Hz 的低通型信号；其中图符0为信号源（单位冲激信号即()t δ，偏移量为0.05）； 图符1为截止频率200Hz ，极点个数为6的模拟低通滤波器； 图符2为采样器，采样频率2000Hz ；图符3为保持电路，Hold Value = Zero ，Gain = 1； 图符4为截止频率250Hz ，极点个数为5的模拟低通滤波器； 频谱选择 |FFT|；2*、自行调整参数，观测并记录混叠失真；3、记录自然采样时样值序列和恢复信号的波形和频谱；用于采样的矩形脉冲序列幅度1V ，频率2000Hz ；4*、调整矩形脉冲序列的占空比，观测并记录样值序列波形和频谱的变化；5、记录瞬时采样时样值序列和恢复信号的波形和频谱；保持电路Hold Value = Last Sample ，增益Gain = 1； 6*、自行调整参数，观测并记录孔径失真。
仿真波形及实验分析：在一系列离散点伤，对模拟信号抽取样值称为抽样，抽样过程可看为用周期性单位冲击脉冲和此模拟信号相乘。
抽样结果得到的是一系列周期性的冲激脉冲，其面积和模拟信号的取值成正比。
设一个连续模拟信号中最高频率小于H f ，则以间隔时间为Hf T 21≤的周期性冲激脉冲对它抽样时，该信号将被这些抽样值完全确定。