Corrigés desexercices
87
e=zeros(M+1,1);
e(M+1)=1;
C=A*(A’*A)^(-1)*e;
n=1000;
[x,yi,yb]=signal2(n);
%F iltrage
yl=yb;
for i=NL+1:n-NR
yl(i)=C’*yb(i-NL:i+NR)’;
end;
%G raphes
subplot(2,2,1)
plot(x,yi)
axis([0 1- 0.1 1.5])
title(’signalinitial’)
subplot(2,2,2)
plot(x,yb)
axis([0 1- 0.1 1.5])
title(’signalbruite’)
subplot (2,2,3)
plot(x,yl)
axis([0 1- 0.1 1.5])
title(’signallisse’)
subplot(2,2,4)
plot(x,yi,x,yl)
axis([0 1- 0.1 1.5])
title(’signalinitial+signallisse’)
✝
✆
Ce qu’il faut retenir de cetexercice
Cettefois-ci une boucle est utilisée pour construire A .Elle évite l’instruction puissance.
6.5 Équation différentielle
EQD1234.m
✞
☎
M=4;
n=6;
h=1/(n+1);
x=(h:h:1-h)’;
A1=x*ones(1,M);
for j=1:M
A2(:,j)=-(j-1)*(j-2)*A1(:,j).^(j-3)+A1(:,j).^(j-1);
end
A=[ zeros(1,M);A2;ones(1,M)];
A(1,1)=1;
b=[0;f(x);0];
alpha=(A’*A)\ (A’*b);
t=0:1/1000:1;
e=exp(1);
uex=e/(2*(e-1/e))*( exp(t)- exp(-t))-t/2.*exp(t);
uapp= polyval(alpha(M:-1:1),t);
plot(t,uex,t,uapp)
erreur= norm(uex-uapp,inf)
✝
✆
 Dunod –Laphotocopie non autorisée est un délit