5阶WENO有限差分法求解Sod型激波管问题

利用5阶WENO有限差分法求解一维Euler方程组，求解Sod型激波管问题的详细程序。

function [LRLambdaPLambdaNA]=Count1DLR（Ur）

for i=1:size（U2）
    rho=U（1i）; E=U（3i）; u=U（2i）/U（1i）;
    p=（r-1）*（E-rho.*u^2/2）; c=sqrt（r*p/rho）;
    L{i}=（r-1）/（2*c^2）*[u^2/2+u*c/（r-1） -u-c/（r-1） 1; -u^2+2*c^2/（r-1） 2*u -2; u^2/2-u*c/（r-1） -u+c/（r-1） 1];
    R{i}=[1 1 1; u-c u u+c; c^2/（r-1）+u^2/2-u*c u^2/2 c^2/（r-1）+u^2/2+u*c];
    Lambda01=u-c; Lambda02=u; Lambda03=u+c; Lambda=1.1*max（abs（[Lambda01 Lambda02 Lambda03]））;
    LambdaP（:i）=[Lambda01+Lambda; Lambda02+Lambda; Lambda03+Lambda]/2;
    LambdaN（:i）=[Lambda01-Lambda; Lambda02-Lambda; Lambda03-Lambda]/2;
    A{i}=[0 1 0; （r-3）*u^2/2 （3-r）*u r-1; -u*c^2/（r-1）+（r-2）*u^3/2 c^2/（r-1）-（r-3/2）*u^2 r*u];
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        689  2014-02-06 09:57  Sod\Count1DLR.m

     文件       3199  2014-02-06 13:48  Sod\SodProgram.m

     文件       1498  2014-02-06 09:59  Sod\WENO51DEulerLu.m

     目录          0  2014-02-07 10:59  Sod

----------- ---------  ---------- -----  ----

                 5386                    4