博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[傅里叶变换及其应用学习笔记] 十七. Ш函数
阅读量:5884 次
发布时间:2019-06-19

本文共 2971 字,大约阅读时间需要 9 分钟。

 

Ш函数的三个性质

上节课我们学习了$Ш_p$函数,其定义如下

$Ш_p = \displaystyle{ \sum_{k=-\infty}^{\infty}\delta(x-kp) }$

 
 
 

$Ш_p$函数有以下三个性质,

1) 采样性质,继承了$\delta$函数的采样性质

$f(x)Ш_p(x) = \displaystyle{ \sum_{k=-\infty}^{\infty}f(kp)\delta(x-kp) }$

 

2) 周期性质,继承了$\delta$函数的移位性质

$(f*Ш_p)(x) = \displaystyle{ \sum_{k=-\infty}^{\infty}f(x-kp) }$

 

3) 傅里叶变换

$\mathcal{F}Ш_p = \frac{1}{p}Ш_{\frac{1}{p}}$

$\mathcal{F}^{-1}Ш_p = \frac{1}{p}Ш_{\frac{1}{p}}$

 

$Ш_p$函数的这三个性质,是本节课后面推导的基础。

 
 
 

内插问题

内插问题是我们接下来要解决的问题,我们需要用可靠的方法对离散的测量或采样值进行内插,通过内插,我们能得到被采样信号的所有的值。(The problem here is and what we're actually gonna solve in a quite remarkable way is the exact interplation of value of a function from a discrete set of measurement or a discrete set of samples. We'll be able to interpolate all values of a signal or a function from a discrete set of samples.)

假设有一个随时间变化的过程,以相等的时间间隔(如几分之一秒)对该过程进行测量,得到一组测量数据$(t_0,y_0),(t_1,y_1),(t_2,y_2)$等等,我们可以把这些数据当成一系列的点

 
 

我们可以用曲线来拟合到采样数据上(fit a curve to the data),或者根据测量数据在中点进行内插(interpolate values of the process at the intermediate points based on the measurements)。

 
 

曲线拟合与内插

1) 曲线拟合是根据采样数据得出一个曲线函数,该函数上的值作为原信号的近似值。
2) 内插是通过公式的代入,求出两采样点中间点的具体数值。内插不一定是线性的,具体数值要根据内插公式来确定。
 
 
两者虽然方法不同,但是目的都是为了从有限的采样点得到原始信号。
 
要求被采样点以外的点数值,这并没有固定的方法,需要根据实际情况来选择合适的方法。当然,更多的测量值,可以提供更高的曲线拟合度或更准确的内插。
 
 

采样点间的不确定性

内插与拟合的不确定性,从极端来看,可以看作振荡,即函数从一个点到另一个点的变化有多快。函数拐弯越频繁,曲线拟合或者内插的不确定性越大,我们需要了解、控制(regulate)这种不确定性。

我们对于信号会从时域,频域两个方面去分析,而在频域的傅里叶变换反映了信号的频率成分,我们能借此分析函数振荡的快慢。傅里叶变换的高频与快速振荡相关,傅里叶变换的低频与低速振荡相关。我们要了解信号的振荡速度,就要对其傅里叶变换进行分析。

 

解决这种不确定性的方法是:规定函数允许振荡的最高频率。如果我们在傅里叶变换后把高频去掉,则相当于把快速振荡去除。

解决方法可以总结为以下定义:

对于一个有限带宽的函数f(x),如果它的傅里叶变换在某频带以外的值恒为零,就是说傅里叶变换$\mathcal{F}f(s)\equiv 0 \ for \ |s|\geqslant \frac{p}{2}$,最小的$p$值就称为带宽。

 
 

对于有限带宽信号,可以完全解决不确定性问题,即可以根据离散的采样值得到该信号的函数表达式$f(x)$

推导过程如下:

 

利用$Ш_p$的周期化性质,用$Ш_p$对$\mathcal{F}f(s)$进行周期化

$\mathcal{F}f * Ш_p$ 

 
 

从周期傅里叶变换恢复为原来的傅里叶变换

$\mathcal{F}f = \Pi_p(\mathcal{F}f * Ш_p)$

 
 

然后求傅里叶逆变换,得到时域信号

$\begin{align*}

f(t)
&=\mathcal{F}^{-1}(\Pi_p(\mathcal{F}f*Ш_p))\\
&=(\mathcal{F}^{-1}\Pi_p)*(\mathcal{F}^{-1}(\mathcal{F}f*Ш_p)) \qquad(Fourier\ Convolution\ Theorem)\\
&=(psinc(pt))*((\mathcal{F}^{-1}\mathcal{F}f(t))(\mathcal{F}^{-1}Ш_p(t)))\\
&=(psinc(pt))*(f(t)\cdot \frac{1}{p}Ш_{\frac{1}{p}}(t))\qquad(Fourier\ Transform\ of\ Ш_p)\\
&=(psinc(pt))*(\frac{1}{p}f(t)\sum_{k=-\infty}^{\infty}\delta(x-\frac{k}{p}))\\
&=(psinc(pt))*(\frac{1}{p}\sum_{k=-\infty}^{\infty}f(\frac{k}{p})\delta(x-\frac{k}{p})) \qquad(Ш_p\ Sampling\ Property)\\
&=\sum_{k=-\infty}^{\infty}f(\frac{k}{p})sinc(pt)*\delta(x-\frac{k}{p})\\
&=\sum_{k=-\infty}^{\infty}f(\frac{k}{p})sinc(p(t-\frac{k}{p}))\qquad (\delta\ shift\ property)
\end{align*}$

 

因此,对于有限带宽函数$f(t)$可写成如下形式,

$f(t) = \displaystyle{ \sum_{k=-\infty}^{\infty}f(\frac{k}{p})sinc(p(t-\frac{k}{p})) }\quad ,\quad \mathcal{F}f(s)\equiv 0 \ for\ |s|\geqslant\frac{p}{2}$

 

结论是:

  • 对于带宽为$P$的函数$f(t)$,如果采样间隔为$\frac{1}{p}$,而且已知所有采样点$f(\frac{k}{p}) ,k=0,\pm1,\pm2 ...$的值,那么我们就能通过该公式内插得到原本的函数$f(t)$的所有的值。

这叫做采样定理(sampling theorem),这可以说是整个课程最重要的公式。

转载地址:http://vgoix.baihongyu.com/

你可能感兴趣的文章
Linux下安装Tomcat 8.0.53
查看>>
10月12日云栖精选夜读 | 程序员工资那么高,为什么会找不到女朋友?
查看>>
轻松搞定RabbitMQ3:消息应答与消息持久化
查看>>
快速下载android源码
查看>>
Linux基础命令---cal
查看>>
Linux下文件重命名、创建、删除、修改及保存文件
查看>>
跟我学习php文件和目录常用函数-上篇
查看>>
关于存量和增量的杂谈
查看>>
JavaScript 单例模式
查看>>
PHP项目学习2
查看>>
Linux Shell常用技巧(目录)
查看>>
在微信网页版下用Chrome控制台发送消息
查看>>
CUBA CLI 2.1.0 发布,企业级应用开发平台
查看>>
lecture 3.1 递归和对象
查看>>
GPS轨迹数据集免费下载资源整理
查看>>
bulk更新mongodb的脚本
查看>>
通用组件示例
查看>>
Flutter 03: 图解第一个程序 Hello World!
查看>>
Js获取当前页面URL各种参数
查看>>
删除链表的倒数第N个节点
查看>>