MATLAB編輯一維熱傳導(dǎo)方程的模擬程序_第1頁
MATLAB編輯一維熱傳導(dǎo)方程的模擬程序_第2頁
MATLAB編輯一維熱傳導(dǎo)方程的模擬程序_第3頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明：本文檔由用戶提供并上傳，收益歸屬內(nèi)容提供方，若內(nèi)容存在侵權(quán)，請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、求解下列熱傳導(dǎo)問題：1=0(0EZEL)g2豆境_2T(z,0)=1-zT(0,t)=1,T(L,t)=0L=1,：-1程序：functionheat_conduction()%一維齊次熱傳導(dǎo)方程options=,空間桿長L,空間點(diǎn)數(shù)N,時(shí)間點(diǎn)數(shù)M,擴(kuò)散系數(shù)alfa,穩(wěn)定條件的值lambda(取值必須小于0.5),;topic=seting;lines=1;def=1,100,1000,1,0.5;h=inputdlg(options,topic,lines,def);L=eval(h1);N=eval(h2);M=eval(h3);alfa=eval(h4);lambda=eval(h5);

2、%lambda的值必須小于0.5%*h=L/N;%空間步長z=0:h:L;z=z;tao=lambda*hA2/alfa;%時(shí)間步長tm=M*tao;%熱傳導(dǎo)的總時(shí)間tmt=0:tao:tm;t=t;%計(jì)算初值和邊值T=zeros(N+1,M+1);Ti=init_fun(z);To=border_funo(t);Te=border_fune(t);T(:,1)=Ti;T(1,:)=To;T(N+1,:)=Te;%用差分法求出溫度T與才f長L、時(shí)間t的關(guān)系fork=1:Mm=2;whilem=NT(m,k+1)=lambda*(T(m+1,k)+T(m-1,k)+(-2*lambda+1)*T(m,k);m=m+1;end;end;%設(shè)置立體網(wǎng)格fori=1:M+1X(:,i)=z;end;forj=1:N+1Y(j,：)=t；endmesh(X,Y,T);view(1-11);xlabel(Z);ylabel(t);zlabel(T);functiony=init_fun(z)%初值條件y=1-z.A2;returnfunctiony=border_funo(t)%z=0的邊界條件y=1+t.*0;returnfunctiony=border_fune(t)%z=L的邊界條件y=t*.0;return運(yùn)行情況：按“run”運(yùn)行時(shí)，彈出窗

溫馨提示

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

評(píng)論

0/150

提交評(píng)論