etude_stabilite_1d.m

来自「１Ｄ非连续性力学分析」· M 代码 · 共 35 行

clear all
close all

% resultats

load('erreur_1d_upwind')
erreur_upwind=erreur_x;
t_up=tmp;

load('erreur_1d_lax')
erreur_lax=erreur_x;
t_lax=tmp;

load('erreur_1d_leap_frog')
erreur_lp=erreur_x;
t_lp=tmp;

load('erreur_1d_upwind_galerkin')
erreur_upwind_galerk=erreur_x;
t_up=tmp;

for n=1:t_up
    tt_up(n)=n;
end
for n=1:t_lax
    tt_lax(n)=n;
end
for n=1:t_lp
    tt_lp(n)=n;
end
for n=1:t_lp
    tt_galerk(n)=n;
end
 hold on
 %title('stabilit