Discussion :

[cozo22]

Bonjour,
je tente de visualiser en 3d cette fonction z suivant les deux parametres x et y mais jái un probleme au niveau de la fonction erreur (ERFC).il m indique

??? Error using ==> mtimes
Inner matrix dimensions must agree.

est ce que quelqu'un pourrait m'aider s'il vous plait....

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
 
clear all;
close all;
Q=188.64;
K=0.802;
w=0.0227861753262656;
tp=600;
Y0=20;
 
[x,y] = meshgrid(1 : 0.1 : 2, 1 : 1 : 600);
A=(x./(2*sqrt(w'*y)));
z= (2*Q./K)'*sqrt(w*y)'*[((exp(-A.^2))./(sqrt(pi)) - ((A)*erfc(A'))];
mesh(z)
figure(5)
surf(z)

[Jerome Briot]

Code :

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

1
2
 
z= (2*Q./K)'*sqrt(w*y)'*[((exp(-A.^2))./(sqrt(pi)) - ((A)*erfc(A'))];

Déjà, on peut simplifier l'expression :

Code :

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

z= (2*Q./K)'*sqrt(w*y)'*(exp(-A.^2)./sqrt(pi) - A*erfc(A'));

Et l'erreur est :

??? Error using ==> minus
Matrix dimensions must agree.

Error in ==> Untitled at 12
z= (2*Q./K)'*sqrt(w*y)'*(exp(-A.^2)./sqrt(pi) - A*erfc(A'));

Le problème vient du dernier terme :

Code :

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

(exp(-A.^2)./sqrt(pi) - A*erfc(A'))

Si on regarde la taille de chaque terme :

Taille de A : 600x11
Taille de erfc(A') : 11x600
Taille de A*erfc(A') : 600x600
Taille de exp(-A.^2)./sqrt(pi) : 600x11

Taille de (exp(-A.^2)./sqrt(pi) - A*erfc(A')) : ??? [600x11] - [600x600] ???

[cozo22]

Ok mais alors pourquoi est ce que lorsaue je tape ce code ca fonctione parce que les matrices sont les memes....

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
 
clear all;
close all;
Q=188.64;
K=0.802;
w=0.0227861753262656;
tp=600;
Y0=20;
 
x=00.0000000010987;
y=zeros(1,1200);
 
for t=1:1:600;
   A=(x./(2*sqrt(w*t)));
   y(t)= Y0 + (2*Q*(1)./K)'*sqrt(w*t)'*[((exp(-A.^2))./sqrt(pi))' - A'*erfc(A)'];
 end 
 
plot(1:600,y(1:600));
title('Temperature evolution With DL 25 A 11.25W- Fe           depth')

[Jerome Briot]

Les dimensions de A ne sont pas les mêmes entre les deux codes.
Dans le premier, A est une tableau 600x11 (une matrice), dans le second, A est une tableau 1x1 (un scalaire)

[cozo22]

il suffit donc que je donne une dimension 600x600 a A et normalement ca marche?

[Jerome Briot]

Et comment serions-nous censé connaitre la taille de A ?

Plutôt que de bidouiller A, tu ferais mieux de contrôler la formule qui calcul z.

Une fois tu fais :

Code :

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

A*erfc(A)'

L'autre fois tu fais :

Code :

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

A'*erfc(A)'

Ceci n'a pas de sens