实验一：低通采样定理和内插与抽取实现a
南理工通信原理实验报告
目录实验一抽样定理实验 (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≤的周期性冲激脉冲对它抽样时，该信号将被这些抽样值完全确定。
低通型采样定理
低通型采样定理低通型采样定理是信号处理领域中的重要理论，它描述了在数字信号处理中对连续信号进行采样的方法和限制条件。
本文将详细介绍低通型采样定理的原理、应用以及一些相关概念。
一、低通型采样定理的原理低通型采样定理是由著名数学家香农（Claude Shannon）在1949年提出的。
它的基本原理是：如果一个连续信号的最高频率为f，则将其进行采样时，采样频率应该大于2f才能完全恢复原始信号。
也就是说，在采样过程中，采样频率必须大于信号的最高频率的两倍，才能保证采样后的信号不发生混叠现象。
二、低通型采样定理的应用低通型采样定理在实际应用中有着广泛的应用。
在音频和视频领域，低通型采样定理被广泛应用于数字音频、数字视频的采样和处理过程中。
通过合理的采样频率选择，可以在不损失信息的情况下，将连续信号转换为数字信号，从而实现信号的存储、传输和处理。
在通信领域，低通型采样定理也起着至关重要的作用。
在无线通信系统中，天线接收到的连续信号首先需要经过模数转换器（ADC）进行采样，然后才能进行数字信号处理和解调。
根据低通型采样定理，合理选择采样频率可以避免信号混叠，保证信号的完整性和准确性。
在生物医学领域，低通型采样定理也被广泛应用于生理信号的采样和处理过程中。
例如，心电图（ECG）信号和脑电图（EEG）信号都是连续信号，为了实现对这些信号的准确分析和诊断，需要首先对其进行采样，然后进行数字信号处理。
三、低通型采样定理的相关概念1. 采样频率：指每秒钟对连续信号进行采样的次数，用赫兹（Hz）表示。
根据低通型采样定理，采样频率应大于信号最高频率的两倍。
2. 采样定理：也称为奈奎斯特采样定理，是信号处理领域中的基本理论，指出连续信号在进行采样时，采样频率应大于信号最高频率的两倍，以避免信号混叠。
3. 混叠现象：也称为折叠现象，是指在采样过程中由于采样频率不满足低通型采样定理的要求，导致高频部分的信号频谱被折叠到低频区域，从而引起信号失真。
低通型采样定理
成绩《通信原理》软件仿真实验报告实验名称：低通型采样定理院系：通信与信息工程学院专业班级： 1学生姓名：学号：（班内序号）0指导教师：报告日期：2012年11月22日●实验目的：1、掌握低通型采样定理；2、掌握理想采样、自然采样和瞬时采样的特点；●知识要点：1、低通型采样定理；2、理想采样及其特点；3、自然采样及其特点；4、瞬时采样及其特点；●仿真模拟图：1、理想采样时信源、样值序列和恢复信号的波形和频谱：理想采样时信源：样值序列：恢复信号：分析:理想采样时的波形与原波形一样频谱也与原波形一致。
2、自然采样时样值序列和恢复信号的波形和频谱：恢复信号：分析:自然采样时的波形是与矩形脉冲相乘 但还是呈原波形的形状只是中间有了间隔而频谱形状会出现某段的频谱衰减或消失。
3、瞬时采样时样值序列和恢复信号的波形和频谱；瞬时采样时样值序列：瞬时采样恢复信号：分析:瞬时采样的波形与自然采样的波形比较像但与自然采样不同的是波形的顶部不是与原波形一样,而是水平直线频谱的顶部形状也会有变化也会出现衰减和消失的现象。
实验成绩评定一览表系统设计与模块布局系统设计合理，模块布局合理，线迹美观清楚系统设计合理，模块布局较合理，线迹清楚系统设计、模块布局较合理，线迹较清楚系统设计基本合理，模块布局较合理，线迹较清楚系统设计不够合理，模块布局较合理，线迹较清楚参数设置与仿真波形参数设置合理，仿真波形丰富、准确参数设置合理，仿真波形较丰富、较准确参数设置较合理，仿真波形较丰富参数设置较合理，仿真波形无缺失、无重大错误参数设置较合理，仿真波形有缺失参数设置不够合理，仿真波形有缺失或重大错误实验分析实验分析全面、准确、表达流畅实验分析较全面、基本无误、表述清楚实验分析基本正确、个别地方表述不清实验分析无原则性错误、表述不清楚实验分析有缺失或存在严重错误实验成绩。
低通采样原理及应用分析
低通采样原理及应用分析低通采样是一种信号处理技术，用于从一个连续时间信号中提取出其频率范围内的低频分量，并将其转换为离散时间信号。
该技术在信号处理、通信系统、音频处理、图像处理等领域有着广泛的应用。
低通采样的原理是基于奈奎斯特采样定理。
奈奎斯特采样定理指出，如果一个连续时间信号的最高频率分量为fmax，那么就需要以不小于2fmax的采样频率对其进行采样，才能完全还原该连续时间信号。
然而，在实际应用中，信号的最高频率通常不会很高，因此可以通过低通采样来降低采样频率，减少存储和处理的复杂度。
低通采样的基本步骤是首先对原始连续时间信号进行采样，然后通过低通滤波器将采样信号的高频成分去除，最后将滤波后的信号转换为离散时间信号。
通常使用的低通滤波器有抗混叠滤波器、Butterworth滤波器、卡尔曼滤波器等。
低通采样在实际应用中有着广泛的应用。
以下是一些低通采样的应用案例：1. 数字音频处理：在音频录制和音频压缩中，通常需要对音频信号进行低通采样。
通过降低采样频率，可以减少存储和传输带宽，同时保留音频信号的主要特征。
2. 图像处理：在数字相机和图像压缩中，通常需要对图像信号进行低通采样。
通过降低采样频率，可以减少图像数据的存储和传输量，同时保留图像的主要特征。
3. 通信系统：在无线通信系统中，为了适应带宽有限的信道，通常需要对调制后的信号进行低通采样。
通过降低采样频率，可以减少信号处理的复杂度和功耗。
4. 传感器信号处理：在传感器信号处理中，通常需要对传感器输出信号进行低通采样。
通过降低采样频率，可以减少信号处理的复杂度和功耗，同时保留传感器信号的主要特征。
总之，低通采样是一种常用的信号处理技术，可以从连续时间信号中提取出其频率范围内的低频分量，并将其转换为离散时间信号。
在音频处理、图像处理、通信系统、传感器信号处理等领域都有着广泛的应用。
通过降低采样频率，低通采样可以减少存储和处理的复杂度，同时保留信号的主要特征。
信号抽样与内插
3）分别在正常采样与欠采样条件下，配置各模块的参数（如信号源的频率，抽样脉冲的间隔，低通滤波器的截止频率等）。
4）在模型文件的菜单中选择Simulation－>Start，运行在正常采样、与欠采样条件下的仿真模型；
5）仿真结束后，打开示波器，观察在正常采样与欠采样条件下的仿真结果。
F=fft(z3,N); FF=F(1:m+1); F13=abs(FF);
F=fft(z4,N); FF=F(1:m+1); F14=abs(FF);
subplot(221);
plot(W,F11,'b',-W,F11,'b');
title('输入信号的幅频特性');
xlabel('频率（Rad/s）');
subplot(222);
plot(W,F12,'b',-W,F12,'b');
title('滤波后信号的幅频特性');
xlabel('频率（Rad/s）');
subplot(223);
plot(W,F13,'b',-W,F13,'b');
title('抽样后信号的幅频特性');
xlabel('频率（Rad/s）');
2、因输入信号是周期的，所以频谱都是离散的，以方波过抽样为例，基波为2*pi，频谱间隔为4*pi，抽样后将频谱搬移200*pi的整数倍。
3、抽样后信号包络是Sa函数，是因为实验采用方波抽样；另外，占空比越大，抽样后信号幅度越大。
通信原理抽样定理及其应用实验报告
当DDS信号源处于《PDM波1》状态，旋转SS01可改变取样脉冲的频率。示波器接在32P02上，可观察取样脉冲波形。
6．取样信号观察：
示波器接在32TP01上，可观察PAM取样信号，示波器接在32P03上，调节“PAM脉冲幅度调制”上的32W01可改变PAM信号传输信道的特性，PAM取样信号波形会发生改变。
实验1 抽样定理及其应用实验
一、实验目的
1．通过对模拟信号抽样的实验,加深对抽样定理的理解；
2．通过PAM调制实验，使学生能加深理解脉冲幅度调制的特点；
3．学习PAM调制硬件实现电路，掌握调整测试方法。
二、实验仪器
1．PAM脉冲调幅模块，位号：H（实物图片如下）
2．时钟与基带数据发生模块，位号：G（实物图片见第3页）
PAM实验原理：它采用模拟开关CD4066实现脉冲幅度调制。抽样脉冲序列为高电平时，模拟开关导通，有调制信号输出；抽样脉冲序列为低电平，模拟开关断开，无信号输出
图1-2 PAM信道仿真电路示意图
四、可调元件及测量点的作用
32P01：模拟信号输入连接铆孔。
32P02：抽样脉冲信号输入连接铆孔。
32TP01：输出的抽样后信号测试点。
32P03：经仿真信道传输后信号的输出连接铆孔。
32W01：仿真信道的特性调节电位器。
五、实验内容及步骤
1．插入有关实验模块：
在关闭系统电源的条件下，将“时钟与基带数据发生模块”、“PAM脉冲幅度调制模块”，插到底板“G、H”号的位置插座上（具体位置可见底板右下角的“实验模块位置分布表”）。注意模块插头与底板插座的防呆口一致，模块位号与底板位号的一致。
3．20M双踪示波器1台
4．频率计1台
5．小平口螺丝刀1只
抽取和内插
多速率信号处理及抽取和内插一：多速率信号处理1、在信号处理系统中有时需要不同的抽样率，这样做的目的有时是为了适应不同系统之间的级联，以利于信号的处理、编码、传输和存储，有时则是为了节省计算工作量。
数据速率的转换两种途径：1）数字信号→数模转换→模拟信号→模数转换→另一抽样率抽样2）数字信号处理→数字信号处理基本方法→抽样率转换目的：改变原有数字信号的频率方法：抽取和内插，低通滤波。
低通滤波：抽取和内插的前提条件是信号频带内没有频谱混叠，实现这一点需要用到低通滤波。
2、多速率滤波器-->具有线性相位的FIR滤波器。
常用的多速率滤波器：多速率FIR滤波器，积分梳状滤波器（CIC）和半带滤波器(HB);3、常用多速率信号处理结构第一级：CIC滤波器。
用于实现抽取和低通滤波第二级：fir实现的半带滤波器优点：工作在较低频率下，且滤波器参数得到优化，更容易以较低阶数实现，达到节省资源，降低功耗的目的。
二：抽取概念：使抽样率降低的转换。
1、整数倍抽取当信号的抽取数据量太大时，为了减少数据量以便于处理和计算，我们把抽样数据每隔（D-1）个取一个，这里D是一个整数。
这样的抽取称为整数抽取，D称为抽取因子。
2、抽取后结果：信号的频谱：信号的频谱周期降低1/D;信号的时域：信号的时域每D个少了（D-1）信号。
3、抗混叠滤波：在抽取前，对信号进行低通滤波，把信号的频带限制在抽样后频率的一半以下，这样，整数倍抽取的的问题就变成了一个低通滤波的问题。
信号时域图信号频域图程序运行后所得到的滤波前后信号的时域图，滤波器的频率响应图如上图。
从图中可以看出，经半带滤波器滤波后的信号，与原信号相比，波形没有改变，但抽样速率降低了一半；半带滤波器通阻带容限相同，具有严格线性相位。
三：内插概念：使抽样率升高的转换。
1、整数倍内插：在已知的相邻抽样点之间等间隔插入（I-1）个零值点。
然后进行低通滤波，即可求得I倍内插的结果。
2、内插后结果：信号的时域：已知抽样序列的两相邻抽样点之间等间隔多了I-1个值信号的频谱：信号的频谱周期增加了I倍。
通信原理实验-抽样定理
学生实验报告）实际上，考虑到低通滤波器特性不可能理想，对最高频率为3400Hz的语言信号，通常采用8KHz 抽样频率，这样可以留出1200Hz的防卫带。
见图4。
如果fs＜fH，就会出现频谱混迭的现象，如图5所示。
在验证抽样定理的实验中，我们用单一频率fH的正弦波来代替实际的语音信号。
采用标准抽样频率fs＝8KHZ。
改变音频信号的频率fH，分别观察不同频率时，抽样序列和低通滤波器的输出信号，体会抽样定理的正确性。
验证抽样定理的实验方框图如图6所示。
在图8中，连接（8）和（14），就构成了抽样定理实验电路。
由图6可知。
用一低通滤波器即可实现对模拟信号的恢复。
为了便于观察，解调电路由射随、低通滤波器和放大器组成，低通滤波器的截止频率为3400HZ2、多路脉冲调幅系统中的路际串话~多路脉冲调幅的实验方框图如图7所示。
在图8中，连接（8）和（11）、（13）和（14）就构成了多路脉冲调幅实验电路。
分路抽样电路的作用是：将在时间上连续的语音信号经脉冲抽样形成时间上离散的脉冲调幅信号。
N路抽样脉冲在时间上是互不交叉、顺序排列的。
各路的抽样信号在多路汇接的公共负载上相加便形成合路的脉冲调幅信号。
本实验设置了两路分路抽样电路。
多路脉冲调幅信号进入接收端后，由分路选通脉冲分离成n路，亦即还原出单路PAM信号。
图7 多路脉冲调幅实验框图冲通过话路低通滤波器后，低通滤波器输出信号的幅度很小。
这样大的衰减带来的后果是严重的。
但是，在分路选通后加入保持电容，可使分路后的PAM信号展宽到100％的占空比，从而解决信号幅度衰减大的问题。
但我们知道平顶抽样将引起固有的频率失真。
PAM信号在时间上是离散的，但是幅度上趋势连续的。
而在PAM系统里，PAM信只有在被量化和编码后才有传输的可能。
本实验仅提供一个PAM系统的简单模式。
3、多路脉冲调幅系统中的路标串话路际串话是衡量多路系统的重要指标之一。
路际串话是指在同一时分多路系统中，某一路或某几路的通话信号串扰到其它话路上去，这样就产生了同一端机中各路通话之间的串话。
低通抽样定理的理解
低通抽样定理的理解带通信号抽样速率的一统性研究（篇一）:低通抽样定理的理解写在前：这是笔者第一篇CSDN博客，希望未来保持更新和分享，督促自己学习和思考，也希望能有助于小伙伴们的解惑，内容是关于通信和信号处理方面的。
我现在基本还是通信和信息的小白，会不断学习。
如果内容有出错或者不够严谨的地方希望指正，相互交流，共同进步。
本篇是《带通信号抽样速率的一统性研究》的篇一，基础知识部分。
从频域和时域两个角度对低通抽样定理进行了探讨和理解。
在通信与信息系统中需要将一个时间连续信号（模拟信号，时间和幅度上均连续）通过抽样（或称为采样）来转换为时间离散信号（时间离散，幅度不离散），继而通过量化、编码得到数字信号（时间离散，幅度离散）。
数字信号相对于模拟信号具有便于存储、处理、传输等优点。
我们只探讨抽样过程和其对应的信号恢复过程。
经过抽样得到的离散序列需要能够还原出原来的模拟信号，这也就意味着抽样过程完整地保留了原信号的信息。
抽样过程追求高效，离散序列的长度被期待尽量短（对应抽样速率尽量低）。
抽样定理包含两个内容：低通抽样定理和带通抽样定理。
我们先从低通抽样定理入手。
低通抽样定理：给定最高非零频率为fH f_HfH的带限信号m(t) m(t)m(t)，如果取抽样间隔Ts<1/(2fH) T_s<1/(2f_H)Ts<1/(2fH )（即抽样速率fs>2fH f_s>2f_Hfs>2fH，有些资料的表达为fs≥2f H f_s\ge2f_Hfs≥2fH），则m(t) m(t)m(t)由其样值序列{mn=m(nTs)，n为整数} \big\{m_n=m(nT_s)，n为整数\big\}{mn=m(nTs)，n为整数}唯一确定，即m(t)⟷只要fs>2fH{mn，n=0,±1,±2,…} m(t)\stackrel{只要f_s>2f_H}{\longleftrightarrow}\big\{m_n，n=0,\pm1,\pm2,…\big\}m(t)⟷只要fs>2fH{mn，n=0,±1,±2,…}低通抽样定理中，抽样速率必须大于2fH 2f_H2fH，该频率2fH 2f_H2fH通常称为奈奎斯特频率。
低通抽样定理
(c) 抽样频率为200Hz (b) 抽样频率为128Hz (a) 抽样频率为100Hz
图4 抽样实验
5. 实际应用中要注意的问题 (1) 抽样前，加截止频率为 F m 的低通滤波器，滤除 F m 赫以上的 频谱成分，从而消除混叠现象和避免由此引起的失真。
(2) 抽样时，抽样速率 f s 要比 2 Fm 大，一般取 f s (2.5 3) Fm 。 因为实际应用中使信号恢复的低通滤波器不可能是理想的。
(3) 抽样时，采用的抽样脉冲序列一般都是高度有限，宽度很窄 的脉冲序列。因为在实际应用中，理想抽样所需的周期性单位冲 激脉冲是不可能实现的。
请多指正 谢谢各位老师
1 2 Fm 秒的间隔对其进行等间隔抽样，则所得的样值可以完全确
定原信号。由于抽样时间间隔相等，此定理又称均匀抽样定理。
定理表明：若 f (t ) 的频谱在某一角频率 m 以上为零，则 f (t )中 的全部信息完全包含在其间隔不大于 1 2 Fm 的均匀抽样序列里。 换句话说，在信号最高频率分量的每一个周期内起码应抽样两次 1 (即 Tm 2Ts )。 ，Tm Fm
3. 低通抽样定理证明 (1) 理论分析 假设抽样脉冲序列是周期为 Ts 的单位冲激函数，其表示式为
T (t )
s
n
(t nT )
s
抽样过程实际上是用信号 f (t ) 与抽样脉冲序列相乘，因此已抽样 信号 f s (t ) 表示式为
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验一：低通采样定理和内插与抽取实现一．实验目的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= 信号组分。