版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
數(shù)值分析常微分方程邊值問題第一頁,共十六頁,編輯于2023年,星期三初值問題歐拉公式:yn+1=yn+hf(xn,yn)修改的歐拉公式(Heun方法):k1=f(xn,yn),k2=f(xn+1,yn+hk1)(n=0,1,2,·······,N)(n=0,1,2,·······,N)2/16第二頁,共十六頁,編輯于2023年,星期三二階Range-Kutta公式一般形式y(tǒng)n+1=yn+h[c1k1+c2k2]k1=f(xn,yn)k2=f(xn+λ2h,yn+μ21hk1)(n=0,1,····,N)c1=0.5,c2=0.5,λ2=1,
μ21=1c1=0,c2=1,λ2=0.5,
μ21=0.5
yn+1=yn+0.5h[k1+k2]k1=f(xn,yn)k2=f(xn+h,yn+hk1)
yn+1=yn+h
k2k1=f(xn,yn)k2=f(xn+.5h,yn+.5hk1)3/16第三頁,共十六頁,編輯于2023年,星期三三階Range-Kutta公式一般形式y(tǒng)n+1=yn+h[k1+4k2+k3]/6k1=f(xn,yn),k2=f(xn+0.5h,yn+0.5hk1)k3=f(xn+h,yn–hk1+2hk2)四階Range-Kutta公式一般形式y(tǒng)n+1=yn+h[k1+2k2+2k3+k4]/6k1=f(xn,yn),k2=f(xn+0.5h,yn+0.5hk1)k3=f(xn+0.5h,yn+0.5hk2),k4=f(xn+h,yn+hk3)4/16第四頁,共十六頁,編輯于2023年,星期三例1數(shù)值實驗:幾種不同求數(shù)值解公式的誤差比較
n10203040h0.20.10.06670.05RK46.862e-0053.747e-0067.071e-0072.186e-007RK30.00121.529e-0044.517e-0051.906e-005RK20.01230.00260.00115.9612e-004Euler0.10590.05210.03420.02565/16第五頁,共十六頁,編輯于2023年,星期三n=input('inputn:=');f=inline('y-x.*y.^2');h=2/n;x=h:h:2;y0=1;k1=f(0,y0);k2=f(.5*h,y0+0.5*h*k1);k3=f(.5*h,y0+.5*h*k2);k4=f(h,y0+h*k3);y(1)=y0+h*(k1+2*k2+2*k3+k4)/6;fork=1:n-1xk=x(k);yk=y(k);k1=f(xk,yk);k2=f(xk+h/2,yk+0.5*h*k1);k3=f(xk+h/2,yk+0.5*h*k2);k4=f(xk+h,yk+h*k3);y(k+1)=yk+h*(k1+2*k2+2*k3+k4)/6;endt=-1:h/2:5;yy=1./(t-1+2*exp(-t));plot(0,1,'ro',x,y,'ro',t,yy)y1=1./(x-1+2*exp(-x));error=max(abs(y-y1))inputn:=20error=3.7475e-0066/16第六頁,共十六頁,編輯于2023年,星期三例2.捕食者與被捕食者模型0≤x≤15
y1(0)=20,y2(0)=20記7/16第七頁,共十六頁,編輯于2023年,星期三x(1)=20;y(1)=20;forn=1:length(t)-1tn=t(n);xn=x(n);yn=y(n);k11=f1(xn,yn);k21=f2(xn,yn);k12=f1(xn+.5*h*k11,yn+.5*h*k21);k22=f2(xn+.5*h*k11,yn+.5*h*k21);k13=f1(xn+.5*h*k12,yn+.5*h*k22);k23=f2(xn+.5*h*k12,yn+.5*h*k22);k14=f1(xn+h*k13,yn+h*k23);k24=f2(xn+h*k13,yn+h*k23);x(n+1)=xn+h*(k11+2*k12+2*k13+k14)/6;y(n+1)=yn+h*(k21+2*k22+2*k23+k24)/6;endf1=inline('x-0.01*x*y');f2=inline('-y+0.02*x*y');h=.15;t=0:h:15;----y1----y2
y1—y2
相位圖8/16第八頁,共十六頁,編輯于2023年,星期三例1.單擺的數(shù)學(xué)模型其中,a=g/L初值條件:(0)=0.4,
’(0)=0第一步:轉(zhuǎn)化為一階方程組令:y1=,y2=
’
初值條件:y1(0)=0.4,y2(0)=0第二步:求解方程組functionf=dan(x,y)f(1,:)=y(2);f(2,:)=-9.8*sin(y(1))/3.2;L=3.2ode23('dan',[0,2],[0.4,0]);9/16第九頁,共十六頁,編輯于2023年,星期三求解邊值問題的數(shù)值方法算例解:取正整數(shù)n,令h=1/(n+1),xj=jh,(j=0,1,···,n+1
).
將常微分方程離散化
整理,得:
–yj-1+(2–h2)yj–yj+1=xjh2
(j=1,2,···,n)y0=0,yn=01.高斯-賽德爾迭代;2.打靶法;3.追趕法10/16第十頁,共十六頁,編輯于2023年,星期三精確解symsxydsolve('D2y+y=-x','y(0)=0','y(1)=0','x')符號系統(tǒng)計算結(jié)果為:ans=-x+1/sin(1)*sin(x)
–yj-1+(2–h2)yj–yj+1=xjh2
(j=1,2,···,n)方程組AY=f
系數(shù)矩陣:11/16第十一頁,共十六頁,編輯于2023年,星期三n=input('inputn:=');y=inline('sin(x)/sin(1)-x');h=1/(n+1);x=0:h:1;hh=h*h;r=2-hh;u(1)=0;u(n+2)=0;er=1;k=0;whileer>0.0001er=0;k=k+1;forj=2:n+1s=(hh*x(j)+u(j-1)+u(j+1))/r;er=max(er,abs(s-u(j)));u(j)=s;endendv=y(x);plot(x,v,x,u,'o')er=max(abs(u-v))—----y(xn);O----yn迭代格式:yj
=[xjh2+yj-1+yj+1]/(2–h2)12/16輸入:n=10計算結(jié)果:
er=0.0012第十二頁,共十六頁,編輯于2023年,星期三追趕法解三對角方程組n=input('inputn:=');h=1/(n+1);hh=h*h;a=-1;b=2-hh;x=0:h:1;f=hh*x;f(1)=[];f(n+1)=[];f(1)=f(1)/b;d=b;fori=2:nc(i-1)=a/d;d=b-a*c(i-1);f(i)=(f(i)-a*f(i-1))/d;endfori=n-1:-1:1f(i)=f(i)-c(i)*f(i+1);endy=sin(x)/sin(1)-x;yk=[0f0];plot(x,y,x,yk,'o')error=max(abs(yk-y))輸入:n=10計算結(jié)果:error=5.4566e-005O--------數(shù)值解;—--------精確解13/16第十三頁,共十六頁,編輯于2023年,星期三邊值問題差分方法:u(a)=
,u(b)=
一般二階方程:指數(shù)變換
1.高精度差分格式;2.組合外推方法;3.大型稀疏方程組解算子參考:
ElsevierSDOS(電子期刊)Appliedmathematicsandcomputation14/16第十四頁,共十六頁,編輯于2023年,星期三線性多步法(n=0,1,···)其中,xn+i=x0+(n+i)h,fn+i=f(xn+i,yn+i)局部載斷誤差A(yù)damas顯格式:yn+2=yn+1+h(3fn+1-fn)/2yn+3=yn+2+h(23fn+
溫馨提示
- 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)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 安徽省皖南地區(qū)2021-2022學(xué)年物理高一第二學(xué)期期末經(jīng)典試題含解析
- 安徽省合肥八中2021-2022學(xué)年物理高一下期末調(diào)研試題含解析
- 故事謎課件圖片
- 2022年新疆兵團(tuán)二師華山中學(xué)物理高一第二學(xué)期期末調(diào)研試題含解析
- 手工油紙傘課件
- 2024年黑色冶金項目提案報告
- 2024年工業(yè)物聯(lián)網(wǎng)項目規(guī)劃申請報告
- 關(guān)于毛巾的選擇需求調(diào)查問卷
- 2024年坐墊項目立項申請報告模稿
- 2024年鉭電解電容器項目申請報告模板
- 正心頌缽-頌缽療愈師培訓(xùn)
- 課后服務(wù)教師聘任協(xié)議書
- 王化成財務(wù)管理學(xué)講義
- 對建筑裝飾工程質(zhì)量檢驗與檢測這門課的心得
- DB31∕T 1205-2020 醫(yī)務(wù)社會工作基本服務(wù)規(guī)范
- 學(xué)生語文課堂學(xué)生學(xué)習(xí)評價量表
- 赤峰市松山區(qū)赤峰國維礦業(yè)有限公司車戶溝礦區(qū)銅鉬礦二〇二二年度礦山地質(zhì)環(huán)境治理計劃書
- GB/T 1972.2-2023碟形彈簧第2部分:技術(shù)條件
- 高中心理健康《拒絕精神內(nèi)耗》主題班會課件
- 摩托車科目一題庫-共400題(附答案)
- 俄羅斯探親邀請函范本3篇
評論
0/150
提交評論