kuaipaimv.m

宽带MVDR波束形成源代码

clear all
close all
clc
N=5
fs=200000;
%设采样频率
fai=-30
%入射角
t=(1:20000)/fs;
n=t*fs;
% 采样点数n=20000

%信号频率为3KHz
c=1500;
d=0.2;
%阵元间距为0.5米
L=20;    %快拍数
T=1;    %数据块长度
td=0.001;
nn=td*fs%滑动步长
Nbeam=180;
Pmv=zeros(1,Nbeam+1);
for f0=1000:100:3000;
X=2*pi*f0*t;
% X=2*pi*t*n/fs
Y=exp(j*X)
A=(0:(N-1));
t0=d*(sin(fai*pi/180))/c;
a=(-j)*2*pi*f0*t0
A=a*A
for ii=1:1:N
   A(:,ii)=exp(A(:,ii))
end
xx=(A')*Y;
[N,M]=size(xx);
x=xx+normrnd(0,10,N,M);
for l=1:1:L
    for i=1:1:N
        X1(i,:)=x(i,(l-1)*nn+1:(T+(l-1)*nn))
    end
%     Y=X_f+Y;
    YY(N*(l-1)+1:N*l,:)=X1;
end
 
        
       Rs=zeros(N,N);
        for l=1:1:l;
            R1=YY(N*(l-1)+1:N*l,:)*YY(N*(l-1)+1:N*l,:)';
            Rs=R1+Rs; 
          
        end
        
          RI=inv(Rs);
pos=d*(0:N-1)';

w=2*pi*f0/c;
theta=linspace(-90,90,Nbeam/2);
P=zeros(1,Nbeam);
mm=1
for k=-90:1:90
V=exp(j*w*pos*sin(k*pi/180));
P(mm)=1/abs(V'*RI *V);
mm=mm+1
end
Pmv=Pmv+P
end
RR=10*log10(Pmv)
m=-90:1:90;
plot(m,RR,'.-');
    grid,xlabel('方位 (度)')
    ylabel('功率(dB)');