clear all;close all;clc;
%---------------------------------------------
jj=0;
N=10900;
% for ii=50:50:500
% jj=jj+1;
% k=ii;
k=300;
u=0.0022;
d=300;
%---------------------------------------------
t1=0:10900-1;
load('loran-c+CWI.mat');
% I=In; %with noise
% I=I_c_n;%tadakhol+noise
%-------------------------
U=[0.025 0.009 u];
%-------------------------
y=zeros(1,N);
xd=zeros(1,N);
xd((d+1):N)=I(1:(N-d));
w=zeros(1,k);
e=zeros(1,N);
e(1:d)=I(1:d);
q=0;
b=1;
p=0;
if b==1
for oo=1:3
p=0;
w=zeros(1,k);
y=zeros(1,N);
xd=zeros(1,N);
xd((d+1):N)=I(1:(N-d));
e=zeros(1,N);
e(1:d)=I(1:d);
for o=1:5
p=0;
for i=(k):N
XN=xd((i-k+1):(i));
y(i)=w*XN';
r(i)=XN*XN';
e(i)=I(i)-y(i);
fai=0.0001;
w=w+U(oo)*XN*e(i)/(r(i)+fai);
p=p+1;
for j=1:p
tavan(oo,j)=(abs(e(j)-i10(j)))^2;
end
mse(oo,p)=sum(tavan(oo,:))/p;
if o==1
mse1(oo,p)=sum(tavan(oo,:))/p;
end
end
end
end
end
if b==2
%
% for o=1:10
for h=1:10
for i=(h*1000-600)*q+(k)*(1-q):1000*h
u=0.0009;
XN=xd(i-k+1:(i));
y(i)=w*XN';
r(i)=XN*XN';
e(i)=I(i)-y(i);
fai=0.0001;
w=w+u*XN*e(i)/(r(i)+fai);
p=p+1;
for j=1:p
tavan(j)=(abs(y(j)-xd(j)))^2;
end
mse1(p)=mean(tavan);
end
for i=1000*h+1:1000*h+400 %N-k+1
u=0.005;
XN=xd((i):(i+k-1));
y(i)=w*XN';
e(i)=I(i)-y(i);
end
q=1;
end
% end
end
% for i=1000*10+400+1:N-k+1
% XN=xd((i):(i+k-1));
% y(i)=w*XN';
% e(i)=I(i)-y(i);
% end
% for i=(k):N-k+1
% XN=xd((i-k+1):(i));
% y(i)=w*XN';
% r(i)=XN*XN';
% e(i)=I(i)-y(i);
% fai=.0001;
% w=w+u*XN*e(i)/(r(i)+fai);
% end
% figure
% e_f=fft(e);
% e_abs=abs(e_f);
% plot(e_abs)
% for j=1:9600
% tavan(j)=(abs(e(j)-i10(j)))^2;
% end
% mse=mean(tavan);
% % mse(jj)=mean(tavan);
% % end
% mse'
plot(t1,e)
xlim([0 9900])
xlabel('t(us)')
ylabel('amplitude')
hold on
plot(t1,i10,'r')
hold off
save('error','e','i10')
%
%
figure
plot(1:p,mse1(1,:),'r')
% semilogy(mse1(1,:),'p','-',0.2)
hold on
plot(1:p,mse1(2,:),'g')
hold on
plot(1:p,mse1(3,:))
xlim([0 9900])
xlabel('iteration')
ylabel('MSE')
hleg1 = legend('u=0.0.25','u=0.009','u=0.003');