資源簡介
在matlab中,利用迎風(fēng)離散格式求解對流方程的代碼。
給定初值和邊界條件,直接求出數(shù)值解。
代碼片段和文件信息
function?u=peYF(adtnminxmaxxM)?
%%方程中的常數(shù);時(shí)間步長;空間節(jié)點(diǎn)數(shù);求解區(qū)間的左端;求解區(qū)間的右端;迭代次數(shù)
format?long;
h=(maxx-minx)/(n-1);
if?a>0
???for?j=1:(n+M)
??????u0(j)=IniU(minx+(j-M-1)*h);%向左延拓M個(gè)節(jié)點(diǎn)的函數(shù)值
???end
else
???for?j=1:(n+M)
??????u0(j)?=?IniU(minx+(j-1)*h);%向左延拓M個(gè)節(jié)點(diǎn)的函數(shù)值
???end
end
u
評論
共有 條評論