MVDRJamming.rar_different angle_jammer _jamming_mvdr

%% problem 1 clear all; clc; %% parameters of simulation N=10; M=1; SNR=20; INR=70; doa_tg=0*pi/180; % target angle doa_jam=-30*pi/180; % jammer angle K=300; % number samples (snapshots) d_lambda=0.5; % distance between elements in wavelength %% Construct noise noise=(sqrt(1/2))*(randn(N,K)+1j*randn(N ,K)); %% Construct jammer jam=((10^(INR/20))/(sqrt(2)))*(randn(M,K)+1j*randn(M,K)); %% streeing vector steering_1=exp(-1j*2*pi*d_lambda*(0:N-1)'*sin(doa_tg)); steering_2=exp(-1j*2*pi*d_lambda*(0:N-1)'*sin(doa_jam)); %% Construct SOI soi=zeros(N,K); soi(:,50)=(10^(SNR/20))*steering_1; %% received signal angles=[-90:0.1:90]; a1=exp(-1j*2*pi*d_lambda*(0:N-1)'*sin([angles(:).']*pi/180)); z=zeros(N,K); U=zeros(N,N); z=steering_2*jam; receiv=soi+z+noise; Q=(10^(INR/10))*(steering_2*steering_2')+eye(N); IQ=inv(Q); for k=1:length(angles) Wopt=steering_1'*IQ*steering_1; mvdr=(1/Wopt)*IQ*steering_1; pmvdr(k)=mvdr'*a1(:,k); R=((10^(SNR/20))/sqrt(M))*(a1(:,k)*a1(:,k)'); snr(k)=(mvdr'*R*mvdr)/(mvdr'*Q*mvdr); end %% plot pattern figure(1) plot(angles,20*log10(abs(pmvdr)),'r-','linewidth',3) grid on xlabel('Angle in degrees') title('Pattern MVDR for one jammer') legend('Normalized Pattern') figure(2) plot(angles,10*log10(abs(snr)),'b-','linewidth',3) grid on xlabel('Angle in degrees') title('SNR LOSS MVDR') legend('Snr Loss')