azi_estimate_max.m

来自「利用频率细化方位估计,水声信号处理中的一项技术」· M 代码 · 共 25 行

%估计目标方位函数azi_estimate_max 
%符合加权
%有线谱情况：估计最大线谱的方位

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);      %能量最大的方位信息