实验三 信号采样与重建(参考答案)
信号处理实验三报告
信号处理实验三报告实验三：时域信号的采样与重构一、实验目的1.学习使用示波器进行时域信号采样；2.学习时域信号重构的方法。
二、实验器材1.数字示波器；2.函数发生器；3.电缆。
三、实验原理1.时域信号的采样时域信号的采样是将连续时间的信号转换为离散时间的信号。
采样过程可以理解为在时间轴上以一定的时间间隔取样，得到采样点的幅值。
采样后的信号可以用离散时间信号表示。
2. Nyquist采样定理Nyquist采样定理指出，要恢复一个最高频率为f的连续时间信号，采样频率必须大于2f，即采样定理为Fs > 2f。
这是由于频谱中的高频分量蕴含着较大的信息量，必须以足够高的采样频率进行采样，否则会出现混叠现象。
3.时域信号的重构时域信号的重构是将采样得到的离散时间信号重新转化为连续时间信号的过程。
重构的方法主要有零阶保持插值、线性插值和插值滤波器等。
实验步骤1.连接示波器和函数发生器。
将函数发生器的输出端通过电缆与示波器的输入端连接。
2.设置函数发生器的频率为1kHz，并选择一个适当的幅度。
3.设置示波器的水平和垂直缩放，使信号在示波器的屏幕上能够完整显示。
4.调节示波器的触发方式和触发电平，使信号的波形稳定。
5.通过示波器的采样功能，进行信号的采样。
选择适当的采样率，观察采样得到的离散时间信号。
6. 根据Nyquist采样定理，选择适当的采样率进行采样，并进行离散时间信号的重构。
选择不同的重构方法，如零阶保持插值和线性插值，观察重构后的信号与原信号的差异。
实验结果1.通过示波器的采样功能，得到了采样频率为1kHz的离散时间信号。
2.通过零阶保持插值和线性插值的方法进行重构，观察到重构后的信号与原信号的差异。
可以发现，零阶保持插值会导致信号的平滑度降低，而线性插值能够更好地重构原信号。
实验分析1. 通过实验结果可以验证Nyquist采样定理的正确性。
当采样频率小于2f时，会出现混叠现象，无法正确恢复原信号。
信号的采样与恢复
实验报告课程名称：信号分析与处理 指导老师： 成绩： 实验名称：信号的采样与恢复 实验类型： 同组学生姓名： 一、实验目的和要求（必填） 二、实验内容和原理（必填） 三、主要仪器设备（必填） 四、操作方法和实验步骤 五、实验数据记录和处理 六、实验结果与分析（必填） 七、讨论、心得一、实验目的和要求1. 了解信号的采样方法与过程以及信号恢复的方法。
2. 验证采样定理。
二、实验内容和原理 2.1信号的自然采样采样信号为周期Ts ，宽度τ的矩形脉冲信号S(t)。
s(t)的傅里叶变换为： 2(t)Sa()()2s s sn S n T ωτπτδωω+∞-∞=-∑ 采样的过程可以视为两个信号相乘：()()()s f t f t s t =在频域中，1()()()2Sa()()2s s s s F F S n F n T ωωωπωττωω+∞-∞=*=-∑可以看到自然采样后的频谱除了左右平移采样信号的角频率ωs 外，还按取样函数Sa(x)的规律衰减。
时域采样定理：如果采样信号的频率为fs ，原信号的最大频率为f m ，为了采样后信号的频谱不混叠，需要有fs ≥2f m 。
2.2信号的恢复在不发生频谱混叠的时候，将信号通过的低通滤波器，理论上可以完全恢复原信号。
低通滤波器的截止频率略大于fm，即“频谱加窗”的方法。
如果发生了频谱混叠，则原信号的频谱不能完全被恢复，通过低通滤波器后输出的信号将产生失真。
本实验分别用500Hz三角波和正弦波作为输入信号，占空比50%和10%的0.4kHz、1kHz、2kHz、5kHz、10kHz的矩形脉冲作为采样信号，使用截止频率1kHz以及2kHz的低通滤波器，观察输出波形，验证采样定理。
实验中，受自然采样、实验滤波器效果的限制，恢复后的波形难免都会有失真。
三、主要仪器设备PC一台、myDAQ设备一套、面包板一块、导线、电容、电阻若干。
四、操作方法和实验步骤1.编辑波形文件：正弦波峰峰值4V、频率500Hz，与10kHz、幅值1V、占空比50%的方波相乘，保存波形文件。
信号的采样与恢复
信号的采样与恢复实验一、任务与目的1. 熟悉信号的采样与恢复的过程。
2. 学习和掌握采样定理。
3. 了解采样频率对信号恢复的影响。
二、原理（条件）PC机一台，TD－SAS系列教学实验系统一套。
1. 采样定理采样定理论述了在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值表示。
这些值包含了该连续信号全部信息，利用这些值可以恢复原信号。
采样定理是连续时间信号与离散时间信号之间的桥梁。
采样定理：对于一个具有有限频谱，且最高频率为ωmax的连续信号进行采样，当采样频率ωs满足ωs>=ωmax时，采样信号能够无失真地恢复出原信号。
三角波信号的采样如图4-1-1所示。
图4-1-1信号的采样2. 采样信号的频谱连续周期信号经过周期矩形脉冲抽样后，抽样信号的频谱为它包含了原信号频谱以及重复周期为的原信号频谱的搬移，且幅度按规律变化。
所以抽样信号的频谱便是原信号频谱的周期性拓延。
某频带有限信号被采样前后频谱如图4-1-2。
图4-1-2 限带信号采样前后频谱从图中可以看出，当ωs ≥2Bf 时拓延的频谱不会与原信号的频谱发生重叠。
这样只需要利用截止频率适当的滤波器便可以恢复出原信号。
3. 采样信号的恢复将采样信号恢复成原信号，可以用低通滤波器。
低通滤波器的截止频率f c 应当满足f max ≤f c ≤f x -f max 。
实验中采用的低通滤波器原理图如图4-1-3所示，其截止频率固定为1802f Hz RCπ=≈图4-1-3 滤波器电路4. 单元构成本实验电路由脉冲采样电路和滤波器两个部分构成，滤波器部分不再赘述。
其中的采样保持部分电路由一片CD4052完成。
此电路由两个输入端，其中IN1端输入被采样信号，Pu 端输入采样脉冲，经过采样后的信号如图4-1-1所示。
三、内容与步骤本实验在脉冲采样与恢复单元完成。
1. 信号的采样（1）使信号发生器第一路输出幅值3V、频率10Hz的三角波信号；第二路输出幅值5V，频率100Hz、占空比50％的脉冲信号。
信号采样与重建
1.软件介绍MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩形计算、视化以线性动态线性系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多领域一面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。
MATLAB是矩阵实验室（Matrix Laboratory）之意。
除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。
经过不断完善MATLAB已经发展成为适合多学科，多种工作平台的功能强大大大型软件。
成为线性代数,自动控制理论，数理统计，数字信号处理，时间序列分析，动态系统仿真等高级课程的基本教学工具。
MTLAB的语言特点：(1)语言简洁紧凑，使用方便灵活，库函数极其丰富。
(2)运算符丰富。
(3)MATLAB既具有结构化的控制语句（如for循环，while循环，break语句和if语句），又有面向对象编程的特性。
(4)程序限制不严格，程序设计自由度大。
(5)MATLAB的图形功能强大。
(6)MATLAB的缺点是，它和其他高级程序相比，程序的执行速度较慢。
由于MATLAB的程序不用编译等预处理，也不生成可执行文件，程序为解释执行，所以速度较慢。
2.课程设计的方案2.1课程设计的原理2.1.1连续信号的采样定理模拟信号经过 (A/D) 变换转换为数字信号的过程称为采样，信号采样后其频谱产生了周期延拓，每隔一个采样频率 fs ，重复出现一次。
为保证采样后信号的频谱形状不失真，采样频率必须大于信号中最高频率成分的两倍，这称之为采样定理。
时域采样定理从采样信号恢复原信号必需满足两个条件:(1)必须是带限信号，其频谱函数在＞各处为零；（对信号的要求，即只有带限信号才能适用采样定理。
抽样定理与信号恢复
抽样定理与信号恢复一、实验目的1、观察离散信号频谱，了解其频谱特点。
2、验证抽样定理并恢复原信号。
二、实验仪器1、双踪示波器 1台 2块 31块 4块三、实验原理1、离散信号不仅可从离散信号源获得，而且也可从连续信号抽样获得。
抽样信号Fs(t)= F(t)·S(t)。
其中F(t)为连续信号（例如三角波），S(t)是周期为Ts 的矩形窄脉冲。
Ts 又称抽样间隔，Fs=1Ts称抽样频率，Fs(t)为抽样信号波形。
F(t)、S(t)、Fs(t)波形如图5-1。
图5-1 连续信号抽样过程将连续信号用周期性矩形脉冲抽样而得到抽样信号，可通过抽样器来实现。
2、连续周期信号经周期矩形脉冲抽样后，抽样信号的频谱()()()2s S a s m m A F j S F m T ωττωωω+∞=-∞=⋅-∑它包含了原信号频谱以及重复周期为fs （f s =ωs /2л）、幅度按A τT Sa （m ωs τ/2）规律变化的原信号频谱，即抽样信号的频谱是原信号频谱的周期性延拓。
因此，抽样信号占有的频带比原信号频带宽得多。
以三角波被矩形脉冲抽样为例。
三角波的频谱：1124()()()S k k k E F j A k k k ωπσωωσωωπ∞∞=-∞=-∞=-=-∑∑ 这里我们取三角波的有效带宽为31ω。
频谱图如5-2所示。
F(f)图5-2 三角波频谱抽样信号的频谱：12ω1()4()(ω)2s S as k m m A F j E S k m T kττωσωωπ∞=-∞=-∞=⋅--∑ 取三角波的有效带宽为31ω，其抽样信号频谱如图5-3所示。
BfFs(f)图5-3 抽样信号频谱图如果离散信号是由周期连续信号抽样而得，则其频谱的测量与周期连续信号方法相同，但应注意频谱的周期性延拓。
3、抽样信号在一定条件下可以恢复出原信号，其条件是2s f f B ≥，其中s f 为抽样频率，f B 为原信号占有频带宽度。
信号的采样与恢复
深圳大学实验报告课程名称：信号与系统实验项目名称：信号的采样与恢复学院：信息工程专业：电子信息指导教师：报告人：学号：班级：实验时间：实验报告提交时间：教务部制一、实验目的和要求1、了解信号的采样方法与过程以及信号恢复的方法。
2、验证采样定理。
二、实验内容和原理实验原理1、离散时间信号可以从离散信号源获得，也可以从连续时间信号采样而得。
采样信号x s （t ）可以看成连续信号x （t ）和一组开关函数s （t ）的乘积。
s （t ）是一组周期性窄脉冲，如图2－5－1，T s 称为采样周期，其倒数f s ＝1/T s 称采样频率。
图2－5－1 矩形采样信号对采样信号进行傅里叶分析可知，采样信号的频率包括了原连续信号以及无限个经过平移的原信号频率。
平移的频率等于采样频率f s 及其谐波频率2f s 、3f s ……。
当采样信号是周期性窄脉冲时，平移后的频率幅度按sinx/x 规律衰减。
采样信号的频谱是原信号频谱周期的延拓，它占有的频带要比原信号频谱宽得多。
2、采样信号在一定条件下可以恢复到原信号。
只要用一截止频率等于原信号频谱中最高频率f n 的低通滤波器，滤除高频分量，经滤波后得到的信号包含了原信号频谱的全部内容，故在低通滤波器输出端可以得到恢复后的原信号。
3、原信号得以恢复的条件是f s ≥2f max ，f s 为采样频率，f max 为原信号的最高频率。
当fs ＜2 f max 时，采样信号的频谱会发生混迭，从发生混迭后的频谱中无法用低通滤波器获得原信号频谱的全部内容。
在实际使用中，仅包含有限频率的信号是极少的，因此即使f s ＝2 f max ，恢复后的信号失真还是难免的。
实验中选用f s ＜2 f max 、f s ＝2 f max 、f s ＞2 f max 三种采样频率对连续信号进行采样，以验证采样定理：要使信号采样后能不失真地还原，采样频率f s 必须大于信号最高频率的两倍。
4、连续信号的采样和采样信号的复原原理框图如图2－5－2所示。
信号与系统课程设计--信号的采样与恢复
一、设计目的与要求1、设计目的通过本课程设计，主要训练和培养学生综合应用所学过的信号及信息处理等课程的相关知识，独立完成信号仿真及信号处理的能力。
包括：查阅资料、合理性的设计、分析和解决实际问题的能力，数学仿真软件Matlab和C语言程序设计的学习和应用，培养规范化书写说明书的能力。
2、设计要求设有一信号Xa(t)=EXP-1000|t|,计算傅立叶变换，分析其频谱，并在精度为1/1000的条件下，分别取采样频率为F=5000Hz,F=1000Hz,绘出对应的采样信号的时域信号波形频谱图。
（1）实现信号时域分析和频谱分析以及滤波器等有关Matlab函数。
（2）写好总结、程序、图表、原理、结果分析。
二、设计原理框图三、设计原理本次课程设计主要涉及采样定理、傅里叶变换、信号时域分析和频谱分析的相关内容的相关知识。
1.采样定理设连续信号)(t x a 属带限信号，最高截止频率为c Ω，如果采样角频率c s Ω≥Ω2，那么让采样性信号)(t x a ∧通过一个增益为T 、截止频率为2/s Ω的理想低通滤波器，可以唯一地恢复出原连续信号)(t x a 。
否则，c s Ω<Ω2会造成采样信号中的频谱混叠现象，不可能无失真地恢复原连续信号。
对连续信号进行等间隔采样形成采样信号，对其进行傅里叶变换可以发现采样信号的频谱是原连续信号的频谱以采样频率s Ω为周期进行周期性的延拓形成的。
对模拟信号进行采样可以看做一个模拟信号通过一个电子开关S ，设电子开关每隔周期T 和上一次，每次和上的时间为τ，在电子开关的输出端得到采样信号x^a(t)。
图1 对模拟信号进行采样2、傅里叶变换(1)对于一个非周期函数f（t），如果在（－∞，＋∞）满足下列条件：①、f（t）在任一有限区间上满足狄利克雷条件；②、f（t）在（－∞，＋∞）上绝对可积（如下积分收敛），即：（1）则有下式的傅立叶积分成立：（2）(2)f(t)满足傅立叶积分定理条件时，下图①式的积分运算称为f(t)的傅立叶变换，②式的积分运算叫做F(ω)的傅立叶逆变换。
连续信号的采样与恢复实验报告
连续信号的采样与恢复实验报告实验报告：连续信号的采样与恢复一、实验目的：1.了解连续信号的采样原理和采样定理；2.理解采样后信号的频谱特性；3.掌握信号恢复的方法。
二、实验原理：采样定理：对于频谱带宽有限的信号，为了保证采样信号不发生混叠现象，必须满足采样频率大于信号频谱的最高分量频率的两倍。
三、实验器材：1.信号发生器；2.示波器；3.编码器；4.数字示波器；5.连接线。
四、实验步骤及结果：1.首先使用信号发生器产生频率为1kHz、幅值为5V的正弦信号作为待采样信号；2.将信号发生器输出的信号连接至示波器进行观察；3.将示波器输出信号连接至编码器进行信号的采样；4.将编码器的输出信号连接至数字示波器，观察离散采样值；5.对离散采样值进行信号恢复，使用零阶保持、线性插值和兰特尔-曼豪姆插值三种恢复方法；6.将恢复后的信号与原信号进行比较，观察恢复的效果。
实验结果：在示波器上观察到频率为1kHz、幅值为5V的正弦信号。
数字示波器上显示出了一系列离散的采样值。
通过零阶保持、线性插值和兰特尔-曼豪姆插值三种方法进行信号恢复后，观察到恢复的信号与原信号基本一致。
五、实验分析：1.信号恢复的效果受到采样频率和采样幅值的影响，采样频率过低或采样幅值过小都会造成信号失真；2.零阶保持方法可以保持离散信号的幅值不变，但是无法恢复信号的高频分量；3.线性插值可以恢复少量的高频分量，但是如果信号存在高频噪声或非线性失真，会导致恢复后信号的质量下降；4.兰特尔-曼豪姆插值是一种高阶插值方法，能够更好地恢复信号的高频分量，但是计算量较大。
六、实验总结：通过本次实验，我了解了连续信号的采样原理和恢复方法，掌握了采样频率的要求和恢复过程中常用的插值方法。
实验中，我观察到了采样信号和恢复信号的特性，并进行了比较分析。
实验结果表明，在合适的采样条件和恢复方法下，可以有效地采样和恢复信号。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2：设连续信号 （－1000|t|）时 ：设连续信号x(t)=exp（－ （－长 ） t=-0.005 : 0.00005 : 0.005; %10mS长，201点密度 点密度 A、求连续时间傅利叶变换 、求连续时间傅利叶变换X(jw) x=exp(－1000*abs(t))；％共有 ；％共有 － ；％共有201点 点 Wmax=2*pi*2000; ％只要算到频率 只要算到频率2000Hz止 止 M=500; m=0:1:M; W=m*Wmax/M; ％得到 得到W(1)~W(501)共501个频谱点值。 个频谱点值。 共 个频谱点值 Xa=x*exp(－j* t’ *W)*0.00005; － 计算谱501点Xa(1)~ Xa(501). ％计算谱 点 Xa=real(Xa); Xa=[fliplr(Xa),Xa(2:501)];%左折但 左折但Xa(1)不动 左折但 不动 W=[－fliplr(W),W(2:501)]; ％频点也左折。 频点也左折。 － figure(1); plot(t*1000,x, ’r’); grid; %毫秒 毫秒 figure(2); plot(W/(2*pi*1000),Xa*1000, ’r’); ％幅度放大 幅度放大1000
两个频率分量 采样频率200 采样频率
2：设连续信号 ：设连续信号x(t)=exp（－ （－1000|t|）时 （－ ）
A、求连续傅利叶变换X(jw)（用很密的点） 、求连续傅利叶变换 （用很密的点）
Xa ( jω ) = ∫ x ( t )e dt = ∫ e e dt + ∫ e − 1000 t e − jω t dt 0 −∞ −∞ 1 1 0 ∞ e j (1000 − jω ) t e − j (1000 + jω ) t = − 0 −∞ 1000 − jω 1000 + jω 0相位，相对 相位， 相位 1 1 0.002 = + = 幅度直流处最 ω 2 1000 − jω 1000 + jω 1+ ( ) 大0.002 1000
频谱
波形
幅度放大1000倍后 倍后 幅度放大
2：设连续信号x(t)=exp（－ ：设连续信号 （－1000|t|）时 （－ ） t=-0.005 : 0.00005 先采样，在求离散 : 0.005; 原信号条件不变 B、Fs=5000,和1000Hz先采样 在求离散DFT频谱。 先采样， %原信号条件不变 频谱。 、 和 频谱
二、实验内容
复合时看不出频率
1、对连续时间信号 ＝sin(24πt)+ sin(40πt)，它 、对连续时间信号y(t)＝ ，
两个等幅度分量。 MATLAB作图 有 12Hz 和 20Hz 两个等幅度分量。用MATLAB作图 法求出Nyquist频率2fmax Nyquist频率2fmax。 法求出Nyquist频率2fmax。
减低采样率
频谱， 频谱，注意两边衰减 没有接近0。 没有接近 。
时域采样序列。 时域采样序列。
D 信号的重构，恢复。 信号的重构，恢复。
Ts=0.0002; n=－25:1:25; Fs=5000； ； －
t=-0.005 : 0.00005 : 0.005; nTs＝n*Ts； ＝ ；
x=exp(－1000*abs(nTs))；％序列 － ；％序列 ；％ xa＝spline（nTs，x，t）；三次样条函数重构。 ）；三次样条函数重构 ＝ （ ， ， ）；三次样条函数重构。 error1＝max（abs（xa－exp（－ （－1000*abs(t)))); ＝ （ （ － （－ figure(1); plot(t*1000,xa,’r’); text(2,0.5,strcat(‘error=’,num2str(error1))); hold on; ％已绘制的曲线保留。 已绘制的曲线保留。 stem (n*Ts*1000,x); hold off;
实验三 信号采样与重建
Lecture Script 2009.5.5
一、实验目的
1、进一步学习MATLAB的函数及其表示。 、进一步学习 的函数及其表示。 的函数及其表示 2、掌握及验证信号的 采样定理。 、掌握及验证信号的SHANNON采样定理。 采样定理 3、由采样序列重构恢复原信号。 、由采样序列重构恢复原信号。
最大误差
ቤተ መጻሕፍቲ ባይዱ
有较大失真
三、思考题： 思考题： 连续时间信号的傅利叶变换matlab求法，这 求法， ①连续时间信号的傅利叶变换 求法 里采用的近似公式是什么？ 里采用的近似公式是什么？ 从序列重构连续信号所采用的matlab函数是什 ②从序列重构连续信号所采用的 函数是什 么？ 采样定理中的信号Nyquist频率是指什 ③shannon采样定理中的信号 采样定理中的信号 频率是指什 与采样频率有什么不同？ 么？与采样频率有什么不同？
最大误差 逐点进行连续的减重 构的
D 信号的重构时，不能恢复。 信号的重构时，不能恢复。
Ts=0.001; n=－5:1:5; Fs=1000； － ；
t=-0.005 : 0.00005 : 0.005; nTs＝n*Ts； ＝ ；
x=exp(－1000*abs(nTs))；％序列 － ；％序列 ；％ xa＝spline（nTs，x，t）；三次样条函数法重构。 ）；三次样条函数法重构 ＝ （ ， ， ）；三次样条函数法重构。 error2＝max（abs（xa－exp（－ （－1000*abs(t)))); ＝ （ （ － （－ figure(1); plot(t*1000,xa,’r’); text(2,0.5,strcat(‘error=’,num2str(error2))); hold on; stem (n*Ts*1000,x); hold off;
所以：满足抽样定理，从而可以用 代替x(t), 所以：满足抽样定理，从而可以用xa[n]代替 代替 ， 范围[－ 之间每间隔5*10－5，20KHz，共200点 ]， 范围 －0.005到0.005之间每间隔 到 之间每间隔 ， 点 ， 到起点0.005时，e－5＝0.0067，已经足够小。 从0到起点 到起点 时 ，已经足够小。
四、实验报告要求： 实验报告要求： ①简述实验目的 ②按实验项目手工写出程序 打印实验图形结果，并粘贴于程序对应处。 ③打印实验图形结果，并粘贴于程序对应处。 ④简要回答3个思考题 简要回答 个思考题
%Analog Signal Dt=0.005; t=-1:Dt:1; % sample frequency 200 Hz y=sin(24*pi*t)+sin(40*pi*t); subplot(3,1,1); plot(t*4,y); xlabel('t in 1/4sec.'); ylabel('y(t)'); 分离成两个分量 title('Analog Signal'); subplot(3,1,2); plot(t*12,sin(24*pi*t)); xlabel('t in 1/12sec.'); ylabel('sin(24*pi*t)'); %f=12Hz 一般的分离成很多个正余弦， 一般的分离成很多个正余弦， subplot(3,1,3); 最高频率的2倍就是 倍就是。 最高频率的 倍就是。 plot(t*20,sin(40*pi*t)); xlabel('t in 1/20sec.'); ylabel('sin(40*pi*t)'); %f=20Hz
另一解法频谱法（适合很多个正弦分量情况） 另一解法频谱法（适合很多个正弦分量情况）
t=0:0.01:1; % T＝0.01 Sec，sample frequency 100 Hz ＝ ， y=sin(24*pi*t)+sin(40*pi*t); Y=fft(y); f=(0:length(y)-1)/length(y)*100; ％每点对应频率 figure(1); plot(f,abs(Y)); grid ;%幅度 幅度 figure(2); plot(f,unwrap(angle(Y))*180/pi); grid ; 相位曲线不作折叠，而是全部展开的画法。 相位曲线不作折叠，而是全部展开的画法。
x1=exp(－1000*abs(t))； － ； Ts=0.0002; n=－25:1:25; %Fs=5000, － x=exp(－1000*abs(n*Ts))； － ； M=500;m=0:1:M; W=pi*m/M; ％频率点 Xa=x*exp(－j* n’ *W）; － ） ％计算离散谱 Xa=real(Xa); Xa=[fliplr(Xa),Xa(2:M+1)]; W=[－fliplr(W),W(2:M+1)]; － figure(1); plot(t*1000,x1,’r’); hold on; stem (n*Ts*1000,x); hold off; figure(2); plot(W/pi,Xa);
− jω t 1000 t − jω t ∞ 0 ∞
When
f 0 = 2000 Hz; ω = 4000π , Xa( jω ) = 1.3*10−5 → 0 1 −5 −5 ; if f = 20kHz >> 2 f 0 = 4000 Hz; ∆t = 5*10 < 25*10 = 2* 2000
频谱。 频谱。
序列
C、减小Fs=1000Hz，先采样，再求离散 、减小 频谱。 ，先采样，再求离散DFT频谱。 频谱 t=-0.005 : 0.00005 : 0.005; %连续的条件不变 连续的条件不变 x1=exp(－1000*abs(t))； － ； Ts=0.001; n=－5:1:5; %Fs=1000,Ts＝1mS － ＝ x=exp(－1000*abs(n*Ts))；％ ；％n*Ts=－0.005～0.005 － ；％ － ～ M=500;m=0:1:M; W=pi*m/M; ％频率点 Xa=x*exp(－j* n’ *W）; － ） ％计算离散谱 Xa=real(Xa); Xa=[ fliplr(Xa), Xa(2:M+1)]; W=[－fliplr(W),W(2:M+1)]; － figure(1); plot(t*1000,x1,’r’); hold on; stem (n*Ts*1000,x); hold off; figure(2); plot(W/pi,Xa);