pehypbmc.m

来自「Matlab语言常用算法程序集.rar 里面有很多有用的Matlab数值运」· M 代码 · 共 21 行

function u = peMC(a,dt,n,minx,maxx,M)
format long;
h = (maxx-minx)/(n-1);
for j=1:(n+2*M)
    u0(j) = IniU(minx+(j-M-1)*h);
end 
u1 = u0;

for k=1:M
    for i=k+1:n+2*M-k
        tmpU1 = -dt*a*(u0(i+1)-u0(i))/h+u0(i);
        tmpU2 = -dt*a*(u0(i)-u0(i-1))/h+u0(i-1);        
        u1(i) = -dt*a*(tmpU1-tmpU2)/h/2+(u0(i)+tmpU1)/2;
    end
    u0 = u1;
end

u = u1((M+1):(M+n));
format short;