1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38
| clear all; clc;
borne =10000;
%% temps
Temps=1:1:1000;
Time=[];
for i=1:length(Temps)
Time=[Time;Temps];
end
%% epaisseur
L=linspace(0.00018,0.5,1000);
epaisseur=[];
for i=1:length(L)
epaisseur=[epaisseur;L];
end
%%
alpha=5e-07;
S=zeros(length(L),length(Temps));
for e=1:length(L)
Somme=zeros(1,length(Temps));
ep=L(e);
for t=1:length(Temps)
interv=-borne:1:borne;
vect=zeros(length(interv),1);
for k=1:length(interv)
n=interv(k);
vect(k)= exp((-1)*(((n-1)*ep)^2)/(alpha*Temps(t)));
end
Somme(t)=sum(vect);
end
S(e,:)=Somme;
end
%%
surf(Time,epaisseur',S)
colorbar |
Partager