MT94 - Introduction ...

Bonjour,

avec un peu de retard, voici la correction. Décommentez les définitions de f,m,a pour chaque fonction

N=200;
k=1:N;
x=-%pi:%pi/50:%pi;

f = 0+(abs(x)<1)
m = 1/%pi*ones(x);
a = 2*sin(k)./k/%pi;

//f = (abs(x)<1).*(1-abs(x))
//m = 1/2/%pi*ones(x);
//a = 2*(1-cos(k))./k.^2/%pi;

//f = (abs(x)<1).*(1-3*x.^2+2*x.^2.*abs(x))
//m=1/2/%pi;
//a = -12*(k.*sin(k)+2*cos(k)-2)./k.^4/%pi;

h=1e-1;
t=0:h:4;
u = zeros(length(x),length(t));
u(:,1) = f;

clf
subplot(1,2,1)
hdl = plot(x,f);
gca().data_bounds(3:4)=[-0.2 1.2];
messagebox("Cliquer sur OK pour demarrer","OK","","","modal")
for i = 2:length(t)
     s = m;
     for k=1:N
            s = s+a(k)*cos(k*x)*exp(-k^2*t(i));
     end
     u(:,i) = s;
     hdl.data(:,2)=s;
     sleep(100)
end

subplot(1,2,2)
mesh(t,x,u)
xlabel("t")
ylabel("x")

Cordialement,

S.