科學(xué)與工程計(jì)算兩點(diǎn)邊值問(wèn)題.ppt

《科學(xué)與工程計(jì)算兩點(diǎn)邊值問(wèn)題.ppt》由會(huì)員分享，可在線閱讀，更多相關(guān)《科學(xué)與工程計(jì)算兩點(diǎn)邊值問(wèn)題.ppt（22頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

,科學(xué)與工程計(jì)算方法,主要內(nèi)容常微分?jǐn)?shù)值計(jì)算方法簡(jiǎn)介偏微分方程數(shù)值計(jì)算方法（差分法、有限元法、有限體積法）統(tǒng)計(jì)方法簡(jiǎn)介選講內(nèi)容計(jì)算流體力學(xué)方法簡(jiǎn)介計(jì)算電磁學(xué)方法簡(jiǎn)介,主要參考書(shū)科學(xué)與工程計(jì)算方法北京理工大學(xué)出版社數(shù)值計(jì)算方法北京理工大學(xué)出版社等偏微分方程數(shù)值解法清華大學(xué)出版社、華中理工大學(xué)出版社等微分方程的數(shù)值方法(英文)springer出版社科學(xué)計(jì)算中的PDE數(shù)值解法科學(xué)出版社微分方程的數(shù)值解法科學(xué)出版社,第一章兩點(diǎn)邊值問(wèn)題數(shù)值解法,1.數(shù)學(xué)模型,例1.電線上的小鳥(niǎo)假設(shè)一根兩端固定的電線上面每個(gè)點(diǎn)都停留一只小鳥(niǎo)，描述此問(wèn)題的數(shù)學(xué)模型,例2.化學(xué)反應(yīng)的動(dòng)力學(xué)模型某種化學(xué)化合物的反應(yīng)可以通過(guò)下面的問(wèn)題描述,,2：線性方程邊值問(wèn)題數(shù)值解法,：導(dǎo)數(shù)逼近方法,第三步：使用適當(dāng)?shù)挠邢薏钌檀鎸?dǎo)數(shù)，如中心差商,則,第五步：將（8）式改寫(xiě)成矩陣的形式，引入向量,,則（8）式可以寫(xiě)為,二：基函數(shù)法,2.1多項(xiàng)式逼近,2.2B-樣條逼近,為了解決高次插值的Runge現(xiàn)象，在多項(xiàng)式插值中會(huì)采用分段Lagrange插值，但是在這里解決二階ODE，對(duì)插值函數(shù)需要一定的光滑性，即二階連續(xù)可導(dǎo)。三次B-樣條插值可以滿(mǎn)足光滑性要求。,,,將（10）代入（4）可得,,2.3Fourier逼近法,三：配置法,例子：B-樣條函數(shù)選取B-樣條中的節(jié)點(diǎn)作為配置點(diǎn)，由此可得,,由邊界條件可得,矩陣形式Ba=W,,四：最小二乘法,由多元函數(shù)求極值,,,五．打靶法,基本思想：利用一階初值問(wèn)題解法來(lái)求解二階二點(diǎn)邊值問(wèn)題。,,求解非線性方程：迭代法求解，比如Newton迭代等,割線迭代法,,,,和,,的選取,如何迭代計(jì)算,例1：用打靶法求解非線性邊值問(wèn)題,,解：線性方程組的初值問(wèn)題,,,,,起源：一些邊值問(wèn)題的解是某些泛函的極小值點(diǎn)。,考慮問(wèn)題,,其中,,為了使討論問(wèn)題簡(jiǎn)單化，齊次化邊界條件,,,六．變分法,準(zhǔn)備工作：,如何求泛函的極小值,基本思想：將在無(wú)窮維的函數(shù)空間上求極值的問(wèn)題變更為在有限維的子空間上求極值,第二步：,,3．非線性邊值問(wèn)題的數(shù)值解法,,考慮問(wèn)題,差商法，即用中心差商作為導(dǎo)數(shù)近似，代入（12）中可得,（13）是非線性代數(shù)方程組，一般來(lái)說(shuō)，需要使用Newton迭代法求近似解,將（13）改寫(xiě)矩陣形式,,Newton迭代法為,,,,4．其他邊界條件的處理,