notes
  • 概述
  • 🐶MATLAB
    • 🏋️1. MATLAB介绍
    • 🤼2. MATLAB基础
    • 🤸3. 分支语句和编程设计
    • ⛹️4. 循环结构
    • 🤺5. 自定义函数
    • 🤾6. 复数数据、字符数据和附加画图类
    • 🏌️7. 稀疏矩阵、单元阵列和结构
    • 🏇8. 输入/输入函数
    • 🧘9. 图形句柄
  • 🐱通信原理
    • 🍔1. 绪论
    • 🥯2. 确知信号
    • 🌮3. 随机过程
    • 🧆4. 信道
    • 🫔5. 模拟调制系统
  • 🐯信息论与编码
    • 📓1. 熵和互信息
    • 📔2. 信道容量
  • 🐷数字信号处理
    • ⚽1. 离散时间信号与系统
  • 🐰电磁场与电磁波
    • 🚗1. 矢量分析
  • 🐼信息安全技术
    • 🖊️1. 信息安全概述
    • ✒️2. 信息安全与密码学
Powered by GitBook
On this page
  • 1.1 离散时间信号——序列
  • 1.1.1 序列
  • 1.1.2 序列的运算
  • 1.1.3 序列的卷积和
  • 1.1.4 序列的相关性
  • 1.1.5. 几种常用的典型序列
  • 1.1.6 序列的周期性
  • 1.1.7 用单位抽样序列表示任意序列
  • 1.2 线性移不变系统
  • 1.2.1 离散时间线性系统
  • 1.2.2 离散时间移不变系统
  • 1.2.3 离散时间线性移不变系统(LSI系统)
  • 1.2.4 因果系统
  • 1.2.5 稳定系统
  • 1.3 常系数线性差分方程
  • 1.4 连续时间信号的采样
  • 1.4.1 模拟信号的连续采样
  • 1.4.2 时域采用定理
  • 1.4.3 模拟信号的实际采样
  • 1.4.4 带通信号的抽样
  • 1.4.5 连续时间信号 、理想抽样信号以及抽样序列
  • 1.4.6 时域信号的插值重构
  • 1.4.7 正弦信号的抽样
  1. 数字信号处理

1. 离散时间信号与系统

Previous数字信号处理Next电磁场与电磁波

Last updated 6 months ago

1.1 离散时间信号——序列

1.1.1 序列

离散时间信号是时间上不连续的序列,用x(n)x(n)x(n)来表示。x(n)x(n)x(n)可以看成对模拟信号xa(t)x_a(t)xa​(t)的等间隔时间抽样,即

x(n)=xa(t)∣t=nT=xa(nT)x(n)=x_a(t)\Big|_{t=nT}=x_a(nT)x(n)=xa​(t)​t=nT​=xa​(nT)

序列的三种表示方法:

  1. 函数表示法,例如x(n)=anu(n)x(n)=a^nu(n)x(n)=anu(n)

  2. 数列表示法,例如x(n)=...,−5,−3,−1‾,0,2,7,9x(n)={...,-5,-3,\underline{-1},0,2,7,9}x(n)=...,−5,−3,−1​,0,2,7,9,n=0n=0n=0时,x(0)x(0)x(0)的值要用下划线标注

  3. 图形表示法,如图

1.1.2 序列的运算

3个基本运算单元:加法器、乘法器和延时单位。

序列的三类运算:

  1. 对幅度进行运算

  2. 基于变量的运算

    1. 时间尺度变换

      1. 插值(下抽样变换),为了增加抽样频率,例如插入零值,

  3. 既对幅度运算又对变量运算

    1. 差分运算

    2. 卷积和运算

    3. 相关运算

1.1.3 序列的卷积和

  1. 卷积和的运算步骤:

  2. 常用的卷积和求解方法:

    1. 图解加上卷积方法,一个结论:

    2. 列表方法,适用于两个有限长序列的卷积和

1.1.4 序列的相关性

  1. 互相关函数序列

    1. 定义

      互相关函数没有翻褶这一步骤

    2. 性质

    3. 用卷积运算表示相关运算,即

  2. 自相关函数序列

    1. 定义

    2. 性质

  3. 对于功率信号,互相关函数及自相关函数定义为

    若为周期信号,可写成

1.1.5. 几种常用的典型序列

  1. 单位抽样(单位冲激,单位脉冲)序列

  2. 单位阶跃序列

  3. 矩形序列

  4. 实指数序列

  5. 复指数序列

  6. 正弦型序列

注意:

1.1.6 序列的周期性

周期性序列:

1.1.7 用单位抽样序列表示任意序列

1.2 线性移不变系统

一个离散时间系统是将输入序列按照所需要的目的变换成输出序列的一种运算,若以$T[\cdot]$表示这种运算,则有

如图

下面讨论线性移不变系统

1.2.1 离散时间线性系统

  1. 满足叠加原理

  2. 证明一个系统是线性系统需对所有常系数(包括复数)以及所有输入(包括复数)满足叠加性(可加性和比例性)的两个条件。

  3. 在线性系统中,在全部时间上,零输入一定产生零输出。

1.2.2 离散时间移不变系统

  1. 若系统的输人输出关系不随时间而变化,则称它为移(时)不变系统,即若

    则

  2. 移不变系统的输出序列随输入序列的移位而作相同的移位,且保持输出序列的形状是不变的

  3. 线性和移不变是系统的两个独立特性

1.2.3 离散时间线性移不变系统(LSI系统)

  1. 单位冲激响应

  2. LSI系统输出序列与输入序列在时域中的关系——卷积和关系

    即LSI系统的输出序列是输入序列与系统单位抽样响应的卷积和。

  3. LSI系统卷积和运算的性质

    1. 交换律

    2. 结合律

    3. 分配律

1.2.4 因果系统

因果系统是指系统的输出不发生在输入之前的系统,即

1.2.5 稳定系统

若系统满足稳定输入产生稳定输出,就是因果系统,即

LSI系统是稳定系统的充要条件是单位抽样响应绝对可和,即

1.3 常系数线性差分方程

离散时间系统的输入和输出关系可用常系数线性差分方程表示,即

求解线性常系数差分方程的三种方法:

  1. 经典解法

  2. 迭代法

  3. 卷积和计算法

常系数差分方程表示的系统,只有当边界条件是使系统松弛的,则该系统是线性,移不变系统,所谓松弛状态是指

则边界条件必须满足:

1.4 连续时间信号的采样

1.4.1 模拟信号的连续采样

  1. 理想抽样信号

    因而有

  2. 理想抽样信号的频谱

1.4.2 时域采用定理

1.4.3 模拟信号的实际采样

1.4.4 带通信号的抽样

带通信号的频谱存在某一频段范围,而不是零频范围,如图

后面的情况同上述分析。

1.4.6 时域信号的插值重构

就可以得到原信号频谱,即

所以输出端即为原模拟信号。

信号重构的抽样内插公式:

1.4.7 正弦信号的抽样

处理周期性正弦序列时,应该注意:

  1. 对抽样后的离散周期性的正弦序列作截断时,其截断长度必须为序列周期的整数倍,才不会产生频域的泄露。

  2. 离散正弦序列不宜补零后作频谱分析,否则会产生频域泄露。

  1. 同一个模拟正弦型信号,如果用两个不同的抽样频率抽样后,所得到的序列仍可能是相同的,我们没法确定其原抽样频率。

加法:z(n)=x(n)+y(n)z(n)=x(n)+y(n)z(n)=x(n)+y(n)

乘法:w(n)=x(n)y(n)w(n)=x(n)y(n)w(n)=x(n)y(n),当x(n)x(n)x(n)或y(n)y(n)y(n)是常数ccc时,称为标度运算,即w(n)=cx(n)w(n)=cx(n)w(n)=cx(n)

累加:y(n)=∑k=−∞nx(k)y(n)=\sum_{k=-\infty}^{n}x(k)y(n)=∑k=−∞n​x(k)

序列的绝对可和:S=∑n=−∞∞∣x(n)∣S=\sum_{n=-\infty}^{\infty}\mid x(n)\mid S=∑n=−∞∞​∣x(n)∣,当S=B<∞S=B<\inftyS=B<∞时,称序列x(n)x(n)x(n)为绝对可和序列

序列的能量:E[x(n)]=∑n=−∞∞∣x(n)∣2E[x(n)]=\sum_{n=-\infty}^\infty\mid x(n)\mid^2E[x(n)]=∑n=−∞∞​∣x(n)∣2,若E[x(n)]=A<∞E[x(n)]=A<\inftyE[x(n)]=A<∞,则称x(n)x(n)x(n)为能量信号有限信号,简称能量信号

序列的平均功率:P[x(n)]=lim⁡N→∞12N+1∑n=−NN∣x(n)∣2P[x(n)]=\lim_{N\to\infty}\frac{1}{2N+1}\sum_{n=-N}^{N}\mid x(n)\mid^2P[x(n)]=limN→∞​2N+11​∑n=−NN​∣x(n)∣2,若极限存在,即P[x(n)]=C<∞P[x(n)]=C<\inftyP[x(n)]=C<∞,则称x(n)x(n)x(n)为功率有限信号,简称功率信号,对于周期信号,只需取一个周期NNN的平均功率即可,即P[x(n)]=1N∑n=0N−1∣x(n)∣2P[x(n)]=\frac{1}{N}\sum_{n=0}^{N-1}\mid x(n)\mid^2P[x(n)]=N1​∑n=0N−1​∣x(n)∣2

移位,序列x(n)x(n)x(n)的移位序列是x(n−m)x(n-m)x(n−m),mmm为正数时,表示向右移位(延时),反之向左移位(超前)

翻褶,序列x(n)x(n)x(n)的翻褶序列为x(−n)x(-n)x(−n),x(n−m)x(n-m)x(n−m)的翻褶序列为x(−n−m)x(-n-m)x(−n−m)

抽取(上抽样变换),为了减少抽样频率,xd(n)=x(Dn),Dx_d(n)=x(Dn),Dxd​(n)=x(Dn),D为整数

xI′(n)={x(n/I),n=mI,I 为整数,m=0,±1,±2,...0,其他 nx_I^{'}(n)=\begin{cases}x(n/I) ,&n=mI ,I\text{ 为整数},m=0,\pm1,\pm2,...\\\\0 ,&\text{其他 }n\end{cases}xI′​(n)=⎩⎨⎧​x(n/I),0,​n=mI,I 为整数,m=0,±1,±2,...其他 n​

前向差分:Δx(n)=x(n+1)−x(n)\Delta x(n)=x(n+1)-x(n)Δx(n)=x(n+1)−x(n)

后向差分:∇x(n)=x(n)−x(n−1)\nabla x(n)=x(n)-x(n-1)∇x(n)=x(n)−x(n−1)

其中:∇x(n)=Δx(n−1)\nabla x(n)=\Delta x(n-1)∇x(n)=Δx(n−1)

y(n)=x(n)∗h(n)=∑m=−∞∞x(m)h(n−m)=∑m=−∞∞x(n−m)h(m)y(n)=x(n)*h(n)=\sum_{m=-\infty}^{\infty}x(m)h(n-m)=\sum_{m=-\infty}^{\infty}x(n-m)h(m)y(n)=x(n)∗h(n)=m=−∞∑∞​x(m)h(n−m)=m=−∞∑∞​x(n−m)h(m)
rxy(m)=∑n=−∞∞x(n)y(n−m)r_{xy}(m)=\sum_{n=-\infty}^{\infty}x(n)y(n-m)rxy​(m)=n=−∞∑∞​x(n)y(n−m)
y(n)=x(n)∗h(n)=∑m=−∞∞x(m)h(n−m)y\left(n\right)=x\left(n\right)*h\left(n\right)=\sum_{m=-\infty}^{\infty}x\left(m\right)h\left(n-m\right)y(n)=x(n)∗h(n)=m=−∞∑∞​x(m)h(n−m)

翻褶:选哑变量为 mmm,作 x(m),h(m)x(m),h(m)x(m),h(m),将 h(m)h(m)h(m)以 m=0m=0m=0 的垂直轴为对称轴翻褶成h(−m)h\left(-m\right)h(−m)。

移位:将h(−m)h(-m)h(−m)移位nnn,得h(n−m),n>0h(n-m),n>0h(n−m),n>0 时,右移nnn 位,n<0,n<0,n<0 时,左移∣n∣|n|∣n∣位。

相乘:将h(n−m)h(n-m)h(n−m)与x(m)x(m)x(m)在相同mmm 处的对应值相乘。

相加:将以上所有mmm处的乘积值叠加,就得到这一个nnn值下的 y(n)y(n)y(n)值。依上法取 n=⋅⋅⋅,−2,−1,0,1,2,⋅⋅⋅n=\cdotp\cdotp\cdotp,-2,-1,0,1,2,\cdotp\cdotp\cdotpn=⋅⋅⋅,−2,−1,0,1,2,⋅⋅⋅各值,即可得到全部y(n)y(n)y(n)值。

若x(n)x(n)x(n)在N3⩽n⩽N4N_3\leqslant n\leqslant N_4N3​⩽n⩽N4​范围有非零值,h(n)h(n)h(n)在N1⩽n⩽N2N_1\leqslant n\leqslant N_2N1​⩽n⩽N2​范围有非零值。 y(n)=x(n)∗h(n)y(n)=x(n)*h(n)y(n)=x(n)∗h(n)存在的范围为:N1+N3⩽n⩽N2+N4N_{1}+N_{3}\leqslant n\leqslant N_{2}+N_{4}N1​+N3​⩽n⩽N2​+N4​

已知:x(n)=1,2,4‾,3x(n)={1,2,\underline{4},3}x(n)=1,2,4​,3,h(n)=2‾,3,5h(n)={\underline{2},3,5}h(n)=2​,3,5,取 y(n)=∑m=02h(m)x(n−m)y(n)=\sum_{m=0}^2h(m)x(n-m)y(n)=∑m=02​h(m)x(n−m),则mmm应为0⩽m⩽20\leqslant m\leqslant 20⩽m⩽2,列表如下

故y(n)=2,7,19‾,28,29,15y(n)={2,7,\underline{19},28,29,15}y(n)=2,7,19​,28,29,15

对位相乘相加法:将序列排成两行,且将其各自nnn最大的序列值对其(即按右端对其),然后作乘法运算,但是不要进位,最后将同一列的乘积值相加即得到卷积和结果,如下图

卷积和的序列长度:x(n)x(n)x(n)为NNN点长序列,h(n)h(n)h(n)为MMM点长序列,则$$y(n)=x(n)*h(n)$$为L=N+M−1L=N+M-1L=N+M−1点长列。

rxy(m)=∑n=−∞∞x(n)y(n−m)r_{xy}\left(m\right)=\sum_{n=-\infty}^{\infty}x\left(n\right)y\left(n-m\right)rxy​(m)=n=−∞∑∞​x(n)y(n−m)

rxy(m)r_{xy}(m)rxy​(m)与ryx(m)r_{yx}(m)ryx​(m)互为偶对称关系,即rxy(m)=ryx(−m)r_{xy}(m)=r_{yx}(-m)rxy​(m)=ryx​(−m)

rxy(m)r_{xy(m)}rxy(m)​不是偶对称函数,即rxy(−m)≠rxy(m)r_{xy}(-m)\neq r_{xy}(m)rxy​(−m)=rxy​(m)

当x(n)x(n)x(n),y(n)y(n)y(n)是绝对可和信号时,则有lim⁡m→∞rxy(m)=0\lim_{m\to\infty}r_{xy}\left(m\right)=0limm→∞​rxy​(m)=0

rxyr_{xy}rxy​中mmm的有值范围

设有限长序列x(n),y(n)x(n),y(n)x(n),y(n)有值的范围分别为

x(n):Nx1⩽n⩽Nx2;y(n):Ny1⩽n⩽Ny2x(n):N_{x1}\leqslant n\leqslant N_{x2} ; y(n):N_{y1}\leqslant n\leqslant N_{y2}x(n):Nx1​⩽n⩽Nx2​;y(n):Ny1​⩽n⩽Ny2​

则rxy(m)r_{xy}(m)rxy​(m)的取值范围为

−(Ny2−Nx1)⩽m⩽(Nx2−Ny1)-\left(N_{y2}-N_{x1}\right)\leqslant m\leqslant\left(N_{x2}-N_{y1}\right)−(Ny2​−Nx1​)⩽m⩽(Nx2​−Ny1​)

若 Ny2⩽Nx1N_{y2}\leqslant N_{x1}Ny2​⩽Nx1​,则 rxy(m)r_{xy}(m)rxy​(m)只在 m⩾0m\geqslant0m⩾0 时有值。

若 Ny1⩾Nx2N_{y1}\geqslant N_{x2}Ny1​⩾Nx2​,则 rxy(m)r_xy(m)rx​y(m)只在 m⩽0m\leqslant0m⩽0 时有值。

rxy(m)=∑n=−∞∞x(n)y(n−m)=∑n=−∞∞x(n)y(−(m−n))=x(m)∗y(−m)\begin{aligned} r_{xy}\left(m\right)& =\sum_{n=-\infty}^{\infty}x(n)y(n-m)=\sum_{n=-\infty}^{\infty}x(n)y(-(m-n)) \\ &=x(m)* y(-m) \end{aligned}rxy​(m)​=n=−∞∑∞​x(n)y(n−m)=n=−∞∑∞​x(n)y(−(m−n))=x(m)∗y(−m)​
rxx(m)=∑x(n)x(n−m)=x(m)∗x(−m)r_{xx}\left(m\right)=\sum x\left(n\right)x\left(n-m\right)=x\left(m\right)*x\left(-m\right)rxx​(m)=∑x(n)x(n−m)=x(m)∗x(−m)

rxx(m)r_{xx}(m)rxx​(m)是实偶序列,即rxx(m)=rxx(−m)r_{xx}(m)=r_{xx}(-m)rxx​(m)=rxx​(−m)

m=0m=0m=0时自相关序列取最大值,即rxx(0)=∑n=−∞∞x2(n)>∣rxx(m)∣r_{xx}(0)=\sum_{n=-\infty}^{\infty}x^2(n)>\mid r_{xx}(m)\mid rxx​(0)=∑n=−∞∞​x2(n)>∣rxx​(m)∣

当x(n)x(n)x(n)是绝对可和信号时,则有lim⁡m→∞rxx(m)=0\lim_{m\to\infty}r_{xx}\left(m\right)=0limm→∞​rxx​(m)=0

rxy(m)=lim⁡M→∞12M+1∑n=−MMx(n)y(n−m)rxx(m)=lim⁡M→∞12M+1∑n=−MMx(n)x(n−m)r_{xy}\left(m\right)=\lim_{M\to\infty}\frac1{2M+1}\sum_{n=-M}^Mx\left(n\right)y(n-m)\\r_{xx}\left(m\right)=\lim_{M\to\infty}\frac1{2M+1}\sum_{n=-M}^Mx\left(n\right)x\left(n-m\right)rxy​(m)=M→∞lim​2M+11​n=−M∑M​x(n)y(n−m)rxx​(m)=M→∞lim​2M+11​n=−M∑M​x(n)x(n−m)
rxy(m)=1N∑n=0N−1x(n)y(n−m)rxx(m)=1N∑n=0N−1x(n)x(n−m)r_{xy}\left(m\right)=\frac{1}{N}\sum_{n=0}^{N-1}x\left(n\right)y(n-m)\\r_{xx}\left(m\right)=\frac{1}{N}\sum_{n=0}^{N-1}x\left(n\right)x\left(n-m\right)rxy​(m)=N1​n=0∑N−1​x(n)y(n−m)rxx​(m)=N1​n=0∑N−1​x(n)x(n−m)
δ(n)={1,n=00,n≠0\delta(n) = \begin{cases}1 ,\quad n = 0\\\\0 ,\quad n \neq 0\end{cases}δ(n)=⎩⎨⎧​1,n=00,n=0​
u(n)={1,n⩾00,n<0u(n)=\begin{cases}1,&n\geqslant0\\[2ex]0,&n<0\end{cases}u(n)=⎩⎨⎧​1,0,​n⩾0n<0​
RN(n)={1,0⩽n⩽N−10,其他 nR_N(n)=\begin{cases}1,&0\leqslant n\leqslant N-1\\\\0,&\text{其他 }n\end{cases}RN​(n)=⎩⎨⎧​1,0,​0⩽n⩽N−1其他 n​
x(n)=anu(n),a 为实数x(n)=a^nu(n),a\text{ 为实数}x(n)=anu(n),a 为实数
x(n)=e(σ+jω0)n=eσn[cos⁡ω0n+jsin⁡(ω0n)]x(n)=\mathrm{e}^{(\sigma+\mathrm{j}\omega_0)n}=\mathrm{e}^{\sigma n}\begin{bmatrix}\cos\omega_0n+\mathrm{j}\sin(\omega_0n)\end{bmatrix}x(n)=e(σ+jω0​)n=eσn[cosω0​n+jsin(ω0​n)​]

当指数为纯虚数时,复指数序列为x(n)=ejω0nx(n)=\mathrm{e}^{\mathrm{j}\omega_0n}x(n)=ejω0​n

x(n)=Asin⁡(ω0n+φ)x(n)=A\sin(\omega_0n+\varphi)x(n)=Asin(ω0​n+φ)

其中 A 为幅度,ω0,\omega_0,ω0​为数字频率,φ,\varphi,φ为起始相位

δ(n)\delta(n)δ(n)与u(n)u(n)u(n)的关系为

δ(n)=u(n)−u(n−1)u(n)=∑k=−∞nδ(k)={1,n⩾00,n<0\delta(n)=u(n)-u(n-1)\\ u(n)=\sum_{k=-\infty}^n\delta(k)=\begin{cases}1,&n\geqslant0\\\\0,&n<0\end{cases}δ(n)=u(n)−u(n−1)u(n)=k=−∞∑n​δ(k)=⎩⎨⎧​1,0,​n⩾0n<0​

δ(n)\delta(n)δ(n)和RN(n)R_N(n)RN​(n)的关系为

RN(n)=∑m=0N−1δ(n−m)=δ(n)+δ(n−1)+⋯+δ(n−N+R_N(n)=\sum_{m=0}^{N-1}\delta(n-m)=\delta(n)+\delta(n-1)+\cdots+\delta(n-N+RN​(n)=m=0∑N−1​δ(n−m)=δ(n)+δ(n−1)+⋯+δ(n−N+

u(n)u(n)u(n)与RN(n)R_N(n)RN​(n)的关系为

RN(n)=u(n)−u(n−N)R_N(n)=u(n)-u(n-N)RN​(n)=u(n)−u(n−N)

因果序列x(n)x(n)x(n)与单位阶跃序列u(n)u(n)u(n)卷积,得到的是该因果序列的累加序列,即

u(n)∗x(n)=∑m=0nx(m)u(n−m)=∑m=0nx(m)u(n)*x(n)=\sum_{m=0}^nx(m)u(n-m)=\sum_{m=0}^nx(m)u(n)∗x(n)=m=0∑n​x(m)u(n−m)=m=0∑n​x(m)

对于anu(n)a^nu(n)anu(n),∣a∣<1|a|<1∣a∣<1时序列收敛,∣a∣>1|a|>1∣a∣>1时序列发散

对于e(σ+jω0)n\mathrm{e}^{(\sigma+\mathrm{j}\omega_0)n}e(σ+jω0​)n,可以将它分解为实部序列与虚部序列或分解为模序列与相角序列

x(n)=e(σ+jω0)n=eσnejω0nx(n)=\mathrm{e}^{(\sigma+\mathrm{j}\omega_0)n}=\mathrm{e}^{\sigma n}\mathrm{e}^{\mathrm{j}\omega_0n}x(n)=e(σ+jω0​)n=eσnejω0​n

数字角频率(ω)(\omega)(ω)和模拟角频率(Ω)(\Omega)(Ω)之间的关系为

ω=ΩT=Ω/fs=2πf/fs(rad)\omega=\Omega T=\Omega/f_{s}=2\pi f/f_{s}\left(rad\right)ω=ΩT=Ω/fs​=2πf/fs​(rad)

其中fs=1Tf_s=\frac{1}{T}fs​=T1​表示抽样频率,所以数字角频率事模拟角频率被抽样频率fsf_sfs​归一化后的弧度,关系可用下表表示

x(n)=x(n+Nr),r=0,±1,±2,⋅⋅⋅x(n)=x(n+Nr),\quad r=0,\pm1,\pm2,\cdotp\cdotp\cdotpx(n)=x(n+Nr),r=0,±1,±2,⋅⋅⋅

正弦型模拟信号xa(t)=Asin(Ω0t+φ)x_a( t) =Asin( \Omega _0t+ \varphi )xa​(t)=Asin(Ω0​t+φ)中,Ω0\Omega_0Ω0​越大,则xa(t)x_a(t)xa​(t)变化越快;但由于x(n)=sin⁡(ω0n+φ)=sin⁡[(ω0+2πm)n+φ]x(n)=\sin(\omega_0n+\varphi)=\sin[(\omega_0+2\pi m)n+\varphi]x(n)=sin(ω0​n+φ)=sin[(ω0​+2πm)n+φ],当ω0\omega_0ω0​变化时,x(n)x(n)x(n)是以 2π2\pi2π 为周期的,并不是 ω0\omega_0ω0​ 越大$,x(n)x(n)x(n)变化越快。

正弦型模拟信号xa(t)=Asin(Ω0t+φ)x_a( t) =Asin( \Omega _0t+ \varphi )xa​(t)=Asin(Ω0​t+φ)一定是周期性信号,周期为T0=2π/Ω0T_0=2\pi /\Omega_0T0​=2π/Ω0​

正弦序列x(n)=Asin⁡(nω0+φ)x(n)=A\sin(n\omega_0+\varphi)x(n)=Asin(nω0​+φ)成为周期性序列的条件:

2π/ω2\pi/\omega2π/ω为整数时,周期为2π/ω2\pi/\omega2π/ω

2π/ω2\pi/\omega2π/ω不是整数但是有理数时,周期将大于2π/ω2\pi/\omega2π/ω

2π/ω2\pi/\omega2π/ω是无理数时,不存在周期

如果正弦型序列时有一个连续正弦信号抽样得到的,抽样频率为fsf_sfs​,抽样时间间隔为TTT,连续信号周期为TTT,则有

2πω0=2πΩ0T=1f0T=T0T=fsf0\frac{2\pi}{\omega_0}=\frac{2\pi}{\Omega_0T}=\frac1{f_0T}=\frac{T_0}T=\frac{f_s}{f_0}ω0​2π​=Ω0​T2π​=f0​T1​=TT0​​=f0​fs​​

抽样时间间隔TTT和连续时间信号周期T0T_0T0​之间的关系如下

当 2π/ω0\pi/\omega_0π/ω0​ 是整数时,例如 2π/ω0=N2\pi/\omega_0=N2π/ω0​=N,则正弦序列的周期为 NNN,有

T0=NT或fs=Nf0T_0=NT\quad\text{或}\quad f_s=Nf_0T0​=NT或fs​=Nf0​

即一个正弦信号周期(T0T_0T0​)中有NNN个抽样周期(TTT),或说抽样频率(fsf_sfs​)是正弦信号频率(f0f_0f0​)的整数(N)倍。

当2π/ω02\pi/\omega_02π/ω0​是有理数时,即2π/ω0=N/M2\pi/\omega_0=N/M2π/ω0​=N/M,其中,NNN、MMM是互为素数的正整数,则正弦序列的周期为 NNN,有

MT0=NT或Mfs=Nf0MT_0=NT\quad\text{或}\quad Mf_s=Nf_0MT0​=NT或Mfs​=Nf0​

即MMM个正弦信号周期(T0T_0T0​)中应有 N 个抽样周期(TTT),或说NNN倍正弦信号频率(f0f_0f0​)应等于MMM倍抽样频率。

δ(n)\delta(n)δ(n)的选择性

x(m)⋅δ(n−m)={x(n),当 m=n0,其他 mx(m)\cdot \delta(n-m)=\begin{cases}x(n),&\text{当 }m=n\\\\0,&\text{其他 }m\end{cases}x(m)⋅δ(n−m)=⎩⎨⎧​x(n),0,​当 m=n其他 m​

任意序列x(n)x(n)x(n)与单位抽样序列δ(n)\delta(n)δ(n)的卷积和等于该序列x(n)x(n)x(n)本身

x(n)=∑m=−∞∞x(m)δ(n−m)=x(n)∗δ(n)x(n)=\sum_{m=-\infty}^\infty x(m)\delta(n-m)=x(n)*\delta(n)x(n)=m=−∞∑∞​x(m)δ(n−m)=x(n)∗δ(n)

任 意 序 列x(n)x( n)x(n)与单位抽样序列的移位序列 δ(n−n0)\delta(n-n_0)δ(n−n0​)的卷积和就得到此序列作相同移位的序列x(n−n0)x(n-n_0)x(n−n0​)

x(n)∗δ(n−n0)=∑m=−∞∞x(m)δ(n−n0−m)=x(n−n0)x(n)*\delta(n-n_0)=\sum_{m=-\infty}^\infty x(m)\delta(n-n_0-m)=x(n-n_0)x(n)∗δ(n−n0​)=m=−∞∑∞​x(m)δ(n−n0​−m)=x(n−n0​)
y(n)=T[x(n)]y(n)=T[x(n)]y(n)=T[x(n)]
a1y1(n)+a2y2(n)=a1T[x1(n)]+a2T[x2(n)]=T[a1x1(n)+a2x2(n)]a_1y_1(n)+a_2y_2(n)=a_1T[x_1(n)]+a_2T[x_2(n)]=T[a_1x_1(n)+a_2x_2(n)]a1​y1​(n)+a2​y2​(n)=a1​T[x1​(n)]+a2​T[x2​(n)]=T[a1​x1​(n)+a2​x2​(n)]

增量线性系统(例如y(n)=4x(n)+6y(n)=4x(n)+6y(n)=4x(n)+6),任意两个输入的响应之差与两个输入之差呈线性关系

T[x(n)]=y(n)T[x(n)]=y(n)T[x(n)]=y(n)
T[x(n−n0)]=y(n−n0)T[x(n-n_0)]=y(n-n_0)T[x(n−n0​)]=y(n−n0​)

若系统有一个移变增益,例如y(n)=nx(n)y(n)=nx(n)y(n)=nx(n),则一定是移变系统

若一个系统在时间轴上有压缩或扩展,例如y(n)=x(Mn)y(n)=x(Mn)y(n)=x(Mn),则一定是移变系统

单位抽样响应(也称单位冲激响应或单位脉冲响应)是指输入为单位抽样序列δ(n)\delta(n)δ(n)时,LSI 系统的输出序列(或称输出响应),一般用h(n)h(n)h(n)表示,即

h(n)=T[δ(n)]h(n)=T[\delta(n)]h(n)=T[δ(n)]
y(n)=∑m=−∞∞x(m)h(n−m)=x(n)∗h(n)y\left(n\right)=\sum_{m=-\infty}^{\infty}x\left(m\right)h\left(n-m\right)=x\left(n\right)*h\left(n\right)y(n)=m=−∞∑∞​x(m)h(n−m)=x(n)∗h(n)
x(n)∗h(n)=∑m=−∞∞x(m)h(n−m)=∑m=−∞∞x(n−m)h(m)=h(n)∗x(n)x(n)* h(n)=\sum_{m=-\infty}^\infty x(m)h(n-m)=\sum_{m=-\infty}^\infty x(n-m)h(m)=h(n)* x(n)x(n)∗h(n)=m=−∞∑∞​x(m)h(n−m)=m=−∞∑∞​x(n−m)h(m)=h(n)∗x(n)
x(n)∗h1(n)∗h2(n)=[x(n)∗h1(n)]∗h2(n)=x(n)∗[h1(n)∗h2(n)]→交换律[x(n)∗h2(n)]∗h1(n)\begin{aligned}x(n)*h_{1}(n)*h_{2}(n)&=\begin{bmatrix}x(n)*h_1(n)\end{bmatrix}*h_2(n)=x(n)*\begin{bmatrix}h_1(n)*h_2(n)\end{bmatrix}\\&\xrightarrow{\text{交换律}}\begin{bmatrix}x(n)*h_2(n)\end{bmatrix}*h_1(n)\end{aligned}x(n)∗h1​(n)∗h2​(n)​=[x(n)∗h1​(n)​]∗h2​(n)=x(n)∗[h1​(n)∗h2​(n)​]交换律​[x(n)∗h2​(n)​]∗h1​(n)​
x(n)∗[h1(n)+h2(n)]=x(n)∗h1(n)+x(n)∗h2(n)x(n)*\begin{bmatrix}h_1(n)+h_2(n)\end{bmatrix}=x(n)*h_1(n)+x(n)*h_2(n)x(n)∗[h1​(n)+h2​(n)​]=x(n)∗h1​(n)+x(n)∗h2​(n)
y(n0) 只取决于 x(n)∣n⩽n0y(n_0)\text{ 只取决于 }x(n)\mid_{n\leqslant n_0}y(n0​) 只取决于 x(n)∣n⩽n0​​

对于因果系统,n<n0n<n_0n<n0​的输入相同则n<n0n<n_0n<n0​的输出也相同

考察系统因果性时,只看输入x(n)x(n)x(n)和输出y(n)y(n)y(n)的关系,而不讨论其他以nnn为变量函数的影响,例

y(n)=(n+2)x(n) 是因果系统y(n)=x(n)sin⁡(n+4)也是因果系统\begin{aligned}&y(n)=(n+2)x(n)\text{ 是因果系统}\\&y(n)=x(n)\sin(n+4)\text{也是因果系统}\end{aligned}​y(n)=(n+2)x(n) 是因果系统y(n)=x(n)sin(n+4)也是因果系统​

LSI系统是因果系统的充要条件是单位冲激响应h(n)h(n)h(n)是因果序列,即

h(n)=0,n<0h(n)=0,n<0h(n)=0,n<0
若∣x(n)∣⩽M<∞则∣y(n)∣⩽P<∞\begin{aligned} &若\mid x(n)\mid\leqslant M<\infty\\ &则\mid{y}(n)\mid\leqslant P<\infty\end{aligned}​若∣x(n)∣⩽M<∞则∣y(n)∣⩽P<∞​
∑n=−∞∞∣h(n)∣=P<∞\sum_{n=-\infty}^{\infty}\mid h\left(n\right)\mid=P<\inftyn=−∞∑∞​∣h(n)∣=P<∞
∑k=0Naky(n−k)=∑m=0Mbmx(n−m)\sum_{k=0}^Na_ky\left(n-k\right)=\sum_{m=0}^Mb_mx\left(n-m\right)k=0∑N​ak​y(n−k)=m=0∑M​bm​x(n−m)
x(n)∣n<n0=0x\left(n\right)|_{n<n_{0}}=0x(n)∣n<n0​​=0
y(n)∣n<n0=0y\left(n\right)|_{n<n_{0}}=0y(n)∣n<n0​​=0

设 xu(t)x_u(t)xu​(t)为模拟信号,x^a(t)\hat{x}_a(t)x^a​(t)为理想抽样信号,则有

x^a(t)=xa(t)⋅p(t)=xa(t)⋅δT(t)\hat{x}_a(t)=x_a(t)\cdot p(t)=x_a(t)\cdot\delta_T(t)x^a​(t)=xa​(t)⋅p(t)=xa​(t)⋅δT​(t)

其中,抽样信号 p(t)p(t)p(t)就是周期性的单位冲激信号

δT(t)=∑m=−∞∞δ(t−mT)\delta_T(t)=\sum_{m=-\infty}^\infty\delta(t-mT)δT​(t)=m=−∞∑∞​δ(t−mT)
x^a(t)=xa(t)∑m=−∞∞δ(t−mT)=∑m=−∞∞xa(t)δ(t−mT)=∑m=−∞∞xa(mT)δ(t−mT)\begin{aligned}\hat{x}_{a}(t)&=x_a(t)\sum_{m=-\infty}^\infty\delta(t-mT)=\sum_{m=-\infty}^\infty x_a(t)\delta(t-mT)\\&=\sum_{m=-\infty}^\infty x_a(mT)\delta(t-mT)\end{aligned}x^a​(t)​=xa​(t)m=−∞∑∞​δ(t−mT)=m=−∞∑∞​xa​(t)δ(t−mT)=m=−∞∑∞​xa​(mT)δ(t−mT)​

TTT为采样周期,fs=1/Tf_s=1/Tfs​=1/T为抽样频率,Ωs=2πfs=2π/T\Omega_s=2\pi f_s=2\pi/TΩs​=2πfs​=2π/T是抽样角频率

X^a(jΩ)=1T∑k=−∞∞Xa[j⁡(Ω−k2πT)]\hat{X}_a\left(\mathrm{j}\Omega\right)=\frac1T\sum_{k=-\infty}^\infty X_a\left[\operatorname{j}\left(\Omega-k\frac{2\pi}T\right)\right]X^a​(jΩ)=T1​k=−∞∑∞​Xa​[j(Ω−kT2π​)]

理想抽样信号的频谱X^a(jΩ)\hat{X} _a(j\Omega)X^a​(jΩ)是被抽样的模拟信号的频谱 Xa(jΩ)X_a(j\Omega)Xa​(jΩ)的周期延拓,在角频率Ω\OmegaΩ轴上其延拓周期为

Ωs=2πT=2πfs\Omega_s=\frac{2\pi}{T}=2\pi f_sΩs​=T2π​=2πfs​

即频率轴fff上以抽样频率fsf_sfs​为周期而周期延拓。

如果模拟信号是带限信号(频带有限信号),信号的最高频率分量为 fhf_hfh​,设抽样频率为 fsf_sfs​,由于抽样后信号的频谱等于模拟信号频谱按抽样频率 fsf_sfs​作周期延拓,只有当

fh<fs2f_h<\frac{f_s}2fh​<2fs​​

时,周期延拓的频谱分量才不会产生交叠,称抽样频率之半(fs/2f_s/2fs​/2)为折叠频率。即

fs2=12T或Ωs2=πT\frac{f_s}2=\frac1{2T}\quad\text{或}\quad\frac{\Omega_s}2=\frac\pi T2fs​​=2T1​或2Ωs​​=Tπ​

奈奎斯特抽样定理:若 xa(t)x_a(t)xa​(t)是频带宽度有限的信号(称带限信号),要想抽样后的信号能够不失真地还原出原信号,则必须抽样频率fsf_sfs​大于信号最高频率分量fhf_hfh​ 的两倍,即

fs>2fhf_{s}>2f_{h}fs​>2fh​

数字频率ω\omegaω是模拟角频率对抽样频率fsf_sfs​的归一化频率,因而抽样频率fsf_sfs​对于的数字频率为

ωs=ΩsT=2πfs/fs=2π\omega_{s}=\Omega_{s}T=2\pi f_{s}/f_{s}=2\piωs​=Ωs​T=2πfs​/fs​=2π

因此数字域的抽样频率ωs\omega_sωs​等于2π2\pi2π,也就是延拓周期等于2π2\pi2π,折叠频率fs/2f_s/2fs​/2所对应的数字域频率为

ωs2=π\frac{\omega_s}{2}=\pi2ωs​​=π

按照采样定理,信号数字域最高频率ωh\omega_hωh​应小于π\piπ。

在实际中,抽样脉冲串不是冲激函数串δT(t)\delta_T(t)δT​(t),而是一点宽度τ\tauτ的矩形周期脉冲串,这时奈奎斯特抽样定理仍然有效

其最高频率为fhf_hfh​,带宽为Δfo\Delta f_oΔfo​,算术中心频率为fa=fh−Δfo2f_a=f_h-\frac{\Delta f_o}{2}fa​=fh​−2Δfo​​,一般来说fa≫Δfof_a\gg\Delta f_ofa​≫Δfo​,即通带中心频率远大于通带带宽。

当fh=rΔfof_h=r\Delta f_ofh​=rΔfo​,即带通信号的最高频率是其通带宽度的整数倍时,则选抽样频率fsf_sfs​为

fs=2Δfof_s=2\Delta f_ofs​=2Δfo​

即所取抽样频率为带通信号通带宽度的两倍,其抽样后的频谱是带通信号频谱以此fsf_sfs​的整数倍而周期延拓后的频谱,如图

图中 r=5r=5r=5,显然没有频谱混叠现象。

当 fh=r′Δfof_h=r^{\prime}\Delta f_ofh​=r′Δfo​,r′≠r^{\prime}\neqr′=整数。这是最一般的情况,即带通信号最高频率不等于其带宽的整数倍时,这时,可保持fhf_hfh​不变,将通带下端延伸到使其带宽为Δfo′\Delta f_o^{\prime}Δfo′​,使得满足

fh=rΔf0′,r=⌊r′⌋(⌊ ⌋表示取整数部分)f_h= r\Delta f_0^{\prime },r= \lfloor r^{\prime }\rfloor( \lfloor \:\rfloor表示取整数部分)fh​=rΔf0′​,r=⌊r′⌋(⌊⌋表示取整数部分)

此时有 Δfo′>Δfo\Delta f_o^{\prime}>\Delta f_oΔfo′​>Δfo​,这时选抽样频率为

fs=2Δfo′=2fh/⌊r′⌋=2fh/⌊fh/Δfo⌋f_s=2\Delta f'_o=2f_h/\lfloor r'\rfloor=2f_h/\lfloor f_h/\Delta f_o\rfloorfs​=2Δfo′​=2fh​/⌊r′⌋=2fh​/⌊fh​/Δfo​⌋

带通信号抽样频率fsf_sfs​的取值范围为(其中Δfo\Delta f_oΔfo​为带宽)

2Δfo⩽fs<4Δfo2\Delta f_{o}\leqslant f_{s}<4\Delta f_{o}2Δfo​⩽fs​<4Δfo​

fsf_{s}fs​的下限是fhΔfo=\frac {f_h}{\Delta f_o}=Δfo​fh​​=整数时的情况,而其上限则对应于fhΔfo≠\frac {f_h}{\Delta f_o}\neqΔfo​fh​​=整数,且是最不利的情况。由于上面讨论的带通信号的抽样频率fsf_sfs​不满足fs>2fhf_s>2f_hfs​>2fh​的要求,因而可称其为亚奈奎斯特抽样频率。

1.4.5 连续时间信号 xa(t)x_a(t)xa​(t)、理想抽样信号x^a(t)\hat{x}_a(t)x^a​(t)以及抽样序列 x(n)x(n)x(n)

x^a(t)\hat{x}_a(t)x^a​(t)在工程应用中无法实现,故工程中采样抽样序列x(n)x(n)x(n)。

x^a(t)\hat{x}_a(t)x^a​(t)和 x(n)x(n)x(n)的本质差别:x^a(t)\hat{x}_a(t)x^a​(t)本质上是连续时间信号,是在 t=nTt=nTt=nT 时的冲激串, 在t≠nTt\neq nTt=nT 时,x^a(t)=0,\hat{x}_a(t)=0,x^a​(t)=0,其每个冲激的幅度都是无穷大,存在时间(宽度)为无穷小,x^a(t)\hat{x}_a(t)x^a​(t)的大小是以冲激的积分面积表示的;抽样序列 x(n)x(n)x(n)则是整数变量 nnn 的函数,这里时间已归一化,x(n)x(n)x(n)本身已没有抽样率的信息,当 n=n=n=整数时,x(n)=xa(nT)x(n)=x_a(nT)x(n)=xa​(nT),即抽样点上抽样序列的幅值是确定数值,在 n≠n\neqn=整数时,x(n)x(n)x(n)无定义,不是零值。

连续时间信号 xa(t)x_a(t)xa​(t)、理想抽样信号x^a(t)\hat{x}_a(t)x^a​(t)以及抽样序列 x(n)x(n)x(n)三者的关系

x(n)=xa(t)∣t=nT=xa(nT),−∞<n<∞x\left(n\right)=x_{a}\left(t\right)\Big|_{t=nT}=x_{a}\left(nT\right),-\infty<n<\inftyx(n)=xa​(t)​t=nT​=xa​(nT),−∞<n<∞

如果满足奈奎斯特抽样定理,即信号谱的最高频率小于折叠频率,则抽样后不会产生频谱混叠,可以由信号的抽样值经插值而重构原信号xa(t)x_a(t)xa​(t)。当∣Ω∣<Ωs/2|\Omega|<\Omega_s/2∣Ω∣<Ωs​/2时,只存在k=0k=0k=0一项,即有

X^a(jΩ)=1TXa(jΩ),∣Ω∣<Ωs2\hat{X}_a(\text{j}\Omega)=\frac1TX_a(\text{j}\Omega),\quad\mid\Omega\mid<\frac{\Omega_s}2X^a​(jΩ)=T1​Xa​(jΩ),∣Ω∣<2Ωs​​

将X^a(jΩ)\hat{X}_a(j\Omega)X^a​(jΩ)送人作为重构用的理想低通滤波器 H(jΩ)H(j\Omega)H(jΩ)

H(jΩ)={T,∣Ω∣<Ωs20,∣Ω∣⩾Ωs2H(\text{j}\Omega)=\begin{cases}T,&\mid\Omega\mid<\frac{\Omega_s}2\\\\0,&\mid\Omega\mid\geqslant\frac{\Omega_s}2\end{cases}H(jΩ)=⎩⎨⎧​T,0,​∣Ω∣<2Ωs​​∣Ω∣⩾2Ωs​​​
Ya(jΩ)=X^a(jΩ)H(jΩ)=Xa(jΩ)Y_a(j\Omega)=\hat{X}_a(j\Omega)H(j\Omega)=X_a(j\Omega)Ya​(jΩ)=X^a​(jΩ)H(jΩ)=Xa​(jΩ)
ya(t)=∑m=−∞xa(mT)sin⁡[π(t−mT)/T]π(t−mT)/Ty_{a}\left(t\right)=\sum_{m=-\infty}x_{a}\left(mT\right)\frac{\sin\left[\pi\left(t-mT\right)/T\right]}{\pi\left(t-mT\right)/T}ya​(t)=m=−∞∑​xa​(mT)π(t−mT)/Tsin[π(t−mT)/T]​

xo(mT)x_o(mT)xo​(mT)经此公式而得到连续信号xa(t)x_a(t)xa​(t),而sin⁡[π(t−mT)/T]/[π(t−mT)/T]\sin[\pi(t-mT)/T]/[\pi(t-mT)/T]sin[π(t−mT)/T]/[π(t−mT)/T]称为内插函数,如图

在抽样点mTmTmT上,函数值为 1,在其余抽样点上,函数值为零,不影响其他抽样点。也就是说,xa(t)x_a(t)xa​(t) 等于各xa(mT)x_{a}(mT)xa​(mT)乘以对应的内插函数的总和。在每一个抽样点上,只有该点所对应的内插函数不为零,这使得各抽样点上信号值不变,而抽样点之间的信号则由各加权抽样函数波形的延伸叠加而成,如图

因而内插函数 h(t)h(t)h(t)的作用是在各抽样点之间起连续插值的作用。

正弦型信号 xa(t)=Asin⁡(Ωot+φ)=Asin⁡(2πfat+φ)x_a(t)=A\sin(\Omega_ot+\varphi)=A\sin(2\pi f_at+\varphi)xa​(t)=Asin(Ωo​t+φ)=Asin(2πfa​t+φ)的频谱在 f=fof=f_of=fo​处为δ\deltaδ函数,一般来说,正弦信号的抽样频率必须满足fs>2fo.f_s>2f_o.fs​>2fo​. 因为若取 fs=2fof_s=2f_ofs​=2fo​,则有以下几种情况发生:

当φ=0\varphi=0φ=0时,一个周期抽取的两个点为x(0)=x(1)=0x(0)=x(1)=0x(0)=x(1)=0,相当于xα(0)x_{\alpha}(0)xα​(0)和xα(π)x_{\alpha}(\pi)xα​(π) 两个点,故不包含原信号的任何信息。

当φ=π2\varphi=\frac\pi2φ=2π​时,则有x(0)=A,x(1)=−Ax(0)=A,x(1)=-Ax(0)=A,x(1)=−A,此时从x(n)x(n)x(n)可以恢复xa(t)x_a(t)xa​(t)。

当φ\varphiφ为已知,且 0<φ<π20<\varphi<\frac\pi20<φ<2π​时,恢复的不是原信号,但经过变换后,可得到原信号。

当φ\varphiφ为未知数时,抽样后不能恢复出原信号xa(t)x_a(t)xa​(t)。所以,至少要取 fs>2fof_s>2f_ofs​>2fo​,避免产生不确定性。

考虑到作DFTDFTDFT时,当要求数据个数为N=2PN=2^PN=2P时,正弦信号一个周期中最好抽取4个点。

满足∣ω2−ω1∣=2πk|\omega_{2}-\omega_{1}|=2\pi k∣ω2​−ω1​∣=2πk,kkk为整数,的两个正弦序列是相同的,有以下两个结论:

两个不同频率的模拟正弦型信号,如果用同一抽样频率fsf_sfs​对其抽样,得到的序列可能是相同的序列,我们没法从序列中区分出它们分别来源于哪一个模拟正弦型信号。

🐷
⚽
Page cover image