Discussion :

[nem007]

Bonjour,
J'ai réalisé un programme sous Matlab et je voudrais le faire tourner dans LV par l'intermédiaire d'un nœud Mathscript. Or dans mon programme j'utilise la fonction refline, mais Labview ne la comprend pas. La fonction refline crée une droite de référence par rapport a un nuage de points.
Y aurai-t-il un moyen de contourner ceci avec une autre fonction?

Merci d'avance pour vos réponses.

Sinon pour ce que ça intéresse, j'ai réalisé un programme "loi de Weibull" qui calcul la pente béta (critère de forme de la loi) et le point B10 (durée de vie de la production) à partir d'un échantillon donné.
Programme :

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
a=80000;
b=87000;
c=89000;
d=81000;
e=83000;
f=94000;
r=[a,b,c,d,e,f];%Vecteur échantillon
r1 = sort(r);%Vecteur échantillon classé dans l'ordre croissant
l=length(r)%Longueur du vecteur échantillon
for i = 1:l
    A(i)=((i-0.3)/(l+0.4));
    A100(i) = A(i)*100;
    B(i)=1/(1-A(i));
    C(i)=log(log(B(i)));
end
D=log(r1)
hold on
plot(D,C,'or')
E=polyfit(D,C,1)
beta=E(1)%critère de forme
B10=exp(-E(2)/E(1))%durée de vie 
refline(E)

A+
Greg

[Jerome Briot]

Y aurai-t-il un moyen de contourner ceci avec une autre fonction?

Vu la documentation de REFLINE, je ne vois pas où se situe la difficulté de la remplacer tout simplement par LINE ...

[nem007]

Salut Dut, merci pour ta réponse,
Dans mon cas la fonction LINE ne me trace pas une droite qui passe par un max de points de mon nuage, ce que fait REFLINE....
Si vous avez d'autres suggestions à me proposer, je suis preneur

[Jerome Briot]

Tu pourrais nous montrer un aperçu de ton nuage de points et de la ligne tracée par REFLINE sous MATLAB ?

[nem007]

Ci joint mon programme et la figure nuage de point et droite d'interpolation.

[Jerome Briot]

Vu que tu as trouvé POLYFIT... tu devrais être en mesure d'utiliser POLYVAL avec la variable E.

Enfin c'est quand même juste un polynôme de degré 1... une régression linéaire c'est pas le bout du monde quand même :

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
xy = [11.29 -2.2
    11.305 -1.2
    11.33 -0.6
    11.37 -0.25
    11.4 0.25
    11.45 0.75];

A = [xy(:,1) ones(size(xy,1),1)];
B = xy(:,2);
C = A\B;

figure
plot(xy(:,1),xy(:,2),'ro')
hold on

x = xlim;
y = C(1)*x+C(2);
plot(x,y,'b-')