lms自适应时间估计,我的程序哪儿错了
各位大虾帮帮忙。我又两序列信号,一个是发射信号,一个是接收信号,两信号假设被传感器检测到的时间相差D。我想用lms自适应的方法求两信号的时间差D。我模拟出了两信号,并且根据自适应原理编程实现,可结果总是不对,希望懂的给指导一下。下面是我的程序。大家看看哪儿出错了第一个
%% the envelop of the simulation original signal and delay signal without
% noise;
% parameter:tau is the delay factor;
% used for solving the delay time;
% new program;
%% emission signal model
function twosignal(tau)
a=1.2*10.^4;f=40000;phi=0;
close all;
t=0:10.^(-6):3*10.^(-3);
alfa=a*t.^1.*exp(-t/(2*10^(-4)));
beta=cos(2*pi*f*t);
e_s=alfa.*beta;
y=hilbert(e_s); % 或者abs(sn)
e_ss=sqrt(real(y).^2+imag(y).^2);
e_ss=e_ss+0.05*randn(size(t));
subplot(211);plot(t,e_ss);title('Original Sigal');grid on
%% the receipt model
alfa=a*(t-tau).^1.*exp(-(t-tau)/(2*10^(-4))); % @
beta=0.*(t<tau)+cos(2*pi*f*(t-tau)+phi).*(t>=tau);% tau此处不能为0,分段函数作图方法
r_s=alfa.*beta;
y=hilbert(r_s); % 或者abs(sn)
r_ss=sqrt(real(y).^2+imag(y).^2);
r_ss=r_ss+0.05*randn(size(t));
hold on;plot(t,r_ss);title('Delay Sigal');
%% save
save time t;
save e_s e_ss; % the envelop of the original signal
save r_s r_ss; % the envelop of the return signal
%% command:twosignal(0.0003);
第二个程序
clear
close all;
load time; % load the time sequence
load e_s; % load the emission signal
load r_s; % load the receive signal
t;e_ss;r_ss;
e_signal=e_ss;
r_signal=r_ss;
subplot(311);plot(t,e_signal);grid on;
hold on;plot(t,r_signal);
title('Emission & Return signal');
%% begin of algorithm
% initialization
L=385; % 滤波器阶数
points=length(t); % the simple-points
d=r_ss(1:points)';
X=e_ss(1:points)';
y=zeros(1,points)'; % system output signal
w=zeros(L,1); % 每列代表每次迭代的L个权值,行代迭代次数,第几次
D=zeros(1,points);
a=0.001;b=10;
% LMS arithmetic
for n=L:points
e(n)=d(n)-w'*X(n:-1:n-L+1);
u=a*(1-exp(-b*(abs(e(n)).^2))); % 王小川的《基于》
w=w+u*e(n)*X(n:-1:n-L+1);
end
%for tt=1:points
% D(tt)=w(:,n)'*(sinc(t(tt)-(1:L)))'; % 参考王小川的《基于》interpolation
%end
subplot(312);plot(w);grid on;axis tight;% the approximate signal after the lms filter
=max(w);hold;plot(m,wm,'or');
%subplot(313);plot(t,D);grid on;axis tight;% the weighted value---filter coefficient
%=max(D); % find the maximum of the weighted value array
%tt=t(c); % look for the corresponding time spot of the maximum
%hold on;plot(tt,Dmax,'or') % draw the spot 所以没报错? 建议至信号版块逛逛! 这个问题就挺专业了.... 版主还不如撤掉,就知道卖书 本帖最后由 meiyongyuandeze 于 2011-4-7 18:20 编辑
买书?没理解LS的意思!可能是论坛里的人对你的问题不熟悉吧!
本帖最后由 ChaChing 于 2011-4-7 22:38 编辑
wweiguo 发表于 2011-4-7 16:15 http://www.chinavib.com/static/image/common/back.gif
版主还不如撤掉,就知道卖书
每个人的专业不同, 不见得什麼都懂吧!
况且也不是一定有义务得回应, 不是吗!?
版主们帮忙建议, 应心存感谢, 怎还如此态度...!
看看下帖, 相信LZ能够了解
建议提问的网友分清 编程问题 和 专业问题 http://www.chinavib.com/thread-36746-1-1.html
提问的智慧!!!!(发帖前请认真阅读) http://www.chinavib.com/thread-21991-1-1.html
回复 5 # meiyongyuandeze 的帖子
不好意思,发错了
页:
[1]