Discussion :

[vampirerose]

Bonjour,
j'ouvre ce post car j'ai un problème dans un calcul pour l'obtention des coef d'une droite de régression linéaire.

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
rpm=num(:,248);
v=num(:,249);
rap=num(:,250);
v(isnan(v)) = [];
n=length(v); % longueur des vecteurs
% Recomposition du niveau de puissance => échantillonnage grossier
for i=1:n
    Lpneu(i)=Apneu+Bpneu*log10(v(i)/Vref)+25.6; % composante pneu
end

Lpneu(isnan(Lpneu)) = [];
ar=[1 1000]; % initialisation 

options=optimset('display','off','tolfun',.1); % options de la fct fminsearch (pas d'affichage, tolérance sur la fonction)
A=fminsearch(@fct_minimisation,[ar],options,rpm,RPMref,Lpneu)

J'ai créé pour cela une fonction :

Code :

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

1
2
3
4
function b=fct_minimisation(A,rpm,RPMref,Lpneu)
L_roulement=A(2)+A(1)*log10(rpm/RPMref);
b=sum(sum((Lpneu-L_roulement).^2));

Ce que je souhaite, c'est à partir des bruit de roulement Lpneu calculés, retrouver la droite de régression en fonction du régime moteur.

Deux problèmes ?

- comment fixer le paramètre ar

- l'erreur que j'obtiens :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
??? Error using ==> minus
Matrix dimensions must agree.
 
Error in ==> fct_minimisation at 5
b=sum(sum((Lpneu-L_roulement).^2));
 
Error in ==> fminsearch at 195
fv(:,1) = funfcn(x,varargin{:});
 
Error in ==> lecture1 at 1342
A=fminsearch(@fct_minimisation,[ar],options,rpm,RPMref,Lpneu)

Merci d'avance, si vous avez des conseils pour résoudre ce problème.

[Jerome Briot]

Les vecteur Lpneu et L_roulement n'ont pas le même nombre d’éléments.

A priori, si tu supprimes les valeurs NaN de Lpneu, il faut aussi supprimer les valeurs correspondantes dans rpm :

Code :

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

1
2
3
idx = isnan(Lpneu);
Lpneu(idx) = [];
rpm(idx) = [];

Non ?

[vampirerose]

J'avais effectivement trouver cette erreur mais matlab me communique tjs la même erreur

[vampirerose]

J'ai trouvé mon problème.

Code :

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

1
2
3
function b=fct_minimisation(A,rpm,RPMref,Lpneu,i)
L_roulement(i)=A(2)+A(1)*log10(rpm(i)/RPMref);
b=sum(sum((Lpneu-L_roulement).^2));

[vampirerose]

Je re-ouvre ce post. En effet j'essaie de simplifier mon problème.
Ma question est la suivante :

peut-on calculer une courbe de tendance sans passer par une toolbox ?
Et si oui ? comment ?

Voyez vous une différence entre une régression linéaire et une courbe de tendance ?

Je précise que je souhaite avoir une courbe de tendance pour des nuages de points avec un abscisse en log

[Jerome Briot]

peut-on calculer une courbe de tendance sans passer par une toolbox ?

=> Ajustement aux moindres carrés de courbes et de surfaces

Voyez vous une différence entre une régression linéaire et une courbe de tendance ?

Une "courbe de tendance" n'est pas forcément linéaire