IDENTIFICATION par moindre carrés.Matlab
Fermé
rahmalse
Messages postés
5
Date d'inscription
dimanche 11 décembre 2011
Statut
Membre
Dernière intervention
8 avril 2013
-
10 janv. 2012 à 01:59
ouzi sim Messages postés 6 Date d'inscription lundi 31 décembre 2012 Statut Membre Dernière intervention 1 février 2013 - 1 févr. 2013 à 23:53
ouzi sim Messages postés 6 Date d'inscription lundi 31 décembre 2012 Statut Membre Dernière intervention 1 février 2013 - 1 févr. 2013 à 23:53
A voir également:
- Moindre carré matlab
- Au carré iphone ✓ - Forum iPhone
- Point d'interrogation dans un carré ✓ - Forum Samsung
- Parenthese carré ✓ - Forum MacOS
- Carré en c ✓ - Forum Programmation
- Symbole racine carré clavier téléphone ✓ - Forum iPhone
1 réponse
ss10
Messages postés
5
Date d'inscription
lundi 2 janvier 2012
Statut
Membre
Dernière intervention
4 février 2012
3
11 janv. 2012 à 22:22
11 janv. 2012 à 22:22
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.
12 janv. 2012 à 18:22
1 févr. 2013 à 23:51
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(:);
1 févr. 2013 à 23:52
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');
1 févr. 2013 à 23:53