MATLAB中數(shù)學(xué)物理建模與方程分類

1、1數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)第第3章章 數(shù)學(xué)物理建模與方程分類數(shù)學(xué)物理建模與方程分類數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)本章內(nèi)容本章內(nèi)容3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例 波動(dòng)型問(wèn)題的建模 熱傳導(dǎo)型問(wèn)題的建模 穩(wěn)定型問(wèn)題的建模3.3 數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件3.4 二階偏微分方程的分類二階偏微分方程的分類數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)相關(guān)基本概念相關(guān)基本概念 什么是什么是數(shù)學(xué)物理？ 建立與研究物理現(xiàn)象的數(shù)學(xué)模型的理論 連接

2、物理學(xué)與數(shù)學(xué)的一門交叉學(xué)科 什么是什么是數(shù)學(xué)物理建模？ 尋求把物理問(wèn)題和其他自然科學(xué)和技術(shù)科學(xué)轉(zhuǎn)換成數(shù)學(xué)模型的過(guò)程 什么是什么是數(shù)學(xué)物理方程(數(shù)理方程)？ 從物理學(xué)及其它各門自然科學(xué)自然科學(xué)、技術(shù)科學(xué)技術(shù)科學(xué)中所導(dǎo)出的函數(shù)函數(shù)方程方程，主要指偏微分方程偏微分方程和積分方程積分方程 凡是在建立和研究描述物理現(xiàn)象的數(shù)學(xué)模型時(shí)所用到的數(shù)學(xué)方法，都屬于數(shù)學(xué)物理方法的范疇3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理的發(fā)展歷史數(shù)學(xué)物理的發(fā)展歷史 17世紀(jì)末世紀(jì)末創(chuàng)建創(chuàng)建 Newton和Leibniz創(chuàng)立了微積分并用于表述力學(xué)基本定律與

3、萬(wàn)有引力定律 18世紀(jì)世紀(jì)形成和發(fā)展形成和發(fā)展 分析學(xué)穩(wěn)步發(fā)展，并廣泛地應(yīng)用于理論物理學(xué) 19世紀(jì)世紀(jì)長(zhǎng)足進(jìn)展長(zhǎng)足進(jìn)展 數(shù)學(xué)物理方法成功地應(yīng)用于研究與各類物理場(chǎng)和波動(dòng)過(guò)程有關(guān)的物理現(xiàn)象的數(shù)學(xué)模型 20世紀(jì)世紀(jì)現(xiàn)代數(shù)學(xué)物理現(xiàn)代數(shù)學(xué)物理 量子物理和相對(duì)論的理論研究，空氣動(dòng)力學(xué)、粒子轉(zhuǎn)換現(xiàn)象和等離物理中新問(wèn)題的出現(xiàn)，大大擴(kuò)展了數(shù)學(xué)物理所使用的數(shù)學(xué)工具3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理的特點(diǎn)數(shù)學(xué)物理的特點(diǎn) 建立具有共性的數(shù)學(xué)模型建立具有共性的數(shù)學(xué)模型 建立具有歸納和演繹雙重功能的的數(shù)學(xué)模型建立具有歸納和演繹雙重功能的的數(shù)

4、學(xué)模型 結(jié)合已有的物理模型開發(fā)新的數(shù)學(xué)方法結(jié)合已有的物理模型開發(fā)新的數(shù)學(xué)方法 結(jié)合計(jì)算技術(shù)推動(dòng)物理學(xué)的發(fā)展結(jié)合計(jì)算技術(shù)推動(dòng)物理學(xué)的發(fā)展3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模數(shù)學(xué)物理建模 數(shù)學(xué)物理方法解決問(wèn)題的立足點(diǎn)是某物理現(xiàn)象數(shù)學(xué)物理方法解決問(wèn)題的立足點(diǎn)是某物理現(xiàn)象的內(nèi)在物理規(guī)律的內(nèi)在物理規(guī)律vgF=ma- -+ +U=Ed3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模數(shù)學(xué)物理建模 在獲得物理規(guī)律的情況下可以根據(jù)現(xiàn)象歸納其在獲得物理規(guī)律的情

5、況下可以根據(jù)現(xiàn)象歸納其具體的物理過(guò)程具體的物理過(guò)程 也可以根據(jù)具體的物理過(guò)程演繹將出現(xiàn)的現(xiàn)象也可以根據(jù)具體的物理過(guò)程演繹將出現(xiàn)的現(xiàn)象 求解定解問(wèn)題就是通過(guò)建立數(shù)學(xué)物理模型、確定定解條件現(xiàn)象現(xiàn)象現(xiàn)象現(xiàn)象現(xiàn)象現(xiàn)象物理模型物理模型3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模的主要步驟數(shù)學(xué)物理建模的主要步驟 將物理問(wèn)題轉(zhuǎn)換成數(shù)學(xué)問(wèn)題將物理問(wèn)題轉(zhuǎn)換成數(shù)學(xué)問(wèn)題 對(duì)物理問(wèn)題根據(jù)相關(guān)的物理定律建立相應(yīng)的數(shù)學(xué)模型，也就是將物理問(wèn)題歸結(jié)成數(shù)學(xué)上的定解問(wèn)題 求解定解問(wèn)題求解定解問(wèn)題 用數(shù)學(xué)方法求出滿足方程和定解條件的解 驗(yàn)證模型的正確性并理

6、解模型的物理意義驗(yàn)證模型的正確性并理解模型的物理意義 對(duì)所得的解通過(guò)輸血的論證和客觀實(shí)踐的檢驗(yàn)鑒定其正確性，并將所得的解作適當(dāng)?shù)奈锢硪饬x解釋，從而理解遵循同一類方程的普遍物理模型3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)本課程的數(shù)學(xué)物理建模的重點(diǎn)本課程的數(shù)學(xué)物理建模的重點(diǎn) 二階線性偏微分方程規(guī)律的物理問(wèn)題的建模二階線性偏微分方程規(guī)律的物理問(wèn)題的建模 靜電勢(shì)和引力勢(shì)滿足的Laplace方程或Poisson方程 波的傳播所滿足的波動(dòng)方程 熱傳導(dǎo)問(wèn)題和擴(kuò)散問(wèn)題中的熱傳導(dǎo)方程 連續(xù)介質(zhì)力學(xué)中的Navier-Stockes方程組合Eule

7、r方程組 描寫電磁場(chǎng)運(yùn)動(dòng)變化的MAxwell方程組 作為微觀物質(zhì)運(yùn)動(dòng)基本規(guī)律的Schrodinger方程和Dirac方程 彈性力學(xué)中的Saint-Venant方程3.1 數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)典型數(shù)學(xué)物理方程的大致分類典型數(shù)學(xué)物理方程的大致分類三類典型的數(shù)學(xué)物理方程三類典型的數(shù)學(xué)物理方程雙曲型方程雙曲型方程波動(dòng)方程為代表波動(dòng)方程為代表拋物型方程拋物型方程熱傳導(dǎo)方程為代表熱傳導(dǎo)方程為代表橢圓型方程橢圓型方程泊松方程為代表泊松方程為代表2,ttuauft r2,tuauft r uf r0Laplaceu 方程3.1 數(shù)

8、學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模的相關(guān)概念數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型弦的微小橫振動(dòng)方程 問(wèn)題：有一條完全柔軟的均勻弦，沿水平直線繃緊，而后用某種方式激發(fā)，使弦在同一個(gè)平面上作小振動(dòng)，列出弦的微小橫振動(dòng)方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型弦的微小橫振動(dòng)方程根據(jù)牛頓第二定律根據(jù)牛頓第二定律F=ma，微元橫向運(yùn)動(dòng)方程為，微元橫向運(yùn)動(dòng)方程為微元的縱向運(yùn)動(dòng)方程為微元的縱向運(yùn)動(dòng)方程為僅考慮微小的橫振動(dòng)，可忽略高階小量，即僅考慮微小的橫振動(dòng)，可忽略

9、高階小量，即2211sinsind( d )ttTTg ss u2211coscos0TT2112cos11, cos12! 311111222sintan, sintan3!222d(d )(d )1 () ddxsxuuxx3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型弦的微小橫振動(dòng)方程注意到注意到因此有因此有這樣弦的橫向和縱向振動(dòng)方程便化簡(jiǎn)成這樣弦的橫向和縱向振動(dòng)方程便化簡(jiǎn)成tansinxuux1122dtansin,tansinxxxxxuu21d21dd0 xxttxxxT uT ug xu

10、xTT3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型弦的微小橫振動(dòng)方程令令T=T1=T2，進(jìn)而可以化簡(jiǎn)成，進(jìn)而可以化簡(jiǎn)成由于由于dx很小，因此有很小，因此有因此可以得到運(yùn)動(dòng)方程因此可以得到運(yùn)動(dòng)方程dddxxttxxxT uug xuxddxxxxxxxuuux0ttxxuTug22ttxxua ugaT3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型弦的微小橫振動(dòng)方程討論：討論：(1)若設(shè)弦的重量遠(yuǎn)小于弦的張力

11、，則重力加速若設(shè)弦的重量遠(yuǎn)小于弦的張力，則重力加速度項(xiàng)可以忽略，可得到齊次偏微分方程度項(xiàng)可以忽略，可得到齊次偏微分方程(弦的自由振動(dòng)方程)(2)若在弦的單位長(zhǎng)度上還有橫向外力若在弦的單位長(zhǎng)度上還有橫向外力F(x, t)的作的作用，則可以得到弦的受迫振動(dòng)方程用，則可以得到弦的受迫振動(dòng)方程2ttxxua u2,ttxxF x tua uf x tf x t3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型薄膜的微小橫振動(dòng)方程 問(wèn)題：有一條完全柔軟的均勻薄膜，沿水平面上繃緊，薄膜的重量和張力相比可以忽略，而后用

12、某種方式激發(fā)，使弦在垂直于水平面方向上作小振動(dòng)3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型薄膜的微小橫振動(dòng)方程注意到二維情況與一維情況的相似之處，可得到注意到二維情況與一維情況的相似之處，可得到微元在微元在x和和x+dx兩邊受力為兩邊受力為微元在微元在y和和y+dy兩邊受力為兩邊受力為那么運(yùn)動(dòng)方程為那么運(yùn)動(dòng)方程為ddd dxxxxxxxT uT uyTux yddd dyyyyyyyT uT uxTux yd dd dd dxxxxttTux yTux yux y3.2 物理模型的數(shù)學(xué)物理建模舉例物

13、理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型薄膜的微小橫振動(dòng)方程經(jīng)過(guò)化簡(jiǎn)可以得到二維模微小橫振動(dòng)方程經(jīng)過(guò)化簡(jiǎn)可以得到二維模微小橫振動(dòng)方程用符號(hào)化的語(yǔ)言可以寫成用符號(hào)化的語(yǔ)言可以寫成如果單位面積上受橫向外力為如果單位面積上受橫向外力為F(x, y, t)，則模的，則模的受迫振動(dòng)方程為受迫振動(dòng)方程為0ttxxyyuT uu20ttuu2, ,ttuuf x y t3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型桿的縱向振動(dòng)問(wèn)題 問(wèn)題：考慮一均勻細(xì)桿

14、，沿桿長(zhǎng)方向作小振動(dòng)，假設(shè)在垂直桿長(zhǎng)方向上的任一橫截面上各點(diǎn)的振動(dòng)情況完全相同3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型桿的縱向振動(dòng)問(wèn)題微元的運(yùn)動(dòng)方程為微元的運(yùn)動(dòng)方程為可以得到桿的縱向振動(dòng)方程可以得到桿的縱向振動(dòng)方程更一般地，在三維空間中的波動(dòng)方程為更一般地，在三維空間中的波動(dòng)方程為dddxxxttxxxuYSuYSuYSxS x ux0ttxxuYu220ttuau桿的縱振動(dòng)和弦的橫振動(dòng)機(jī)桿的縱振動(dòng)和弦的橫振動(dòng)機(jī)理并不完全相同，但滿足的理并不完全相同，但滿足的偏微分方程形式完全一樣，偏微分方程形

15、式完全一樣，這類方程統(tǒng)稱波動(dòng)方程。這類方程統(tǒng)稱波動(dòng)方程。3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程) 問(wèn)題：兩條長(zhǎng)的平行傳輸線在加上交流電壓的時(shí)候會(huì)出現(xiàn)電感、電阻、靜電容和漏電導(dǎo)等效應(yīng)，建立傳輸線上的電壓和電流滿足的方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程)分析微元電路，應(yīng)用基爾霍夫第二定律可得到分析微元電路，應(yīng)用基爾霍夫第二定律可得到考慮到考慮到dx0，

16、則方程可以化簡(jiǎn)成，則方程可以化簡(jiǎn)成同理利用基爾霍夫第一定律可以得到同理利用基爾霍夫第一定律可以得到類似地可以化簡(jiǎn)得到類似地可以化簡(jiǎn)得到,dd ,0tRi x tLix txv xx tv x t0txRiLiv,dd ,0tCvx tGv x txi xx ti x t0txCvGvi3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程)注意到電流和電壓所服從的兩個(gè)微分方程形式相注意到電流和電壓所服從的兩個(gè)微分方程形式相似，對(duì)他們進(jìn)行變量代換，即可得以下兩個(gè)方程似，對(duì)他們進(jìn)行變量代換，即

17、可得以下兩個(gè)方程這兩個(gè)方程就是一般的傳輸線方程這兩個(gè)方程就是一般的傳輸線方程(電報(bào)方程電報(bào)方程)。xxtttvLCvRCGL vGRvxxtttiLCiRCGL iGRi3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程)討論：討論：(1) 無(wú)失真線：無(wú)失真線：RC=LG，信號(hào)無(wú)失真。方程化為，信號(hào)無(wú)失真。方程化為2212xxtttvvvv2212xxtttiiii221,RGLC3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算

18、機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程)(2) 無(wú)損耗線：無(wú)損耗線：RG 0，無(wú)損耗，如高頻條件下，無(wú)損耗，如高頻條件下有有LR，CG。方程化簡(jiǎn)為。方程化簡(jiǎn)為具有與振動(dòng)方程類似的數(shù)學(xué)形式，盡管它們的物具有與振動(dòng)方程類似的數(shù)學(xué)形式，盡管它們的物理本質(zhì)根本不同理本質(zhì)根本不同xxttvLCvxxttiLCi3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型傳輸線方程(電報(bào)方程)(3) 無(wú)漏導(dǎo)，無(wú)電感線：無(wú)漏導(dǎo)，無(wú)電感線： RL0，如同軸電纜。，如同軸電纜。方程可以化簡(jiǎn)為方程可以化簡(jiǎn)為這與后面將講

19、到的一維熱傳導(dǎo)方程具有類似的數(shù)這與后面將講到的一維熱傳導(dǎo)方程具有類似的數(shù)學(xué)形式，盡管它們的物理本質(zhì)根本不同學(xué)形式，盡管它們的物理本質(zhì)根本不同xxtvRCvxxtiRCi3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型電磁波傳播方程 問(wèn)題：設(shè)空間沒(méi)有電荷，且E和H分別表示電場(chǎng)強(qiáng)度和磁感應(yīng)強(qiáng)度。建立空間電場(chǎng)強(qiáng)度和磁感應(yīng)強(qiáng)度的數(shù)學(xué)物理模型3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型電磁波傳播方程由電磁理論，描述介質(zhì)

20、中電磁場(chǎng)運(yùn)動(dòng)的方程為麥由電磁理論，描述介質(zhì)中電磁場(chǎng)運(yùn)動(dòng)的方程為麥克斯韋方程組，其微分形式為克斯韋方程組，其微分形式為為了分析方便，考慮各向同性的均勻介質(zhì)，則介為了分析方便，考慮各向同性的均勻介質(zhì)，則介電常數(shù)、磁導(dǎo)率和電導(dǎo)率都是常數(shù)電常數(shù)、磁導(dǎo)率和電導(dǎo)率都是常數(shù)00tt EHHEEEH3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型電磁波傳播方程對(duì)方程進(jìn)行化簡(jiǎn)可以得到對(duì)方程進(jìn)行化簡(jiǎn)可以得到這兩個(gè)方程是矢量方程，將其在特定方向上這兩個(gè)方程是矢量方程，將其在特定方向上(如如三個(gè)正交坐標(biāo)軸方向三個(gè)正交坐標(biāo)軸方向

21、)進(jìn)行投影，則可以得到普進(jìn)行投影，則可以得到普通的標(biāo)量方程通的標(biāo)量方程這兩個(gè)方程就是電磁波傳播方程這兩個(gè)方程就是電磁波傳播方程ttt EEEttt HHH3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)物理模型波動(dòng)物理模型電磁波傳播方程討論：討論：(1) 考慮電導(dǎo)率很小的情況，即考慮電導(dǎo)率很小的情況，即0，如真空中。，如真空中。方程化簡(jiǎn)為方程化簡(jiǎn)為(2) 考慮介質(zhì)具有很高的導(dǎo)電性的情況，即考慮介質(zhì)具有很高的導(dǎo)電性的情況，即 ，如金屬材料中，則方程可以寫成，如金屬材料中，則方程可以寫成21ttu = aua21tu = au

22、a3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型熱傳導(dǎo)方程 問(wèn)題：由于溫度不均勻，熱量將從溫度高的地方向溫度低的地方轉(zhuǎn)移，這種現(xiàn)象叫熱傳導(dǎo)。試建立熱傳導(dǎo)方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型熱傳導(dǎo)方程推導(dǎo)熱傳導(dǎo)方程的時(shí)候主要根據(jù)熱傳導(dǎo)的推導(dǎo)熱傳導(dǎo)方程的時(shí)候主要根據(jù)熱傳導(dǎo)的Foriour定律在在dt時(shí)間內(nèi)從時(shí)間內(nèi)從ABCD進(jìn)入小體積元的熱量為進(jìn)入小體積元的熱量為同理從同理從EFGH進(jìn)入體積元

23、的熱量為進(jìn)入體積元的熱量為3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例dd duQkS tndd d dd d dxxxuuQkt y zkt y znx ddddd d dd d dxxxxxxuuQkt y zkt y znx數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型熱傳導(dǎo)方程那么通過(guò)這兩個(gè)垂直于那么通過(guò)這兩個(gè)垂直于x軸方向流入小體積元的軸方向流入小體積元的熱量為熱量為類似地有類似地有3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例ddd+dd d dd d d dxxxxxxxuuQQkkt y zxxukt x y

24、zxxdd+dd d d dyyyyuQQkt x y zyydd+dd d d dzzzzuQQkt x y zzz數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型熱傳導(dǎo)方程那么經(jīng)過(guò)六個(gè)面進(jìn)入小體積元的熱量為那么經(jīng)過(guò)六個(gè)面進(jìn)入小體積元的熱量為根據(jù)物體吸熱和溫度變化的關(guān)系可以得到根據(jù)物體吸熱和溫度變化的關(guān)系可以得到3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例dd d d duuuQkkkt x y zxxyyzz0d d d dd d d duuukkkt x y zCu x y zxxyyzz0uuuukkkCxxyyzzt數(shù)學(xué)物理建模與計(jì)

25、算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型熱傳導(dǎo)方程討論：討論：(1) 對(duì)各向同性的均勻物體，對(duì)各向同性的均勻物體，k為常數(shù)，方程變?yōu)闉槌?shù)，方程變?yōu)?2) 若物體內(nèi)有熱源，單位時(shí)間單位體積內(nèi)發(fā)出若物體內(nèi)有熱源，單位時(shí)間單位體積內(nèi)發(fā)出熱量為熱量為F(x, y, z, t)，則有非齊次熱傳導(dǎo)方程，則有非齊次熱傳導(dǎo)方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例0 xxyyzztk uuuC u2200tkuauaC 20, , ,tFuauf x y z tfC 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型擴(kuò)散方程

26、擴(kuò)散方程 問(wèn)題：設(shè)有一塊橫截面積為S的半導(dǎo)體材料，把所需的雜質(zhì)涂敷在材料的表面上，雜質(zhì)就會(huì)向材料里擴(kuò)散。試建立雜質(zhì)的擴(kuò)散方程。3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型擴(kuò)散方程擴(kuò)散方程用用u(x, t)表示表示t時(shí)刻在時(shí)刻在x處的雜質(zhì)濃度。雜質(zhì)的擴(kuò)處的雜質(zhì)濃度。雜質(zhì)的擴(kuò)散服從的物理定律為擴(kuò)散定律散服從的物理定律為擴(kuò)散定律單位時(shí)間內(nèi)，從單位時(shí)間內(nèi)，從x+dx處橫截面流過(guò)的雜質(zhì)量是處橫截面流過(guò)的雜質(zhì)量是流入流入x到到x+dx一個(gè)薄層內(nèi)的雜質(zhì)將引起濃度增加一個(gè)薄層內(nèi)的雜質(zhì)將引起濃度增加,u x tm

27、 x tDSx d ,d ,u xx tm xx tDSx d ,dMum xx tm x tSxtt3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳導(dǎo)物理模型擴(kuò)散方程擴(kuò)散方程因此可以得到因此可以得到即即此即一維形式的擴(kuò)散方程。類似有三維擴(kuò)散方程此即一維形式的擴(kuò)散方程。類似有三維擴(kuò)散方程d ,du xx tu x tu x tDSSxxxt22txxua uaD22tuauaD3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)物理模型熱傳

28、導(dǎo)物理模型擴(kuò)散方程擴(kuò)散方程若外界有擴(kuò)散源，擴(kuò)散源的強(qiáng)度若外界有擴(kuò)散源，擴(kuò)散源的強(qiáng)度(單位時(shí)間、單單位時(shí)間、單位體積內(nèi)產(chǎn)生的雜質(zhì)粒子數(shù)位體積內(nèi)產(chǎn)生的雜質(zhì)粒子數(shù))為為f(x, y, z, t)，則擴(kuò)，則擴(kuò)散方程為散方程為我們注意到，熱傳導(dǎo)和擴(kuò)散這兩種不同的物理現(xiàn)我們注意到，熱傳導(dǎo)和擴(kuò)散這兩種不同的物理現(xiàn)象，但可以用同一類方程來(lái)描述。象，但可以用同一類方程來(lái)描述。2, , ,tuaufx y z t 3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定型物理模型穩(wěn)定型物理模型靜電場(chǎng)電勢(shì)方程靜電場(chǎng)電勢(shì)方程 問(wèn)題：靜電場(chǎng)是有源無(wú)旋場(chǎng)，電

29、力線不閉合，始于正電荷，終于負(fù)電荷。反映靜電場(chǎng)的定理是高斯定理和電場(chǎng)強(qiáng)度的無(wú)旋性。試推導(dǎo)靜電場(chǎng)的方程。3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定型物理模型穩(wěn)定型物理模型靜電場(chǎng)電勢(shì)方程靜電場(chǎng)電勢(shì)方程微分形式的高斯定理為微分形式的高斯定理為靜電場(chǎng)中電場(chǎng)和電勢(shì)的關(guān)系為靜電場(chǎng)中電場(chǎng)和電勢(shì)的關(guān)系為因此可以得到因此可以得到泊松(Poisson)方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例0EU E0U 0U數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定型物理模型穩(wěn)定型物理模型穩(wěn)定溫度分布穩(wěn)定溫度

30、分布 問(wèn)題：如果熱源分布和邊界條件都不隨事件發(fā)生變化，則過(guò)相當(dāng)?shù)臅r(shí)間后物體內(nèi)部的溫度分布將達(dá)到穩(wěn)定狀態(tài)。求其溫度分布方程。 在這種情況下熱傳導(dǎo)方程中隨時(shí)間變化的項(xiàng)將消失，類似地可以得到泊松方程和拉普拉斯方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例21, ,uf x y za 0u 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定型物理模型穩(wěn)定型物理模型穩(wěn)定濃度分布穩(wěn)定濃度分布 問(wèn)題：如果擴(kuò)散源強(qiáng)度f(wàn)(x, y, z)不隨時(shí)間和邊界條件都不隨事件發(fā)生變化，則擴(kuò)散運(yùn)動(dòng)將持續(xù)下去，最終達(dá)到穩(wěn)定狀態(tài)，空間中各點(diǎn)濃度不變。求其濃度分布方程。 與溫度穩(wěn)定分布相似有泊松方程和

31、拉普拉斯方程3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例21, ,uf x y za 0u 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定型物理模型穩(wěn)定型物理模型單色波單色波 問(wèn)題：若電磁波是單色的，即只有一個(gè)頻率，那么電磁波將會(huì)在空間形成穩(wěn)定的分布。根據(jù)要求電磁波場(chǎng)以單一頻率周期變化為根據(jù)要求電磁波場(chǎng)以單一頻率周期變化為那么復(fù)振幅將滿足那么復(fù)振幅將滿足亥姆霍茲(Helmhotz)方程其中其中k =/a為波數(shù)。為波數(shù)。, , , ,expiu x y z tv x y zt20vk v 3.2 物理模型的數(shù)學(xué)物理建模舉例物理模型的數(shù)學(xué)物理建模舉例數(shù)學(xué)物理建模與計(jì)算

32、機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)完整描述一個(gè)物理過(guò)程需兩個(gè)條件完整描述一個(gè)物理過(guò)程需兩個(gè)條件 描述過(guò)程變化的規(guī)律描述過(guò)程變化的規(guī)律數(shù)學(xué)物理方程數(shù)學(xué)物理方程 描述過(guò)程的滿足條件描述過(guò)程的滿足條件定解條件定解條件 例：3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)初始條件初始條件 初始條件是物理過(guò)程剛開始的時(shí)候，也就是初初始條件是物理過(guò)程剛開始的時(shí)候，也就是初始時(shí)刻的系統(tǒng)狀態(tài)始時(shí)刻的系統(tǒng)狀態(tài) 對(duì)時(shí)間的二階偏微分方程來(lái)說(shuō)初始條件包括變對(duì)時(shí)間的二階偏微分方程來(lái)說(shuō)初始條件包括變量對(duì)時(shí)間的零階和一階導(dǎo)數(shù)的初始時(shí)刻的值。量對(duì)時(shí)間的零階和一階導(dǎo)

33、數(shù)的初始時(shí)刻的值。3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)初始條件舉例初始條件舉例 例1：一根長(zhǎng)度為l的弦，兩端固定于x = 0和x = l處，然后在距離坐標(biāo)原點(diǎn)b的位置沿橫向拉開距離h，然后放手讓其振動(dòng)，試寫出其初始條件。3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件,00tux0,0hxxbbu xhlxbxllb數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)邊界條件邊界條件 一個(gè)物理系統(tǒng)，除了描述其運(yùn)動(dòng)特征的數(shù)學(xué)物一個(gè)物理系統(tǒng)，除了描述其運(yùn)動(dòng)特征的數(shù)學(xué)物理方程和描述初始狀態(tài)的初始條件外，還需要理方程和描述初始

34、狀態(tài)的初始條件外，還需要考慮系統(tǒng)所處的考慮系統(tǒng)所處的特定的環(huán)境 周邊的環(huán)境對(duì)物理系統(tǒng)的影響通常由邊界上的周邊的環(huán)境對(duì)物理系統(tǒng)的影響通常由邊界上的物理狀況決定，也就是邊界條件物理狀況決定，也就是邊界條件3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件000,Hf xyzunut第二類邊界條件第二類邊界條件第一類邊界條件第一類邊界條件第三類邊界條件第三類邊界條件數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)波動(dòng)方程的邊界條件波動(dòng)方程的邊界條件 例例1： 例例2： 例例3：3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件 0,utf tu l tg t x lF tuxYS,0 x l

35、uku l txYS數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)熱傳導(dǎo)方程的邊界條件熱傳導(dǎo)方程的邊界條件 例例4：邊界：邊界上的溫度分布為上的溫度分布為f(, t) 例例5：?jiǎn)挝粫r(shí)間內(nèi)從單位面積上流入的熱量為：?jiǎn)挝粫r(shí)間內(nèi)從單位面積上流入的熱量為 例例6：物體表面通過(guò)輻射或?qū)α髋c外界交換熱：物體表面通過(guò)輻射或?qū)α髋c外界交換熱量的時(shí)候肯能夠與邊界上的溫度有關(guān)，有牛頓量的時(shí)候肯能夠與邊界上的溫度有關(guān)，有牛頓冷卻定律：冷卻定律：3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件, , ,u x y z tft1,utnk 11duQH uukuhun數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模

36、與計(jì)算機(jī)輔助設(shè)計(jì)穩(wěn)定場(chǎng)方程的邊界條件穩(wěn)定場(chǎng)方程的邊界條件 穩(wěn)定場(chǎng)不存在初始條件穩(wěn)定場(chǎng)不存在初始條件 對(duì)穩(wěn)定場(chǎng)一樣也存在三種邊界條件對(duì)穩(wěn)定場(chǎng)一樣也存在三種邊界條件 第一類邊界條件：狄利克雷問(wèn)題，給出邊界上的變量值 第二類邊界條件：諾依曼問(wèn)題，給出邊界上變量導(dǎo)數(shù)值 第三類邊界條件：洛平問(wèn)題，給出邊界上的變量值和變量導(dǎo)數(shù)值的線性組合3.3數(shù)學(xué)物理方程的定解條件數(shù)學(xué)物理方程的定解條件,uutn數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的基本概念二階線性偏微分方程的基本概念 偏微分方程偏微分方程 方程的階：未知函數(shù)偏導(dǎo)數(shù)的最高階數(shù) 方程的次數(shù)：最高階偏導(dǎo)數(shù)的次數(shù) 線性方程

37、：所有階(組合)偏導(dǎo)數(shù)的次數(shù)都為1 準(zhǔn)線性方程：僅最高階偏導(dǎo)數(shù)是線性的 自由項(xiàng)：不含未知函數(shù)及其偏導(dǎo)數(shù)的項(xiàng) 齊次方程：沒(méi)有自由項(xiàng)的偏微分方程3.4 二階偏微分方程的分類二階偏微分方程的分類2222, , ,0uuuuF x yuxyxy數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的基本概念二階線性偏微分方程的基本概念 方程的解：若某函數(shù)帶入偏微分方程后，使方程化簡(jiǎn)成一個(gè)恒等式，則此函數(shù)為方程的解 通解：包含任意獨(dú)立函數(shù)的方程的解，且獨(dú)立函數(shù)的個(gè)數(shù)等于方程的階數(shù) 特解：從通解中選擇任意函數(shù)而得到的解3.4 二階偏微分方程的分類二階偏微分方程的分類22sin cosx

38、xyyuuxy2220txxyyxyua ub uc u0 xxyyzzuuuyyuuxyz2yyyxxx uyuln0 xxxxyuuulnsinxxxxyuuuxy數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的基本概念二階線性偏微分方程的基本概念 通解與特解舉例：二階線性非齊次偏微分方程通解與特解舉例：二階線性非齊次偏微分方程 的通解為的通解為 其中其中F(x)和和G(y)為任意函數(shù)。若令為任意函數(shù)。若令 則可以得到方程的一個(gè)特解則可以得到方程的一個(gè)特解3.4 二階偏微分方程的分類二階偏微分方程的分類2xyuyx 221,2u x yxyx yF xG y 4

39、25,2sinF xxG yy2241,252sin2u x yxyx yxy 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類 在利用在利用MATLAB的一些特別的命令和工具求解的一些特別的命令和工具求解二階偏微分方程的時(shí)候需要把方程進(jìn)行歸類和二階偏微分方程的時(shí)候需要把方程進(jìn)行歸類和標(biāo)準(zhǔn)化標(biāo)準(zhǔn)化 前面已經(jīng)初步給出了波動(dòng)型方程、熱傳導(dǎo)型方前面已經(jīng)初步給出了波動(dòng)型方程、熱傳導(dǎo)型方程和穩(wěn)定場(chǎng)型方程，代表三種不同類別的物理程和穩(wěn)定場(chǎng)型方程，代表三種不同類別的物理過(guò)程過(guò)程 二階偏微分方程的分類標(biāo)準(zhǔn)：類似于二次曲線二階偏微分方程的分類標(biāo)準(zhǔn)：類似于二

40、次曲線的分類標(biāo)準(zhǔn)的分類標(biāo)準(zhǔn)3.4 二階偏微分方程的分類二階偏微分方程的分類數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類 兩個(gè)自變量的二階偏微分方程的一般形式為：兩個(gè)自變量的二階偏微分方程的一般形式為：3.4 二階偏微分方程的分類二階偏微分方程的分類,xxxyyyxyA x y uB x y uC x y uD x y uE x y uF x y uG x y220axbxycydxeyf數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類 考慮到考慮到A、B和和C不同時(shí)全為不同時(shí)全為

41、0。作變換。作變換 同時(shí)保證變換的雅可比行列式必須滿足同時(shí)保證變換的雅可比行列式必須滿足 利用多元函數(shù)求導(dǎo)法則進(jìn)行化簡(jiǎn)可以得到利用多元函數(shù)求導(dǎo)法則進(jìn)行化簡(jiǎn)可以得到3.4 二階偏微分方程的分類二階偏微分方程的分類,x yx y,0,xyxyJx yx y aubucudueufug數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類 其中的系數(shù)為其中的系數(shù)為3.4 二階偏微分方程的分類二階偏微分方程的分類222222xxyyxxxyyxyyxxyyxxxyyyxyxxxyyyxyaABCbABCcABCdABCDEeABCDEfFgG 數(shù)學(xué)物理

42、建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類 注意到關(guān)系注意到關(guān)系 注意到：注意到： (1) 在作滿足雅可比行列式不為零的變換的時(shí)候仍然有兩個(gè)相互獨(dú)立的變量，并且變化過(guò)程中的“判別式”的符號(hào)不變。 (2) 可以通過(guò)合適的變換，使得a、b和c的一個(gè)或幾個(gè)為零，達(dá)到化簡(jiǎn)方程的目的3.4 二階偏微分方程的分類二階偏微分方程的分類222,44,bacBACx y 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類為了實(shí)現(xiàn)化簡(jiǎn)，我們首先看下面方程如何求解：為了實(shí)現(xiàn)化簡(jiǎn)，我們首先看下面方程如何求解：為

43、了求解此方程，引入定理：為了求解此方程，引入定理：定理定理：如果：如果 是方程是方程的一般積分，則的一般積分，則 是方程是方程的一個(gè)特解。的一個(gè)特解。3.4 二階偏微分方程的分類二階偏微分方程的分類220 xxyyaABC 0, x yC22dd dd0AyB x yCx, x y220 xxyyABC 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類此定理說(shuō)明：可以通過(guò)求解常微分方程此定理說(shuō)明：可以通過(guò)求解常微分方程得到兩個(gè)變換。一般情況下這兩個(gè)解是無(wú)關(guān)的，得到兩個(gè)變換。一般情況下這兩個(gè)解是無(wú)關(guān)的，稱為偏微分方程的特征線方程，構(gòu)成的曲線

44、為特稱為偏微分方程的特征線方程，構(gòu)成的曲線為特征曲線。征曲線。在求解此方程的時(shí)候，需要判斷其判別式在求解此方程的時(shí)候，需要判斷其判別式3.4 二階偏微分方程的分類二階偏微分方程的分類2dd0ddyyABCxx24BAC 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類討論二階線性偏微分方程的分類討論(1) = B2 - 4AC 0。方程可以得到兩個(gè)實(shí)數(shù)解：。方程可以得到兩個(gè)實(shí)數(shù)解：令令則可以得到則可以得到a=c=0，則偏微分方程可以化簡(jiǎn)為，則偏微分方程可以化簡(jiǎn)為再作進(jìn)一步變換再作進(jìn)一步變換偏微分方程可以化簡(jiǎn)成雙曲型方程偏微分方程可以化簡(jiǎn)成雙曲型方程3.4 二階偏

45、微分方程的分類二階偏微分方程的分類12,x yCx yC,x yx y, , ,0uu u u ,1, , ,0uuu uu 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類討論二階線性偏微分方程的分類討論(2) = B2 - 4AC = 0。方程可以得到兩個(gè)重根，因。方程可以得到兩個(gè)重根，因而只能得到一個(gè)解而只能得到一個(gè)解 。此時(shí)作變換。此時(shí)作變換 即可使得即可使得a=0，進(jìn)一步有，進(jìn)一步有b=0。只要任取變換。只要任取變換 使得雅可比行列式不為零，即可保證兩個(gè)變換相使得雅可比行列式不為零，即可保證兩個(gè)變換相互獨(dú)立。若取互獨(dú)立。若取c=1，即可得到，即可得到拋

46、物線方程拋物線方程3.4 二階偏微分方程的分類二階偏微分方程的分類, , ,0uu u u 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類二階線性偏微分方程的分類(3) = B2 - 4AC 0。方程可以得到兩個(gè)互為共軛。方程可以得到兩個(gè)互為共軛復(fù)數(shù)的根。此時(shí)偏微分方程的兩條特征線是一對(duì)復(fù)數(shù)的根。此時(shí)偏微分方程的兩條特征線是一對(duì)共軛復(fù)函數(shù)族。可以得到兩個(gè)互為共軛復(fù)數(shù)的變共軛復(fù)函數(shù)族?？梢缘玫絻蓚€(gè)互為共軛復(fù)數(shù)的變換換進(jìn)一步地作變換進(jìn)一步地作變換因此可以得到因此可以得到橢圓方程橢圓方程3.4 二階偏微分方程的分類二階偏微分方程的分類,x yx y,i2, , ,0

47、uuu uu 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類討論二階線性偏微分方程的分類討論綜上所述，要判斷二階線性偏微分方程屬于何種綜上所述，要判斷二階線性偏微分方程屬于何種類型，只需討論判別式類型，只需討論判別式= B2 - 4AC 的符號(hào)即可。的符號(hào)即可。在作變換的時(shí)候，需要弄清的幾點(diǎn)是：在作變換的時(shí)候，需要弄清的幾點(diǎn)是：(1) 判別式的符號(hào)，直接決定偏微分方程的類型(2) 如何通過(guò)常微分方程求的偏微分方程的特征線方程，進(jìn)而得到所需的變換(3) 對(duì)拋物線型的偏微分方程，如何尋求另一個(gè)變換下面將通過(guò)幾個(gè)例子來(lái)加以說(shuō)明。下面將通過(guò)幾個(gè)例子來(lái)加以說(shuō)明。3.4

48、二階偏微分方程的分類二階偏微分方程的分類數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類舉例例例1：討論方程：討論方程 的類型，并將其的類型，并將其化成標(biāo)準(zhǔn)型。化成標(biāo)準(zhǔn)型。解：解：偏微分方程的特征方程為偏微分方程的特征方程為求解可以得到求解可以得到3.4 二階偏微分方程的分類二階偏微分方程的分類220 xxyyy ux u22222,0,0440AyBCxx yBACx y 222d0dyyxx222212,2222yxyxcc數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類

49、舉例因此特征曲線為因此特征曲線為作變換作變換則偏微分方程化為則偏微分方程化為3.4 二階偏微分方程的分類二階偏微分方程的分類222212,2222yxyxcc2222,2222yxyx222222uuu 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類舉例再作變換再作變換則偏微分方程可以進(jìn)一步化為則偏微分方程可以進(jìn)一步化為3.4 二階偏微分方程的分類二階偏微分方程的分類,1122uuuu 數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類舉例例例2：將偏微分方程：將偏微分方程

50、化為標(biāo)準(zhǔn)型，并求其通解?；癁闃?biāo)準(zhǔn)型，并求其通解。解：解：上述偏微分方程對(duì)應(yīng)的特征方程為上述偏微分方程對(duì)應(yīng)的特征方程為可以解得可以解得3.4 二階偏微分方程的分類二階偏微分方程的分類22200 xxxyyyx uxyuy uy222,2,040AxBxyCyyBAC 222dd20ddyyxxyyxx11d0dyyxyyc xcxx數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類舉例令令則原方程可以化為則原方程可以化為注意，這兒的注意，這兒的 的取法不是唯一的，只要的取法不是唯一的，只要滿足雅可比行列式不為零即可滿足雅可比行列式不為零即

51、可3.4 二階偏微分方程的分類二階偏微分方程的分類,yyx200uy數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階線性偏微分方程的分類舉例例例3：判斷偏微分方程：判斷偏微分方程 的類的類型，并化為標(biāo)準(zhǔn)型。型，并化為標(biāo)準(zhǔn)型。解：解：方程是橢圓型的，特征方程為方程是橢圓型的，特征方程為求解可以得到求解可以得到3.4 二階偏微分方程的分類二階偏微分方程的分類0 xxxyyyuuu21,1,1430ABCBAC 2dd10ddyyxx 1233i,i2222xxyxcyxc數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)二階線性偏微分方程的分類舉例二階

52、線性偏微分方程的分類舉例令令則方程可以化簡(jiǎn)成則方程可以化簡(jiǎn)成需注意的是，這里的變量代換也不是唯一的，如需注意的是，這里的變量代換也不是唯一的，如3.4 二階偏微分方程的分類二階偏微分方程的分類3,22xyx 22 333uuuu322 3,2233xyxuuuu數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)數(shù)學(xué)物理建模與計(jì)算機(jī)輔助設(shè)計(jì)PDE工具箱的方程分類工具箱的方程分類 橢圓型方程：橢圓型方程：Elliptic 拋物型方程：拋物型方程：Parabolic 雙曲型方程：雙曲型方程：hyperbolic 本征型方程：本征型方程：Eigenmodes3.4 二階偏微分方程的分類二階偏微分方程的分類 div*grad*cua uf * div*grad*d ucua uf * div*grad*d ucua uf div*grad*cua ula