IDENTIFICATION par moindre carrés.Matlab
rahmalse
Messages postés
5
Date d'inscription
Statut
Membre
Dernière intervention
-
ouzi sim Messages postés 6 Date d'inscription Statut Membre Dernière intervention -
ouzi sim Messages postés 6 Date d'inscription Statut Membre Dernière intervention -
Bonjour,
Votre aide svp, y a t'il un modèle d'identification de n'importe quel Système simple par exemple système de 1er ordre (RL) par Matlab en se basant sue la méthode des moindres carrés pour avoir une idée de ce genre d'identification et l'adapté à mon système de recherche les machines électriques.
Merciiii
Votre aide svp, y a t'il un modèle d'identification de n'importe quel Système simple par exemple système de 1er ordre (RL) par Matlab en se basant sue la méthode des moindres carrés pour avoir une idée de ce genre d'identification et l'adapté à mon système de recherche les machines électriques.
Merciiii
A voir également:
- Moindre carré matlab
- Signe au carré ✓ - Chromebook
- Parenthese carre ✓ - Forum Word
- Symbole accolade {} ✓ - Forum Clavier
- Metre carré clavier ✓ - Forum Windows
- Faire un ² sur clavier QWERTZ ✓ - Forum Bureautique
1 réponse
bn soir ,voici un exemple que nous avons fait en tp sur la méthode des moindres carrée.
x=[0.1 0.2 0.5 1 1.5 1.9 2 3 4 6];
y=[0.95 0.89 0.79 0.7 0.63 0.58 0.56 0.45 0.36 0.28];
n=length(x)
f=log(y)
phi=ones(n,2)
phi(:,1)=x'
teta=inv(phi'*phi)*phi'*f'
j'éspée que c'est ce que tu cherche
bonne chance.
x=[0.1 0.2 0.5 1 1.5 1.9 2 3 4 6];
y=[0.95 0.89 0.79 0.7 0.63 0.58 0.56 0.45 0.36 0.28];
n=length(x)
f=log(y)
phi=ones(n,2)
phi(:,1)=x'
teta=inv(phi'*phi)*phi'*f'
j'éspée que c'est ce que tu cherche
bonne chance.
function [e,w]=lms1(mu,M,u,d);
w=zeros(M,1);
N=length(u);
u=u(:);
d=d(:);
for n=M:N
uvec = u(n:-1:n-M+1);
e(n)=d(n)-w'*uvec;
w=w+mu*uvec*conj(e(n));
end
e=e(:);
fs=500;
l=256;
n=2000;
t=(1:n)/fs;
a=0.5;
w=(1:n)*4*pi/fs;
x=sawtooth(w,.5)';
subplot(311),plot(t,x,'k');
interfer=sin(w*2.33)';
x=x+interfer;
ref=.45*interfer;
subplot(312),plot(t,x,'k');
px=(1/(n+1))*sum(x.^2);
delta=(1/(10*l*px))*a;
[b,y]=lms1(delta,l,ref,x);
subplot(313),plot(t,b,'k');