一維非穩(wěn)態(tài)導熱方程求解(附Matlab程序)_第1頁
一維非穩(wěn)態(tài)導熱方程求解(附Matlab程序)_第2頁
一維非穩(wěn)態(tài)導熱方程求解(附Matlab程序)_第3頁
全文預覽已結束

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權，請進行舉報或認領

文檔簡介

1、使用差分方法求解下面的熱傳導方程l；(x,t)=a3rIx(x,t)(0 xl,0t0.2,a=1)初值條件：T(x,0)=4x-4x2；邊值條件:T(O,t)=0T(l,t)=O+O(h2)一2幾+1擊1?使用差分公式TOqh)2丁(務)+丁(務+1】占)KKh2上而兩式帶入原熱傳導方程4k令r=,化簡上式的lrT,j-h=(1-2rJRj+r(K-i,)+Xr,j)如下圖：(i(1-2嘰it=meshgrid(x,t),mesh(x,t,T),xlabelCx1)ylabelCOzlabel(T)functionUNt=Techuandao(cgl,g2xspan,tspan,ngnd)%

2、熱傳導方程：%Ut(x,t)=cA2*Uxx(xt)axbtstgl(t)%u(b,tg2(t)%參數(shù)說明%c：方程中的系數(shù)%f：初值條件%gl,g2：邊值條件%xspan=a,b：x的取值范圍%tspan=ts,tf：t的取值范圍%ngnd=n,m：網(wǎng)格數(shù)帚,m為x網(wǎng)格點數(shù)量，n為t的網(wǎng)格點數(shù)量%U：方程的數(shù)值解%x,t：x和t的網(wǎng)格點n=ngiid(l);mF=ngnd(2),hange(xspan)/(m-1),x=hnspace(xspan(l),xspan(2),m);k=ange(tspan)/(n-1),t=linspace(tspan(l),tspan(2),n),心k/h=ifr0.5errorf為了保證算法的收斂，請増大步長h或減小步長k!1)end5=1-2*!；U=zeips(ngnd),%邊界條件u(：,i)=gi(0；U(：)m)=g2(t),%初值條件U(l)=f(x),

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網(wǎng)頁內(nèi)容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間，僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理，對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯，并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容，請與我們聯(lián)系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論