



全文預(yù)覽已結(jié)束
下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
微分方程數(shù)值解報(bào)告三 迎風(fēng)格式和Lax-Friedrichs格式求解微分方程一 、引言雙曲型方程差分格式的性質(zhì)和定解問題解析解的性質(zhì)之間有著密切的聯(lián)系,由于其對初值的局部依賴關(guān)系和特征關(guān)系是其他兩類方程所沒有的,其初值函數(shù)的一些性質(zhì)也會沿特征線傳播,從而使解不具有光滑性質(zhì),在構(gòu)造雙曲型方程的差分格式時應(yīng)充分考慮這些特性,下面江將用兩種常見格式求解一個簡單的偏微分方程。二、方法原理本題求解首先用到一階迎風(fēng)格式(Upwind Scheme)其中定義:另外,為了克服中心差分格式的不穩(wěn)定性,還有Lax-Friedrichs格式:三、求解問題四、編程設(shè)計(jì)和方法分析考慮網(wǎng)格參數(shù):接下來我們按照迎風(fēng)格式和Lax-Friedrichs格式分別對方程進(jìn)行離散處理迎風(fēng)格式:Lax-Fredrichs格式為:迎風(fēng)格式代碼:%迎風(fēng)格式clear allclose allclc%設(shè)置計(jì)算網(wǎng)格的參數(shù)L=1.0;%計(jì)算區(qū)域的長度T=0.5;%計(jì)算的時間長度dx=0.01;dt=0.005;x=-L:dx:L;t=0:dt:T;xn=size(x,2)-1;%網(wǎng)格段的個數(shù)tn=size(t,2)-1;%時間段的個數(shù)u=zeros(tn+1,xn+1);%設(shè)置初始條件for xi=1:xn+1if(x(xi)=0)u(1,xi)=1;%x0時,uo(x)=0endend%在計(jì)算開始前,將中間的網(wǎng)格點(diǎn)的值設(shè)為零for i=2:xnu(i,1)=1;%將中間的網(wǎng)格點(diǎn)的值設(shè)置為0endfor ti=1:tnfor j=2:xnu(ti+1,j)=dt*u(ti,j-1)/dx+(1.0-dt/dx)*u(ti,j);endend%繪制計(jì)算結(jié)果for i=2:30:92subplot(2,2,1+(i-2)/30);plot(x,u(i,:);str=num2str(t(i),T=%f);xlabel(str);str=num2str(i,u(%d,x);ylabel(str);endLax-Friedrichs格式代碼2%Lax-Friedrichs格式clear allclose allclc%設(shè)置計(jì)算網(wǎng)格的參數(shù)L=1.0;%計(jì)算區(qū)域的長度T=0.8;%計(jì)算的時間長度dx=0.01;%網(wǎng)格尺寸dt=0.005;%時間步長miu=1;x=-L:dx:L;t=0:dt:T;xn=size(x,2)-1;%x的分段數(shù)tn=size(t,2)-1;%時間的分段數(shù)u=zeros(tn+1,xn+1);%u的上標(biāo)為時間標(biāo)號,下標(biāo)為坐標(biāo)標(biāo)號%設(shè)置初始條件for xi=1:xn+1if(x(xi)=0)u(1,xi)=1;%x0時,uo(x)=0endendfor i=2:xnu(i,1)=1;endfor ti=1:tnfor j=2:xnu(ti+1,j)=0.5*(1-dt/dx)*u(ti,j+1)+0.5*(1+dt/dx)*u(ti,j-1);endend%繪制計(jì)算結(jié)果for i=2:30:92subplot(2,2,1+(i-2)/30);plot(x,u(i,:);str=nu
溫馨提示
- 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)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 石家莊市第四十中學(xué)2025年高一下化學(xué)期末檢測模擬試題含解析
- 2025屆山東省新泰第一中學(xué)北?;瘜W(xué)高一下期末學(xué)業(yè)質(zhì)量監(jiān)測試題含解析
- 2025屆浙江省浙南聯(lián)盟化學(xué)高一下期末教學(xué)質(zhì)量檢測試題含解析
- 校園讀物日常管理辦法
- 民工衣服庫存管理辦法
- 晚會捐贈收入管理辦法
- 冬季水管防護(hù)管理辦法
- 合肥苗木采伐管理辦法
- 現(xiàn)代語文教育數(shù)字化資源整合平臺建設(shè)研究
- 小企業(yè)會計(jì)準(zhǔn)則的實(shí)務(wù)操作與案例分析
- 2024年遼寧沈陽水務(wù)集團(tuán)有限公司招聘筆試真題
- 濰坊交通發(fā)展集團(tuán)有限公司招聘筆試題庫2025
- 胸痛中心質(zhì)控管理
- 2025時政試題及答案(100題)
- 第七章城市軌道交通屏蔽門設(shè)備接口68課件
- 國家開放大學(xué)漢語言文學(xué)本科《中國現(xiàn)代文學(xué)專題》期末紙質(zhì)考試第三大題分析題庫2025春期版
- 成都大學(xué)附屬中學(xué)英語新初一分班試卷含答案
- 新22J01 工程做法圖集
- 創(chuàng)新創(chuàng)業(yè)大賽項(xiàng)目商業(yè)計(jì)劃書模板
- 2025年1月國家開放大學(xué)漢語言文學(xué)本科《心理學(xué)》期末紙質(zhì)考試試題及答案
- 糖尿病酮癥酸中毒疑難病例護(hù)理
評論
0/150
提交評論