最美情侣中文字幕电影,在线麻豆精品传媒,在线网站高清黄,久久黄色视频

歡迎光臨散文網(wǎng) 會員登陸 & 注冊

分離變數(shù)法||數(shù)理方程

2021-05-08 10:19 作者:湮滅的末影狐  | 我要投稿

//上一章還有一點點內(nèi)容需要收尾。

//這一篇基本上是關(guān)于分離變數(shù)法這一章的內(nèi)容。

//前面有一次專欄做到最后他提示我圖片超過一百張投不了,我這才了解到那些公式都是以特殊圖片形式插入的。因此這一次為了盡量避免這一問題有一些簡短公式就不專門用公式格式打了,略有不規(guī)范請諒解...

我們考慮以下問題,并以此為例介紹一些微分方程解法。

考慮長度為l的弦上的波動方程

u_%7Btt%7D-a%5E2%20u_%7Bxx%7D%3D0

有兩端固定邊界條件u(0%2Ct)%3Du(l%2Ct)%3D0

和初始條件u(x%2C0)%3D%5Cphi(x)%2C%20u_t(x%2C0)%3D%5Cpsi(x).

I 達(dá)朗貝爾公式

我們在上一章中對于相關(guān)方程的分析提示我們,可以考慮變量代換:

%5Cleft%5C%7B%0A%5Cbegin%7Barray%7D%7Bc%7D%0A%5Cxi%3Dx%2Bat%5C%5C%0A%5Ceta%3Dx-at%0A%5Cend%7Barray%7D%5Cright.

并且我們也證明了微分方程有通解u%3Df(x%2Bat)%2Bg(x-at)

這里我們直接代入初始條件得到:

%5Cleft%20%5C%7B%5Cbegin%7Balign%7D%0Af(x)%2Bg(x)%3D%5Cphi%20(x)%5C%5C%0Aaf'(x)-bg'(x)%3D%5Cpsi%20(x)%0A%5Cend%7Balign%7D%5Cright%20.

后一式積分,結(jié)合前式即可證明:

u%3D%5Cfrac12%5B%5Cphi(x%2Bat)%2B%5Cphi(x-at)%5D%2B%5Cfrac1%7B2a%7D%5Cint_%7Bx-at%7D%5E%7Bx%2Bat%7D%5Cpsi(x)%5Cmathrm%20d%20x

上式被稱為達(dá)朗貝爾公式。

II 分離變數(shù)法

事實上,我們知道固定端點會引起波的反射,這讓我們聯(lián)想到駐波。猜駐波解:

u(x%2Ct)%3DX(x)T(t)

代入微分方程得

XT''-a%5E2X''T%3D0

%5CRightarrow%20%5Cfrac%7BT%5E%7B%5Cprime%20%5Cprime%7D%7D%7Ba%5E%7B2%7D%20T%7D%3D%5Cfrac%7BX%5E%7B%5Cprime%20%5Cprime%7D%7D%7BX%7D%3D-%5Clambda

一邊關(guān)于時間,一邊關(guān)于空間,總是相等,因此必為常數(shù)%5Clambda.

%5CRightarrow%20X''%2B%5Clambda%20X%3D0

又有邊界條件:X(0)%3DX(l)%3D0

因為任何指數(shù)發(fā)散的解都是無意義的,這里只能%5Clambda%3E0,得三角函數(shù)解:

X_k(x)%3DC_k%5Csin%20%5Cfrac%7Bk%5Cpi%20x%7D%7Bl%7D%20

%5Clambda%3D%5Cfrac%7Bk%5E2%5Cpi%5E2%7D%7Bl%5E2%7D,這種特定取值稱為本征值,以上對X的求解為本征值問題。

同理不難解得

T(t)%3DA%20%5Ccos%20%5Cfrac%7Bk%20%5Cpi%20a%20t%7D%7Bl%7D%2BB%20%5Csin%20%5Cfrac%7Bk%20%5Cpi%20a%20t%7D%7Bl%7D

u_k(x%2Ct)%3D(A_k%5Ccos%20%5Cfrac%7Bk%5Cpi%20at%7D%7Bl%7D%20%2BB_k%5Csin%20%5Cfrac%7Bk%5Cpi%20at%7D%7Bl%7D)%5Csin%5Cfrac%7Bk%5Cpi%20x%7D%7Bl%7D%20

以上就是所有可能的駐波,每一個k值均對應(yīng)一個一種駐波,或稱本征振動。k=1稱為基波,而其他為k次諧波。

分離變數(shù)法原理

考慮到微分方程是線性的,各個本征值對應(yīng)的解可以疊加起來:

u(x%2Ct)%3D%5Csum_%7Bk%3D1%7D%5E%7B%5Cinfty%20%7D%20(A_k%5Ccos%20%5Cfrac%7Bk%5Cpi%20at%7D%7Bl%7D%20%2BB_k%5Csin%20%5Cfrac%7Bk%5Cpi%20at%7D%7Bl%7D)%5Csin%5Cfrac%7Bk%5Cpi%20x%7D%7Bl%7D%20%20

這個形式讓我們想到傅里葉級數(shù)。事實上,將初值條件:

u(x%2C0)%3D%5Cphi(x)%2Cu_t(x%2C0)%3D%5Cpsi(x)

代入,即發(fā)現(xiàn)它們分別展開為傅里葉正弦級數(shù),即可求得A_k%2CB_k.

%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%0AA_%7Bk%7D%3D%5Cfrac%7B2%7D%7Bl%7D%20%5Cint_%7B0%7D%5E%7Bl%7D%20%5Cphi(%5Cxi)%20%5Csin%20%5Cfrac%7Bk%20%5Cpi%20%5Cxi%7D%7Bl%7D%20%5Cmathrm%7B~d%7D%20%5Cxi%20%5C%5C%0AB_%7Bk%7D%3D%5Cfrac%7B2%7D%7Bk%20%5Cpi%20a%7D%20%5Cint_%7B0%7D%5E%7Bl%7D%20%5Cpsi(%5Cxi)%20%5Csin%20%5Cfrac%7Bk%20%5Cpi%20%5Cxi%7D%7Bl%7D%20%5Cmathrm%7B~d%7D%20%5Cxi%0A%5Cend%7Barray%7D%5Cright.

III 極坐標(biāo)系的處理

勻強(qiáng)電場中,放置垂直導(dǎo)體圓柱。自由空間的靜電場滿足拉普拉斯方程:

%5Cnabla%5E2%5Cphi%3D0%5C%3B%20(r%3ER)

邊界條件:(假設(shè)勻強(qiáng)電場沿x軸)

%5Cphi(R%2C%5Ctheta)%3D0%2C%20%5C%3B%5Clim_%7Br%5Crightarrow%20%5Cinfty%7D%5Cphi%5Csim%20-E_0r%5Ccos%5Cphi

我們?nèi)匀环诸愖兞?img type="latex" class="latex" src="http://api.bilibili.com/x/web-frontend/mathjax/tex?formula=%5Cphi%3DR(r)%5CPhi(%5Ctheta%20%20)" alt="%5Cphi%3DR(r)%5CPhi(%5Ctheta%20%20)">

極坐標(biāo)中的拉普拉斯方程為%5Cfrac1r%5Cfrac%7B%5Cpartial%20%7D%7B%5Cpartial%20r%7D%20(r%5Cfrac%7B%5Cpartial%20%5Cphi%20%7D%7B%5Cpartial%20r%7D)%2B%5Cfrac1%7Br%5E2%7D%5Cfrac%7B%5Cpartial%5E2%20%5Cphi%20%7D%7B%5Cpartial%20%5Ctheta%20%5E2%7D%3D0

%5CRightarrow%20%5Cfrac%7Br%7D%7BR%7D%5Cfrac%7B%5Cmathrm%7Bd%7D%20%7D%7B%5Cmathrm%7Bd%7D%20r%7D(rR')%20%20%3D-%5Cfrac%7B%5CPhi''%7D%7B%5CPhi%7D%3D%5Clambda

%5CPhi''%2B%5Clambda%5CPhi%3D0

我們注意到電勢必然是單值函數(shù),這意味著存在自然邊界條件:%5CPhi(%5Ctheta%2B2%5Cpi)%3D%5CPhi(%5Ctheta)

以上構(gòu)成了本征值問題,從而

%5CPhi(%5Ctheta)%3DA_m%5Ccos%20m%5Ctheta%2BB_m%5Csin%20m%5Ctheta%2C%5C%3B%5Clambda%3Dm%5E2%2C%5C%3Bm%3D0%2C1%2C2%2C...

把本征值m2代入關(guān)于R的常微方程

r%5E2R''%2BrR'-m%5E2R%3D0

這是歐拉型常微方程,令r%3De%5Et可得到%5Cfrac%7B%5Cmathrm%7Bd%7D%5E2R%20%7D%7B%5Cmathrm%7Bd%7D%20t%5E2%7D-m%5E2R%3D0

%5CRightarrow%20...%5CRightarrow%20R(r)%3D%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bc%7D%0AC%2BD%5Cln%20r%20%2C%20%5C%3Bm%3D0%5C%5C%0AC%20r%5Em%2BDr%5E%7B-m%7D%2C%5C%3Bm%3D1%2C2%2C...%0A%5Cend%7Barray%7D%5Cright.

對結(jié)果進(jìn)行整理,我們得到極坐標(biāo)下的靜電場通解:

%5Cphi(r%2C%5Ctheta)%3DC_0%2BD_0%5Cln%20r%20%2B%5Csum_%7Bm%3D1%7D%5E%5Cinfty%20%5B%5Crho%5Em(A_m%5Ccos%20m%5Ctheta%2BB_m%5Csin%20m%5Ctheta)%2B%5Crho%5E%7B-m%7D(C_m%5Ccos%20m%5Ctheta%2BD_m%5Csin%20m%5Ctheta)%5D

非齊次方程的處理

例如,初始條件不變,范定方程改為

u_%7Bt%20t%7D-a%5E%7B2%7D%20u_%7Bx%20x%7D%3DA%20%5Ccos%20%5Cfrac%7B%5Cpi%20x%7D%7Bl%7D%20%5Csin%20%5Comega%20t

而邊界條件改為u_x%7C_%7Bx%3D0%7D%3Du_x%7C_%7Bx%3Dl%7D%3D0?(即將u改為對x的偏導(dǎo))

我們考慮將解展開成傅里葉級數(shù),在當(dāng)前邊界條件下應(yīng)考慮余弦級數(shù):

u(x%2Ct)%3D%5Csum_%7Bk%3D0%7D%5E%5Cinfty%20T_k(t)%5Ccos%7B%5Cfrac%7Bk%5Cpi%20x%20%7D%7Bl%7D%7D

代入范定方程,會得到

%5Csum_%7Bk%3D0%7D%5E%7B%5Cinfty%7D%5Cleft%5BT_%7Bk%7D%5E%7B%5Cprime%20%5Cprime%7D%2B%5Cfrac%7Bk%5E%7B2%7D%20%5Cpi%5E%7B2%7D%20a%5E%7B2%7D%7D%7Bl%5E%7B2%7D%7D%20T_%7Bk%7D%5Cright%5D%20%5Ccos%20%5Cfrac%7Bk%20%5Cpi%20x%7D%7Bl%7D%3DA%20%5Ccos%20%5Cfrac%7B%5Cpi%20x%7D%7Bl%7D%20%5Csin%20%5Comega%20t

而等式右邊恰好也是x的傅里葉余弦級數(shù)。比較系數(shù)得

T_k''%2B%5Cfrac%7Bk%5E2%5Cpi%5E2a%5E2%7D%7Bl%5E2%7DT_k%3D%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%0A0%2Ck%5Cneq1%5C%5C%0AA%5Csin%5Comega%20t%2C%20k%3D1%0A%5Cend%7Barray%7D%5Cright.

可以看到,非齊次項僅影響了T1.

我們再將初始條件也展開為傅里葉級數(shù),就可以解得最終解了。

...先更這么多吧...最近太忙了...


分離變數(shù)法||數(shù)理方程的評論 (共 條)

分享到微博請遵守國家法律
台湾省| 金华市| 兴国县| 疏勒县| 墨玉县| 佛坪县| 岑巩县| 湘阴县| 方山县| 岳池县| 济宁市| 柳江县| 饶河县| 平潭县| 林芝县| 鹿邑县| 榆社县| 贡嘎县| 铜山县| 潞城市| 米易县| 平阴县| 双城市| 西城区| 抚州市| 屏东县| 南岸区| 贡嘎县| 津市市| 新密市| 高平市| 渑池县| 徐水县| 南乐县| 英德市| 西和县| 普兰县| 商南县| 宿迁市| 武汉市| 汉阴县|