%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%PVM算法            PVM.m%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

R=(Xt*Xt')/N;
theta1=[-90:0.1:90];%基陣接受信號方向
i=[0:M-1]';
atheta1=exp(-j*2*pi*f*i*d*sin(theta1/180*pi)/c);
for i=1:size(theta1,2)
% wopt=(R^-1*atheta1(:,i))/(atheta1(:,i)'*(R^-1)*atheta1(:,i));
% Pw1(i)=wopt'*R*wopt;
Pw1(i)=1/(atheta1(:,i)'*(R^-1)*atheta1(:,i));

end

plot(theta1,10*log10(real(Pw1)/max(real(Pw1))));