Corrigés desproblèmes
193
plaque.m
✞
☎
clear, close all
n=10
[x,y,ue,uh]=diffini(n);
er= abs(ue-uh);
disp(’||solex-solapp||’)
err= max( max(er))
disp(’erreur*(n+1)^2’)
err*(n+1)^2
subplot(1,2,1)
surf(x,y,uh)
title(’Solutionapprochee’,’Fontsize’,14)
axis([0 101-1.1 1])
subplot(1,2,2)
surf(x,y,ue)
title(’Solutionexacte’,’Fontsize’,14)
axis([0 101-1.1 1])
figure
surf(x,y,uh-ue)
title(’Solapp-Solexacte’,’Fontsize’,20)
figure
surf(x,y,uh,’LineStyle’,’none’,’FaceColor’,’red’);
camlight left;l ighting phong
title(’Solutionapprochee’,’Fontsize’,16)
✝
✆
Ce qu’il faut retenir de cetexercice
Noterles instructions qui donnent un lissage des solutions :
surf(x,y,uh,’LineStyle’,’none’,’FaceColor’,’red’) ;
camlight left ;l ighting phong
etuderreur.m
✞
☎
tN=[2 5710 15 20 38 75];
for k=1:length(tN)
n=tN(k)
[x,y,ue,uh]=diffini(n);
er= abs(ue-uh);
taberreur(k)=max( max(er));
taberreur1(k)=taberreur(k)*(n+1)^2;
end
tN,taberreur
disp(’erreur*(N+1)^2’)
taberreur1
a=polyfit( log(tN+1), log(taberreur),1);
plot( log(tN+1), log(taberreur));
title([’pentedeelaadroitedeeregression::’...
, num2str(a(1))],’FontSize’,16)
xlabel(’log(N+1)’,’FontSize’,16)
ylabel(’log(Erreur)’,’FontSize’,16)
✝
✆
 Dunod –Laphotocopie non autorisée est un délit