ls.m

自己编的基于最抖下降法

%
%主函数
%XS是周期性信号
%XN是干扰信号
%
clear
t=0:0.1:999.9;
xs=(5*sin(0.05*t))';
figure(1);
subplot(211);
plot(t,xs);grid;
ylabel('fuzhi');
title('\it{输入周期性信号}');
xn=(sqrt(0.5)*randn(1,length(t)))';
subplot(212);
plot(t,xn);grid;
ylabel('幅度');
xlabel('时间');
title('\it{随机噪声信号}');
d=xs;
x=xs+xn;
N=10;
w=(zeros(1,N))';
M=length(x);
p=0.1*eye(N,N);
a=0.1;

y=(zeros(1,M))';
e=(zeros(1,M))';
sum1=zeros(N,N);
sum2=zeros(N,1);


for n=N:M
    x1=x(n:-1:n-N+1);
    juzhen=x1*x1';
    
    k=((1/a)*p*x1)/(1+(1/a)*x1'*p*x1);
    e(n)=d(n)-w'*x1;
    w=w+k*conj(e(n));
    p=(1/a)*p-(1/a)*k*x1'*p;
    y(n)=w'*x1;
end

figure(2);
plot(y);grid;
figure(3);
plot(e);grid;