postdf2d.m

来自「２Ｄ力学模型.有限单元法.巴黎六大实验教材.」· M 代码 · 共 37 行

%*****************************************************************************
%  PostDF2D.m




function []=PostDF2D(Tresult,Nx,Ny,nx,ny,L,H);

%1)Affichage de la solution calculee
%=======================================================================

figure
plot(Nx,Ny,Tresult,'ko');

%2) Affichage de la solution exacte
%=========================================================================

for i=1:nx & j=1:ny
    Tex(i,j)=A*exp((-B)*((r0)^2));
    
end

hold on
plot(Nx,Ny,Tex,'k-');

legend('t^h','t^(ex)')

xlabel('x','Fontsize',12);
ylabel('T^h(x)','Fontsize',12);


%3=Calcul de l'erreur en norme L^(inf)
%=========================================================================

Err=Tex-Tresult';
err=max(abs(Err))