精华内容
下载资源
问答
  • 抽样信号为什么是模拟信号
    千次阅读
    2022-02-16 15:09:50





    一、模拟信号、离散时间信号、数字信号



    时间是 连续 的 , 幅度也是 连续 的 , 该信号是 模拟信号连续时间信号 ;

    时间是 离散 的 , 幅度是 连续 的 , 该信号是 离散时间信号 ;

    时间是 离散 的 , 幅度是 离散 的 , 该信号是 数字信号 ;





    二、采样导致时间离散



    采样导致时间离散 :

    时间离散是因为采样形成的 , 如音频采样 44100 Hz , 一秒钟采样 44100 个音频样本 , 即使是这样 , 时间也是离散的 ;

    时间是离散的 , 如果幅度是连续的 , 说明可以取值 f ( t ) f(t) f(t) 对应的 y y y 轴上任意一点 , 该信号是 离散时间信号 ;





    三、量化导致幅度离散



    量化导致幅度离散 :

    幅度的离散 , 指的是 , 信号的采样值 , 只能是给定的几个值 , 如音频采样位数 8 位 , 则音频采样的取值范围只能是 0 ~ 2 8 2^8 28 之间的值 , 不能取值 0.0001 这个值 ;


    下图中 左侧的信号 是 离散时间信号 , 其时间是离散的 , 时间取值只能是 整数 , 但是 幅度取值可以是 52 ~ 64 之间的任意值 , 如 0.0001 取值 ;

    下图中 右侧的信号 是 数字信号 , 其时间和幅度都是离散的 , 时间取值只能是 整数 , 幅度取值也只能是整数 ;

    将 连续的幅度值 转为 离散的幅度值 的过程就是 量化 ;

    在这里插入图片描述

    更多相关内容
  • 资源名:matlab程序,可实现脉冲编码调制,模拟信号的数字传输:抽样、量化、编码 资源类型:matlab项目全套源码 源码说明: 全部项目源码都是经过测试校正后百分百成功运行的,如果您下载后不能运行可联系我进行...
  • 信号数据可用于表示任何信息,如符号、文字、语音、图像等,从表现形式上可归结两类:模拟信号和数字信号。  模拟信号与数字信号的区别可根据幅度取什是否离散来确定。模拟信号指幅度的取值是连续的(幅值可由...
  • 调制和解调、带通滤波系统应用抽样信号恢复模拟信号.ppt
  • 一、 模拟数据 编码 数字信号、 二、 音频信号 PCM 编码、 三、 抽象、 四、 量化、 五、 编码、 六、 采样定理、





    一、 模拟数据 编码为 数字信号



    模拟数据 编码为 数字信号 :

    ① 计算机数据形式 : 计算机内部的数据都是 二进制数据 ;

    ② 数字音频 : 计算机内的音频都是 数字音频 ;

    ③ 音频数字化 :模拟信号的音频 , 通过 采样 , 量化转换 为有限个 数字表示的 离散序列 ;





    二、 音频信号 PCM 编码



    模拟数据 编码为 数字信号 , 最典型的应用 , 就是将 模拟的 音频信号 进行 脉码调制 ( PCM ) 编码 , 转为 数字信号 ;

    PCM 音频数据 , 就是 高保真 音频 , 没有经过压缩的原始音频数据 ; 其被存储于 WAV 格式的音频中 ; MP3 , OGG 等格式都是被压缩过的 ;


    PCM 编码过程主要有三个步骤 :

    ① 抽象

    ② 量化

    ③ 编码





    三、 抽象



    抽样 :模拟信号 进行 周期性 扫描 , 将 时间上 连续的信号 , 转为 时间上 离散的信号 ;

    采样定理 : 为了使所有的离散信号 , 能够 不失真地代表 被抽样的模拟数据 , 需要使用 采样定理 :

    f 采 样 频 率 ≥ 2 f 信 号 最 高 频 率 f_{采样频率} \geq 2f_{信号最高频率} f2f

    采样定理 规定了 采样频率 必须 大于等于 信号最高频率的 2 2 2 倍 ;






    四、 量化



    量化 :

    ① 抽样结果 : 抽样取得了 一系列的 电平幅值 集合 ;

    ② 分级标度 : 将上述 电平幅值 按照一定 分级标度 , 转为对应的数值 , 这些数值取整 ;

    ③ 离散数值 :连续的 电平幅值 转为 离散的数值 ;


    分级标度 示例 :

    • 256 种分级 , 对应着每个采样需要使用 l o g 2 256 = 8 log_2 256 = 8 log2256=8 比特来表示 , 对应着 8 8 8 位的音频采样 ;
    • 65536 种分级 , 对应着每个采样需要使用 l o g 2 65536 = 16 log_2 65536 = 16 log265536=16 比特来表示 , 对应着 16 16 16 位的音频采样 ;

    如 : 音频格式是 44100 Hz , 单声道 , 16 位采样 , 就意味着 , 每个采样的取值有 65536 种 ;





    五、 编码



    编码 :量化的结果 , 转为 二进制编码 ;





    六、 采样定理



    采样定理 : 为了使所有的离散信号 , 能够 不失真地代表 被抽样的模拟数据 , 需要使用 采样定理 :

    f 采 样 频 率 ≥ 2 f 信 号 最 高 频 率 f_{采样频率} \geq 2f_{信号最高频率} f2f

    采样定理 规定了 采样频率 必须 大于等于 信号最高频率的 2 2 2 倍 ;


    模拟信号都是 正弦波 构成的 , 每个模拟信号都可以过滤出一个正弦波 ;

    正弦波 S S S , 1 1 1 秒有 2 2 2 个完整的波形 , 即 2 2 2 个波峰 , 2 2 2 个波谷 , 其周期是 0.5 0.5 0.5 秒 , 频率是 2 H z 2Hz 2Hz ;

    正弦波 B B B , 1 1 1 秒有 1 1 1 个完整的波形 , 即 1 1 1 个波峰 , 1 1 1 个波谷 , 其周期是 1.0 1.0 1.0 秒 , 频率是 1 H z 1Hz 1Hz ;


    信号频率就是带宽 , 1 1 1 秒钟震荡的次数 , S S S 信号波形 频率 2 H z 2Hz 2Hz , B B B 信号波形 频率 1 H z 1Hz 1Hz ;


    针对 B B B 波形 , 确定采样频率 :

    • 正弦波的公式 : y = A s i n ( ω t + ϕ ) y = A sin( \omega t + \phi ) y=Asin(ωt+ϕ)

    • 已知参数 : 其中的 ω \omega ω 就是频率 , 1 H z 1Hz 1Hz ;

    • 未知参数 : 还剩下未知数 A A A ϕ \phi ϕ ;

    • 未知参数求值条件 : 只要代入两组数据 , 就可以将该正弦函数的公式求出 , 因此采样时 , 采两组数据 , 就能完全还原该正弦函数 ;


    采样定理 结论 :

    • 在一个 正弦波周期内 , 采样 2 2 2 次 , 就可以还原该正弦波 ;

    • 如果 1 1 1 秒钟有 1 1 1 个完整的正弦波 , 那么采样 2 2 2 次即可 ;

    • 如果 1 1 1 秒钟有 2 2 2 个完整的正弦波 , 那么采样 2 × 2 = 4 2 \times 2 = 4 2×2=4 次即可 ;

    • 因此采样定理中规定 , 采样频率 必须 大于等于 2倍信号最高频率 ;

    在这里插入图片描述


    也可以采更多的样本 , 采样频率越高 , 正弦波形恢复的更准确 , 就越不容易失真 ;

    如 : 音频的采样 44100 H z 44100Hz 44100Hz , 48000 H z 48000 Hz 48000Hz , 96000 H z 96000 Hz 96000Hz , 都非常大 ;


    人耳能听到的声音是 20 H z 20Hz 20Hz ~ 20000 H z 20000Hz 20000Hz , 如果让人耳能够听不出来区别 , 必须在 20000 × 2 = 40000 H z 20000 \times 2 = 40000 Hz 20000×2=40000Hz 以上的采样率才能达到最基本 高保真 要求 ;

    高频失真 , 就是高频的波形没有完整的还原出来 , 采样率不足导致的 ;

    展开全文
  • 基于MATLAB软件,制作demo来模拟信号采样,实现由用户输入被采样正弦信号的幅度、相位、频率信息,产生相应信号并显示出该正弦信号的时域波形图和频域波形图;输入抽样信号的频率,产生抽样脉冲;最后作用于正弦信号...
  • 数字信号是对连续变化的模拟信号进行抽样、量化和编码产生的,称为PCM(Pulse-code modulation),即脉冲编码调制。
  • 什么是模拟信号

    千次阅读 2016-12-02 09:41:48
    模拟信号介绍 信号采集 模拟信号和数字信号的区别模拟信号(介绍) 模拟信号(英语:analog signal)是指在时间和数值上均具有连续性,即对应于任意时间值t均有确定的函数值u或i,并且u或i的幅值是连续取值的信号。与...

    模拟信号(介绍)

    模拟信号(英语:analog signal)是指在时间和数值上均具有连续性,即对应于任意时间值t均有确定的函数值u或i,并且u或i的幅值是连续取值的信号。与模拟信号对应的是数字信号,后者采取分立的逻辑值,而前者可以取得连续值。模拟信号分布于自然界的各个角落,如每天温度的变化,模拟信号的概念常常在涉及电的领域中被使用,不过经典力学、气动力学(pneumatic)、水力学等学科有时也会使用模拟信号的概念。


    信号采集

    介绍下面的信号采集先明确几个重要的概念:
    1. V/F转换

    V/F转换器是把电压信号转变为频率信号的器件,有良好的精度、线性和积分输入特点,他的应用电路简单,外围元件性能要求不高,对环境适应能力强,转换速度不低于一般的双积分型A/D器件,且价格较低,因此在一些非快速A/D转换过程中,V/F转换技术倍受青睐。
    优点:
    1. 接口简单,占用计算机硬件资源少。频率信号可输入微机的任一根I/O线或作为中断源及计数输入等。
    2. 抗干扰性能好。V/F转换本身是一个积分过程,且用V/F转换器实现A/D转换,就是频率计数过程,相当于在计数时间内对频率信号进行积分,因而有转强的抗干扰能力。另外可采用光电耦合器连接V/F转换器与计算机之间的通道,实现光电隔离。
    3. 易于远距离传输。可通过调制进行无线传输或光传输。基于以上这些特点由V/F转换器构造的系统可以简化电路,降低成本,提高性价比。

    1. A/D转换:A/D转化电路。 亦称“模拟数字转换器”,简称“模数转换器”。将模拟量或连续变化的量进行量化(离散化),转换为相应的数字量的电路。 A/D变换包含三个部分:抽样、量化和编码。一般情况下,量化和编码是同时完成的。 抽样是将模拟信号在时间上离散化的过程; 量化是将模拟信号在幅度上离散化的过程; 编码是指将每个量化后的样值用一定的二进制代码来表示。


      • 这里的模拟信号是指电压和电流信号,对模拟信号的处理技术主要包括模拟量的选通、模拟量的放大、信号滤波、电流电压的转换、V/F转换、A/D转换等。
      • 模拟通道—-单片机测控系统有时需要进行多路和多参数的采集和控制,除特殊情况下采用多路独立的放大、A/D外,通常采用公共的采样/保持及A/D转换电路。
      • 信号滤波—-从传感器或其它接收设备获得的电信号,往往会有多种频率成分的噪声信号,需采取滤波措施,滤除干扰噪声,提高系统的信噪比(S/N)。
      • 电流电压—-电压信号可以经由A/D转换器件转换成数字信号然后采集,但是电流不能直接由A/D 转换器转换。在应用中,先将电流转变成电压信号,然后进行转换。
      • 电压频率—-V/F转换器适用于一些非快速而需进行远距离信号传输的A/D转换过程。利用V/F变换,还可以减化电路、降低成本、提高性价比。
      • A/D转换—-A/D转换是指将模拟输入信号转换成N位二进制数字输出信号的过程。
      • 积分ADC—-积分型ADC又称为双斜率或多斜率ADC,它的应用也比较广泛。在此时间间隔内利用计数器对时钟脉冲进行计数,从而实现A/D转换。
      • 压频变换—-压频变换型ADC是间接型ADC,它先输入模拟信号的电压转换成频率与其成正比的脉冲信号,然后在固定的时间间隔内对此脉冲信号进行计数,计数结果即为正比于输入模拟电压信号的数字量。
      • 判断采样—-采样/保持器主要用于稳定信号量,实现平顶抽样。对于高频信号的采集,采样/保持器是非常必要的。
      • 选择量程—-模拟信号的动态范围较大,有时还有可能出现负电压。在选择时,待测信号的动态范围最好在A/D器件的量程范围内。以减少额外的硬件付出。
      • 选择线形—-在A/D采集过程中,线形度越高越好。但是线形度越高,器件的价格也越高。当然,也可以通过软件补偿来减少非线性的影响。
      • 输出接口—-A/D器件接口的种类很多,有并行总线接口的,有SPI、I2C、1-Wire等串行总线接口的。它们在原理和精度上相同,但是控制方法和接口电路会有很大差异。在接口上的选择,主要决定于系统要求、已经开发者对于各种接口的熟练程度。
        这里写图片描述

    模拟信号和数字信号的区别

    模拟信号与数字信号的区别在于:不同的数据必须转换为相应的信号才能进行传输:模拟数据一般采用模拟信号,例如用一系列连续变化的电磁波或电压信号来表示;数字数据则采用数字信号,例如用一系列断续变化的电压脉冲或光脉冲来表示。


    当模拟信号采用连续变化的电磁波来表示时,电磁波本身既是信号载体,同时作为传输介质;而当模拟信号采用连续变化的信号电压来表示时,它一般通过传统的模拟信号传输线路来传输。当数字信号采用断续变化的电压或光脉冲来表示时,一般则需要用双绞线、电缆或光纤介质将通信双方连接起来,才能将信号从一个节点传到另一个节点。
    这里写图片描述

    展开全文
  • 展开的关系从傅里叶展开说起傅里叶变换的输入输出复指数形式的傅里叶级数(推导)复指数形式的傅里叶变换(结论)非周期函数傅里叶变换(推导)非周期函数傅里叶变换(结论)冲激函数及其性质推奈奎斯特抽样定理式子原始...

    本文的主要线路是证明采样定理中的频域变换式。

    其中知识框架如下:
    采样定理流程
    以下将从按照拓扑序进行相关模块的设计

    傅里叶级数与傅里叶展开

    傅里叶级数

    傅里叶级数类似于泰勒级数,是一种特殊的描述函数的方式。

    傅里叶级数的定义为:

    f ( x ) = a 0 + ∑ n = 1 + ∞ ( a n cos ⁡ ( n w x ) + b n sin ⁡ ( n w x ) ) f ( x ) = ∑ n = 0 + ∞ ( a n cos ⁡ ( n w x ) + b n sin ⁡ ( n w x ) ) , b 0 = 0 f(x) = a_0 + \sum_{n=1}^{+\infty}(a_n \cos(nwx) + b_n \sin(nwx)) \\ f(x) = \sum_{n=0}^{+\infty}(a_n \cos(nwx) + b_n \sin(nwx)), b_0=0 f(x)=a0+n=1+(ancos(nwx)+bnsin(nwx))f(x)=n=0+(ancos(nwx)+bnsin(nwx)),b0=0

    它能够描述函数的原理为:傅里叶认为所有周期函数都可以展开成为若干三角函数的求和。

    上面的原理我们可以描述为以下的式子:

    f ( x ) = ∑ i = 1 + ∞ k i sin ⁡ ( w i x + b i ) = ∑ i = 1 + ∞ k i ( sin ⁡ ( w i x ) cos ⁡ ( b i ) + cos ⁡ ( w i x ) sin ⁡ ( b i ) ) = ∑ i = 1 + ∞ k i cos ⁡ ( b i ) sin ⁡ ( w i x ) + k i sin ⁡ ( b i ) cos ⁡ ( w i x ) = ∑ i = 1 + ∞ a i sin ⁡ ( w i x ) + b i cos ⁡ ( w i x ) f(x) = \sum_{i=1}^{+\infty}{k_i \sin(w_ix+b_i)} \\ = \sum_{i=1}^{+\infty}k_i (\sin(w_ix)\cos(b_i)+\cos(w_ix)\sin(b_i)) \\ = \sum_{i=1}^{+\infty} k_i \cos(b_i) \sin(w_ix) + k_i \sin(b_i)\cos(w_ix)\\ = \sum_{i=1}^{+\infty} a_i \sin(w_ix) + b_i \cos(w_ix)\\ f(x)=i=1+kisin(wix+bi)=i=1+ki(sin(wix)cos(bi)+cos(wix)sin(bi))=i=1+kicos(bi)sin(wix)+kisin(bi)cos(wix)=i=1+aisin(wix)+bicos(wix)

    这个是不是长得很像上面的那个式子啊 ^ _ ^,唯一不同的就是上面那个式子中 w i w_i wi取的都是整数, 而不是像下面这样随便取值 ( 随便取值这个叫做傅里叶变换, 后面会有的 )

    傅里叶展开

    傅里叶展开是啥? 傅里叶展开就是把一个函数的傅里叶级数给求出来. 可以类比泰勒展开是求泰勒级数一样, 傅里叶展开的目的就是求傅里叶级数.

    这里需要用到的知识有三角函数的正交性.

    具体来说的话就是对于在集合 { 1 , sin ⁡ x , cos ⁡ x , sin ⁡ 2 x , cos ⁡ 2 x , . . . , sin ⁡ n x , cos ⁡ n x } \{1, \sin x, \cos x, \sin 2x, \cos 2x, ... , \sin nx, \cos nx\} {1,sinx,cosx,sin2x,cos2x,...,sinnx,cosnx}中任取两个函数A(x), B(x). 则满足以下性质:

    ∫ − π π A ( x ) B ( x ) d x = { 0 , A ( x ) ≠ B ( x ) π , A ( x ) = B ( x ) ≠ 1 2 π , A ( x ) = B ( x ) = 1 \int_{-\pi}^{\pi} A(x) B(x) dx = \begin{cases} 0, &A(x) \neq B(x)\\ \pi, &A(x) = B(x) \neq 1\\ 2\pi , &A(x)=B(x) = 1 \end{cases} ππA(x)B(x)dx=0,π,2π,A(x)=B(x)A(x)=B(x)=1A(x)=B(x)=1

    于是就有了求解上面的 a n , b n a_n, b_n an,bn的方法:

    a 0 = 1 2 π ∫ − π π f ( x ) d x a n = 2 2 π ∫ − π π f ( x ) cos ⁡ ( n x ) d x b n = 2 2 π ∫ − π π f ( x ) sin ⁡ ( n x ) d x a_0 = {1 \over 2\pi} \int_{-\pi}^{\pi} f(x) dx\\ a_n = {2 \over 2\pi} \int_{-\pi}^{\pi} f(x) \cos(nx) dx\\ b_n = {2 \over 2\pi} \int_{-\pi}^{\pi} f(x) \sin(nx) dx a0=2π1ππf(x)dxan=2π2ππf(x)cos(nx)dxbn=2π2ππf(x)sin(nx)dx

    这里需要注意一个问题, 就是关于傅里叶级数上面的 a 0 a_0 a0的表示方法, 这个有等价形式, 详情请看: 傅里叶级数

    然后这里就完成了傅里叶级数和傅里叶展开的简单回顾.

    欧拉公式

    先上欧拉公式的公式

    e i x = cos ⁡ ( x ) + i sin ⁡ ( x ) e^{i x} = \cos( x) + i \sin( x) eix=cos(x)+isin(x)

    这个公式的证明是欧拉完成的, 但是可以有一些有趣的理解:

    − x = x ∗ i 2 = x ∗ − 1 -x = x*i^2 = x*-1 x=xi2=x1如果放在数轴上看, 是不是可以把-x看成x逆时针旋转了180°, 那么两个 i i i能让它转180°, 一个 i i i自然可以让它转 90°吧? 所以啊, 这个 i i i就可以看成是一个90°旋转角.

    于是任何一个二维平面上的点都可以表示成 ∣ A ∣ e ρ i |A|e^{\rho i} Aeρi的形式. 同样一个点经过了 x ∗ 90 ° x * 90° x90°的旋转, 是不是可以表示成这个点原先的位置乘了 e x i e^{xi} exi

    除此之外, 这个复指数函数还满足普通指数函数满足的一些性质.

    ρ 1 e a i ⋅ ρ 2 e b i = ρ 1 ρ 2 e ( a + b ) i e − i x = cos ⁡ ( x ) − i sin ⁡ ( x ) \rho_1 e^{ai} \cdot \rho_2 e^{bi} = \rho_1 \rho_2 e^{(a+b)i}\\ e^{-ix} = \cos(x) - i \sin(x) ρ1eaiρ2ebi=ρ1ρ2e(a+b)ieix=cos(x)isin(x)

    剩下的一些其他性质可以看复变函数相关的资料.

    时域频域关系

    音符举例

    先来一个看起来很形象的例子: 一个音符.

    时域中的样子可能是这样的

    在这里插入图片描述
    其中横轴是时间, 纵轴是你接收到的模拟信号的强度

    那么它在频域的信号可能是这样的:

    在这里插入图片描述
    它可能表示的信息就是在八分音符的这半拍里面, 它会有一个跟这个音符音高表示的音一样的频率的正弦波出现, 也就是上面那种正弦波

    所以, 上图的时域上的正弦波可以用下图频域上一个音高来表示.

    相互表示

    在演奏会上, 音乐家手里拿的肯定是长得跟下面的"频域"上的谱子, 而不是很多音符演奏出去之后的交错的时域上的谱子(类似于下面这样?)

    在这里插入图片描述
    为啥? 是因为频域在某些时候能更精确简洁反应时域上的变化

    那既然知道了从频域到时域的意义, 可是我们一般能见到的也只有时域上的谱子啊, 我们怎么搞成频域上面的?

    答案是通过傅里叶变换可以将时域上的东西转化到频域上去.

    又或者我们知道了频域上的信息, 怎么去还原这个时域上的信息?

    答案是通过傅里叶逆变换, 可以将频域还原回时域上面去.

    傅里叶变换

    傅里叶变换和级数, 展开的关系

    傅里叶级数假定了函数是连续的而且是周期函数, 傅里叶变换则将所有函数看成了一个周期无限大的函数, 去做这个傅里叶级数的展开.

    个人感觉本质上傅里叶变换就是更强的傅里叶展开成为傅里叶级数.

    从傅里叶展开说起

    我们先不看傅里叶变换, 普通的傅里叶展开就可以将时域信号转化为频域的信号. 那咋做呢?

    既然我们说任何函数都可以看成是多个三角函数的叠加, 那是不是我们也可以把一个函数分解成多个三角函数的和?

    那必然可以, 我们这里取一个只有奇数次的三角函数相加的例子:

    在这里插入图片描述

    好, 你已经得到了离散的频域的分量.

    咋得到的呢? 这就要看傅里叶变换的输入输出是啥了.

    傅里叶变换的输入输出

    傅里叶变换的输入是一个随时间t变化的函数 f ( t ) f(t) f(t)
    傅里叶变换的输出是一个随频率w变化的函数 F ( w ) F(w) F(w)

    其中这个w啊, 就是三角函数里面那个频率啊!

    如果你把上面那个图的第一条线的频率看作是1, 那么因为后面的频率都是它的整数倍, 所以后面的频率也是一系列的整数.

    对于每个w, 你都能知道一个它的具体形式中 A sin ⁡ ( w x ) A\sin(wx) Asin(wx)里面的 A 和 w A和w Aw, 那既然w是自变量, A理所应当的就成为了因变量. 也就是对应的 F ( w ) = A F(w) = A F(w)=A. (当然这只是一个形式化的理解)

    在这里插入图片描述
    等下哈, 你把这三个图放在一起看, 你看看能发现什么 >_<

    在这里插入图片描述
    是不是非常妙?

    复指数形式的傅里叶级数(推导)

    啥? 为啥我一开始这么明白的东西还能扯上复指数?

    别急, 且听我慢慢道来.

    你看那个优美的式子(w是那个最小的频率, n就是对应频谱的那个下标)

    f ( t ) = ∑ n = 0 + ∞ a n cos ⁡ ( n w t ) + b n sin ⁡ ( n w t ) e i n w t = cos ⁡ ( n w t ) + i sin ⁡ ( n w t ) e − i n w t = cos ⁡ ( n w t ) − i sin ⁡ ( n w t ) f(t) = \sum_{n=0}^{+\infty} a_n \cos(nwt) + b_n \sin(nwt)\\ e^{inwt} = \cos(nwt) + i\sin(nwt)\\ e^{-inwt} = \cos(nwt) - i\sin(nwt) f(t)=n=0+ancos(nwt)+bnsin(nwt)einwt=cos(nwt)+isin(nwt)einwt=cos(nwt)isin(nwt)

    是不是可以做做等价变换啊? 那就先看看三角函数怎么变换吧

    cos ⁡ ( n w t ) = e i n w t + e − i n w t 2 sin ⁡ ( n w t ) = e i n w t − e − i n w t 2 i \cos(nwt) = {e^{inwt} + e^{-inwt} \over 2}\\ \sin(nwt) = {e^{inwt} - e^{-inwt} \over 2i} cos(nwt)=2einwt+einwtsin(nwt)=2ieinwteinwt

    然后我们就把它带到上面那个傅里叶展开式中去

    f ( t ) = ∑ n = 0 + ∞ a n e i n w t + e − i n w t 2 + b n e i n w t − e − i n w t 2 i = ∑ n = 0 + ∞ a n ⋅ e i n w t 2 + a n ⋅ e − i n w t 2 + − i b n ⋅ e i n w t 2 + i b n ⋅ e − i n w t 2 = ∑ n = 0 + ∞ e i n w t a n − i b n 2 + ∑ n = 0 + ∞ e − i n w t a n + i b n 2 f(t) = \sum_{n=0}^{+\infty}a_n {e^{inwt} + e^{-inwt} \over 2} + b_n {e^{inwt} - e^{-inwt} \over 2i}\\ =\sum_{n=0}^{+\infty} {a_n \cdot e^{inwt} \over 2} + {a_n \cdot e^{-inwt} \over 2} + {-ib_n \cdot e^{inwt} \over 2} + {ib_n \cdot e^{-inwt}\over 2}\\ =\sum_{n=0}^{+\infty} e^{inwt} {a_n - ib_n \over 2} + \sum_{n=0}^{+\infty} e^{-inwt} {a_n + ib_n\over 2} f(t)=n=0+an2einwt+einwt+bn2ieinwteinwt=n=0+2aneinwt+2aneinwt+2ibneinwt+2ibneinwt=n=0+einwt2anibn+n=0+einwt2an+ibn

    然后我们把后半部分式子的n反过来, 也就是

    ∑ n = 0 + ∞ e − i n w t a n + i b n 2 = ∑ n = 0 − ∞ e i n w t a − n + i b − n 2 \sum_{n=0}^{+\infty}e^{-inwt}{a_n + ib_n\over 2} = \sum_{n=0}^{-\infty} e^{inwt} {a_{-n} + ib_{-n}\over 2} n=0+einwt2an+ibn=n=0einwt2an+ibn

    然后把它带回去原来的式子, 并且单独把n=0给拿出来

    f ( t ) = ∑ n = 1 + ∞ e i n w t a n − i b n 2 + ∑ n = − 1 − ∞ e i n w t a − n + i b − n 2 + ∑ n = 0 0 e i n w t a 0 f(t) = \sum_{n=1}^{+\infty} e^{inwt} {a_n - ib_n \over 2} + \sum_{n=-1}^{-\infty} e^{inwt} {a_{-n} + ib_{-n}\over 2} + \sum_{n=0}^{0} e^{inwt} a_0 f(t)=n=1+einwt2anibn+n=1einwt2an+ibn+n=00einwta0

    然后你看它是不是都有 e i n w t e^{inwt} einwt?那我们不如?

    f ( t ) = ∑ n = − ∞ + ∞ c n e i n w t c n = { a 0 , n = 0 a n − i b n 2 , n > 0 a − n + i b − n 2 , n < 0 f(t) = \sum_{n=-\infty}^{+\infty} c_n e^{inwt}\\ c_n = \begin{cases} a_0 &,n = 0\\ {a_n-ib_n\over 2} &,n>0\\ {a_{-n}+ib_{-n}\over 2} &,n<0 \end{cases} f(t)=n=+cneinwtcn=a02anibn2an+ibn,n=0,n>0,n<0

    展开说说?

    { a 0 = 1 T ∫ 0 T f ( t ) d t a n = 2 T ∫ 0 T f ( t ) cos ⁡ ( n w t ) d t b n = 2 T ∫ 0 T f ( t ) sin ⁡ ( n w t ) d t { a n = { c n , n = 0 c n + c − n , n > 0 b n = i ( c n − c − n ) , n > 0 \begin{cases} a_0 = {1\over T} \int_{0}^{T} f(t) dt\\ a_n = {2 \over T} \int_{0}^{T} f(t) \cos(nwt) dt\\ b_n = {2 \over T} \int_{0}^{T} f(t) \sin(nwt) dt \end{cases} \\ \begin{cases} a_n = \begin{cases} c_n &,n=0\\ c_n + c_{-n} &,n>0 \end{cases}\\ b_n = i(c_n-c_{-n}) ,n>0 \end{cases} a0=T10Tf(t)dtan=T20Tf(t)cos(nwt)dtbn=T20Tf(t)sin(nwt)dtan={cncn+cn,n=0,n>0bn=i(cncn),n>0

    然后我们把这个ab的求法带到c里面去, 看看c长啥样?

    先推导那个 c n , n = 0 c_n, n=0 cn,n=0

    c 0 = a 0 = 1 T ∫ 0 T f ( t ) d t = 1 T ∫ 0 T f ( t ) e i n w t d t , n = 0 c_0 = a_0 = {1 \over T} \int_{0}^{T}f(t) dt = {1 \over T} \int_{0}^{T}f(t) e^{inwt} dt, n=0 c0=a0=T10Tf(t)dt=T10Tf(t)einwtdt,n=0

    然后推导 c n , n > 0 c_n, n>0 cn,n>0

    c n = a n − i b n 2 = 1 2 ( 2 T ∫ 0 T f ( t ) cos ⁡ ( n w t ) d t − i 2 T ∫ 0 T f ( t ) sin ⁡ ( n w t ) d t ) = 1 T ∫ 0 T f ( t ) ( cos ⁡ ( n w t ) − i sin ⁡ ( n w t ) ) d t = 1 T ∫ 0 T f ( t ) e − i n w t d t \begin{aligned} c_n = {a_n-ib_n\over 2} &= {1\over 2}({{2\over T}\int_{0}^{T}f(t)\cos(nwt)dt -i {2 \over T} \int_{0}^{T}f(t)\sin(nwt)dt}) \\ &= {1\over T}\int_{0}^{T} f(t) (\cos(nwt) - i\sin(nwt)) dt\\ &= {1\over T} \int_{0}^{T} f(t) e^{-inwt} dt \end{aligned} cn=2anibn=21(T20Tf(t)cos(nwt)dtiT20Tf(t)sin(nwt)dt)=T10Tf(t)(cos(nwt)isin(nwt))dt=T10Tf(t)einwtdt

    然后推导 c n , n < 0 c_n, n<0 cn,n<0

    c n = a − n + i b − n 2 = 1 2 ( 2 T ∫ 0 T f ( t ) cos ⁡ ( − n w t ) d t + i 2 T ∫ 0 T f ( t ) sin ⁡ ( − n w t ) d t ) = 1 T ∫ 0 T f ( t ) ( cos ⁡ ( − n w t ) + i sin ⁡ ( − n w t ) ) d t = 1 T ∫ 0 T f ( t ) e i ( − n w t ) d t = 1 T ∫ 0 T f ( t ) e − i n w t d t \begin{aligned} c_n = {a_{-n}+ib_{-n}\over 2} &= {1\over 2}({{2\over T}\int_{0}^{T}f(t)\cos(-nwt)dt +i {2 \over T} \int_{0}^{T}f(t)\sin(-nwt)dt}) \\ &= {1\over T}\int_{0}^{T} f(t) (\cos(-nwt) + i\sin(-nwt)) dt\\ &= {1\over T} \int_{0}^{T} f(t) e^{i(-nwt)} dt\\ &= {1\over T} \int_{0}^{T} f(t) e^{-inwt} dt \end{aligned} cn=2an+ibn=21(T20Tf(t)cos(nwt)dt+iT20Tf(t)sin(nwt)dt)=T10Tf(t)(cos(nwt)+isin(nwt))dt=T10Tf(t)ei(nwt)dt=T10Tf(t)einwtdt

    哎? 是不是长得非常一样啊? 那确实啊! 整理一下

    c n = 1 T ∫ 0 T f ( t ) e − i n w t d t , n ∈ Z c_n = {1\over T} \int_{0}^{T} f(t) e^{-inwt} dt, n\in Z cn=T10Tf(t)einwtdt,nZ

    让我们回归本源( w w w是一个常量, 表示最小的频率尺度):

    f ( t ) = ∑ n = 0 + ∞ a n cos ⁡ ( n w t ) + b n sin ⁡ ( n w t ) = ∑ n = − ∞ + ∞ c n e i n w t = 1 T ∑ n = − ∞ + ∞ ( ∫ 0 T f ( t ) e − i n w t d t ) e i n w t \begin{aligned} f(t) &= \sum_{n=0}^{+\infty} a_n \cos(nwt) + b_n \sin(nwt)\\ &= \sum_{n=-\infty}^{+\infty} c_n e^{inwt}\\ &= {1\over T} \sum_{n=-\infty}^{+\infty}( \int_{0}^{T} f(t) e^{-inwt} dt ) e^{inwt} \end{aligned} f(t)=n=0+ancos(nwt)+bnsin(nwt)=n=+cneinwt=T1n=+(0Tf(t)einwtdt)einwt

    哦对了, 这里补充一下, 为啥会从普通的实数域的 a , b a,b a,b扩展到复数域的这个 c c c.
    我个人理解是: 原先通过两个变量 ( a , b ) (a,b) (a,b)来确定这个函数, 现在可以认为是 ( c r , c i ) (c_r, c_i) (cr,ci)实部和虚部来确定的

    复指数形式的傅里叶变换(结论)

    f ( t ) = ∑ n = − ∞ + ∞ c n e i n w t c n = 1 T ∫ 0 T f ( t ) e − i n w t d t , n ∈ Z f(t) = \sum_{n=-\infty}^{+\infty} c_n e^{inwt}\\ c_n = {1\over T} \int_{0}^{T} f(t) e^{-inwt} dt, n\in Z f(t)=n=+cneinwtcn=T10Tf(t)einwtdt,nZ

    非周期函数傅里叶变换(推导)

    你看上面那个东西是不是都有一个最小的频率w, 所有的w取值只能是离散的吧! 那能不能搞成连续的w呢?

    等等, 似乎不用太麻烦, 如果我 w → 0 w\rightarrow 0 w0足够小, 是不是说所有值都能通过 W = n w W=nw W=nw的形式去表示出来?

    那就开始推式子呗, 开始从普通的离散过渡到连续上面去

    F ( n w ) = c n = 1 T ∫ 0 T f ( t ) e − i n w t d t F(nw) = c_n = {1\over T} \int_{0}^{T} f(t) e^{-inwt} dt F(nw)=cn=T10Tf(t)einwtdt

    让我们取 w → 0 w\rightarrow 0 w0, 带进去看看?

    f ( t ) = lim ⁡ w → 0 ∑ n = − ∞ + ∞ 1 T ( ∫ − T 2 T 2 f ( t ) e − i n w t d t ) e i n w t = lim ⁡ w → 0 ∑ n = − ∞ + ∞ w 2 π ( ∫ − T 2 T 2 f ( t ) e − i n w t d t ) e i n w t = lim ⁡ T → + ∞ ∑ n = − ∞ + ∞ w 2 π ( ∫ − ∞ + ∞ f ( t ) e − i n w t d t ) e i n w t \begin{aligned} f(t) &= \lim_{w \rightarrow 0} \sum_{n=-\infty}^{+\infty} {1\over T} (\int_{-T\over 2}^{T\over 2} f(t)e^{-inwt}dt) e^{inwt}\\ &= \lim_{w \rightarrow 0} \sum_{n=-\infty}^{+\infty} {w\over 2\pi} (\int_{-T\over 2}^{T\over 2} f(t)e^{-inwt}dt) e^{inwt}\\ &= \lim_{T\rightarrow +\infty} \sum_{n=-\infty}^{+\infty} {w\over 2\pi} (\int_{-\infty}^{+\infty}f(t) e^{-inwt}dt)e^{inwt} \end{aligned} f(t)=w0limn=+T1(2T2Tf(t)einwtdt)einwt=w0limn=+2πw(2T2Tf(t)einwtdt)einwt=T+limn=+2πw(+f(t)einwtdt)einwt

    我们把这个求和啊给它换成积分好吧

    f ( t ) = lim ⁡ T → + ∞ ∑ n = − ∞ + ∞ w 2 π ( ∫ − ∞ + ∞ f ( t ) e − i n w t d t ) e i n w t = lim ⁡ T → + ∞ ∫ − ∞ + ∞ w 2 π ( ∫ − ∞ + ∞ f ( t ) e − i n w t d t ) e i n w t d n \begin{aligned} f(t) &= \lim_{T\rightarrow +\infty} \sum_{n=-\infty}^{+\infty} {w\over 2\pi} (\int_{-\infty}^{+\infty}f(t) e^{-inwt}dt)e^{inwt}\\ &= \lim_{T\rightarrow +\infty} \int_{-\infty}^{+\infty} {w\over 2\pi} (\int_{-\infty}^{+\infty}f(t) e^{-inwt}dt)e^{inwt} dn \end{aligned} f(t)=T+limn=+2πw(+f(t)einwtdt)einwt=T+lim+2πw(+f(t)einwtdt)einwtdn

    既然 n w nw nw都连续了, 那要不直接用 W W W来代替 n w nw nw呗, 给这个积分换个积分变量!

    f ( t ) = lim ⁡ T → + ∞ ∫ − ∞ + ∞ w 2 π ( ∫ − ∞ + ∞ f ( t ) e − i n w t d t ) e i n w t d n = ∫ − ∞ + ∞ 1 2 π ( ∫ − ∞ + ∞ f ( t ) e − i W t d t ) e i W t d W = 1 2 π ∫ − ∞ + ∞ ( ∫ − ∞ + ∞ f ( t ) e − i W t d t ) e i W t d W \begin{aligned} f(t) &= \lim_{T\rightarrow +\infty} \int_{-\infty}^{+\infty} {w\over 2\pi} (\int_{-\infty}^{+\infty}f(t) e^{-inwt}dt)e^{inwt} dn\\ &= \int_{-\infty}^{+\infty} {1\over 2\pi} (\int_{-\infty}^{+\infty}f(t) e^{-iWt}dt)e^{iWt} dW\\ &= {1 \over 2\pi} \int_{-\infty}^{+\infty} (\int_{-\infty}^{+\infty}f(t) e^{-iWt}dt)e^{iWt} dW \end{aligned} f(t)=T+lim+2πw(+f(t)einwtdt)einwtdn=+2π1(+f(t)eiWtdt)eiWtdW=2π1+(+f(t)eiWtdt)eiWtdW

    想想咱们傅里叶变换的输出是啥来着?

    输出是对于某个频率的强度 F ( W ) F(W) F(W)啊! 这个是不是就是中间括号里面的东西啊!

    这就是傅里叶变换正变换 / 傅里叶变换

    F ( W ) = ∫ − ∞ + ∞ f ( t ) e − i W t d t F(W) = \int_{-\infty}^{+\infty}f(t) e^{-iWt}dt F(W)=+f(t)eiWtdt

    然后还原回去的过程就是傅里叶变换逆变换 / 逆傅里叶变换

    f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( W ) e i W t d W f(t) = {1\over 2\pi} \int_{-\infty}^{+\infty} F(W) e^{iWt} dW f(t)=2π1+F(W)eiWtdW

    非周期函数傅里叶变换(结论)

    F ( W ) = ∫ − ∞ + ∞ f ( t ) e − i W t d t f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( W ) e i W t d W \begin{aligned} &F(W) = \int_{-\infty}^{+\infty}f(t) e^{-iWt}dt\\ &f(t) = {1\over 2\pi} \int_{-\infty}^{+\infty} F(W) e^{iWt} dW \end{aligned} F(W)=+f(t)eiWtdtf(t)=2π1+F(W)eiWtdW

    上面这个叫傅里叶变换正变换 / 傅里叶变换, 下面这个叫傅里叶变换逆变换 / 逆傅里叶变换.

    冲激函数及其性质

    首先, 冲激函数并不是单位冲激函数, 单位冲激函数 n = 0 n=0 n=0的时候为1, 但是冲激函数 n = 0 n=0 n=0的时候为无穷. 而且通过广义函数相关理论得到冲激函数的一些性质, 而且对于单位冲激函数, 除了第一条不同, 后两条性质相同

    δ ( x ) = { + ∞ , x = 0 0 , x ≠ 0 ∫ − ∞ ∞ δ ( x ) d x = 1 δ ( x ) ∗ h ( x ) = h ( x ) ∫ − ∞ ∞ δ ( x ) h ( x ) d x = h ( 0 ) \delta(x) = \begin{cases} +\infty&, x = 0\\ 0&, x\neq 0 \end{cases}\\ \int_{-\infty}^{\infty} \delta(x) dx= 1\\ \delta(x) * h(x) = h(x)\\ \int_{-\infty}^{\infty} \delta(x) h(x) dx = h(0)\\ δ(x)={+0,x=0,x=0δ(x)dx=1δ(x)h(x)=h(x)δ(x)h(x)dx=h(0)

    这些性质从简单大学数学上是不好解释的.记住就好,可以去查阅广义函数、奇异函数相关内容

    推奈奎斯特抽样定理式子

    原始抽样序列的表示

    x a ( t ) x_a(t) xa(t)为抽样信号, x ′ ( t ) x'(t) x(t)为理想抽样信号, p ( t ) = δ T ( t ) = ∑ m = − ∞ + ∞ δ ( t − m T ) p(t) = \delta_T(t)=\sum_{m=-\infty}^{+\infty}\delta(t-mT) p(t)=δT(t)=m=+δ(tmT)为周期性的单位抽样函数. 注意, 以上函数都是在t上取值连续的函数

    x ′ ( t ) = x a ( t ) ⋅ p ( t ) = x a ( t ) ⋅ ∑ m = − ∞ + ∞ δ ( t − m T ) = ∑ m = − ∞ + ∞ x a ( t ) δ ( t − m T ) x'(t) = x_a(t) \cdot p(t) = x_a(t) \cdot \sum_{m=-\infty}^{+\infty} \delta(t-mT)\\ =\sum_{m=-\infty}^{+\infty} x_a(t) \delta(t-mT) x(t)=xa(t)p(t)=xa(t)m=+δ(tmT)=m=+xa(t)δ(tmT)

    因为 δ ( t − m T ) \delta(t-mT) δ(tmT)当且仅当 t = m T t=mT t=mT的时候有值, 因此可以认为上式为:

    x ′ ( t ) = ∑ m = − ∞ + ∞ x a ( t ) δ ( t − m T ) = ∑ m = − ∞ + ∞ x a ( m T ) δ ( t − m T ) x'(t) = \sum_{m=-\infty}^{+\infty} x_a(t) \delta(t-mT) = \sum_{m=-\infty}^{+\infty} x_a(mT) \delta(t-mT) x(t)=m=+xa(t)δ(tmT)=m=+xa(mT)δ(tmT)

    求抽样完成后序列频谱与原频谱关系

    X a ( i Ω ) = F x a ( Ω ) X_a(i\Omega) = F_{x_a}(\Omega) Xa(iΩ)=Fxa(Ω), Δ T ( i Ω ) = F δ T ( Ω ) \Delta_T(i\Omega) = F_{\delta_T}(\Omega) ΔT(iΩ)=FδT(Ω), X a ′ ( i Ω ) = F x a ′ ( Ω ) X'_a(i\Omega) = F_{x'_a}(\Omega) Xa(iΩ)=Fxa(Ω)也就是对应的频谱的输出结果

    x a ′ ( t ) = δ T ( t ) ⋅ x a ( t ) X a ′ ( i Ω ) = 1 2 π [ Δ T ( i Ω ) ∗ X a ( i Ω ) ] x'_a(t) = \delta_T(t) \cdot x_a(t)\\ X'_a(i\Omega) = {1\over 2\pi} [\Delta_T(i\Omega) * X_a(i\Omega)] xa(t)=δT(t)xa(t)Xa(iΩ)=2π1[ΔT(iΩ)Xa(iΩ)]

    (第二个式子还没学会怎么证明, 等学会了补上>_<)

    那么我们就从 Δ T ( t ) \Delta_T(t) ΔT(t)开始求解吧

    求解周期冲激函数傅里叶展开

    首先明确一个事情, 就是我们这个冲激函数求出来之后会近似于原始的冲激函数.(不信你做做逆变换). 为啥要求这个傅里叶展开, 这是为了解决这个函数不好计算的问题

    这个周期冲激函数是个周期函数, 所以直接用周期函数的傅里叶展开就行了 Ω s \Omega_s Ωs为抽样频率

    δ T ( t ) = ∑ k = − ∞ + ∞ c k e i k Ω s t d t \delta_T(t) = \sum_{k=-\infty}^{+\infty} c_k e^{ik\Omega_st}dt δT(t)=k=+ckeikΩstdt

    然后周期函数的 c k c_k ck求法带进去

    c k = 1 T ∫ − T 2 T 2 δ T ( t ) e − i k Ω s t d t ( 展 开 δ T ) = 1 T ∫ − T 2 T 2 ∑ p = − ∞ + ∞ δ ( t − p T ) e − i k Ω s t d t \begin{aligned} c_k &= {1\over T} \int_{-T\over 2}^{T\over 2} \delta_T(t) e^{-ik\Omega_st}dt\\ (展开\delta_T)&= {1\over T} \int_{-T\over 2}^{T\over 2}\sum_{p=-\infty}^{+\infty}\delta(t-pT) e^{-ik\Omega_st}dt \end{aligned} ck(δT)=T12T2TδT(t)eikΩstdt=T12T2Tp=+δ(tpT)eikΩstdt

    发现只有 p = 0 p=0 p=0的时候, t = 0 t=0 t=0在最靠近0的这个周期里面. 所以上面这个式子可以写成(利用一下冲激函数的性质)

    c k = 1 T ∫ − T 2 T 2 ∑ p = − ∞ + ∞ δ ( t − p T ) e − i k Ω s t d t = 1 T e − i k Ω s 0 = 1 T \begin{aligned} c_k &= {1\over T} \int_{-T\over 2}^{T\over 2}\sum_{p=-\infty}^{+\infty}\delta(t-pT) e^{-ik\Omega_st}dt\\ &= {1\over T} e^{-ik\Omega_s0} = {1\over T} \end{aligned} ck=T12T2Tp=+δ(tpT)eikΩstdt=T1eikΩs0=T1

    所以啊, 原来那个周期性冲激函数就可以被展开成

    δ T ( t ) = 1 T ∑ k = − ∞ + ∞ e i k Ω s t \delta_T(t) = {1\over T} \sum_{k=-\infty}^{+\infty} e^{ik\Omega_st} δT(t)=T1k=+eikΩst

    然后这个展开式的图呢, 在最靠近0的这个区间里啊, 大概长这样(我只模拟了80层叠加, 而且是展开之后的叠加)

    在这里插入图片描述
    你看它这个趋势就是在向无穷去叠加, 也就是去模拟这个周期冲激函数.

    然后我们去考虑一下把这个东西给做一下傅里叶变换.

    求解周期冲激函数的频域

    嗯, 列式子吧(看起来很简单!), F T [ ⋅ ] FT[\cdot] FT[] 表示傅里叶变换

    Δ T ( i Ω ) = F T [ 1 T ∑ k = − ∞ + ∞ e i k Ω s t ] = 1 T ∑ k = − ∞ + ∞ F T [ e i k Ω s t ] \begin{aligned} \Delta_T(i\Omega) &= FT[{1\over T} \sum_{k=-\infty}^{+\infty} e^{ik\Omega_st}]\\ &= {1\over T} \sum_{k=-\infty}^{+\infty} FT[e^{ik\Omega_st}] \end{aligned} ΔT(iΩ)=FT[T1k=+eikΩst]=T1k=+FT[eikΩst]

    然后我们看这个 F T [ e i k Ω s t ] FT[e^{ik\Omega_st}] FT[eikΩst], 它是不是很简单做出来

    它有一条性质, 也就是结果, 是 F T [ e i k Ω s t ] = 2 π δ ( Ω − k Ω s ) FT[e^{ik\Omega_st}] = 2\pi\delta(\Omega-k\Omega_s) FT[eikΩst]=2πδ(ΩkΩs)

    然后带进去就能得到:

    Δ T ( i Ω ) = 1 T ∑ k = − ∞ + ∞ F T [ e i k Ω s t ] = 2 π T ∑ k = − ∞ + ∞ δ ( Ω − k Ω s ) = Ω s ∑ k = − ∞ + ∞ δ ( Ω − k Ω s ) \begin{aligned} \Delta_T(i\Omega) &= {1\over T} \sum_{k=-\infty}^{+\infty} FT[e^{ik\Omega_st}]\\ &= {2\pi\over T} \sum_{k=-\infty}^{+\infty} \delta(\Omega - k\Omega_s)\\ &= {\Omega_s} \sum_{k=-\infty}^{+\infty} \delta(\Omega - k\Omega_s)\\ \end{aligned} ΔT(iΩ)=T1k=+FT[eikΩst]=T2πk=+δ(ΩkΩs)=Ωsk=+δ(ΩkΩs)

    什么? 你非要推出上面那个性质来? 那就单独开一节吧!

    求解一个特殊的傅里叶变换(建议跳过)

    求解下面这个式子

    F ( W ) = ∫ − ∞ + ∞ e i k Ω s t ⋅ e − i W t d t F(W) = \int_{-\infty}^{+\infty} e^{ik\Omega_st} \cdot e^{-iWt} dt F(W)=+eikΩsteiWtdt

    然后求解一下这个问题哈:

    F ( W ) = { 0 , W ≠ k Ω s + ∞ , W = k Ω s F(W) = \begin{cases} 0&, W\neq k\Omega_s\\ +\infty &,W = k\Omega_s \end{cases} F(W)={0+,W=kΩs,W=kΩs

    哎你看, 这东西是不是冲激函数啊?那我们给它搞成冲激函数一样的东西呗

    F ( W ) = δ ( W − k Ω s ) F(W) = \delta(W-k\Omega_s) F(W)=δ(WkΩs)

    但是这…真的对吗?

    让我们利用一下逆傅里叶变换, 看看能得到些啥?

    f ( t ) = 1 2 π ∫ − ∞ + ∞ F ( W ) e i W t d W = 1 2 π ∫ − ∞ + ∞ δ ( W − k Ω s ) e i W t d W = 1 2 π e i k Ω s t \begin{aligned} f(t) &= {1\over 2\pi} \int_{-\infty}^{+\infty} F(W) e^{iWt} dW\\ &= {1\over 2\pi} \int_{-\infty}^{+\infty} \delta(W-k\Omega_s) e^{iWt} dW\\ &= {1\over 2\pi} e^{ik\Omega_st} \end{aligned} f(t)=2π1+F(W)eiWtdW=2π1+δ(WkΩs)eiWtdW=2π1eikΩst

    哎?是不是少了个系数? 少了个系数 2 π 2\pi 2π哎, 那就搞进去呗

    F ( W ) = 2 π δ ( W − k Ω s ) F(W) = 2\pi\delta(W-k\Omega_s) F(W)=2πδ(WkΩs)

    然后就得到了那个显然成立的式子.

    理想抽样信号频谱

    X a ′ ( i Ω ) = 1 2 π [ ( 2 π T ∑ k = − ∞ + ∞ δ ( Ω − k Ω s ) ) ∗ X a ( i Ω ) ] = 1 T [ ∑ k = − ∞ + ∞ δ ( Ω − k Ω s ) ) ∗ X a ( i Ω ) ] = 1 T ∑ k = − ∞ + ∞ [ δ ( Ω − k Ω s ) ∗ X a ( i Ω ) ] = 1 T ∑ k = − ∞ + ∞ [ ∫ − ∞ + ∞ X a ( i θ ) ⋅ δ ( Ω − θ − k Ω s ) d θ ] = 1 T ∑ k = − ∞ + ∞ X a ( i ( Ω − k Ω s ) ) = 1 T ∑ k = − ∞ + ∞ X a ( i ( Ω − k 2 π T ) ) \begin{aligned} X'_a(i\Omega) &= {1\over 2\pi} [({2\pi\over T} \sum_{k=-\infty}^{+\infty} \delta(\Omega - k\Omega_s)) * X_a(i\Omega)]\\ &= {1\over T} [\sum_{k=-\infty}^{+\infty} \delta(\Omega - k\Omega_s)) * X_a(i\Omega)]\\ &= {1\over T} \sum_{k=-\infty}^{+\infty} [\delta(\Omega - k\Omega_s) * X_a(i\Omega)]\\ &= {1\over T} \sum_{k=-\infty}^{+\infty} [\int_{-\infty}^{+\infty} X_a(i\theta) \cdot \delta(\Omega-\theta - k\Omega_s)d\theta ]\\ &= {1\over T} \sum_{k=-\infty}^{+\infty} X_a(i(\Omega - k\Omega_s))\\ &= {1\over T} \sum_{k=-\infty}^{+\infty} X_a(i(\Omega - k{2\pi\over T})) \end{aligned} Xa(iΩ)=2π1[(T2πk=+δ(ΩkΩs))Xa(iΩ)]=T1[k=+δ(ΩkΩs))Xa(iΩ)]=T1k=+[δ(ΩkΩs)Xa(iΩ)]=T1k=+[+Xa(iθ)δ(ΩθkΩs)dθ]=T1k=+Xa(i(ΩkΩs))=T1k=+Xa(i(ΩkT2π))

    这不就是把原先的那个 X a X_a Xa给它左右平移 2 π T 2\pi \over T T2π嘛, 那就是说的延拓周期

    Ω s = 2 π T = 2 π f s \Omega_s = {2\pi \over T} = 2\pi f_s Ωs=T2π=2πfs

    就是抽样频率啊!

    X a ′ ( i Ω ) = 1 T ∑ k = − ∞ + ∞ X a ( i ( Ω − k 2 π T ) ) = Ω s 2 π ∑ k = − ∞ + ∞ X a ( i ( Ω − k Ω s ) ) \begin{aligned} X'_a(i\Omega) &= {1\over T} \sum_{k=-\infty}^{+\infty} X_a(i(\Omega - k{2\pi\over T}))\\ &= {\Omega_s \over 2\pi} \sum_{k=-\infty}^{+\infty} X_a(i(\Omega - k\Omega_s)) \end{aligned} Xa(iΩ)=T1k=+Xa(i(ΩkT2π))=2πΩsk=+Xa(i(ΩkΩs))

    展开全文
  • 抽样信号的matlab实现

    2011-06-03 20:50:30
    matlab对信号与系统的实现有很好的帮助
  • 首先我们要知道什么是抽样信号抽样信号对于信号的分析和处理有什么作用?我们能用它解决什么样的问题?这些是我们要学习的重点。 什么是抽样? 抽样就是指利用所谓的抽样脉冲序列从连续信号中“抽取”一系列的...
  • 本设计方案旨在提出一套基于TMS320C6x11系列DSP的图像获取方案,利用模拟视频信号的统一性,实现随意更换带有标准模拟视频信号输出接口的图像设备而无需在图像处理系统的硬件和软件上作修改。
  • 通信原理各种模拟信号的调制仿真,包括AM、FM、2DPSK等调制的仿真。
  • 有关对于信号的抽样和恢复问题 最近在学习数字信号处理这门课程,真的是头都要炸了,第...关于模拟信号抽样抽样之后的信号的频谱是原模拟信号的频谱的以2pifs 周期的周期延拓,幅值是以fs加权之后的幅值 通过...
  • 信号抽样实验报告

    2021-04-20 07:33:57
    大连理工大学实验报告学院(系): 专业: 班级:姓 名: 学号: 组: ___ 实验时间: 实验室: 实验台:指导教师签字: 成绩:实验三 信号抽样一、实验目的1 学会运用MATLAB 完成信号抽样及对抽样信号的频谱进行分析...
  • 里面有低通抽样的代码,有均匀量化的代码,还有非均匀量化的代码,写了简单的例子调用这些函数方便大家查看效果。
  • 第七章 模拟信号数字化理论 7.7 PCM和DM性能比较 7.7-1 不同编码位数n时PCM与DM的性能比较   7.8 时分复用(TDM)
  • matlab 信号抽样与恢复

    2021-04-21 11:50:27
    2.信号重建经采样后得到信号)(t f s 经理想低通)(t h 则可得到重建信号)(t f ,即: )(t f =)(t f s *)(t h其中:)(t f s =)(t f ∑∞∞--)(s nT t δ=∑∞∞--)()(ss nT t nT f δ )()(t Sa T t h c c s ωπω= ...
  • 通信原理各种模拟信号的调制仿真,包括AM、FM、2DPSK等调制的仿真。(Modulation simulation of various analog signals in communication principle, Including AM, FM, 2DPSK modulation simulation.)
  • 信号与系统】Multisim 仿真抽样定理与信号恢复

    千次阅读 多人点赞 2020-07-22 16:46:13
    Multisim 仿真抽样定理与信号恢复,观察离散信号频谱,了解其频谱特点, 验证抽样定理并恢复原信号。要恢复原信号,低通滤波器的截止频率满足的条件。
  • 最近在准备考研复试,查询了很久连续时间信号、离散时间信号、模拟信号和数字信号的区别,发现网上各种说法不一,在这里以熊庆旭老师编写的《信号与系统》和樊昌信老师编写的《通信原理》依据来给大家捋一捋这四种...
  • \cfrac{\pi }{4} ω=4π​ , Ω = 2 π f 0 \Omega = 2\pi f_0 Ω=2πf0​ , 其中 f 0 f_0 f0​ 是模拟频率 , 没有单位 , f 0 = T T 0 f_0 = \cfrac{T}{T_0} f0​=T0​T​ , 其中 T 0 T_0 T0​ 是模拟信号 ...
  • 为了让MATLAB数字信号处理的相关博文能够得到一个梳理,我开通了一个专栏:数字信号处理的MATLAB实现模拟信号经过采样后得到x(n),从x(n)中重建模拟信号在数学上可用公式来描述: 式中,是一种内插函数。 对于样本...
  • 离散时间信号:模拟信号经过采样得到的时间上离散、幅值上连续的信号,简称离散信号,也称为抽样信号或序列。 量化:采用一组数码(如二进制数码)来逼近离散信号的幅值,将其转换数字信号。 数字信号:离散信号...
  • 一、DSP 知识领域、 二、抽样定理、 三、多抽样
  • 什么是模拟通信呢?比如在电话通信中,用户线上传送的电信号是随着用户声音大小的变化而变化的。这个变化的电信号无论在时间上或是在幅度上都是连续的,这种信号称为模拟信号。在用户线上传输模拟信号的通信方式称为...
  • 模拟信号采样过程

    2022-01-02 14:48:51
  • 如题所示几种信号均计算机网络中出现的信号,且易搞混,在这里进行简单梳理。 模拟信号:形状连续,可以设想形如sin...宽带信号(对应模拟信号):将基带信号调制后形成模拟信号,并在模拟信道上传输,而这种传输...
  • 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 f =)(t f s *)(t h ...
  • 摘要:介绍了TI公司的信号处理器TMS320VC5402串行口的主要特点,及其和音频模拟接口芯片TLC320AD50C的结构及其使用注意事项。详细讲述了两者的硬件连接及软件实现。  关键词:数字信号处理器 TMS320CV5402 音频模拟...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 9,167
精华内容 3,666
关键字:

抽样信号为什么是模拟信号

友情链接: f45.rar