etude_diffusion_1d.m

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

clear all
close all

% resultats

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

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

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

load('erreur_1d_upwind_galerkin')
erreur_upwind_galerk=erreur_dif;
t_lp=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('diffusion des trois sch閙as')
 title('diffusion avec upwind')
 xlabel('temps');
 ylabel('erreur diffusion');
a=plot(tt_up,erreur_upwind,'r*');
b=plot(tt_lax,erreur_lax,'b*');
c=plot(tt_lp,erreur_lp,'g*');
d=plot(tt_galerk,erreur_upwind_galerk,'k*');

%legend([a,d],'diffusion upwind','diffusion up galerkin')

legend([a,b,c,d],'diffusion upwind','diffusion lax','diffusion lp','diffusion up galerkin')