%function main()
close  all

% 周期信號的產生 
t=0:99;
xs=10*sin(0.5*t);
figure;
subplot(2,1,1);
plot(t,xs);grid;
ylabel('幅值');
title('it{輸入周期性信號}');

% 噪聲信號的產生
randn('state',sum(100*clock));
xn=randn(1,100);
subplot(2,1,2);
plot(t,xn);grid;
ylabel('幅值');
xlabel('時間');
title('it{隨機噪聲信號}');

% 信號濾波
xn = xs+xn;
xn = xn.' ;   % 輸入信號序列
dn = xs.' ;   % 預期結果序列
M  = 20   ;   % 濾波器的階數

rho_max = max(eig(xn*xn.'));   % 輸入信號相關矩陣的最大特征值
mu = rand()*(1/rho_max)   ;    % 收斂因子 0 < mu < 1/rho

[yn,W,en] = LMS(xn,dn,M,mu);

% 繪制濾波器輸入信號
figure;
subplot(2,1,1);
plot(t,xn);grid;
ylabel('幅值');
xlabel('時間');
title('it{濾波器輸入信號}');

% 繪制自適應濾波器輸出信號
subplot(2,1,2);
plot(t,yn);grid;
ylabel('幅值');
xlabel('時間');
title('it{自適應濾波器輸出信號}');

% 繪制自適應濾波器輸出信號,預期輸出信號和兩者的誤差
figure 
plot(t,yn,'b',t,dn,'g',t,dn-yn,'r');grid;
legend('自適應濾波器輸出','預期輸出','誤差');
ylabel('幅值');
xlabel('時間');
title('it{自適應濾波器}');