Znn3bq.jpeg
²é¿´: 1435  |  »Ø¸´: 5
µ±Ç°Ö÷ÌâÒѾ­´æµµ¡£

zzbei

ľ³æ (ÕýʽдÊÖ)

[½»Á÷] ×ÔÊÊÓ¦matlab³ÌÐò£¨LMSºÍRLS£©

%lmsËã·¨Ô´³ÌÐò

clear all
close all
%channel system order
sysorder = 5 ;
% Number of system points
N=2000;
inp = randn(N,1);
n = randn(N,1);
[b,a] = butter(2,0.25);
Gz = tf(b,a,-1);
%This function is submitted to make inverse Z-transform (Matlab central file exchange)
%The first sysorder weight value
%h=ldiv(b,a,sysorder)';
% if you use ldiv this will give h :filter weights to be
h= [0.0976;
0.2873;
0.3360;
0.2210;
0.0964;];
y = lsim(Gz,inp);
%add some noise
n = n * std(y)/(10*std(n));
d = y + n;
totallength=size(d,1);
%Take 60 points for training
N=60 ;  
%begin of algorithm
w = zeros ( sysorder , 1 ) ;
for n = sysorder : N
  u = inp(n:-1:n-sysorder+1) ;
y(n)= w' * u;
e(n) = d(n) - y(n) ;
% Start with big mu for speeding the convergence then slow down to reach the correct weights
if n < 20
mu=0.32;
else
mu=0.15;
end
  w = w + mu * u * e(n) ;
end
%check of results
for n = N+1 : totallength
  u = inp(n:-1:n-sysorder+1) ;
y(n) = w' * u ;
e(n) = d(n) - y(n) ;
end
hold on
plot(d)
plot(y,'r');
title('System output') ;
xlabel('Samples')
ylabel('True and estimated output')
figure
semilogy((abs(e))) ;
title('Error curve') ;
xlabel('Samples')
ylabel('Error value')
figure
plot(h, 'k+')
hold on
plot(w, 'r*')
legend('Actual weights','Estimated weights')
title('Comparison of the actual weights and the estimated weights') ;
axis([0 6 0.05 0.35])

% RLS Ëã·¨
randn('seed', 0) ;
rand('seed', 0) ;

NoOfData = 8000 ;  % Set no of data points used for training
Order = 32 ;    % Set the adaptive filter order

Lambda = 0.98 ;    % Set the forgetting factor
Delta = 0.001 ;    % R initialized to Delta*I

x = randn(NoOfData, 1) ;% Input assumed to be white
h = rand(Order, 1) ;  % System picked randomly
d = filter(h, 1, x) ;  % Generate output (desired signal)

% Initialize RLS

P = Delta * eye ( Order, Order ) ;
w = zeros ( Order, 1 ) ;

% RLS Adaptation

for n = Order : NoOfData ;
   
  u = x(n:-1:n-Order+1) ;
  pi_ = u' * P ;
  k = Lambda + pi_ * u ;
  K = pi_'/k;
  e(n) = d(n) - w' * u ;
  w = w + K * e(n) ;
  PPrime = K * pi_ ;
  P = ( P - PPrime ) / Lambda ;
  w_err(n) = norm(h - w) ;

end ;

% Plot results

figure ;
plot(20*log10(abs(e))) ;
title('Learning Curve') ;
xlabel('Iteration Number') ;
ylabel('Output Estimation Error in dB') ;

figure ;
semilogy(w_err) ;
title('Weight Estimation Error') ;
xlabel('Iteration Number') ;
ylabel('Weight Error in dB') ;

[ Last edited by »ÃÓ°ÎÞºÛ on 2006-10-26 at 07:47 ]
»Ø¸´´ËÂ¥

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

source

» ²ÂÄãϲ»¶

¾­Àú¾ÍÊÇÒ»ÖÖÏíÊÜ
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wwfjhx863

0.5

Ì«ºÃÁË
2Â¥2006-05-09 07:56:42
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

0.5

3Â¥2006-05-09 17:21:55
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wanglong

Ìú¸Ëľ³æ (ÕýʽдÊÖ)

·Ç³£¸Ðл£¡£¡£¡£¡ÀÏ´ó£¬ÓÐɶºÃ¶«¶«£¬¶àÉÏ´©Ò»µã
4Â¥2006-05-11 12:51:08
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

wanglong

Ìú¸Ëľ³æ (ÕýʽдÊÖ)

0.5

·Ç³£¸Ðл,¸Õ²ÅûÆÀ¼Û,ÖØÀ´
5Â¥2006-05-11 12:52:37
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

sanheiliu

ľ³æ (СÓÐÃûÆø)

ºÃÓÃ!лл!!!
6Â¥2006-05-16 23:02:49
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ zzbei µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0856ר˶Çóµ÷¼Á Ï£ÍûÊÇaÇøÔºÐ£ +22 ºÃºÃÐÝÏ¢ºÃ²»ºÃ 2026-04-09 25/1250 2026-04-12 17:09 by ajpv·çÀ×
[¿¼ÑÐ] 085410 273Çóµ÷¼Á +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by ÄæË®³Ë·ç
[¿¼ÑÐ] 347Çóµ÷¼Á +4 mhyqyy 2026-04-06 4/200 2026-04-12 02:27 by Çï¶¹²ËÑ¿
[¿¼ÑÐ] 332Çóµ÷¼Á +14 ½¶½¶123 2026-04-10 14/700 2026-04-12 00:27 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +18 MAXâê㯠2026-04-09 20/1000 2026-04-11 23:31 by zhen¡«
[¿¼ÑÐ] 267Çóµ÷¼Á +8 ÔÙæҲҪ³Ô·¹°¡ 2026-04-09 8/400 2026-04-11 21:42 by cfdbai
[¿¼ÑÐ] µ÷¼ÁÇóÖú +6 ¹ûÈ»ÓÐÎÒ 2026-04-11 7/350 2026-04-11 16:22 by Ã÷Ô´ËʱÓÐ
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤µ÷¼Á 10+11 ÏÂÒ»Õ¾Éϰ¶@ 2026-04-10 36/1800 2026-04-11 10:26 by 89436494
[¿¼ÑÐ] 080100Á¦Ñ§316Çóµ÷¼Á +8 L_Hairui 2026-04-07 8/400 2026-04-11 10:00 by zhq0425
[¿¼ÑÐ] 346£¬¹¤¿Æ0854Çóµ÷¼Á£¬×¨Ë¶ +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by Öí»á·É
[¿¼ÑÐ] ²ÄÁÏÀà284µ÷¼Á +40 Ïë»»ÊÖ»ú²»Ïë½âÊ 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[¿¼ÑÐ] Çൺ¿Æ¼¼´óѧ²ÄÁÏѧԺ£¬»·¾³Ñ§Ôºµ÷¼Á²¹Â¼4ÔÂ10ÈÕÒÔǰ¶¼¿ÉÒÔ +3 1Çà¿Æ´ó¡£ 2026-04-09 5/250 2026-04-10 09:58 by ôæôæÒ»ÊéÉú
[¿¼ÑÐ] Ò»Ö¾Ô¸Öпƴó070300»¯Ñ§£¬314·ÖÇóµ÷¼Á +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[¿¼ÑÐ] 0703»¯Ñ§Çóµ÷¼Á +21 ²»ÖªÃûµÄСئ 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[¿¼ÑÐ] ±¾¿ÆÖ£ÖÝ´óѧ£¬Ò»Ö¾Ô¸»ª¶«Ê¦·¶´óѧ282Çóµ÷¼Á +23 Ðܸçxtk 2026-04-07 26/1300 2026-04-09 17:17 by 18446523
[¿¼ÑÐ] һ־Ը³¶«´óѧ071000ÉúÎïѧѧ˶³õÊÔ·ÖÊý276Çóµ÷¼Á +3 Ľ¾øcc 2026-04-09 3/150 2026-04-09 09:57 by liuhuiying09
[¿¼ÑÐ] ÉúÎïѧ328·ÖÇóµ÷¼Á +9 ÉÁµçkkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[¿¼ÑÐ] 293·ÖÇóµ÷¼Á£¬ÍâÓïΪ¶íÓï +7 ¼ÓÒ»Ò»¾Å 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[¿¼ÑÐ] ²ÄÁϹ¤³ÌרҵÈÕÓïÉúÇóµ÷¼Á +9 111623 2026-04-07 9/450 2026-04-07 23:31 by Ò»Ö»ºÃ¹û×Ó?
[¿¼ÑÐ] 344Çóµ÷¼Á +11 κ×Óper 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û