十九世纪的偏微分方程(一)

1 引言

偏微分方程诞生于18世纪,发展于19世纪,随着物理研究在深度和广度上的进展,微分方程在数量和类型上增加了,过去已知的方程如波动方程和位势方程也应用到新的物理领域了。偏微分方程变成数学的中心,一是因为在物理中应用广泛,二是因为偏微分方程的求解促进了函数论、变分法、级数展开、常微分方程、代数、微分几何等方面的发展,本章只能讨论其中的一小部分。

今天人们习惯按类型对偏微分方程分类,但19世纪初对偏微分方程的了解还不足以进行分类,由物理指导应讨论哪种方程,数学家则随意从一种类型的问题转向另一种,忽略了其中的差别(对今天来说是最基本的差别),毕竟物理从过去到现在都不关心数学家的分类。

2 热方程与傅里叶级数

傅里叶(Joseph Fourier,1768-1830)迈出了19世纪第一也是最重要的一步(还记得18世纪偏微分方程提到的欧拉、伯努利、达朗贝尔之争么)。傅里叶年轻时数学很好,但他想当士官,因为出身裁缝家庭,不让当,后来军校让他留校任教,他接受了,然后搞了一辈子数学。

当时流行搞热流动研究,这个研究有很多用途,比如冶炼金属等工业应用,以及确定地球内部温度、研究温度随时间变化等科研课题,1807年,傅里叶向巴黎科学院提交了一篇关于热传导的基本论文,被拉格朗日、拉普拉斯和勒让德评审拒绝了(感觉这仨人喜欢搞天体),不过科学院想鼓励他接着研究,宣布1812年给搞得好这个课题的发高额奖金,于是1811年傅里叶提交了修改过的论文,获得了奖金,但受到了缺乏严密性的批评,未能发表在当时科学院的《报告》里。傅里叶很生气,继续搞热学,在1822年发表了数学的经典文献《热的解析理论》(傅里叶思想的主要出处),把1811年论文的第一部分直接放进去了,两年后他成为科学院秘书,又把1811年论文原封不动地发表在《报告》中(傅里叶:莫欺中年穷ok?话说看了半天数学史难得有个人到中年出成果的,可能因此学界比较歧视中年人的智慧……)。

在吸放热的物体内部,温度分布一般是不均匀的,在任何点上都随时间变化,所以温度T是时间和空间的函数。函数的准确形式依赖于物体形状、密度、材料的比热、T的初始分布(t=0时的温度分布)以及物体表面的边界条件。傅里叶在书中考虑的第一个主要问题是在均匀和各向同性的物体内确定作为x,y,z,t函数的温度T, 根据物理原理T必须满足偏微分方程:(\frac{\partial ^2T}{\partial x^2}+\frac{\partial ^2T}{\partial y^2}+\frac{\partial ^2T}{\partial z^2})=k^2\frac{\partial T}{\partial t} ,称为三维空间的热方程,其中k^2是一个常数,其值依赖于物体的质料。傅里叶解决了特殊的热传导问题,对两端T保持在0度,侧面绝热因而无热流通过的柱轴,求解热方程,这根轴只涉及一维空间,即\frac{\partial ^2T}{\partial x^2}=k^2\frac{\partial T}{\partial t},边界条件T(0,t)=0,T(l,t)=0,t>0,初始条件T(x,0)=f(x),0<x<1,傅里叶使用分离变量法,令T(x,t)=Φ(x)ψ(t),代入微分方程得到Φ''(x)/(k^2*Φ(x))=ψ''(t)/ψ(t),因为x变化时ψ''(t)/ψ(t)不变(是常数),同理t变化时Φ''(x)/Φ(x)也是常数,设Φ''(x)/(k^2*Φ(x))=ψ''(t)/ψ(t)=-λ,有Φ''(x)+λk^2Φ(x)=0和ψ''(t)+λψ(t)=0,把变量分离代入边界条件,有Φ(0)=0和Φ(l)=0,Φ(x)的通解是\phi(x)=bsin(\sqrt{\lambda}kx+c),得到λ的取值(称为本征值or特征值)\lambda_{v}=(\frac{v\pi}{kl})^2,v\in Z ,ψ的通解是指数函数,现在为了满足初始条件,对t=0必须有:f(x)=\sum_{v=1}^{\infty}b_{v}sin\frac{v\pi x}{l}。于是傅里叶面临的问题是f(x)能否表示成三角级数,特别地,bv能否确定。

傅里叶进而回答这些问题(虽然略有严密性问题)。为简单起见,设l=π,他把每个正弦函数按麦克劳林定理展开为幂级数,最终把f(x)表示成x的幂级数。

这隐含了一个假设,即这个幂级数必须是f(x)的麦克劳林级数,因此f(x)=\sum_{k=0}^{\infty} \frac{1}{k!}f^{(k)}(0) x^k,令两个f(x)中x的同次幂的系数相等,傅里叶发现对偶数k,f^{(k)}(0)=0,此外bv是无穷线性代数方程组里的未知数的一个无穷集合。

对f(x)能否表示成三角级数,他面临同类的方程组,取前k项和前k个方程的右端常数,解前k个方程得bv,k表示bv的近似值,得到bv,k的一般表达式后,他大胆地下结论称:b_{v}=\lim_{k\to\infty}b_{v,k} ,但这次确定bv有很多困难,他用很复杂的程序说明几个不同的f(x)如何确定bv,再用这些特殊情形作为指导,得到了bv的含有无穷乘积和无穷和的一个表达式,他觉得没啥用,经过了一些大胆和创造性的步骤得到了公式b_v=\frac{2}{\pi}\int_{0}^{\pi }f(s)\sin vs ds,在某种程度上这不算新结论,之前克莱罗和欧拉已经知道怎么把某些函数展开为傅里叶级数并得到公式:

此外,傅里叶得到的结果是很局限的,因为他假定f(x)有麦克劳林展开,即有无穷阶导数,最后,他的方法不够严密,且比欧拉的方法复杂。傅里叶不得不用无穷线性方程组,而欧拉使用三角函数的性质,更为简单。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,793评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,567评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,342评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,825评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,814评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,680评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,033评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,687评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 42,175评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,668评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,775评论 1 332
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,419评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,020评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,978评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,206评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,092评论 2 351
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,510评论 2 343

推荐阅读更多精彩内容