qr104.m

自己编写的程序包括 主要是天线的自适应算法仿真

clc;
clear;
close all;
tic
M=4;%阵元数目

snap=1000;%快拍数目
thetai=0;
thetaj=0;
fs=1000;
ts=1/fs;
t=(0:snap-1)*ts;
u0=5;
u1=10;
C=3e8;
f0=10e6;
lamda=C/f0;
d=0.5*lamda;
a=[0:M-1]';%阵列矢量
si=exp(j*2*pi*(f0*t+0.5*u0*t.^2));%信号
sj=2^0.5*exp(j*2*pi*(f0*t+0.5*u1*t.^2));%干扰
a_thetai=exp(j*2*pi*d/lamda*a*sin(thetai/180*pi));
a_thetaj=exp(j*2*pi*d/lamda*a*sin(thetaj/180*pi));
A=[a_thetai a_thetaj];
S=[si;sj];
X0=A*S;
randn('state',0);
real_noise=randn(size(X0));
randn('state',3);
imag_noise=randn(size(X0));
noise=(real_noise+j*imag_noise)/2^0.5;
X=X0+noise;%信号加噪声


R_x=X*X'/length(t);
[V,D]=eig(R_x)
conj(V)'

Q=V*conj(V)'
    
RX_1=V*inv(D)*conj(V)'
w=RX_1*a_thetai/(conj(a_thetai)'*RX_1*a_thetai)
theta=-90:0.1:90;
F=conj(w)'*exp(j*2*pi*d/lamda*a*sin(theta/180*pi));
plot(theta,10*log(abs(F)/max(abs(F))))