波动方程的分离变量法

波动方程初边值问题

首先给出波动方程初边值问题的格式：（记为方程A）
$u_{tt} - a^2u_{xx} = f(x,t)\;\;(1)$ $u(x,0) = \varphi (x),u_t(x,0) = \psi(x)\;\;(2)$ $u(0,t) = 0,u(l,t) = 0\;\;(3)$ 其中，（1）为运动方程，（2）为初始条件，（3）为边界条件（齐次形式）。

利用叠加原理分成两个问题

方程B
$u_{tt} - a^2u_{xx} = 0\;\;(4)$ $u(x,0) = \varphi (x),u_t(x,0) = \psi(x)\;\;(5)$ $u(0,t) = 0,u(l,t) = 0\;\;(6)$ 方程C
$u_{tt} - a^2u_{xx} = f(x,t)\;\;(7)$ $u(x,0) = 0,u_t(x,0) = 0\;\;(8)$ $u(0,t) = 0,u(l,t) = 0\;\;(9)$ 由叠加原理，若方程B和方程C分别有解 $u_1,u_2$ 那么方程A的解 $u=u_1+u_2$ 。过程略。由齐次化原理我们可以将方程C转化为方程B的形式求解，因此在此我们只讨论方程B的解。

分离变量法（驻波法）

简要的说明一下驻波法的数学推导步骤

构造满足运动方程（4）和边值条件（6）特解（单音） $u(x,t)=X(x)T(t)$
先代入运动方程（4）得到T和X满足的带系数 $\lambda$ 的常微分方程
将X的通解代入边值条件（6），得到 $\lambda$ 的取值（一列 $\{\lambda(k)\}$ ）
将 $\{\lambda(k)\}$ 代入T的常微分方程，求出T的解
将X和T的解（有系数 $A_kB_k$ ）代入 $u_k(x,t)=X_k(x)T_k(t)$
通解 $u$ 为 $u_k$ 的线性组合（此处需要边值条件（6）必须是齐次的！也就是说驻波法要求齐次边值条件）
要求通解 $u$ 满足初值条件（5），刚好满足傅里叶展开的形式。由傅里叶展开“正交系”的性质求出系数 $A_kB_k$
至此，通解 $u$ 已求出。

step1变量分离的构造（驻波法的物理背景）

想法：由物理知识（还有傅里叶展开），我们想到

物理背景：弦振动所发出的声音可以分解为各种不同频率的单音的叠加，相对于每种单音，弦振动时波形保持不变，从而当时间变化时各点的振幅作同步的变化，也就是说，每种单音都是具有形式为 $u(x,t) =X(x)T(t)\;\;(10)$ 形式的特殊解，而整个复杂振动过程可以通过这种特殊解的叠加得到。

step2代入运动方程（4）

将（10）代入（4），得到
$X(x)T''(t) - a^2X''(x)T(t) = 0$ 变量分离，得到
$\frac{T''(t)}{a^2T(t)}=\frac{X''(x)}{X(x)}=-\lambda$ 其中 $\lambda$ 是特征值。
$\left\{\begin{matrix} T''(t) + \lambda a^2T(t)=0\;\;(11) \\ X''(x) + \lambda X(x)=0\;\;(12) \end{matrix}\right.$

step3将X的通解代入边值条件（6）

（ $\lambda <=0$ 的计算过程与下面类似，但只能得到0解，不是我们想要的。在此省略过程）
当 $\lambda >0$ ，（12）的通解为
$X(x) = C_1 cos(\sqrt{\lambda} x)+ C_2 sin(\sqrt{\lambda} x)$ 由边值条件（6） $X(0)=0$ ，有 $C_1=0$
再由 $X(l) = 0$ ，有 $X(l) = C_2 sin(\sqrt{\lambda} l) = 0$
若要 $C_2$ 不为零，必须有 $\sqrt{\lambda}l = k\pi,(k = 1,2, \dots )$
即 $\lambda = k^2\pi^2/l^2,(k = 1,2, \dots )$

step4将 $\{\lambda(k)\}$ 代入T的常微分方程，求出T的解

$X_k = C_k sin(k\pi x/l)$ $T_k = A_k cos(k\pi at/l) + B_k sin(k\pi at/l)$

step5将X和T的解（有系数 $A_kB_k$ ）代入 $u_k(x,t)=X_k(x)T_k(t)$ $u_k(x,t)=X_k(x)T_k(t) = (A_k cos(k\pi at/l) + B_k sin(k\pi at/l))sin(k\pi x/l)$

由齐次边值条件（6），通解 $u$ 为
$u(x,t) = \sum_{k=1}^\infty (A_k cos(k\pi at/l) + B_k sin(k\pi at/l))sin(k\pi x/l) \;\;(13)$ 对t求导，得
$u_t(x,t) = \sum_{k=1}^\infty k\pi a/l(-A_k sin(k\pi at/l) + B_k cos(k\pi at/l))sin(k\pi x/l) \;\;(14)$

step7要求通解 $u$ 满足初值条件（5）

初值条件（5）： $u(x,0) = \varphi (x),u_t(x,0) = \psi(x)$
由（13）（14），有
$u(x,0) = \sum_{k=1}^\infty A_ksin(k\pi x/l) = \varphi(x)\;\;(15)$ $u_t(x,0) = \sum_{k=1}^\infty (k\pi a/l) B_k sin(k\pi x/l) = \psi(x)\;\;(16)$ 以（15）为例计算：
由 $sin(k\pi x)$ 的“正交系”性质，即
$\left\{\begin{matrix} \int_0^1 sin^2(k\pi x)dx = 1/2 \\ \int_0^1 sin(k_1\pi x)sin(k_2\pi x)dx = 0,\forall k_1 \neq k_2 \end{matrix}\right.$ 因此，将（15）左右同乘 $sin(k\pi x/l)$ 并在 $[0,l]$ 上积分，可得
$A_k l/2= \int_0^l \varphi(x)sin(k\pi x/l)dx$，即$A_k =(2/l) \int_0^l \varphi(x)sin(k\pi x/l)dx$ 同理可得 $B_k =(2/k\pi a) \int_0^l \psi(x)sin(k\pi x/l)dx$ 。
将 $A_kB_k$ 代入（13）可得方程B的解
$u_t(x,t) = \sum_{k=1}^\infty k\pi a/l(-(2/l) \int_0^l \varphi(x)sin(k\pi x/l)dx*sin(k\pi at/l) + (2/k\pi a) \int_0^l \psi(x)sin(k\pi x/l)dx* cos(k\pi at/l))sin(k\pi x/l)$

最后编辑于：2018.12.10 16:16:19

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 215,539评论 6赞 497
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 91,911评论 3赞 391
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 161,337评论 0赞 351
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 57,723评论 1赞 290
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 66,795评论 6赞 388
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 50,762评论 1赞 294
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 39,742评论 3赞 416
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 38,508评论 0赞 271
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 44,954评论 1赞 308
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 37,247评论 2赞 331
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 39,404评论 1赞 345
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 35,104评论 5赞 340
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 40,736评论 3赞 324
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 31,352评论 0赞 21
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 32,557评论 1赞 268
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 47,371评论 2赞 368
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 44,292评论 2赞 352