Probleme avec une courbe 3D sur Matlab

Version imprimable

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

Citation:

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

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

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)

06/03/2008, 12h35
Jerome Briot
Citation:
Envoyé par cozo22

Code:

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:

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

Et l'erreur est :

Citation:

??? 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:

(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] ???

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

Code:

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')

06/03/2008, 12h47
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)
06/03/2008, 12h48
cozo22

il suffit donc que je donne une dimension 600x600 a A et normalement ca marche?
06/03/2008, 12h54
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:

A*erfc(A)'

L'autre fois tu fais :

Code:

A'*erfc(A)'

Ceci n'a pas de sens :?