trynlms.m

来自「本次的算法是相位周期展开聚类」· M 代码 · 共 15 行

x=0.1*randn(1,10000);
fidG1=fopen('G1.bin','r');  
G1=fread(fidG1,2500,'double');
fclose(fidG1);
G1=G1';
%b=fir1(31,0.5);
d=filter(G1,1,x);
w0=zeros(1,2500);
u=0.8;
for i=2500:10000
    x1=x(i:-1:i-2500+1);
    e=d(i)-x1*w0';
    w0=w0+u*e*x1/(x1*x1');
end
w0