Fit sur un système d'équations linéaires avec paramètres fixes

Bonjour. J'aimerais faire un fit linéaire, mais en gardant des paramètres fixes. Voilà comment je m'y prenais auparavant, avec la méthode classique utilisant l'opérateur \ :

Code:

---

1 2 3 4 5

x=rand(3); b=[1;2;3]; y=x*b; a=x\y

---

On voit que a est bien égal à b. Ce que je voudrais faire maintenant, c'est lui dire que je veux que le deuxième paramètre de a soit égale à 2.5 par exemple et qu'il fasse le fit sur les autres paramètres. Problème, je ne sais pas comment m'y prendre. Quelqu'un à une idée? Merci!!!

Salut,

essayes ceci :

Code:

---

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

%---------------------------------- % Fit   % Normal x=rand(3); b=[1;2;3]; y=x*b; a1=x\y   % New a2_2 = 2.5; ix = [1 3]; nix = 2;   x2 = x(:,ix); y2 = y - x(:,nix)*a2_2; % devrit marche pour a2_2 vecteur ...   atmp = x2\y2;   a2 = zeros(size(a1)); a2(ix) = atmp; a2(nix) = a2_2;   plot(b,a1) hold on plot(b,a2,'r')

---

Normalement le même principe devrait marche si tu veux fixer plusieurs valeurs

++

Merci beaucoup mr_samurai!

Bon, maintentant, il faut que je l'incorpore dans le gros programme et que je vois si ça marche comme ce qui est esconté. Je vous tiens au courant.