籍工场分享 http://blog.sciencenet.cn/u/jlpemail

博文

卷积公式的介绍

已有 10837 次阅读 2016-9-21 11:30 |个人分类:读书笔记|系统分类:科研笔记| 公开课, 麻省理工

http://open.163.com/movie/2006/1/C/M/M6TUO44DQ_M6TUPUBCM.html

麻省理工的公开课。


不同于一般常微分方程课程千篇一律地从分离变量和一阶线性方程讲起,MIT《微分方程》第一讲就以独特的视角从全局的角度诠释了微分方程的内涵。课程从方向场和积分曲线入手,深入透彻地剖析了微分方程的实质。一上来,撇开那些有解的特殊的微分方程不谈,却从几何方向通俗易懂,而又全面深入地告诉我们什么是微分方程,解微分方程其实是什么。

[第2课]欧拉数值方法及推广

老头爽约了,他没有按之前说的,讲线性方程的解法,而是开始讲数值方法。按他自己的话说:“线性方程还是推迟到下一讲吧,多数微分方程都是通过数值方法解出来的,先讲这个更好”。他还说:“现在已经是二十一世纪了,计算机都能帮你搞定”。听了他的课才领略,数学不只是那几个臭公式,更重要的是应用。听了他的课,让人深刻地意识到,计算机和数学之间的联系如此紧密。

[第3课]一阶线性常微分方程解法

这一讲的主要内容是一阶线性ODE:y'+p(x)y=q(x),及其解法积分因子法。这一讲通过两个实际问题——“热传导问题”和“溶液浓度扩散问题”,引出了ODE中“最重要”的一节线性微分方程,并透彻讲解。

[第4课]一阶方程代换法

这一讲介绍换元法(或译作代换法,substitution method),并以此为思想将某些特定形式的一阶方程转化为可分离变量方程或线性方程。本讲用换元法解决了两类特定的一阶方程,即伯努利方程和齐次方程。伯努利方程y'=p(x)y+q(x)yⁿ,通过换元化为可分离变量方程。齐次方程y'=F(y/x),令z=y/x可化为线性方程。

[第5课]一阶自治微分方程

这一讲的主题是一阶自治方程y'=f(y)。这一讲不涉及到此类方程的解法,转而考虑在不求解方程的前提下,进行定性分析,直观地获得方程的相关信息,从而避免了由于积分复杂造成不必要的无用功。这一讲还详细讲解了自治方程的一些实用模型:银行存款模型、人口增长模型。

[第6课]复数及复指数

复数在ODE中应用相当广泛。这一讲从复数的运算着手,落脚于复数的极坐标形式。围绕欧拉公式e^iθ=cosθ+isinθ展开,从各个方面详细介绍了这种美妙形式的由来。这一讲还利用复指数巧妙地解决了∫e^x(sinx)dx这种指数、三角函数混合型积分,方法效率远大于常规的分部积分法。

[第7课]一阶常系数线性方程

这一讲特别介绍了一阶常系数线性方程y'+ky=q(t),并解释了k>0时稳态和暂态的内涵。特别地,这一讲强调了y'+ky=kq(t)形式的方程及在相应模型中的应用,并引入输入-响应的概念。最后以正弦波输入作为例子,讲解了分析和求解此类方程的复方法。

[第8课]一阶常系数线性方程(续)

这一讲继续强调一阶常系数线性方程和复数思想。特别强调了正弦输入的情况,并巧妙地通过向量法和复数法给出了三角恒等式acosθ+bsinθ=Ccos(θ-φ)的证明。这一讲的最后,用温度、混合、RC电路、衰变和增长等多个模型为一阶常系数线性方程画上了完美的句号。

[第9课]二阶常系数线性方程

这一讲的主题是二阶常系数齐次线性ODE:y''+Ay'+By=0。这种方程在实际中对应弹簧-质量-阻尼系统,其一般性解法是代入e^(rt),然后通过特征方程r²+Ar+B=0求出r。根据特征方程根的性质,分为两个不同实根、二重实根和复根三种情况,分别对应过阻尼、临界阻尼和欠阻尼三种情况。

[第10课]二阶常系数线性方程(续)

这一讲首先深入讲解了二阶常系数齐次线性常微分方程y''+Ay'+By=0的解如何在实解和复解之间进行转换。然后将方程化为具有物理意义的形式的振动方程y''+2py'+ω²y=0,分别讨论了无阻尼情形(p=0)时解的性质和意义,以及阻尼情况下解的性质和振动的情况。

[第11课]二阶齐次线性方程

这一讲的讨论对象是二阶齐次线性方程y''+p(x)y'+q(x)y=0,讨论了其通解的性质,为何用两个线性独立的解就能表示所有解,而且所有解都在通解的集合内。并解释了叠加原理、唯一性定理、朗斯基行列式等概念。

[第12课]二阶非齐次方程

这一讲的重点是二阶非齐次线性方程y''+p(x)y'+q(x)y=f(x)。首先是将f(x)看成输入或驱动,用弹簧和电路两个例子强调方程的重要性。然后用线性算子,描述了解的一般形式和结构。这一讲的另一个重点是暂态和稳态,在什么条件下对二阶线性方程成立,教授用一句精辟的结论总结了这个问题。

[第13课]非齐次方程特解求法

本讲用算子方法求解高阶非齐次线性方程p(D)y=e^(αx),α为复数,p(D)为D的多项式。考虑p(α)≠0时,特解为e^(αx)/p(α)[用到了代换法则];p(α)=0时,需要分情况讨论,其中单根时,特解xe^(αx)/p'(α)[用到指数位移法则]。

[第14课]案例解读:共振

这一讲是关于共振的。为什么输入频率等于固有频率时,振幅会达到最大?教授从微分方程和数学的角度解释了这个问题。之后教授讲解了带阻尼情况下的"共振",考虑了输入频率和阻尼伪频率之间什么关系时,才能实现这种"共振"。

[第15课]傅里叶级数简介

傅里叶级数在数论、组合数学、信号处理、概率论、统计学、密码学、声学、光学等领域都有着广泛的应用。这一讲首先介绍以2π为周期的函数f(t)可以写作c0+∑(ancosnt+bnsinnt)的傅里叶无穷级数形式。教授通过三角函数正交关系的证明,给出了an和bn的表达式。

[第16课]傅里叶级数简介(续)

这一讲是上一讲的续集,首先考虑了奇函数和偶函数两种情况,讲解了傅里叶级数在这些情况下如何简化运算(以及如果将积分简化到半个周期内)。然后将2π周期延伸到了任意周期2L的情况。最后课程介绍了非周期函数的延伸,任意有限区间都可以用到傅里叶级数。特别地,教授还讲到了傅里叶级数和泰勒级数着眼点的异同。

[第17课]通过傅里叶级数求特解

这一讲主题是利用傅里叶级数求x''+ω0²x=f(t)的特解,其中f(t)化为傅里叶级数,通过sin和cos的可解性来求特解。这一讲采用了方波的例子,告诉我们方程的输入响应系统是如何自然选出与固有频率最接近的共振项的,并以此简单介绍了人耳识别乐音的机理。

[第19课]拉普拉斯变换简介

记得幂级数吧,如1/(1-x)=Σ(x^n)、e^x=∑(x^n/n!),考虑某种变换,让两个幂级数的系数1和1/n!分别对应于f(x)=1/(1-x)或f(x)=e^x,这很容易。其实拉普拉斯变换与这是对应的。教授用这种深入浅出的讲解,让我们了解了拉普拉斯变换的由来。然后分别计算了1、e^at、cos(at)等几种常见函数的变换,并讲解了指数位移的重要公式。大名鼎鼎的拉普拉斯变换,其实并不难。

[第20课]利用拉普拉斯变换求解线性常微分方程

这一讲的主要目标是用拉氏变换求解线性ODE,特别的,解y''+py'+qy=f(t)形式方程。为此,教授首先引入导数的拉氏变换公式,即已知y(t)经过拉氏变换得到Y(t),那么y'及y''如何用Y(t)来表示。拉氏变换解法也就是方程两边同时进行拉氏变换,然后求解得到的代数方程,之后运用部分分式,最后用拉氏逆变换求出解y(t)。

播放中

[第21课]卷积公式

这一讲引入了卷积公式f(t)*g(t)=∫f(u)g(t-u)du。教授从两个方面介绍了卷积的由来和用途:理论方面,卷积和拉氏变换密切相关,L(f)L(g)=L(f*g),卷积由拉氏变换乘积关系的自然产生;实践方面,卷积最普遍的例子是用作放射物质倾泻的积累量问题。教授另外还举了三个实际例子。这一讲全面剖析了卷积公式,并做到了真正的深入浅出。

[第22课]利用拉普拉斯变换求解非连续输入ODE

这一讲主要是讲跳跃式不连续函数u(t)=1(t>0); 0(t<0)的情况,重新定义拉普拉斯逆变换的唯一性,即L(u(t))=1/s。之后教授讲到了函数平移之后的拉普拉斯变换如何进行,之后推广到更一般的不连续输入问题。最后教授以几个实用的例题作结。

[第24课]-一阶常微分方程组简介

这是一阶方程组的第一讲,首先引入了形如x'=f(x,y,t);y'=g(x,y,t)的一阶方程组。教授讲了一些实际用到一阶方程组的例子,然后利用煮鸡蛋的例子,演示了如何用比较直观的消元法来求解。最后教授给出了速度场的几何解释。

[第26课]常系数齐次线性方程组(续)

这一讲继续以矩阵形式x'=Ax讨论常系数齐次线性方程组。课堂上引入了重复实特征值和复特征值两种特殊情况,即特征方程解出重根或复根的情况,两种情况教授分别举出一个实际例子进行讨论。一个是鱼缸温度传递的例子,一个是苏飞传中的爱情例子,引起满堂哄笑。

[第27课]2x2齐次线性方程组作图

这一讲教授讲到了2x2常系数齐次线性方程组各种情况的图像,以此希望给学生一个比较直观的感受,此类方程组解是什么样子。为此,教授引入了两州旅游竞争模型,分别就特征方程中存在两负实根、一正一负实根、以及复根的三种情况给出了方程组解的草图。

[第28课]非齐次方程组矩阵方法

这一讲过渡到非齐次方程组,还是以2x2常系数方程组为例,以矩阵形式x'=Ax+r进行讲解。首先,教授介绍了两个相关定理,为求解做了铺垫。然后介绍了x'=Ax的基本矩阵X。最后通过参数变分的方法,给出了非齐次方程组的特解xp=X∫X^(-1)rdt。

[第29课]矩阵指数

这一讲给出了齐次微分方程组x'=Ax的解的一般公式,即用矩阵指数e^(At)表示基本矩阵X。同单个微分方程x'=ax中,a可以看作是1x1矩阵,其解是e^at。这里就是方程组在nxn矩阵上的推广,以此引入矩阵指数及其在解方程组中的应用。

[第30课]常系数解耦线性方程组

这一讲给出了齐次线性微分方程组x'=Ax的解耦解法,这是第三种方法。由于在自科和工程领域,方程组通常具有物理意义,解耦解法能偶提供对解更为本质的认识,因此教授将其作为这一讲的主题。首先是一个实际例子,然后是一般方程组的解法。

[第31课]非线性自治系统

这一讲介绍非线性的情况,主要是通过轻微阻尼的非线性摆的例子,介绍了该情况下如何求临界点,并作轨迹草图。简谐振动中,摆使用的是小角近似为线性情况,这一讲是一个推广,摆使用的不一定是小角,不过仍然通过线性化得到解释。

[第32课]极限环

这一讲的主题是极限环,首先教授给出了极限环的定义,它首先是方程组的解形成的一条闭合轨迹,另外它不同于一般闭合轨迹,它必须是附近轨迹在t趋于无穷时逼近的轨迹。然后教授介绍了极限环何时不存在的两个准则,分别是本迪克松准则和临界点准则,证明本迪克松准则时,证明过程中涉及了反证法,以及逆否命题逻辑。最后教授介绍了极限环的一些历史,并用他经历的一个有趣故事结束了本课,与某位中国教授有关。




https://blog.sciencenet.cn/blog-255-1004170.html

上一篇:看油画写诗句
下一篇:约稿之痛
收藏 IP: 124.207.244.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (1 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-27 03:08

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部