Discussion :

[Newenda]

Bonjour,

Question: j'ai deux jeux de données que je trace avec plotyy

Code :

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

1
2
3
4
5
6
7
8
 
GWC = [54 58 62 69] ;
G = [11000 7000 3600 1600] ;
dens = [1715 1710 1616 1601] ;
Vs = sqrt(G./dens) 
figure
hold on
[AX,H1,H2] = plotyy(GWC, Vs, GWC, G)

et j'aimerais tracer également leurs régressions mais je ne vois pas comment en tracer d'autre dans la même figure.

Merci d'avance pour votre aide.

[Jerome Briot]

De quelle régression parles-tu ?
De régression linéaire ?
Comme 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
GWC = [54 58 62 69] ;
G = [11000 7000 3600 1600] ;
dens = [1715 1710 1616 1601] ;
Vs = sqrt(G./dens); 
figure
 
[AX,H1,H2] = plotyy(GWC, Vs, GWC, G);
 
hold(AX(1),'on')
hold(AX(2),'on')
 
x = xlim;
 
M = [ones(numel(GWC),1) GWC(:)]; 
 
T = M\Vs(:);
plot(AX(1),x,T(2)*x+T(1),'r--')
 
T = M\G(:);
plot(AX(2),x,T(2)*x+T(1),'m--')

[Newenda]

Regression c'est pas vraiement le problème, mais les syntaxes de hold et du plot ! merci !

[Newenda]

Hum, je pensais avoir pigé mais pourriez vous me dire pourquoi la ligne ne fonctionne pas (j'essaye de faire passer une fonction d'ordre2) :

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
GWC = [54 58 62 69] ;
G = [110000 70000 36000 16000] ;
dens = [1715 1710 1616 1601] ;
Vs = sqrt(G./dens) 
 
%% FIT
GWCfit=[54:0.01:70] ;
Gfit = 399.45*GWCfit.^2 - 55474*GWCfit + 1.9416e+006 ; %399.45*x^{2} - 55474*x + 1.9416e+006
Vsfit = 0.01092*GWCfit.^2 - 1.6709*GWCfit + 66.442 ;
 
figure
hold on
[AX,H1,H2] = plotyy(GWC, Vs, GWC, G)
hold(AX(1),'on')
hold(AX(2),'on')
 
plot(AX(1),GWCfit,Gfit,'g--')
 
plot(AX(2),GWCfit,Vsfit,'b--')

>> il trace mal la 2eme regression et même pas la 1ere... et en changeant les plots comme suis :

Code :

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

1
2
3
plot(GWCfit,Gfit,'g--')
 
plot(GWCfit,Vsfit,'b--')

> il me trace bien mais juste la dernière (en bleu)...

Comprends pas, merci de votre aide.

[Jerome Briot]

Tu as inversé les deux objet Axes :

Code :

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

1
2
3
4
plot(AX(2),GWCfit,Gfit,'g--')

plot(AX(1),GWCfit,Vsfit,'b--')

[Newenda]

hum...merci..