Een vraag over Lattice Filteren op LMS-algoritme

[elvaman2004]

Ik weet dat de Matlab de diatropic Filter realiseren op LMS-algoritme,
maar Hoe gebruik Matlab om de Lattice Filter realiseren op LMS algoritme? Please help me resove!

Het programma van Matlab te realiseren de diatropic Filteren op LMS-algoritme:

Alles wissen
sluit alle
HOLD uit
% channel systeem om
sysorder = 5;
% Aantal punten systeem
N = 2000;
INP = randn (N, 1);
n = randn (N, 1);
[b, a] = boter (2,0.25);
Gz = tf (b, a, -1);
% Deze functie wordt ingediend om inverse Z-transformatie (Matlab centraal bestand uitwisseling)
% De eerste sysorder gewicht waarde
% h = ldiv (b, a, sysorder) ';
% Als u ldiv gebruik van deze zullen geven h: filter gewichten te
h = [0.0976;
0,2873;
0,3360;
0,2210;
0,0964;];
y = lsim (GZ, INP);
% voeg wat lawaai
n = n * std / (10 * STD );
d = y n;
totallength = size (d, 1);
Neem% 60 punten voor de opleiding
N = 60;
% begin van algoritme
w = nullen (sysorder, 1);
voor n = sysorder: N
U = INP (n: -1: n-sysorder 1);
y = w '* u;
e = d - y ;
% Begin met grote mu voor snelheidsovertredingen de convergentie dan vertragen de juiste gewichten te bereiken
Als n <20
mu = 0.32;
anders
mu = 0.15;
eindigen
w = w mu * u * e ;
eindigen
% controle van de resultaten
voor n = N 1: totallength
U = INP (n: -1: n-sysorder 1);
y = w '* u;
e = d - y ;
eindigen
wacht even
plot (d)
plot (y, 'r');
titel ( 'System output');
xlabel ( 'Monsters')
ylabel ( 'Ware en geschatte output')
figuur
semilogy ((abs (e)));
titel ( 'Fout curve');
xlabel ( 'Monsters')
ylabel ( 'Fout waarde')
figuur
plot (h, 'k ')
wacht even
plot (w, 'r *')
legende ( 'Actual gewichten' Geschatte gewichten ')
titel ( 'Vergelijking van de werkelijke gewichten en de geschatte gewichten');
as ([0 6 0,05 0,35])

 

[mro83]

Hoi
MATALB (mijn ver.: 7) heeft een functie:
LMS
past het om te werken ur?

 

[elvaman2004]

HI
wo hier geen gebruik lms fuction!
realiseren we de LMS fuction.Ben je er een hebt?