E
elvaman2004
Guest
Žinau Matlab realizuoti diatropic Filtruoti LMS algoritmas,
bet, kaip naudotis Matlab realizuoti grotelių Filtruoti LMS algoritmas? Prašome padėti man resove!
Nuo Matlab programą įgyvendinti diatropic Filtruoti LMS algoritmas:
Išvalyti visus
uždaryti visus
Atidedamas
% kanalų sistemą siekiant
sysorder = 5;
% Taškų sistema taškai
N = 2000;
INP = randn (N, 1);
n = randn (N, 1);
= sviesto (2,0.25);
Gz = TF (b,, -1);
% Ši funkcija yra pateikta padaryti atvirkštinę z-transformacija (Matlab Centrinis failų mainų)
% Pirmojo sysorder masės vertė
% h = ldiv (B, sysorder) ";
%, Jeigu jūs naudojate ldiv tai duos h: filtro masė turi būti
h = [0,0976;
0,2873;
0,3360;
0,2210;
0,0964;]
y = lsim (GZ, INP);
% pridėti triukšmo
n = n * STD
/ (10 * STD );
d = y n;
totallength = dydis (d, 1);
60% visų taškų mokymas
N = 60;
% pradėti algoritmus
w = nuliai (sysorder, 1);
, n = sysorder: N
u = INP (n: -1: N-sysorder 1);
y = w '* u;
e = d - y ;
% Pradėti didelis mu už greičio viršijimą konvergencijos tada sulėtinti pasiekti teisingą svorio
Jei n <20
mu = 0,32;
kitas
mu = 0,15;
pabaiga
w = w mu * u * e;
pabaiga
% patikrinti rezultatus
, n = N 1: totallength
u = INP (n: -1: N-sysorder 1);
y = w '* u;
e = d - y ;
pabaiga
laikykis
sklypo plotas (d)
plot (y, 'r');
title ( 'Sistemos išėjimo ");
xlabel (fragmentai)
ylabel ( "Tiesa ir apskaičiuotas produkcija)
pav
semilogy ((ABS (e)));
title ( 'Klaida kreivė');
xlabel (fragmentai)
ylabel ( 'Klaida vertė)
pav
sklypo plotas (H "k ")
laikykis
sklypo plotas (W, R *)
legenda (Tikrasis svorius "," Numatomas svorio)
Pavadinimas ( "palyginimas faktinė masė ir apskaičiuotas svoriai");
axis ([0 6 0.05 0.35])
bet, kaip naudotis Matlab realizuoti grotelių Filtruoti LMS algoritmas? Prašome padėti man resove!
Nuo Matlab programą įgyvendinti diatropic Filtruoti LMS algoritmas:
Išvalyti visus
uždaryti visus
Atidedamas
% kanalų sistemą siekiant
sysorder = 5;
% Taškų sistema taškai
N = 2000;
INP = randn (N, 1);
n = randn (N, 1);
= sviesto (2,0.25);
Gz = TF (b,, -1);
% Ši funkcija yra pateikta padaryti atvirkštinę z-transformacija (Matlab Centrinis failų mainų)
% Pirmojo sysorder masės vertė
% h = ldiv (B, sysorder) ";
%, Jeigu jūs naudojate ldiv tai duos h: filtro masė turi būti
h = [0,0976;
0,2873;
0,3360;
0,2210;
0,0964;]
y = lsim (GZ, INP);
% pridėti triukšmo
n = n * STD
/ (10 * STD );d = y n;
totallength = dydis (d, 1);
60% visų taškų mokymas
N = 60;
% pradėti algoritmus
w = nuliai (sysorder, 1);
, n = sysorder: N
u = INP (n: -1: N-sysorder 1);
y
= w '* u;e
= d - y ;% Pradėti didelis mu už greičio viršijimą konvergencijos tada sulėtinti pasiekti teisingą svorio
Jei n <20
mu = 0,32;
kitas
mu = 0,15;
pabaiga
w = w mu * u * e
;pabaiga
% patikrinti rezultatus
, n = N 1: totallength
u = INP (n: -1: N-sysorder 1);
y
= w '* u;e
= d - y ;pabaiga
laikykis
sklypo plotas (d)
plot (y, 'r');
title ( 'Sistemos išėjimo ");
xlabel (fragmentai)
ylabel ( "Tiesa ir apskaičiuotas produkcija)
pav
semilogy ((ABS (e)));
title ( 'Klaida kreivė');
xlabel (fragmentai)
ylabel ( 'Klaida vertė)
pav
sklypo plotas (H "k ")
laikykis
sklypo plotas (W, R *)
legenda (Tikrasis svorius "," Numatomas svorio)
Pavadinimas ( "palyginimas faktinė masė ir apskaičiuotas svoriai");
axis ([0 6 0.05 0.35])