【压缩感知合集1】（背景知识）香农奈奎斯特采样定理的数学推导和图解分析

最新推荐文章于 2024-05-08 10:05:31 发布

呆呆象呆呆

最新推荐文章于 2024-05-08 10:05:31 发布

阅读量4.5k

点赞数 9

分类专栏：理论知识学习

本文链接：https://blog.csdn.net/qq_41554005/article/details/118616089

版权

理论知识学习专栏收录该内容

50 篇文章 115 订阅

订阅专栏

【压缩感知合集1】（背景知识）香农奈奎斯特采样定理的数学推导和图解分析
 【压缩感知合集2】（背景知识）信号稀疏表示的数学推导和解释理解
 【压缩感知合集3】压缩感知的背景与意义
 【压缩感知合集4】（背景知识）理想采样信号和随机采样信号两种采样信号的频谱分析，以及采样效果比较

香农奈奎斯特采样定理

1 评价

支配先阶段几乎所有信号的获取、处理、存储、传输过程。
采样后再进行压缩的方式浪费了大量的采样资源，如果采样后的信号长度仍然很长，那么变换会消耗很长时间。
由于需要保留的 $K$ 个重要分量的位置是随着信号的不同而不同，所以这种编解码方式是自适应的，需要分配多余的存储空间以保留 $K$ 个重要分量的位置。
$K$ 个重要分量有可能在传输过程中丢失其中的某几个分量从而造成较差的抗干扰能力。

2 数学解释

傅里叶分析示意图

数学模型证明

输入连续信号 $x_a(t)$

采样输出信号 $x_s(t)$

通常意义上采样的过程是一段周期的矩形脉冲载波信号 $S (t)$ ，信号周期为 $T$ ，宽度为 $\tau$ ， $\tau$ 越小采样越精准（理解：保留的值越接近时序点的对应值），当 $\tau<<T$ 的时候，采样脉冲信号接近于 $\delta$ 函数信号 $S(t)=\sum_{-\infty}^{\infty} \delta(t-n T)$

理想采样

$x_{s}(t)=x_{a}(t) s(t)$

信号 $x_a(t)$ 的傅里叶变换

$X_{a}(j \omega)=\int_{-\infty}^{\infty} x_{a}(t) e^{-j \omega t} d t$

信号 $S (t)$ 的傅里叶变换

$\omega)=\int_{-\infty}^{\infty} S(t) e^{-j \omega t} d t$

信号 $x_s(t)$ 的傅里叶变换

$X_{s}(j \omega)=\int_{-\infty}^{\infty} x_{s}(t) e^{-j \omega t} d t$

对于理想采样 $x_{s}(t)=x_{a}(t) s(t)$ 有时域相乘对应频域卷积

$X_{s}(j \omega)=\frac{1}{2 \pi} X_{a}(j \omega) * S(j \omega)$

根据补充证明可以知道

$\omega)=\frac{2 \pi}{T} \sum_{k=-\infty}^{\infty} \delta\left(\omega-k \omega_{s}\right)$

所以
$\begin{aligned} &X_{s}(j \omega)=\frac{1}{2 \pi}\left[\frac{2 \pi}{T} \sum_{k=-\infty}^{\infty} \delta\left(\omega-k \omega_{s}\right) * X_{a}(j \omega)\right] \\ &=\frac{1}{T} \int_{-\infty}^{\infty} X_{a}(j \theta) \sum_{k=-\infty}^{\infty} \delta\left(\omega-k \omega_{s}-\theta\right) d \theta \end{aligned}$
根据冲击函数的性质（见下图的冲激函数示意图）

可以得到下面的公式推导
$X_{s}(j \omega)=\frac{1}{T} \sum_{k=-\infty}^{\infty} X_{a}\left(j \omega-j k \omega_{s}\right)$
这也就相当于说针对原始被采样信号的频谱做了一个周期延拓，具体可以看如下的示意图

补充说明

每一个周期函数都可以写成完全用傅里叶级数展开的求和形式，将周期采样函数 $S (t)$ 写成如下形式
$S(t)=\sum_{k=-\infty}^{\infty} a_{k} e^{j k \omega_{s} t}$
其中 $\omega_{s}=\frac{2 \pi}{T}$

然后我们求一下傅里叶级数的系数
$a_{k}=\frac{1}{T} \int_{-T / 2}^{T / 2} \sum_{-\infty}^{\infty} \delta(t-n T) e^{-j k \omega_{s} t} d t=\frac{1}{T}$
发现以采样频率 $\omega_{s}$ 作为基波的任意整数倍频带内都有相同的能量也就是 $\frac{1}{T}$ 在此提一下，正是因为有这个均匀的能量分布才有了后面的周期延拓

傅里叶分解后采样信号的每一个级数分量（在不考虑系数情况下），都是单独频率点唯一的正弦类函数，又因为线性时不变系统的原因，采样信号 $S (t)$ 和被采样信号 $x_{s}(t)$ 相乘可以理解为被采样信号的每一个傅里叶级数分量进行分别的调制，而被调制搬移到了不同的频谱。

现在考虑系数理想采样信号，拥有性质：任意傅里叶级数能量相同，所以搬移调制后会有着很好的频谱能量不变性。（即在每一个频率分量段内都保持同等放大缩小幅度）。

继续推导，有了刚才的傅里叶级数的系数，带入傅里叶展开之后再看一次
$S(t)=\frac{1}{T} \sum_{k=-\infty}^{\infty} e^{j k \omega_{s} t}$
有了这个级数展开之后我们去看频域的公式
$\begin{aligned} S(j \omega)=\int_{-\infty}^{\infty} S(t) e^{-j \omega t} d t \\=\int_{-\infty}^{\infty} (\frac{1}{T} \sum_{k=-\infty}^{\infty} e^{j k \omega_{s} t}) e^{-j \omega t} d t \\=\frac{1}{T} \sum_{k=-\infty}^{\infty} F\left[e^{j k \omega_{s} t}\right] \end{aligned}$
由于（其实也就是因为 $S (t)$ 只有固定频率点 $k\omega_{s}$ 有能量，其他任何与整数倍 $\omega_{s}$ 不相同的 $\omega$ 积分的时候都会变成0）
$F\left[e^{j k \omega_{s} t}\right]=2 \pi \delta\left(\omega-k \omega_{s}\right)$
所以
$\omega)=\frac{2 \pi}{T} \sum_{k=-\infty}^{\infty} \delta\left(\omega-k \omega_{s}\right)$