Discussion :

[danaee]

Bonjour,

Je viens à vous car cela fait 2 jours que je tente de résoudre mon problème...sans résultat.
Je découvre le logiciel matlab depuis peu. Votre aide me serait bien précieuse pour pouvoir avancer dans mes recherches.
Voici mon code:

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
22
23
24
25
26
27
28
29
30
31
32
33
34
%debit du fluide à la paroi en fonction du temps
function[Q,t]=d(K,phi,mu,bf,br,b,psi,a,g,tx)

K=10e-15;
phi=0.1;
mu=10e-3;
bf=5*10e-10;
br=1*10e-10;
b=1.;
psi=10e+6;
a=0.1524;
g=0.57722;
% K,phi, mu....que des constantes

t=(0:.01:60);
%on fait varier t de 0 a 60 par saut de 0.1

tx=(K*t)/(mu*(a^2)*((b-phi)*br+(phi*bf)));
% tx ne depends que de t

Q=((2*pi*K*psi)/mu)*((1./(log(4.*tx)-(2*g))-(g/((log(4.*tx)-(2*g))^2))));

plot(t,Q);
xlabel('Temps (min)');
ylabel('Debit (litre/min)');

%if tx<1.275
%Q=((2*pi*K*psi)/mu)*((1/sqrt(pi*tx))+(1/2)-((1/4)*sqrt(tx/pi)+(tx/8)))
%elseif 1.275<tx<3
%Q=((2*pi*K*psi)/mu)*((1/(log(4*tx)-2*g)-(g/((log(4*tx)-(2*g))^2))));
%else
%Q=((2*pi*K*psi)/mu)*(6.97323-(0.54196*tx));
%end;

J'ai d'abord essayé d'obtenir le résultat pour une des équations avant de faire tourner la boucle en if, mais matlab m'affiche la reponse suivante:

Code :

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

1
2
3
4
5
??? Error using ==> mpower
Matrix must be square.

Error in ==> d at 24
Q=((2*pi*K*psi)/mu)*((1./(log(4.*tx)-(2*g))-(g/((log(4.*tx)-(2*g))^2))));

voila,
merci bien

[tug83]

Code :

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

 Q=((2*pi*K*psi)/mu)*((1./(log(4.*tx)-(2*g))-(g./((log(4.*tx)-(2*g)).^2))))

Le plus simple c'est de décomposer cette ligne , exemple:

Code :

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

1
2
3
4
5
6
A=((2*pi*K*psi)/mu);
B1 =1./(log(4.*tx)-(2*g));
B2 = log(4.*tx);
B3 =B2-(2*g);
B =B1-(g./((B3).^2));
Q=A*B

[danaee]

Merci bcp!!
je n'avais pas essayé le point après les parenthèses