clear all;close all;

m=6 ;% Array number陣元
p=4; %  signal number信號數
N=2000;% The number of sample采樣數或迭代次數
A=zeros(m,p); % array pattern陣列流型
theta=[30 0 -45 60]*pi/180;% the direction for each signal and the signal from the direction of 30 degree is desired. DOA 30為期望信號方向
j=sqrt(-1);
w=[0.1 0.2 0.3 0.4]; % frequency for each signal.各個信號的數字頻率
s=zeros(p,N); %signal matrix.信號矩陣
s_rec=zeros(m,N);%matrix for received signal ,assume channel gain is 1.接收到的信號矩陣,假設信道增益為 1
%%
s=to_get_s(w,N,p); % degital complex siginal 復指數字信號exp(j*wi*n)
s_rec=get_s_rec(s,m,p,theta);
S=s_rec; % mxn array date. m*N 的陣元數據

X=S;
d=s(1,1:N); %desired siginal期望信號
%%%%—————————————LMS
%[e,wop]=lms(m,N,X,d)
%%%%--------------------------RLS
[e,wop]=rls(m,N,X,d)

%%-------------------------畫圖
drawpp(m,e,wop);