mvdrb.m

波束形成中MVDR算法 matlab程序实现

clc;
clear all;
N=48;        %% N
k=1000;          %% sample number
thetas=10;%信号入射角度
thetai=[-30 30];%干扰入射角度
n=[0:N-1]';%n
sita=45*[-1:0.001:1];                                            %% scan vector
vs=exp(-j*pi*n*sin(thetas/180*pi));%信号阵列流型
vi=exp(-j*pi*n*sin(thetai/180*pi));%干扰流型
v=exp(-j*pi*n*sin(sita/180*pi));  %% 扫描流型 
f=15000;
t=[0:1:k-1]/1000;
snr=10;%信噪比
inr=10;%干噪比

xs=sqrt(10^(snr/10))*vs*exp(j*2*pi*f*t);%构造有用信号
xi=sqrt(10^(inr/10)/2)*vi*[randn(length(thetai),k)+j*randn(length(thetai),k)];%构造干扰信号(2个)
noise=[randn(N,k)+j*randn(N,k)]/sqrt(2);
% w=1/N*vs';
% figure(1),plot(sita,20*log10(abs(w*v)));
% axis([-90 90 -50 0]);
% hold on

%% mvdr 
NI=xi+noise;
Rni=NI*NI'/N;
wop1=Rni\vs/(vs'/Rni*vs);
B=abs(wop1'*v);
plot(sita,20*log10(B/max(B)),'r');
title('mvdr自适应的波束图(db)');xlabel('角度');ylabel('波束图');
grid on
axis([-45 45 -50 0]);
hold off