版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
第五章偏微分方程和特殊函數(shù).5.1引言5.2偏微分方程的分類5.3典型方程的建立5.4定解條件和定解問題5.5線性迭加原理5.6分離變量法5.7非齊次邊界條件的處理第五章偏微分方程和特殊函數(shù).5.1引言5.1引言5.1.1偏微分方程的定義5.1.1.1方程定義:描述物理量在時、空域中變化規(guī)律的方程,若含有未知函數(shù)的偏導(dǎo)數(shù),則稱之為偏微分方程。5.1引言5.1.1偏微分方程的定義5.1.1.2方程的規(guī)定:(1)方程中出現(xiàn)的偏導(dǎo)數(shù)的最高階數(shù)稱為方程的階數(shù)。(2)若方程中沒有未知函數(shù)及其偏導(dǎo)數(shù)的乘積或冪等非線性項稱方程為線性的,反之統(tǒng)稱成為非線性的。在非線性方程中若僅對未知函數(shù)的所有最高階偏導(dǎo)是線性的——擬線性的。(3)不含未知函數(shù)及其偏導(dǎo)數(shù)的項稱之為自由項。自由項為零的方程稱為齊次方程,否則稱為非齊次方程。5.1.1.2方程的規(guī)定:例如:例如:5.1.2偏微分方程的定解問題5.1.2.1n階線性常微分方程對于n階線性常微分方程其中線性無關(guān);
——任意常數(shù)當存在n個邊界條件可以來確定系數(shù):→特解的求解過程。5.1.2偏微分方程的定解問題5.1.2.1n階線性常5.1.2.2偏微分方程在偏微分方程中特解是具有特定形式的任意函數(shù)。如方程:的通解是,特解:
而滿足二維拉普拉斯方程滿足一維熱傳導(dǎo)方程5.1.2.2偏微分方程
由此可以得出兩點結(jié)論:①偏微分方程的通解包含有任意函數(shù),因此解偏微分方程,一般都不是先解通解,后由定解條件確定特解,而是直接求特解。②一個特定形式的偏微分方程可以描述許多物理現(xiàn)象的共性規(guī)律,可以有很多不同形式的特解。所以可稱為“泛定方程”。由此可以得出兩點結(jié)論:
確定地描述某個系統(tǒng)的運動過程,除了反映運動一般規(guī)律的偏微分方程(泛定方程)外,還必須根據(jù)實際問題的模型提出定解條件。泛定方程加定解條件構(gòu)成一個確定的物理過程的“定解問題”,由此求得特定的解。定解條件包括初始條件(當方程含有時間變量時)和邊界條件(關(guān)于空間變量的約束條件)。確定地描述某個系統(tǒng)的運動過程,除了反映運動一般規(guī)律的5.1.2.3偏微分方程的求解方法1、方程的建立建模。2、求解(1)解析解:①分離變量法包括貝塞爾函數(shù)、勒讓德多項式(球柱坐標);②拉普拉斯變換法等。(2)數(shù)值解:尤拉法、龍格-庫塔法等。5.1.2.3偏微分方程的求解方法5.2二階偏微分方程分類.限兩個自變量的二階線性方程,未知函數(shù)u(x,y)一般形式:
F(x,y,u,ux,uy,uxx,uyy,uxy)=0(7-4)線性形式:
Auxx+2Buxy+Cuyy+Dux+Euy+Gu+f=0
(7-5)A,B,C,D,E,f,G是x,y的函數(shù)當A,B,C,D,E,G是常數(shù)時,式(7-5)是二階常系數(shù)線性偏微分方程,f=f(x,y)為已知函數(shù)是自由項。5.2二階偏微分方程分類.限兩個自變量的二階線性方程,未知函由參數(shù)A,B,C判斷二階線性方程的分類設(shè)M=M(x,y)為自變量域內(nèi)的某一點,若在該點處有:(1)B2-AC>0,則方程在該點處為雙曲線型的,如:
uxx-uyy=0
(7-6)(2)B2-AC=0,則方程在該點處為拋物線型,如:
ut-uxx=0
(7-7)(3)B2-AC<0,則方程在該點處為橢圓型的,如:
uxx+uyy=0
(7-8)由參數(shù)A,B,C判斷二階線性方程的分類方程的類型在域內(nèi)不一定是唯一的。如:
xuxx+yuyy+2yux-xuy=0B2-AC=0-xy,①xy<0,M在二,四象限雙曲線型②xy>0,M在一,三象限橢圓型③xy=0(x或y=0),M在兩數(shù)軸拋物線型三類方程:(最典型的物理含義)雙曲線型:utt=a2(uxx+uyy+uzz)波動方程拋物線型:ut=a2(uxx+uyy+uzz)熱傳導(dǎo)方程橢圓型:uxx+uyy+uzz=0
拉普拉斯方程方程的類型在域內(nèi)不一定是唯一的。如:5.3典型方程的建立.
5.3典型方程的建立.
例題例題青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件[例]半無限介質(zhì)中線性熱傳導(dǎo)方程(本無界桿的熱傳導(dǎo)方程),端點溫度恒定為T1,初始時刻溫度均為T0,求溫度分布解:定解問題首先確定對哪個變量做拉氏變換,就本例題而言,對x與t都可進行拉氏變換,但由于缺少u
'(0,t),故選t作L變[例]半無限介質(zhì)中線性熱傳導(dǎo)方程(本無界桿的熱傳導(dǎo)方程)設(shè)即對方程(1)兩端作拉氏變換,有(1)變換為對(3)的邊界條件做拉氏變換設(shè)(4)是關(guān)于x的二階線性常微分方程方程的通解為因為,B=0;因為,所以(4)是關(guān)于x的二階線性常微分方程代入(6)得做逆變換誤差函數(shù)
余誤差函數(shù)代入(6)得用拉氏變換解偏微分方程的要點是:(1)首先確定對哪個自變量作拉氏變換。要求該自變量變化范圍(0,∞),而且根據(jù)拉氏變換的微分性質(zhì)該變量必需具備上式有關(guān)的初值條件。如若有兩個自變量都滿足要求,那應(yīng)取決于對哪個自變量求解過程最簡單為準。(2)除對方程作拉氏變換外,還要對凡在方程變化中沒有用到的定解條件都要作拉氏變換,使其作為變換后新方程的定解條件。(3)最后得到定解問題的解的關(guān)鍵是對新方程之解作拉氏逆變換。當象函數(shù)較復(fù)雜時,運用查表和拉氏變換一章介紹的幾種求逆變換方法也不得其解時,就只能運用拉氏變換的反演公式,通常用復(fù)變函數(shù)的圍道積分法求解。用拉氏變換解偏微分方程的要點是:5.4定解條件和定解問題.
上節(jié)建立的數(shù)學(xué)物理方程是具有某類共性的物理現(xiàn)象的泛定方程。在引言中我們也看到了,同一個泛定方程可以有多個不同函數(shù)的解。因此對實際的物理現(xiàn)象特性的討論還需對其特定的“環(huán)境”和起始狀態(tài)加以描述和限定——定解條件,結(jié)合泛定方程,便可確定定解問題的特解。5.4定解條件和定解問題.上節(jié)建立的數(shù)學(xué)物5.4.1初始條件(初值條件)對于隨著時間而發(fā)生變化的問題,必須考慮研究對象的初始時刻的狀態(tài),即初始條件。凡泛定方程中只含t的一階偏導(dǎo)數(shù)的只需要一個初始條件,u的初始分布。(5-37)泛定方程中含有t的二階偏導(dǎo)數(shù)的則需要兩個初始條件,初始分布和初始速度。
(5-38)初始條件給出了整個系統(tǒng)的狀態(tài)(t=0)。穩(wěn)態(tài)過程因與t無關(guān),則不存在初始條件5.4.1初始條件(初值條件)5.4.2邊界條件(1)第一邊界條件——已知函數(shù)直接給出在邊界上的值(s——Γ上的動點)(7-39)如弦振動,長為的弦兩端固定,則邊界條件為,(7-40)5.4.2邊界條件(1)第一邊界條件——已知函數(shù)又高h,半徑r0的圓柱體的穩(wěn)態(tài)導(dǎo)熱問題
(7-41)第一邊界條件稱為Dirichlet(狄利克利)條件,問題稱為Dirichlet問題。又高h,半徑r0的圓柱體的穩(wěn)態(tài)導(dǎo)熱問題(2)第二類邊界條件——已知導(dǎo)數(shù)一維熱傳導(dǎo)(桿的導(dǎo)熱),設(shè)桿的一端x=a絕熱,則由外到內(nèi)經(jīng)過桿端的熱量流速為零因K,S是常數(shù),故(7-42)對于二維、三維應(yīng)以邊界的外法向?qū)?shù)表述
(7-43)第二類邊界條件稱為Noumann(牛曼)條件,問題稱為Noumann問題。(2)第二類邊界條件——已知導(dǎo)數(shù)(3)第三類邊界條件——混合邊界條件給出邊界上函數(shù)值與其法向?qū)?shù)構(gòu)成的線性關(guān)系。如一維導(dǎo)熱,桿端x=a處自由冷卻,環(huán)境介質(zhì)溫度為u0
,則——桿端散發(fā)出的熱流效率與端點溫度與介質(zhì)溫度之差成正比??筛膶懀?-44)(3)第三類邊界條件——混合邊界條件對于長為的桿兩端自由冷卻(7-45)第三類邊界條件的一般形式(7-46)對于長為的桿兩端自由冷卻5.5線性迭加原理.
在講如何用分離變量法求解偏微分方程的定解問題前,先介紹一下線性偏微分方程解的迭加原理。
所謂線性疊加就是幾種不同因素綜合作用于系統(tǒng),產(chǎn)生的效果等于各因素獨立作用產(chǎn)生的效果的總和。5.5線性迭加原理.迭加原理:設(shè)函數(shù)是齊次線性偏微分方程的特解,若級數(shù)可逐項求偏微分,則該級數(shù)也是齊次線性偏微分方程的解。迭加原理:5.6分離變量法.對于多個自變量的偏微分方程定解問題的求解,在可能的情況下,我們總設(shè)法使自變量的個數(shù)減少。分離變量法就是基于這種想法產(chǎn)生的。分離變量法也叫傅立葉方法,它利用變量分離形式的解法,將求解偏微分方程的定解問題化為求解常微分方程的固有值問題,再利用定解條件及有關(guān)數(shù)學(xué)方法,求得定解問題的解。分離變量法對定解條件尤其是邊界條件的要求比較苛刻,一般只涉及較為規(guī)則的邊界問題。我們主要是通過各種例題來介紹分離變量法的具體應(yīng)用。5.6分離變量法.對于多個自變量的偏微分方程定解問題的求[例]有界弦的自由振動解:弦長為兩端張緊固定且無外力作用的弦振動問題,可用下述定解問題表述
這是齊次方程,齊次邊界條件的定解問題。[例]有界弦的自由振動(1)分離變量設(shè)其解函數(shù)可以表示為兩個單自變量函數(shù)的乘積。令u(x,t)=X(x)T(t),X(x)——x的函數(shù),T(t)——t的函數(shù)代入(7-52)得分離變量改寫為(1)分離變量不妨令其等于常數(shù)λ得到兩個常微分方程(7-56)(7-57)由邊界條件(7-53)
(7-58)不妨令其等于常數(shù)λ(2)求本征值(固有值)λ(i)設(shè)λ>0,則(5-56)的通解為代入條件(5-58)
解得A=-B=0,即X(x)≡0,不合題意舍去。(u≡0)(ii)設(shè)λ=0得,
由條件(5-58)得A=B=0,X(x)≡0,不合題意舍去。(2)求本征值(固有值)λ(iii)設(shè)λ<0,不妨令λ=-β2
,式(5-56)的通解為
由條件(7-58)
因為β≠0
,所以
(iii)設(shè)λ<0,不妨令λ=-β2,式(5-56)的λn——固有值或本征值.Xn(x)——固有函數(shù),常寫成不帶系數(shù)的形式:
λn——固有值或本征值.(3)求(5-57)的通解Tn(t)及定解問題的一組特解un(x,t)對,有通解為(5-61)(5-60)×(5-61)得一組特解(5-62)(3)求(5-57)的通解Tn(t)及定解問題的一組特解un(4)由傅立葉級數(shù)確定系數(shù)Cn,Dn,求由解的迭加原理(5-63)代入初始條件(5-55)、(5-56)
(4)由傅立葉級數(shù)確定系數(shù)Cn,Dn,求
式(5-64)和式(5-65)分別是φ(x),ψ(x)
的傅立葉正弦級數(shù)的展開式,而φ(x),ψ(x)是由初始條件給出的定義在[0,l]上的連續(xù)函數(shù)(或只有有限個第一類間斷點,且至多有有限個極值點),所以只要選取即代入(5-63)即得定解問題的完整特解。式(5-64)和式(5-65)分別是φ(x),ψ(x)的[例]
這是齊次方程,齊次邊界條件的定解問題。[例]解:u(x,t)是其一個解函數(shù)假設(shè)函數(shù)可以表示為各個自變量單元函數(shù)的乘積,代入方程后可分離為各自變量的常微分方程。設(shè)u(x,t)=X(x)T(t),X(x)-x的函數(shù),T(t)-t的函數(shù)解:u(x,t)是其一個解函數(shù)代入原方程中:代入原方程中:將邊界條件代入:將邊界條件代入:運用疊加原理運用正交函數(shù):兩邊同乘以,并從積分。(為需要求的待定系數(shù))運用正交函數(shù):兩邊同乘以,并從運用正交函數(shù):兩邊同乘以,并從運用正交函數(shù):兩邊同乘以,并從積分。運用正交函數(shù):兩邊同乘以,并從運用疊加原理運用正交函數(shù):兩邊同乘以,并從積分。(為需要求∴∴u(x?t)=總結(jié):①通過假設(shè),將變量分離②確定待定系數(shù)(通過已知條件)③利用疊加原理得到解函數(shù)∴∴u(x?t)=總結(jié):①通過假設(shè),將變量分離青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件
得到F(詳見P223頁)得到F(詳見P223頁)【習(xí)題】用分離變量法解下列振動問題初始條件:(2)兩端固定,解:(2)的定解問題為:【習(xí)題】用分離變量法解下列振動問題1°分離變量令u=X?T,代入方程(1)得(λ為常量)(4)(5)由邊界條件(2)得1°分離變量2°求固有值λ①當λ>0時,(6)的通解為代入邊界條件(7)得不合題意,舍去;②當λ=0時,(6)的通解為代入邊界條件(7)得不合題意,舍去;2°求固有值λ③當λ<0時,令λ=β2
,方程(6)的通解為由邊界條件(7)得③當λ<0時,令λ=β2,方程(6)的通解為3°求將代入方程(5)得通解為3°求4°確定系數(shù)Cn,Dn
,求u(x,t)由解的疊加原理由初始條件(3)4°確定系數(shù)Cn,Dn,求u(x,t)青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件【例(補充)】有限桿的熱傳導(dǎo)長為l的均勻細桿放置x軸上,其側(cè)面及兩端面絕熱,桿內(nèi)無熱源且桿內(nèi)初始溫度分布為φ(x)。定解問題:【例(補充)】有限桿的熱傳導(dǎo)解:(1)分離變量令代入方程式(1)得解:(2)求λ(i)設(shè)λ>0,通解為代入邊界條件得:
即X(x)≡0,不合題意。(ii)λ=0,通解為代入邊界條件得所以(2)求λ(iii)設(shè)λ<0,取,通解
代入邊界條件得:所以固有值固有函數(shù)(iii)設(shè)λ<0,?。?)求(i)當λ0=0時
化為得解
(ii)當(3)求(4)確定cn及u(x,t),由疊加原理代入初始條件代入u(x,t)的表達式即得定解問題的解。(4)確定cn及u(x,t),由疊加原理5.7非齊次邊界條件的處理.
前面分離變量涉及的定解條件是齊次邊界條件,對非齊次邊界條件則要進行齊次化的處理。以一維熱傳導(dǎo)定解問題為例介紹非齊次邊界條件的處理方法——邊界條件齊次化。一維熱傳導(dǎo)定解問題:5.7非齊次邊界條件的處理.前面分離變量涉作代換:(7-98)選取適當?shù)腤(x,t),使V(x,t)滿足齊次邊界條件(7-99)則作代換:W(x,t)的形式力求簡單,達到V(x,t)邊界條件齊次化的目的即可,不妨設(shè)W(x,t)為x的一次函數(shù)。令
由邊界條件式(7-100)可得于是
W(x,t)的形式力求簡單,達到V(x,t)邊界條件齊次化的將式(5-101)及式(5-98)代入(A)得問題(A)轉(zhuǎn)化為問題(B)
將式(5-101)及式(5-98)代入(A)得
經(jīng)邊界條件齊次化處理后,定解問題(A)的求解轉(zhuǎn)化為定解問題(B)的求解。由(B)解得V(x,t),則(A)的解為。一般情況下(B)是非齊次方程的定解問題。在這一節(jié)我們只看形如的齊次方程的求解,而形如的非齊次方程的求解放到下一節(jié)討論。經(jīng)邊界條件齊次化處理后,定解問題(A)的求解轉(zhuǎn)類似于定解問題(A)的其它類型的非齊次邊界條件有如相應(yīng)的W(x,t)有不同的結(jié)果。類似于定解問題(A)的其它類型的非齊次邊界條件有如[例5-13]一維熱傳導(dǎo)問題
設(shè)一均勻細桿,初始時全桿有一均勻溫度,然后使其一端保持不變溫度,另一端則有恒定的熱流輸入,試求溫度分布規(guī)律。此定解問題可描述如下[例5-13]一維熱傳導(dǎo)問題設(shè)一均勻細桿,解:(1)邊界條件齊次化設(shè)為使則解:(1)邊界條件齊次化設(shè)
代入(7-104)得所以代入原問題(A)得到定解問題(B)
(2)求解定解問題(B)這是關(guān)于V(x,t)的齊次邊界問題①分離變量令代入(7-105)并分離變量,并結(jié)合齊次邊界條件得(7-108)
(2)求解定解問題(B)這是關(guān)于V(x,t)的齊次邊②求λ由前面的例題分析可知,λ<0時(B)有解,不妨設(shè),式(7-109)的通解為由邊界條件式(7-110),得固有函數(shù)
②求λ由前面的例題分析可知,λ<0時(B)有解,不妨設(shè)③求λ代入(7-108),得③求④確定系數(shù)由初始條件式(7-107)④確定系數(shù)利用Fourier級數(shù)展開法求系數(shù)從而利用Fourier級數(shù)展開法求系數(shù)(3)求
解式(5-111)的結(jié)構(gòu)可以看出隨時間的推延,此級數(shù)解迅速收斂。(3)求第五章偏微分方程和特殊函數(shù).5.1引言5.2偏微分方程的分類5.3典型方程的建立5.4定解條件和定解問題5.5線性迭加原理5.6分離變量法5.7非齊次邊界條件的處理第五章偏微分方程和特殊函數(shù).5.1引言5.1引言5.1.1偏微分方程的定義5.1.1.1方程定義:描述物理量在時、空域中變化規(guī)律的方程,若含有未知函數(shù)的偏導(dǎo)數(shù),則稱之為偏微分方程。5.1引言5.1.1偏微分方程的定義5.1.1.2方程的規(guī)定:(1)方程中出現(xiàn)的偏導(dǎo)數(shù)的最高階數(shù)稱為方程的階數(shù)。(2)若方程中沒有未知函數(shù)及其偏導(dǎo)數(shù)的乘積或冪等非線性項稱方程為線性的,反之統(tǒng)稱成為非線性的。在非線性方程中若僅對未知函數(shù)的所有最高階偏導(dǎo)是線性的——擬線性的。(3)不含未知函數(shù)及其偏導(dǎo)數(shù)的項稱之為自由項。自由項為零的方程稱為齊次方程,否則稱為非齊次方程。5.1.1.2方程的規(guī)定:例如:例如:5.1.2偏微分方程的定解問題5.1.2.1n階線性常微分方程對于n階線性常微分方程其中線性無關(guān);
——任意常數(shù)當存在n個邊界條件可以來確定系數(shù):→特解的求解過程。5.1.2偏微分方程的定解問題5.1.2.1n階線性常5.1.2.2偏微分方程在偏微分方程中特解是具有特定形式的任意函數(shù)。如方程:的通解是,特解:
而滿足二維拉普拉斯方程滿足一維熱傳導(dǎo)方程5.1.2.2偏微分方程
由此可以得出兩點結(jié)論:①偏微分方程的通解包含有任意函數(shù),因此解偏微分方程,一般都不是先解通解,后由定解條件確定特解,而是直接求特解。②一個特定形式的偏微分方程可以描述許多物理現(xiàn)象的共性規(guī)律,可以有很多不同形式的特解。所以可稱為“泛定方程”。由此可以得出兩點結(jié)論:
確定地描述某個系統(tǒng)的運動過程,除了反映運動一般規(guī)律的偏微分方程(泛定方程)外,還必須根據(jù)實際問題的模型提出定解條件。泛定方程加定解條件構(gòu)成一個確定的物理過程的“定解問題”,由此求得特定的解。定解條件包括初始條件(當方程含有時間變量時)和邊界條件(關(guān)于空間變量的約束條件)。確定地描述某個系統(tǒng)的運動過程,除了反映運動一般規(guī)律的5.1.2.3偏微分方程的求解方法1、方程的建立建模。2、求解(1)解析解:①分離變量法包括貝塞爾函數(shù)、勒讓德多項式(球柱坐標);②拉普拉斯變換法等。(2)數(shù)值解:尤拉法、龍格-庫塔法等。5.1.2.3偏微分方程的求解方法5.2二階偏微分方程分類.限兩個自變量的二階線性方程,未知函數(shù)u(x,y)一般形式:
F(x,y,u,ux,uy,uxx,uyy,uxy)=0(7-4)線性形式:
Auxx+2Buxy+Cuyy+Dux+Euy+Gu+f=0
(7-5)A,B,C,D,E,f,G是x,y的函數(shù)當A,B,C,D,E,G是常數(shù)時,式(7-5)是二階常系數(shù)線性偏微分方程,f=f(x,y)為已知函數(shù)是自由項。5.2二階偏微分方程分類.限兩個自變量的二階線性方程,未知函由參數(shù)A,B,C判斷二階線性方程的分類設(shè)M=M(x,y)為自變量域內(nèi)的某一點,若在該點處有:(1)B2-AC>0,則方程在該點處為雙曲線型的,如:
uxx-uyy=0
(7-6)(2)B2-AC=0,則方程在該點處為拋物線型,如:
ut-uxx=0
(7-7)(3)B2-AC<0,則方程在該點處為橢圓型的,如:
uxx+uyy=0
(7-8)由參數(shù)A,B,C判斷二階線性方程的分類方程的類型在域內(nèi)不一定是唯一的。如:
xuxx+yuyy+2yux-xuy=0B2-AC=0-xy,①xy<0,M在二,四象限雙曲線型②xy>0,M在一,三象限橢圓型③xy=0(x或y=0),M在兩數(shù)軸拋物線型三類方程:(最典型的物理含義)雙曲線型:utt=a2(uxx+uyy+uzz)波動方程拋物線型:ut=a2(uxx+uyy+uzz)熱傳導(dǎo)方程橢圓型:uxx+uyy+uzz=0
拉普拉斯方程方程的類型在域內(nèi)不一定是唯一的。如:5.3典型方程的建立.
5.3典型方程的建立.
例題例題青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件[例]半無限介質(zhì)中線性熱傳導(dǎo)方程(本無界桿的熱傳導(dǎo)方程),端點溫度恒定為T1,初始時刻溫度均為T0,求溫度分布解:定解問題首先確定對哪個變量做拉氏變換,就本例題而言,對x與t都可進行拉氏變換,但由于缺少u
'(0,t),故選t作L變[例]半無限介質(zhì)中線性熱傳導(dǎo)方程(本無界桿的熱傳導(dǎo)方程)設(shè)即對方程(1)兩端作拉氏變換,有(1)變換為對(3)的邊界條件做拉氏變換設(shè)(4)是關(guān)于x的二階線性常微分方程方程的通解為因為,B=0;因為,所以(4)是關(guān)于x的二階線性常微分方程代入(6)得做逆變換誤差函數(shù)
余誤差函數(shù)代入(6)得用拉氏變換解偏微分方程的要點是:(1)首先確定對哪個自變量作拉氏變換。要求該自變量變化范圍(0,∞),而且根據(jù)拉氏變換的微分性質(zhì)該變量必需具備上式有關(guān)的初值條件。如若有兩個自變量都滿足要求,那應(yīng)取決于對哪個自變量求解過程最簡單為準。(2)除對方程作拉氏變換外,還要對凡在方程變化中沒有用到的定解條件都要作拉氏變換,使其作為變換后新方程的定解條件。(3)最后得到定解問題的解的關(guān)鍵是對新方程之解作拉氏逆變換。當象函數(shù)較復(fù)雜時,運用查表和拉氏變換一章介紹的幾種求逆變換方法也不得其解時,就只能運用拉氏變換的反演公式,通常用復(fù)變函數(shù)的圍道積分法求解。用拉氏變換解偏微分方程的要點是:5.4定解條件和定解問題.
上節(jié)建立的數(shù)學(xué)物理方程是具有某類共性的物理現(xiàn)象的泛定方程。在引言中我們也看到了,同一個泛定方程可以有多個不同函數(shù)的解。因此對實際的物理現(xiàn)象特性的討論還需對其特定的“環(huán)境”和起始狀態(tài)加以描述和限定——定解條件,結(jié)合泛定方程,便可確定定解問題的特解。5.4定解條件和定解問題.上節(jié)建立的數(shù)學(xué)物5.4.1初始條件(初值條件)對于隨著時間而發(fā)生變化的問題,必須考慮研究對象的初始時刻的狀態(tài),即初始條件。凡泛定方程中只含t的一階偏導(dǎo)數(shù)的只需要一個初始條件,u的初始分布。(5-37)泛定方程中含有t的二階偏導(dǎo)數(shù)的則需要兩個初始條件,初始分布和初始速度。
(5-38)初始條件給出了整個系統(tǒng)的狀態(tài)(t=0)。穩(wěn)態(tài)過程因與t無關(guān),則不存在初始條件5.4.1初始條件(初值條件)5.4.2邊界條件(1)第一邊界條件——已知函數(shù)直接給出在邊界上的值(s——Γ上的動點)(7-39)如弦振動,長為的弦兩端固定,則邊界條件為,(7-40)5.4.2邊界條件(1)第一邊界條件——已知函數(shù)又高h,半徑r0的圓柱體的穩(wěn)態(tài)導(dǎo)熱問題
(7-41)第一邊界條件稱為Dirichlet(狄利克利)條件,問題稱為Dirichlet問題。又高h,半徑r0的圓柱體的穩(wěn)態(tài)導(dǎo)熱問題(2)第二類邊界條件——已知導(dǎo)數(shù)一維熱傳導(dǎo)(桿的導(dǎo)熱),設(shè)桿的一端x=a絕熱,則由外到內(nèi)經(jīng)過桿端的熱量流速為零因K,S是常數(shù),故(7-42)對于二維、三維應(yīng)以邊界的外法向?qū)?shù)表述
(7-43)第二類邊界條件稱為Noumann(牛曼)條件,問題稱為Noumann問題。(2)第二類邊界條件——已知導(dǎo)數(shù)(3)第三類邊界條件——混合邊界條件給出邊界上函數(shù)值與其法向?qū)?shù)構(gòu)成的線性關(guān)系。如一維導(dǎo)熱,桿端x=a處自由冷卻,環(huán)境介質(zhì)溫度為u0
,則——桿端散發(fā)出的熱流效率與端點溫度與介質(zhì)溫度之差成正比??筛膶懀?-44)(3)第三類邊界條件——混合邊界條件對于長為的桿兩端自由冷卻(7-45)第三類邊界條件的一般形式(7-46)對于長為的桿兩端自由冷卻5.5線性迭加原理.
在講如何用分離變量法求解偏微分方程的定解問題前,先介紹一下線性偏微分方程解的迭加原理。
所謂線性疊加就是幾種不同因素綜合作用于系統(tǒng),產(chǎn)生的效果等于各因素獨立作用產(chǎn)生的效果的總和。5.5線性迭加原理.迭加原理:設(shè)函數(shù)是齊次線性偏微分方程的特解,若級數(shù)可逐項求偏微分,則該級數(shù)也是齊次線性偏微分方程的解。迭加原理:5.6分離變量法.對于多個自變量的偏微分方程定解問題的求解,在可能的情況下,我們總設(shè)法使自變量的個數(shù)減少。分離變量法就是基于這種想法產(chǎn)生的。分離變量法也叫傅立葉方法,它利用變量分離形式的解法,將求解偏微分方程的定解問題化為求解常微分方程的固有值問題,再利用定解條件及有關(guān)數(shù)學(xué)方法,求得定解問題的解。分離變量法對定解條件尤其是邊界條件的要求比較苛刻,一般只涉及較為規(guī)則的邊界問題。我們主要是通過各種例題來介紹分離變量法的具體應(yīng)用。5.6分離變量法.對于多個自變量的偏微分方程定解問題的求[例]有界弦的自由振動解:弦長為兩端張緊固定且無外力作用的弦振動問題,可用下述定解問題表述
這是齊次方程,齊次邊界條件的定解問題。[例]有界弦的自由振動(1)分離變量設(shè)其解函數(shù)可以表示為兩個單自變量函數(shù)的乘積。令u(x,t)=X(x)T(t),X(x)——x的函數(shù),T(t)——t的函數(shù)代入(7-52)得分離變量改寫為(1)分離變量不妨令其等于常數(shù)λ得到兩個常微分方程(7-56)(7-57)由邊界條件(7-53)
(7-58)不妨令其等于常數(shù)λ(2)求本征值(固有值)λ(i)設(shè)λ>0,則(5-56)的通解為代入條件(5-58)
解得A=-B=0,即X(x)≡0,不合題意舍去。(u≡0)(ii)設(shè)λ=0得,
由條件(5-58)得A=B=0,X(x)≡0,不合題意舍去。(2)求本征值(固有值)λ(iii)設(shè)λ<0,不妨令λ=-β2
,式(5-56)的通解為
由條件(7-58)
因為β≠0
,所以
(iii)設(shè)λ<0,不妨令λ=-β2,式(5-56)的λn——固有值或本征值.Xn(x)——固有函數(shù),常寫成不帶系數(shù)的形式:
λn——固有值或本征值.(3)求(5-57)的通解Tn(t)及定解問題的一組特解un(x,t)對,有通解為(5-61)(5-60)×(5-61)得一組特解(5-62)(3)求(5-57)的通解Tn(t)及定解問題的一組特解un(4)由傅立葉級數(shù)確定系數(shù)Cn,Dn,求由解的迭加原理(5-63)代入初始條件(5-55)、(5-56)
(4)由傅立葉級數(shù)確定系數(shù)Cn,Dn,求
式(5-64)和式(5-65)分別是φ(x),ψ(x)
的傅立葉正弦級數(shù)的展開式,而φ(x),ψ(x)是由初始條件給出的定義在[0,l]上的連續(xù)函數(shù)(或只有有限個第一類間斷點,且至多有有限個極值點),所以只要選取即代入(5-63)即得定解問題的完整特解。式(5-64)和式(5-65)分別是φ(x),ψ(x)的[例]
這是齊次方程,齊次邊界條件的定解問題。[例]解:u(x,t)是其一個解函數(shù)假設(shè)函數(shù)可以表示為各個自變量單元函數(shù)的乘積,代入方程后可分離為各自變量的常微分方程。設(shè)u(x,t)=X(x)T(t),X(x)-x的函數(shù),T(t)-t的函數(shù)解:u(x,t)是其一個解函數(shù)代入原方程中:代入原方程中:將邊界條件代入:將邊界條件代入:運用疊加原理運用正交函數(shù):兩邊同乘以,并從積分。(為需要求的待定系數(shù))運用正交函數(shù):兩邊同乘以,并從運用正交函數(shù):兩邊同乘以,并從運用正交函數(shù):兩邊同乘以,并從積分。運用正交函數(shù):兩邊同乘以,并從運用疊加原理運用正交函數(shù):兩邊同乘以,并從積分。(為需要求∴∴u(x?t)=總結(jié):①通過假設(shè),將變量分離②確定待定系數(shù)(通過已知條件)③利用疊加原理得到解函數(shù)∴∴u(x?t)=總結(jié):①通過假設(shè),將變量分離青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件
得到F(詳見P223頁)得到F(詳見P223頁)【習(xí)題】用分離變量法解下列振動問題初始條件:(2)兩端固定,解:(2)的定解問題為:【習(xí)題】用分離變量法解下列振動問題1°分離變量令u=X?T,代入方程(1)得(λ為常量)(4)(5)由邊界條件(2)得1°分離變量2°求固有值λ①當λ>0時,(6)的通解為代入邊界條件(7)得不合題意,舍去;②當λ=0時,(6)的通解為代入邊界條件(7)得不合題意,舍去;2°求固有值λ③當λ<0時,令λ=β2
,方程(6)的通解為由邊界條件(7)得③當λ<0時,令λ=β2,方程(6)的通解為3°求將代入方程(5)得通解為3°求4°確定系數(shù)Cn,Dn
,求u(x,t)由解的疊加原理由初始條件(3)4°確定系數(shù)Cn,Dn,求u(x,t)青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課件青海大學(xué)-化工應(yīng)用數(shù)學(xué)-第五章-偏微分方程及特殊課
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 酒店總經(jīng)理助理崗位職責
- 電力線路雨季冬季特殊現(xiàn)場施工措施
- 資產(chǎn)收購協(xié)議
- 證券期貨交易委托合同
- 酒店承包合同范本
- 電力施工臨時用工合同范本
- 2025年一級醫(yī)院放射防護規(guī)章制度(三篇)
- 醫(yī)院創(chuàng)建二甲醫(yī)院工作實施方案樣本(4篇)
- 辦公環(huán)境管理制度范例(2篇)
- 2025年開學(xué)國旗下的講話演講稿樣本(3篇)
- 湖南省湘西州吉首市2023屆九年級上學(xué)期期末素質(zhì)監(jiān)測數(shù)學(xué)試卷(含解析)
- 2023-2024學(xué)年湖北省武漢市東西湖區(qū)三年級(上)期末數(shù)學(xué)試卷
- GB/T 31771-2024家政服務(wù)母嬰護理服務(wù)質(zhì)量規(guī)范
- 2023-建筑施工技02課件講解
- 期末試卷:福建省廈門市集美區(qū)2021-2022學(xué)年八年級上學(xué)期期末歷史試題(原卷版)
- 美容院2024年度規(guī)劃
- 裝飾裝修巡查記錄表
- 公司安全生產(chǎn)事故隱患內(nèi)部報告獎勵工作制度
- 艾滋病預(yù)防知識講座
- 零售服務(wù)質(zhì)量提升
- 《4 平平安安回家來》 說課稿-2024-2025學(xué)年道德與法治一年級上冊統(tǒng)編版
評論
0/150
提交評論