資源簡介
用迎風格式、Lax-Wendroff格式的差分方法、Lax-Friedrichs格式的差分方法來解對流方程,并且有詳細的matlab程序。
代碼片段和文件信息
%輸入a=1N=151M=100dt=0.5
x=-0.5:0.01:1;
for?i=1:151
????if?i<=101
????????w(i)=1;
????else
????????w(i)=0;
????end
end
%以上程序的操作是精確解。
a=input(‘請輸入a=‘)
N=input(‘請輸入N=‘)
M=input(‘請輸入M=‘)
dt=input(‘請輸入dt=‘)
for?j=1:N
????u0(j)=-0.5+0.01*(j-1)
????if?u0(j)<=0
?????????u0(j)=1;
評論
共有 條評論