[傅里叶变换及其应用学习笔记] 六. 热方程讨论

这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用。

上节课讲到,在对非周期函数进行傅里叶分析时,有

$C_k = \displaystyle{\frac{1}{T}\int_{-\frac{T}{2}}^{\frac{T}{2}}f(t)e^{-2\pi i\frac{k}{T}t}dt }$

$f(t) = \displaystyle{\sum_{k=-\infty}^{\infty}C_ke^{2\pi i\frac{k}{T}t} }$

我们希望仅让$T\to \infty$就能得到我们希望的结果:傅里叶变换适用于非周期函数。但结果证明了这样还不可行,最后得出:对任意$C_k$,都有$C_k \leqslant \frac{M}{T}$,当$T \to \infty \ ,C_k \to 0$。$C_k$跟$T$是成反比例的。

按照这种关系,我们是否能把$T$引入到$C_k$这边?

新符号$\eta$

$\displaystyle{\eta f(\frac{k}{T}) =C_k \times T = \int_{-\frac{T}{2}}^{\frac{T}{2}}e^{-2\pi i \frac{k}{T}t}f(t)dt }$

即有,

$f(t) = \displaystyle{\sum_{k=-\infty}^{\infty}\eta f(\frac{k}{T})e^{2\pi i\frac{k}{T}t} \frac{1}{T} }$

现在令$T \to \infty$,那么$\frac{k}{T}$的取值范围为$(\frac{k=-\infty}{T\to \infty},\frac{k=+\infty}{T\to \infty})$,即$(-\infty, +\infty)$。取值间隔为$\frac{1}{T} \to 0$,趋于连续变量。现在用连续变量$s$来表示$\frac{k}{T}$:

$s = \frac{k}{T} \ , –\infty < s < \infty$

$\eta f(s) = \displaystyle{\int_{-\infty}^{\infty}e^{-2\pi ist}f(t)dt }$

由于$\frac{k}{T}$被替换成了连续变量$s$,那么傅里叶级数的多项式会被替换成积分,其中$\frac{1}{T}$为$\bigtriangleup s$,即$ds$

$f(t) = \displaystyle{\int_{-\infty}^{\infty}\eta f(s)e^{2\pi ist}ds }$

结论(定义)

如果$f(t)$的周期被定义在整个实数域中,即$-\infty < T < \infty$,那么其

傅里叶变换:

$\eta f(s) = \displaystyle{\int_{-\infty}^{\infty} e^{-2\pi ist}f(t)dt }$

傅里叶逆变换:

$f(t) = \displaystyle{\int_{-\infty}^{infty}e^{2\pi is0} \eta f(s)ds }$

也可以写作如下形式:

$\eta f(s) = \displaystyle{\int_{-\infty}^{\infty} e^{-2\pi ist}f(t)dt }$

$\eta^{-1}g(t) = \displaystyle{\int_{-\infty}^{infty}e^{2\pi ist}g(s)ds }$

符号$\eta$代表傅里叶正变换,$\eta^{-1}$代表傅里叶逆变换。

傅里叶正变换吧函数分解成连续复指数;傅里叶逆变换把这些连续复指数组合成原函数。

零点的值

$\eta f(0) = \displaystyle{\int_{-\infty}^{\infty} e^{-2\pi i0t}f(t)dt  = \int_{-\infty}^{\infty}f(t)dt } $

$\eta^{-1}g(0) = \displaystyle{\int_{-\infty}^{infty}e^{2\pi is0} g(s)ds  = \int_{-\infty}^{infty}g(s)ds }$

傅里叶变换例子

1. 矩形函数

$\pi (t) =
\left\{\begin{matrix}
1 & \left| t \right| < \frac{1}{2}\\
0 & \left| t \right| \geqslant \frac{1}{2}
\end{matrix}\right.$

傅里叶变换如下:

$\begin{align*}
\eta \pi(s)
&= \displaystyle{\int_{-\infty}^{\infty}e^{-2\pi ist}\pi(t)dt } \\
&= \displaystyle{\int_{-\frac{1}{2}}^{\frac{1}{2}}e^{-2\pi ist}dt } \\
&= \left . -\frac{1}{2\pi is}e^{-2\pi ist} \right |_{\frac{1}{2}}^{\frac{1}{2}} \\
&= -\frac{1}{2\pi is}e^{-2\pi is\frac{1}{2}} - (-\frac{1}{2\pi is}e^{-2\pi is(-\frac{1}{2})}) \\
&= -\frac{1}{2\pi is}e^{-\pi is} + \frac{1}{2\pi is}e^{\pi is} \\
&= \frac{1}{\pi s}(\frac{e^{\pi is} - e^{-\pi is}}{2i}) \\
&= \frac{1}{\pi s}(\frac{cos(\pi s)+isin(\pi s) - cos(-\pi s) - isin(-\pi s)}{2i}) \\
&= \frac{1}{\pi s}(\frac{2isin(\pi s)}{2i}) \\
&= \frac{sin(\pi s)}{\pi s} \\
&= sinc \ s
\end{align*}$

该函数被称为$sinc$函数

2. 三角形函数

$\Lambda (t) =
\left\{\begin{matrix}
1 - \left|t\right| & \left|t\right|<1 \\
0 & \left|t\right| \geqslant 1
\end{matrix}\right.$

傅里叶变换如下:

\begin{align*}\eta\Lambda(s)
&= \displaystyle{\int_{-\infty}^{\infty}e^{-2\pi ist}\Lambda(t)dt}\\
&=\displaystyle{\int_{-1}^{0}e^{-2\pi ist}(1+t)dt + \int_{0}^{1}e^{-2\pi ist}(1-t)dt}\\
&=\left(\left.(1+t)(-\frac{2\pi is}{1}e^{-2\pi ist})\right|_{-1}^0-\displaystyle{\int_{-1}^{0}-\frac{1}{2\pi is}e^{-2\pi ist}dt }\right)+\left(\left.(1-t)(-\frac{2\pi is}{1}e^{-2\pi ist})\right|_{0}^1-\displaystyle{\int_{0}^{1}-\frac{1}{2\pi is}e^{-2\pi ist}dt }\right)\\
&=\left(-\frac{1}{2\pi is}-\left. \frac{1}{4\pi^2i^2s^2}e^{-2\pi ist}\right|_{-1}^{0}\right)+\left(\frac{1}{2\pi is}+\left. \frac{1}{4\pi^2i^2s^2}e^{-2\pi ist}\right|_{0}^{1}\right )\\
&=-\left(\frac{1}{-4\pi^2s^2}-\frac{1}{-4\pi^2s^2}e^{2\pi is}\right)+\left(\frac{1}{-4\pi^2s^2}e^{-2\pi is} -\frac{1}{-4\pi^2s^2}\right)\\
&=\frac{-2+cos(2\pi s)+isin(2\pi s)+cos(-2\pi s)+isin(-2\pi s)}{-4\pi^2s^2}\\
&=\frac{-2+2cos(2\pi s)}{-4\pi^2s^2}\\
&=\frac{-4cos^2(\pi s)}{-4\pi^2s^2}\\
&=\frac{cos^2(\pi s)}{(\pi s)^2}\\
&=sinc^2s
\end{align*}

时间: 2024-10-13 01:33:03

[傅里叶变换及其应用学习笔记] 六. 热方程讨论的相关文章

[傅里叶变换及其应用学习笔记] 五. 傅里叶级数连续性讨论,热方程

这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用. 热方程后续 上节课推导出热方程的傅里叶系数: $C_k(t) = C_k(0)e^{-2\pi ^2 k^2t}$ 那么$C_k(0)$是什么? 上节课有提到温度有如下关系式: $U(x,t) = \displaystyle{\sum_{k=-\infty}^{\infty}C_k(t)e^{2\pi ikx} }$ 当$t=0$,代表初始时刻圆环上的温度分布 $f(x) = U(x,0) = \displayst

python之raw_input()(学习笔记六)

python之raw_input()(学习笔记六) 我们经常使用raw_input()读取用户的输入,如下例子所示: >>> name = raw_input('please input your name:'),截图如下: 下面简单说下,raw_input()与if搭配使用,脚本如下: #!/usr/bin/env python # -*- coding:utf-8 -*- birth = raw_input('birth:') if birth < 2000: print '0

swift学习笔记(六)析构过程和使用闭包对属性进行默认值赋值

一.通过闭包和函数实现属性的默认值 当某个存储属性的默认值需要定制时,可以通过闭包或全局函数来为其提供定制的默认值. 注:全局函数结构体和枚举使用关键字static标注    函数则使用class关键字标注 当对一个属性使用闭包函数进行赋值时,每当此属性所述的类型被创建实例时,对应的闭包或函数会被调用,而他们的返回值会被作为属性的默认值. ESC: Class SomeCLass{ let someProperty:SomeType={ //给someProperty赋一个默认值 //返回一个与

java之jvm学习笔记六-十二(实践写自己的安全管理器)(jar包的代码认证和签名) (实践对jar包的代码签名) (策略文件)(策略和保护域) (访问控制器) (访问控制器的栈校验机制) (jvm基本结构)

java之jvm学习笔记六(实践写自己的安全管理器) 安全管理器SecurityManager里设计的内容实在是非常的庞大,它的核心方法就是checkPerssiom这个方法里又调用 AccessController的checkPerssiom方法,访问控制器AccessController的栈检查机制又遍历整个 PerssiomCollection来判断具体拥有什么权限一旦发现栈中一个权限不允许的时候抛出异常否则简单的返回,这个过程实际上比我的描述要复杂 得多,这里我只是简单的一句带过,因为这

初探swift语言的学习笔记六(ARC-自动引用计数,内存管理)

Swift使用自动引用计数(ARC)来管理应用程序的内存使用.这表示内存管理已经是Swift的一部分,在大多数情况下,你并不需要考虑内存的管理.当实例并不再被需要时,ARC会自动释放这些实例所使用的内存. 另外需要注意的: 引用计数仅仅作用于类实例上.结构和枚举是值类型,而非引用类型,所以不能被引用存储和传递. swift的ARC工作过程 每当创建一个类的实例,ARC分配一个内存块来存储这个实例的信息,包含了类型信息和实例的属性值信息. 另外当实例不再被使用时,ARC会释放实例所占用的内存,这些

Linux System Programming 学习笔记(六) 进程调度

1. 进程调度 the process scheduler is the component of a kernel that selects which process to run next. 进程调度器需要使 处理器使用率最大化,并且提供 使多个进程并发执行的虚拟 Deciding which processes run, when, and for how long is the process scheduler's fundamental responsibility. 时间片:th

Lua学习笔记(六):函数-续

Lua中的函数是带有词法定界(lexical scoping)的第一类值(first-class values).第一类值指:在Lua中函数和其他值(数值.字符串)一样,函数可以被存放在变量中,也可以存放在表中,可以作为函数的参数,还可以作为函数的返回值.词法定界指:嵌套的函数可以访问他外部函数中的变量.这一特性给Lua提供了强大的编程能力. Lua中关于函数稍微难以理解的是函数也可以没有名字,匿名的.当我们提到函数名(比如print),实际上是说一个指向函数的变量,像持有其他类型的变量一样:

laravel3学习笔记(六)

原作者博客:ieqi.net ==================================================================================================== ORM Laravel3中MVC体系中Model里最重要的组成部分无疑是ORM了,ORM — object-relational mapper — 将数据操作面向对象化,使得整个web框架的核心风格统一,降低整体复杂度,为开发者提供便利. Laravel3中的ORM叫

[傅里叶变换及其应用学习笔记] 九. 继续卷积的讨论

这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用. 卷积在滤波中的应用 浑浊度(Turbidity)研究是关于测量水的清澈度的研究.大致方法是把光传感器放置到深水区域,然后测量光线的昏暗程度,测量出来的值将随时间变化. (由于没有真实数据,下面用mathematica比较粗糙地模拟水域的浑浊度数据)         能看到信号主要集中在低频,我们需要把毛刺去除,也就是把高频去除,在频域进行低通滤波(Low Pass Filtering)         滤波后的波形