%估計目標方位函數(shù)azi_estimate_max 
%符合加權(quán)
%有線譜情況：估計最大線譜的方位

function     [amplitude_theta,theta,result]=azi_estimate_max(sp,svx,svy,spvx,spvy,nfft,fs,ind_up,ind_down);

%限頻帶范圍

sp=sp(ind_down:ind_up);
spvx=spvx(ind_down:ind_up);
spvy=spvy(ind_down:ind_up);
svx=svx(ind_down:ind_up);
svy=svy(ind_down:ind_up);

%估計方位
svr=svx+svy;
energy=real(spvy).^2+real(spvx).^2;			%求出各個頻率能量

theta=atan2(real(spvy),real(spvx))*180/pi;				%求出各個頻率指定的方位
amplitude_theta=energy./(sp.*svr).*(sp+svr);            %各譜線幅度信息
ind_max=find(amplitude_theta==max(amplitude_theta));

result=theta(ind_max);      %能量最大的方位信息