clc;
clear?all;
%********************************4?array************************************
M=4;%4?antenna
N=[0:M-1]‘;%line
ratio=0.5;
Fd=1;%bit?rate
Fs=4;
Lx=2000;
delay=3;
r=0.5;
SNR=[103020];
SNR=sqrt（10.^（SNR/10））*sqrt（1/2）;
theta=[0;
??????????0.7;
??????????-0.4];?%jam2_theta
rand（‘seed‘sum（1000*clock））;
u=rand（1Lx）;????????
v=rand（1Lx）;
%?源信號
signal=sqrt（-2*log（1-u））.*cos（2*pi*v）*SNR（1）;
%?signal=signal（1:Lx）;
u=rand（1Lx）;????????
v=rand（1Lx）;
%?干擾信號（高斯白噪聲，方差為1均值為0）
jam_1=sqrt（-2*log（1-u））.*cos（2*pi*v）*SNR（2）;
jam_2=sqrt（-2*log（1-v））.*sin（2*pi*u）*SNR（3）;
%?jam_1=sqrt（1/2）*SNR（1）*（randn+j*randn）*signal;
%?jam_2=sqrt（1/2）*SNR（2）*（randn+j*randn）*signal;
rec_sig=[signal;jam_1;jam_2];
a_sig=exp（-j*2*pi*ratio*N*sin（theta（1）））;
a_jam1=exp（-j*2*pi*ratio*N*sin（theta（2）））;
a_jam2=exp（-j*2*pi*ratio*N*sin（theta（3）））;
a=[a_siga_jam1a_jam2];
noise=crandn（4Lx）;
ss=a*rec_sig+noise;
%******************************lms?processing*****************
Tr_lx=1000;
mu=1e-5;
%?d=awgn（signal_t20‘measured‘）;
d=signal‘;
w=zeros（41）;
w（11）=1;
en=[];
for?i=1:Lx
????en（i）=d（i）-w‘*ss（:i）;
????w=w+mu*ss（:i）*conj（en（i））;
?end
g2=w‘*a_g;
g_theta2=abs（g2）;
g_theta2=g_theta2./max（g_theta2）;
plot（k20*log10（g_theta2）‘g‘）;
%******************************RLS?processing*****************
M=4;
lamda=0.995;
delta=0.005;
P=delta^（-1）*eye（4）;
C=zeros（M1）;
g=zeros（M1）;????
%?wopt5=zeros（M1）;
%?wopt5（11）=randn;
wopt5=[];
for?j=1:M
????wopt5（j1）=rand;
end
for?i=1:Lx
?????u=ss（:i）‘*P*ss（:i）;
?????g=P*ss（:i）./（lamda+u）;?
?????e_rls?=d（i）?-?ss（:i）‘*wopt5;
?????wopt5=wopt5+e_rls.*g;
?????P=lamda^（-1）*（P-g*ss（:i）‘*P）;?????????????
end
g3=wopt5‘*a_g;
g_theta3=abs（g3）;
g_theta3=g_theta3./max（g_theta3）;
plot（k20*log10（g_theta3）‘c‘）;
hold
grid?
legend（‘lms‘‘rls‘）
title（‘Beam?Nulling?arithmetic?of?three?Arithmetic‘）
ylabel（‘Normalized?Antenna?Gain（db）‘）
axis（[-22-600]）

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

?????文件???????1984??2008-09-21?19:41??lms_4array.m

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

?????????????????1984????????????????????1