數(shù)理方程-淺談齊次化原理

筆者剛學(xué)《數(shù)理方程》不久能曾,就齊次化原理談?wù)勛约旱睦斫舛认酰胁粚?duì)的地方歡迎交流批評(píng)指正。

齊次化原理是數(shù)理方程中非常重要的定理寿冕,在各種情況的方程都會(huì)用到蕊程。在此以最簡(jiǎn)單的一維波動(dòng)方程柯西問題為例:(為方便起見,命名為方程A)
(A)\left\{\begin{array}{} \frac{\partial^2u}{\partial t^2}-a^2 \frac{ \partial^2u}{\partial x^2} = f(x,t) ,&(t>0,-\infty<x<+\infty) \\ u=\varphi(x), \frac{ \partial u}{\partial t}=\psi(x),& (t=0,-\infty<x<+\infty) \end{array}\right.

解決本題的第一步是驼唱,由“疊加原理”藻茂,分成一個(gè)f=0的問題和一個(gè)\varphi=\psi=0的問題。第二個(gè)問題可由達(dá)朗貝爾公式直接求解玫恳,第一個(gè)方程使用“齊次化原理”化成第二個(gè)方程的形式辨赐,再由達(dá)朗貝爾公式求解。

下面具體分析第二個(gè)方程的“齊次化原理”部分:(方程B)
(B)\left\{\begin{array}{} \frac{ \partial^2u}{\partial t^2}-a^2 \frac{ \partial^2u}{\partial x^2} = f(x,t),&(t>0,-\infty <x<+\infty) \\ u=0, \frac{ \partial u}{\partial t}=0,& (t=0,-\infty<x<+\infty) \end{array}\right.

首先給出齊次化原理的定義:
構(gòu)造一個(gè)方程組(方程C)
(C)\left\{\begin{array}{} \frac{ \partial^2W}{\partial t^2}-a^2 \frac{ \partial^2W}{\partial x^2} = 0 &(t>\tau,-\infty<x<+\infty) \\ W=0, \frac{ \partial W}{\partial t}=f(x,\tau)& (t=\tau,-\infty<x<+\infty) \end{array}\right.

若它的解是W(x,t;\tau)京办,則原問題的解為
u(x,t) = \int_0^t W(x,t;\tau)d\tau

下面介紹它的意思:

首先是這個(gè)問題的物理模型掀序,畢竟我們是數(shù)理方程課,幾乎所有的數(shù)學(xué)結(jié)論都有物理解釋惭婿,而反過來(lái)也可以從物理現(xiàn)象為數(shù)學(xué)定理提供思路森枪。

本題(方程B)描述的是,一根弦审孽,無(wú)限長(zhǎng)(即不考慮邊界的影響)县袱,且在初始時(shí)刻(t=0),繩子的位置在正中間的平衡位置(u=0)(不考慮重力佑力,只考慮繩子本身的微小振動(dòng))式散,也沒有初始速度(\frac{ \partial u}{\partial t}=0)。但每時(shí)每刻有外力f(x,t)作用于繩子打颤。

首先暴拄,為了簡(jiǎn)單研究,我們將連續(xù)的時(shí)間分成一段一段的编饺,這也是一個(gè)微元法的思想乖篷,先切割成小段離散的研究,最后再取極限變成連續(xù)的情況透且。打個(gè)比方撕蔼,我們研究t\in [0,T]的繩子運(yùn)動(dòng)豁鲤,那我們就將它看成n個(gè)\Delta\taun\Delta\tau=T),每個(gè)\Delta\tau內(nèi)各個(gè)部分受到的力不變鲸沮,這樣就有\int_t^{t+\Delta\tau}f(x,y)dy=f(x,t)\Delta\tau琳骡。

再由物理中的動(dòng)量守恒定律,簡(jiǎn)單的寫就是ft = m\Delta v讼溺,沖量等于動(dòng)量的增量楣号。由這個(gè)物理定律,我們可以把“在一段時(shí)間內(nèi)沒有初速度但受到外力作用導(dǎo)致的沖量”等同與“在一段時(shí)間內(nèi)有初速度但沒有外力導(dǎo)致的沖量”怒坯。按照這個(gè)思路炫狱,我們能寫出方程C并推導(dǎo)出齊次化原理。

下面使用微元法試圖推導(dǎo):(這個(gè)過程是按自己理解寫的剔猿,不一定嚴(yán)謹(jǐn))

我們可以近似的將每\Delta\tau時(shí)間的力f導(dǎo)致的沖量f\Delta\tau變成一個(gè)速度增量f\Delta\tau(這里忽略質(zhì)量m视译,我物理不好我解釋不清楚m了……),就是繩子每過\Delta\tau就獲得一個(gè)速度增量艳馒。再由疊加原理憎亚,我們可以證明在本問題中(方程B)员寇,力對(duì)繩子的作用近似的等于由好多個(gè)時(shí)刻的速度增量的作用疊加近似弄慰。

也就是說(shuō),在我們的例子中蝶锋,可以給出n個(gè)方程:
(D)\left\{\begin{array}{} \frac{ \partial^2W}{\partial t^2}-a^2 \frac{ \partial^2W}{\partial x^2} = 0 ,&(t>\tau_i,-\infty<x<+\infty) \\ W=0, \frac{ \partial W}{\partial t}=f(x,\tau_i)\Delta\tau.& (t=\tau_i,-\infty<x<+\infty) \end{array}\right.\\ \tau_i=i\Delta\tau,i=0,\cdots,n-1

最后他們的解的和等于原問題的解陆爽。當(dāng)i=0,表示[0,\tau_1]的力轉(zhuǎn)化為\tau_0時(shí)刻的初速度扳缕。當(dāng)i=1慌闭,表示[\tau_1,\tau_2]的力轉(zhuǎn)化為\tau_1時(shí)刻的初速度,這條繩子在t\in [0,\tau_1]內(nèi)靜止不動(dòng)躯舔,在\tau_1時(shí)獲得一個(gè)初速度驴剔,之后由方程的第一個(gè)式子描述的運(yùn)動(dòng)開始運(yùn)動(dòng)。后幾個(gè)類似粥庄。

當(dāng)我們解出這n個(gè)方程的解W(x,t;\tau_i)丧失,(求解的時(shí)候需要用一步換元把t>\tau改成t'>0,之后套用達(dá)朗貝爾公式)
W(x,t;\tau_i)=\frac{1}{2a}\int_{x-a(t-\tau_i)}^{x+a(t-\tau_i)}f(\alpha,\tau_i)\Delta\tau \mathrmyvzkjbr\alpha=\Delta\tau\frac{1}{2a}\int_{x-a(t-\tau_i)}^{x+a(t-\tau_i)}f(\alpha,\tau_i) \mathrmqqqqqlf\alpha


W'(x,t;\tau_i)=\frac{W(x,t;\tau_i)}{\Delta\tau}=\frac{1}{2a}\int_{x-a(t-\tau_i)}^{x+a(t-\tau_i)}f(\alpha,\tau_i) \mathrmvkzdbkh\alpha

則原方程的解為
u(x,T) = \sum_{i=0}^{n-1}W(x,T;\tau_i)=\sum_{i=0}^{n-1}W'(x,T;\tau_i)\Delta\tau

并且可以看出W'是以下方程的解
(E)\left\{\begin{array}{} \frac{ \partial^2W}{\partial t^2}-a^2 \frac{ \partial^2W}{\partial x^2} = 0 ,&(t>\tau_i,-\infty<x<+\infty) \\ W=0, \frac{ \partial W}{\partial t}=f(x,\tau_i).& (t=\tau_i,-\infty<x<+\infty) \end{array}\right. \\ \tau_i=i\Delta\tau,i=0,\cdots,n-1

方程E與方程D只有第二行\frac{\partial W}{\partial t}不同惜互,對(duì)應(yīng)方程C布讹。

\Delta\tau\rightarrow0,就得到了積分的形式训堆,也就是齊次化原理:
u(x,T) = \int_0^T W'(x,T;\tau)d\tau

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末描验,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子坑鱼,更是在濱河造成了極大的恐慌膘流,老刑警劉巖,帶你破解...
    沈念sama閱讀 221,820評(píng)論 6 515
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場(chǎng)離奇詭異睡扬,居然都是意外死亡盟蚣,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 94,648評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門卖怜,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)屎开,“玉大人,你說(shuō)我怎么就攤上這事马靠⊙俪椋” “怎么了?”我有些...
    開封第一講書人閱讀 168,324評(píng)論 0 360
  • 文/不壞的土叔 我叫張陵甩鳄,是天一觀的道長(zhǎng)逞度。 經(jīng)常有香客問我,道長(zhǎng)妙啃,這世上最難降的妖魔是什么档泽? 我笑而不...
    開封第一講書人閱讀 59,714評(píng)論 1 297
  • 正文 為了忘掉前任,我火速辦了婚禮揖赴,結(jié)果婚禮上馆匿,老公的妹妹穿的比我還像新娘。我一直安慰自己燥滑,他們只是感情好渐北,可當(dāng)我...
    茶點(diǎn)故事閱讀 68,724評(píng)論 6 397
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著铭拧,像睡著了一般赃蛛。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上搀菩,一...
    開封第一講書人閱讀 52,328評(píng)論 1 310
  • 那天呕臂,我揣著相機(jī)與錄音,去河邊找鬼肪跋。 笑死歧蒋,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的澎嚣。 我是一名探鬼主播疏尿,決...
    沈念sama閱讀 40,897評(píng)論 3 421
  • 文/蒼蘭香墨 我猛地睜開眼,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼易桃!你這毒婦竟也來(lái)了褥琐?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 39,804評(píng)論 0 276
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤晤郑,失蹤者是張志新(化名)和其女友劉穎敌呈,沒想到半個(gè)月后贸宏,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 46,345評(píng)論 1 318
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡磕洪,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,431評(píng)論 3 340
  • 正文 我和宋清朗相戀三年吭练,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片析显。...
    茶點(diǎn)故事閱讀 40,561評(píng)論 1 352
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡鲫咽,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出谷异,到底是詐尸還是另有隱情分尸,我是刑警寧澤,帶...
    沈念sama閱讀 36,238評(píng)論 5 350
  • 正文 年R本政府宣布歹嘹,位于F島的核電站箩绍,受9級(jí)特大地震影響,放射性物質(zhì)發(fā)生泄漏尺上。R本人自食惡果不足惜材蛛,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,928評(píng)論 3 334
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望怎抛。 院中可真熱鬧卑吭,春花似錦、人聲如沸抽诉。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,417評(píng)論 0 24
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)迹淌。三九已至,卻和暖如春己单,著一層夾襖步出監(jiān)牢的瞬間唉窃,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 33,528評(píng)論 1 272
  • 我被黑心中介騙來(lái)泰國(guó)打工纹笼, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留纹份,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 48,983評(píng)論 3 376
  • 正文 我出身青樓廷痘,卻偏偏與公主長(zhǎng)得像蔓涧,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子笋额,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,573評(píng)論 2 359

推薦閱讀更多精彩內(nèi)容