clc
clear
abar=input('Masukkan tebakan frekuensi= ');
fi=input('Sudut fasa berapa pi?');
dt=input('dt= ');
mu1=input('mu untuk estimasi frekuensi = ');
mu2=input('mu untuk estimasi fasa = ');
ta=input('serta ta=');
a1(1)=abar;
fi1(1)=fi;
a=1;
fi0=0;
t=0:dt:ta;
n=length(t);
for i=1:(n-1)
   home
   k=i
   clc
   ab=a1(i);
   fib=fi1(i);
   x(i)=sin(a*t(i)+fi0*pi);
   x1(i)=sin(ab*t(i)+fib*pi);
   e(i)=abs(x(i))-abs(x1(i));
   a1(i+1)=a1(i)+mu1*e(i);
   fi1(i+1)=fi1(i)+mu2*e(i);
end;
x(n)=sin(a*t(n)+fi0*pi);
ab=a1(n);
fib=fi1(n);
x1(n)=sin(ab*t(n)+fib*pi);
e(n)=abs(x(n))-abs(x1(n));
subplot(511),plot(t,x),xlabel('t'),ylabel('x(t)'),grid;
subplot(512),plot(t,x1),xlabel('t'),ylabel('x^(t)'),grid;
subplot(513),plot(t,e),xlabel('t'),ylabel('e(t)'),grid;
subplot(514),plot(t,a1),xlabel('t'),ylabel('est. frek.'),grid;
subplot(515),plot(t,fi1),xlabel('t'),ylabel('est. fasa'),grid;