注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

RLS 时间延迟估计 Demo  

2012-07-06 23:36:46|  分类: 学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

clear;
N=1024;
n=200;
FS=1000;
t=0:1/FS:(N-1)/FS;
s=exp(-30*t).*sin(80*pi*t);
D=50;
S=[zeros(1,N) s];
figure(1);
plot(s);
hold on;
plot(S,'r');
figure(2);
x1=awgn(s,10,'measured');
subplot(211);
plot(x1);
x2=awgn(S,10,'measured');
subplot(212);
plot(x2);
Cxx=1000*eye(n);
g=zeros(N,n);
h=zeros(N,n);
e=zeros(1,N);
y=zeros(1,N);
tr=zeros(1,N);
for (j=1:N-n+1)
    wm=max(h);
    if h(j,:)==wm
        f(j,:)=wm;
    end
end
figure(3)
plot(f)
for i=n+1:N,
    g(i,:)=(Cxx*x1(i-n+1:i)'/(1+x1(i-n+1:i)*Cxx*x1(i-n+1:i)'))';
    e(1,i)=x2(i)-h(i-1,:)*x1(i-n+1:i)';
    h(i,:)=h(i-1,:)+e(1,i)* g(i,:);
    Cxx= Cxx-g(i,:)'*x1(i-n+1:i)*Cxx;
    y(1,i)= h(i,:)*x1(i-n+1:i)';
    tr(1,i)=trace(Cxx);
end
figure(4)
plot(0:N-n,x1(1,n:N)),grid
title('x1(k) input signal in V')
xlabel('Sanples')
figure(5)
plot(0:N-n,x2(1,n:N),'r')
grid
axis([0 800 -1.2 1.2])
title('x2(k) reference signal in V')
xlabel('Samples')
figure(6)
plot(0:N-n,e(1,n:N)),hold on
plot(0:N-n,y(1,n:N),'r'),hold off
grid
title('e(k) error and y(k) output in V')
xlabel('Samples')
gtext('e(k)'),gtext('y(k)')
figure(7)
plot(0:N-n,h(n:N,1)),hold on
plot(0:N-n,h(n:N,2),'r'),hold off
grid 
title('filter coefficients evolution')
xlabel('Samples')


 

  评论这张
 
阅读(499)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016