clear?all
close?all
clc
j=sqrt（-1）;
M=10;
m=0:M-1;
theta=[1040];?%信號源10度，干擾40度
theta=theta*pi/180;
us=0;
A0=exp（j*pi*m‘*sin（theta（1）））;%信號源的方向矢量?
Ai=exp（j*pi*m‘*sin（theta（2）））;%干擾的方向矢量
SNR=20;??%?信噪比（可變）
INR=1;?
P0=10.^（SNR/10）;?
Pi=10.^（INR/10）;
o=pi/180;
xxl=[-90:1:90];
xx=xxl*o;
a=exp（j*pi*m‘*sin（xx））;???????%?掃描矢量;
%%%%%%%%%%%最優波束形成%%%%%%%%%%%%%%%%%%%%%%%
L=5000;
Xs=A0*sqrt（P0*2）*wgn（1L1‘complex‘）;%期望
Xi=Ai*sqrt（Pi*2）*wgn（1L1‘complex‘）;%干擾
Xn=sqrt（2）*wgn（ML1‘complex‘）;%噪聲
Xx=Xs+Xi+Xn;
?%Xx=Xi+Xw;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure
Rs=P0*A0*A0‘;??????????????????????????????????????%期望信號協方差矩陣
Ri=Pi*Ai*Ai‘;??????????????????????????????????????%干擾信號協方差矩陣
Rn=eye（M）;?????????????????????????????????????????%噪聲信號協方差矩陣
Rx=Rs+Ri+Rn;????
Pout_min=pinv（A0‘*pinv（Rx）*A0）;
Wopt=Pout_min*pinv（Rx）*A0;
%Wop=pinv（Rx）*A0/（A0‘*pinv（Rx）*A0）;
s1=20*log10（abs（Wopt‘*a）/max（abs（Wopt‘*a）））;
plot（xxls1‘--r‘）
grid?on
%%%%%%%%%%%%%%%LMS算法%%%%%%%%%%%%%%%%%%
w=zeros（M1）;
eta=0.01;
miu=0.0001;
e=1;
p=1;?
d=Xs（1:）;
while?abs（e）>eta?
?e=d（p）-w‘*Xx（:p）;
?w=w+miu*conj（e）*Xx（:p）;
?e1（p）=e;
?p=p+1;
end
s2=20*log10（abs（w‘*a）/max（abs（w‘*a）））;
hold?on
plot（xxls2‘-b‘‘linewidth‘1.5）;
legend（‘最優情況‘‘LMS算法‘）;
xlabel（‘到達角度‘）;
ylabel（‘輸出幅度增益/dB‘）;

for?k?=?1:p-1
????en（k）=（abs（e1（k）））.^2;
end
figure
semilogy（en）;
xlim（[0?500]）
%?axis（[0?500?0.001?10]）
xlabel（‘樣本數‘）;
ylabel（‘誤差e^{2}（n）‘）;
%?for?k?=?1:p-1
%?????en（k）=sum（abs（e1（1:k））.^2）/k;
%?end
%?figure
%?semilogy（en）;
%?ylim（[0?1000]）
%?%?axis（[0?500?0.001?10]）
%?xlabel（‘樣本數‘）;
%?ylabel（‘均方誤差E[e^{2}（n）]‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1925??2014-05-08?22:30??SMI.m

?????文件???????1755??2014-05-08?22:30??LMS.m

?????文件???????1808??2014-05-08?22:30??RLS.m

-----------?---------??----------?-----??----

?????????????????5488????????????????????3