弦振动方程的解法（分离变量法+行波法）

弦振动方程的解法

微分方程基础

本文所言的方程基本上是微分方程，而非中学阶段所言的代数方程，微分方程的求解目标一般是一个函数，而代数方程的解一般是一个值。下面给出两例来介绍微分方程和代数方程的区别。

下列方程是一个典型的代数方程：

$x^2-3x+2=0$

根据中学阶段所学的知识，可以解出方程的两个代数解：

$x_1=1,\ x_2=2$

下列方程是一个典型的微分方程：

$\frac{{\rm d}^2y}{{\rm d}x^2}-3\frac{{\rm d}y}{{\rm d}x}+2y=0$

根据本科阶段所学的知识，可以解出方程的通解：

$y=c_1{\rm e}^{x}+c_2{\rm e}^{2x}$

但是上述的微分方程是常微分方程，常微分方程只有一个变量。而对于弦振动来说，位置 $x$

和时间 $t$ 都是变量，故而弦振动方程是一种偏微分方程。

弦振动方程介绍

弦振动方程是波动方程的一种，是局限在一维空间内的波动方程。我们假定弦上在位置 $x$ 处的一点，在时间 $t$ 时刻，偏移平衡位置的偏移量是 $u(x,t)$ ，那么偏移的加速度就是 $u$ 关于 $t$ 求两次偏导数 $u_{tt}$ ，而弦的弯曲程度（凹凸性）是 $u$ 关于 $x$ 的二阶偏导数 $u_{xx}$ 。

易于理解，当 $u_{xx}$ 为正数时，意味着弦在这里向下凸，类似于”U“字型。这个地方在两端的拉力下，拉力的左右分力相互抵消，只剩下向上的分力。这个地方应该向上有加速度，所以 $u_{tt}$ 也为正。同样地，当 $u_{xx}$ 为负数时，类似于”∩“字型，加速度为负。

通过物理背景，我们可以列出弦振动方程：

$u_{tt}=a^2u_{xx}$

偏移量 $u$ 的国际单位是 $\rm m$ ，它关于 $t$ 求两次偏导数，所得加速度的单位是 $\rm m/s^2$ ，而右边 $u_xx$ 的单位是 $\rm m^{-1}$ ，这样，我们可以得到 $a$ 的单位是 ${\rm m/s}$ ，实验现象表明，这里的 $a$ 就是波速。

分离变量法解弦振动方程

上文我们列出了弦振动方程的偏微分方程，对于有限长弦振动问题，我们可以尝试用分离变量法来求解弦振动方程。

如果弦长 $L$ ，并且弦两端固定，那么可以得到偏微分方程的两个边界条件：

$u|_{x=0}=0,\ u|_{x=L}=0$

在物理现象中，代表空间的 $x$ 和代表时间的 $t$ 是可以分离的。

假定 $u(x,t)=X(x)T(t)$ ，代入原偏微分方程，得到：

$X(x)T''(t)=a^2X''(x)T(t)$

变形

$\frac{X''}{X}= \frac{T''}{a^2T}$

令这个数值为 $-\lambda$ ，即 $\frac{X''}{X}= \frac{T''}{a^2T}=-\lambda$ ，我们得到了两个常微分方程：

$X''+\lambda X=0$

$T''+\lambda a^2 T=0$

当 $\lambda \le 0$ 时，考虑到边界条件的话，只能得到零解 $X(x)=0$

当 $\lambda > 0$ 时，根据本科学习的知识，可以解出 $X(x)$ 的通解：

$X(x)=a \sin \sqrt{\lambda }x + b \cos \sqrt{\lambda }x$

考虑到边界条件 $u|_{x=0}=0,\ u|_{x=L}=0$ ，即

$X(0)=b =0,\ X(L)=a \sin \sqrt{\lambda}L=0$

可知 $\sqrt{\lambda} L = n\pi,\ n=1,2,...$

我们可以解出 $\lambda _n=n^2\pi^2/L^2$

下面我们再关注 $T''+\lambda a^2 T=0$ ，我们可以解出 $T(t)$ 的通解：

$T(t)=c \sin \sqrt{\lambda}at + d \cos \sqrt{\lambda}at$

实际上，因为 $\lambda$ 可以取不同的值， $T(t)$ 是一系列波的叠加：

$T(t)=\sum_{n=1}^{\infty} (C_n \sin \sqrt{\lambda_n}at + D_n \cos \sqrt{\lambda_n}at)$

故而，我们得到了有限长弦振动方程的一般解：

$u(x,t)=\sum_{n=1}^{\infty} (C_n \sin n\pi at/L + D_n \cos n\pi at/L)\sin n\pi x/L$

如果能给出初始条件，便能确定 $C_n$ 和 $D_n$ 的值，从而得到了有限长弦振动方程的特解。

行波法求解弦振动方程

在无限长的弦上，一切的波动都可以看做一个左行波和一个右行波的叠加。故而可以假设

$u(x,t)=f_1(x-at)+f_2(x+at)$

上面这个方程恰好满足弦振动方程 $u_{xx}=a^2u_{tt}$ ，是无限长弦振动方程的通解。

如果我们有初始条件（弦的初始形状和初始速度）：

$u|_{t=0}=\varphi (x),\ u_t|_{t=0}=\psi (x)$

将通解代入

$f_1(x)+f_2(x)=\varphi (x),\ -af'_1(x)+af'_2(x)=\psi (x)$

对后者求一次不定积分可得

$-f_1(x)+f_2(x)=\frac{1}{a}\int \psi (x){\rm d}x$

可以分别解出 $f_1(x)$ 和 $f_2(x)$ ：

$f_1(x)=\frac{1}{2} \varphi(x) - \frac{1}{2a}\int \psi (x){\rm d}x$

$f_2(x)=\frac{1}{2} \varphi(x) + \frac{1}{2a}\int \psi (x){\rm d}x$

那么， $u(x, t)=f_1(x-at)+f_2(x+at)$

$u(x,t)=\frac 12 (\varphi (x+at)+\varphi (x-at))+\frac 1{2a}\int_{x-at}^{x+at} \psi(x){\rm d} x$

这一结论又称为d'Alembert公式、

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念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