Discussion :

[conanran]

Bonjour,

Je reviens vers vous encore une fois,
En effet, je tracer les courbes sur la figure 1 en utilisant l’équation à côté.
J’ai essayé mais je n’ai pas obtenu la même forme.
Merci pour votre aide si vous avez des propositions dans le petit code.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
x=-4:0.1:4;
 
for t=1:100
    for n=1:size(x,2)
    T(t,:)=(0.5/sqrt(pi*t))*exp((-x(n)^2)/4*t);
    end
end
plot(x,T(10,:))

[Jerome Briot]

Attention ici :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

exp((-x(n)^2)/(4*t))

Par contre, je ne comprends pas bien pourquoi t varie de 1 à 100 dans ton code ?

Ne serait-ce pas plutôt ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
x = -4:0.1:4;
t = [1/16 1/4 1];
 
for nt = 1:numel(t)
    for nx = 1:numel(x)
        T(nx,nt) = 0.5/sqrt(pi*t(nt))*exp((-x(nx)^2)/(4*t(nt)));
    end
end
 
figure
plot(x, T)

[conanran]

Merci Jerome Briot
oui une erreur de ma part dans la définition de l'intervalle de t.
Merci beaucoup....

[conanran]

En effet, de la même façon j’ai essayé de faire l’équation sur la figure 2, sauf que au lieu de x, je vais raisonner sur un rayon r

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
clear; close all; clc;
%%
[X,Y] = meshgrid(0:0.01:1,0:0.01:1);
R=0.03; 
Qmax=1000;
x_0=0.5;
y_0=0.5;
r=(X-y_0).^2+(Y-x_0).^2;
Q=Qmax*exp(-2*((r^2))/(R^2)); 
a=4.3812e-07;                 
b=2.2975e-05;
%%
t= 0.01:0.01:1;
for nt = 1:numel(t)
Tm=(2.*Q)./((b*sqrt((pi^3)*t(nt)))*((R^2)+8*a*t(nt)));
 T(:,:,nt)=Tm*exp(((-2.*(r.^2))./((R^2)+8.*a.*t(nt)))); 
 
end

normalement j'obtiens des courbes semblables...
Merci pour votre aide

Merci pour votre aide

[Jerome Briot]

Je ne suis pas sûr de bien comprendre mais pourquoi ne pas faire ceci ?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
R = 0.03;
Qmax = 1000;
 
a = 4.3812e-07;
b = 2.2975e-05;
 
r = linspace(0, R, 100);
t = [1/16 1/4 1];
 
for nr = 1:numel(r)
    Q = Qmax*exp(-2*r(nr)^2/R^2);
    for nt = 1:numel(t)
        Tm = 2*Q/((b*sqrt(pi^3*t(nt)))*(R^2+8*a*t(nt)));
        T(nr,1,nt) = Tm*exp((-2*r(nr)^2)/(R^2+8*a*t(nt)));
    end
end
 
T = squeeze(T);
 
figure
plot(r, T)