| 查看: 790 | 回复: 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 |
» 猜你喜欢
女博士高校择业三天之内签合同，求支招
已经有32人回复
-
山东省基金：作为项目负责人承担面上项目数量累计不得超过2项
已经有5人回复
-
又见奇葩学生，可咋办
已经有31人回复
-
又忍不住想吐槽文科的学术圈
已经有11人回复
-
化学口面上
已经有6人回复
-
转发，”朋友说招呼都定点打到他那里了＂
已经有13人回复
-
化学B口多少分能上会呀
已经有15人回复
-
TsCl保护羟基为什么不反应
已经有25人回复
-
苯磺酰氯与醇羟基反应
已经有12人回复
-
海南大学国家级人才田新龙教授团队氢燃料电池与器件方向招收科研助理 (重点支持方向)
已经有11人回复
太好了2楼2006-05-09 07:56:42
0.5
 |
3楼2006-05-09 17:21:55
wanglong
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5612.6
- 红花: 1
- 帖子: 640
- 在线: 236小时
- 虫号: 92438
- 注册: 2005-09-09
- 性别: GG
- 专业: 大气与海洋光学
4楼2006-05-11 12:51:08
wanglong
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5612.6
- 红花: 1
- 帖子: 640
- 在线: 236小时
- 虫号: 92438
- 注册: 2005-09-09
- 性别: GG
- 专业: 大气与海洋光学
5楼2006-05-11 12:52:37
6楼2006-05-16 23:02:49
48