RLS2.zip

clear all; %close all; %参数设置 N = 8;%天线数 n = 1:N; K = 50;%块长 k = 1:K; alpha = 0.9;%遗忘因子，取[0-1] d_lamba = .5; T = 1e-3;%1ms sigma = 0.01; f = 20*1e9; theta0 = 30; theta1 = -60; theta2 = 0; a0 = exp(1j*2*pi*d_lamba*(n-1)*sind(theta0)); n1 = exp(1j*2*pi*d_lamba*(n-1)*sind(theta1)); n2 = exp(1j*2*pi*d_lamba*(n-1)*sind(theta2)); S = cos(2*pi*f*(k-1)/(K-1))+1j*sin(2*pi*f*(k-1)/(K-1));%期望信号 %S = cos(2*pi*f*(k-1)/(K-1));%期望信号 I1 = randn(1,K); I2 = randn(1,K); % I1 = sin(pi*(k-1)/(K-1));%干扰信号1 % I2 = sin(pi*(k-1)/(K-1));%干扰信号2 xN = sqrt(sigma)*randn(N,K); Rn = xN*xN'/K; xS = a0.'*S; xI1 = n1.'*I1; xI2 = n2.'*I2; x = xS+xI1+xI2; Rx = x(:,1)*x(:,1)'+Rn; Rx_inv = inv(Rx); alpha = 1/alpha; w = zeros(N,1); for i = 2:K Rx_inv = alpha*Rx_inv-alpha^2*Rx_inv*x(:,i)*x(:,i)'*Rx_inv/(1+alpha*x(:,i)'*Rx_inv*x(:,i)); g = Rx_inv*x(:,i); w = w+g*(conj(S(i))-x(:,i)'*w); end w = w/w(1); theta = -90:1:90; AF = zeros(1,length(theta)); for t = 1:N AF = AF + w(t)'.*exp(1j*2*pi*d_lamba*(t-1)*sind(theta)); end plot(theta,abs(AF)/max(abs(AF)));hold on