Probleme algorithme des moindres sous matlab
Fermé
zmeryem
Messages postés
3
Date d'inscription
dimanche 13 janvier 2008
Statut
Membre
Dernière intervention
14 janvier 2008
-
13 janv. 2008 à 22:05
Mouhamed Aoun - 16 avril 2012 à 10:16
Mouhamed Aoun - 16 avril 2012 à 10:16
A voir également:
- Probleme algorithme des moindres sous matlab
- Logiciel algorithme gratuit - Télécharger - Édition & Programmation
- Code ascii algorithme - Guide
- Ecrire un algorithme qui permet de resoudre ax²+bx+c=0 - Forum Programmation
- Remplir une matrice algorithme - Forum Pascal
- Ecrire un algorithme qui calcule le produit de deux nombres - Forum Programmation
4 réponses
sfel
Messages postés
1640
Date d'inscription
lundi 18 juin 2007
Statut
Membre
Dernière intervention
15 juillet 2009
430
13 janv. 2008 à 22:13
13 janv. 2008 à 22:13
bonsoir,
déjà ne fait pas
// probleme
mais % probleme
sinon, peut-être le N-L
utilise une variable temporaire juste avant, par exemple temp
déjà ne fait pas
// probleme
mais % probleme
sinon, peut-être le N-L
utilise une variable temporaire juste avant, par exemple temp
temp=N-L; for i=1:temp
zmeryem
Messages postés
3
Date d'inscription
dimanche 13 janvier 2008
Statut
Membre
Dernière intervention
14 janvier 2008
1
13 janv. 2008 à 22:45
13 janv. 2008 à 22:45
ca na pa marché :(
essaie ça:
N= 500;
h=[1 -2.5 1];
RSB=10;
L=length(h)-1;
delta =10^(-9);
x=randn(1,N);
x=x-mean(x);
b=randn(1,N);
b=b-mean(b);
b=(b/std(b))*sqrt(var(x)*10^(RSB/10));
y=filter(h,1,x)+b;
%iinitialisation des parametres de l'algorithme
Q=eye(L+1)/delta; h1=zeros(L+1,1);
for i=1:N-L %probleme
xn= x(i+L:-1:i);
Qx = Q * xn';
Kn= Qx /(1 + xn*Qx);
h1 = h'+ Kn* y(i+L)-xn*h1;
Q = Q-Kn*(xn*Q);
e(:,i)=h1'-h;
end
N= 500;
h=[1 -2.5 1];
RSB=10;
L=length(h)-1;
delta =10^(-9);
x=randn(1,N);
x=x-mean(x);
b=randn(1,N);
b=b-mean(b);
b=(b/std(b))*sqrt(var(x)*10^(RSB/10));
y=filter(h,1,x)+b;
%iinitialisation des parametres de l'algorithme
Q=eye(L+1)/delta; h1=zeros(L+1,1);
for i=1:N-L %probleme
xn= x(i+L:-1:i);
Qx = Q * xn';
Kn= Qx /(1 + xn*Qx);
h1 = h'+ Kn* y(i+L)-xn*h1;
Q = Q-Kn*(xn*Q);
e(:,i)=h1'-h;
end
16 avril 2012 à 10:16