%%%%%%%%%%%高分辨脉组误差成像子函数
function [minvs]=h_profile(n,deltaf,prf,f1,scat_range,v1,snr,scat_rcs,vc,winid,vs,nscat)
f0=f1;
num_pulses=n;
freq_step=deltaf;
PRI=1./prf;
signal1=zeros(num_pulses,1);
signal2=zeros(num_pulses,1);
nfft=512;
phase1_con=zeros(num_pulses,1);
phase2_con=zeros(num_pulses,1);
signal1_com=zeros(num_pulses,1);
signal2_com=zeros(num_pulses,1);
window=zeros(num_pulses,1);
for jscat=1:nscat
for i=1:num_pulses
res_freq=((i-1)*freq_step);
noise=exp(-4*j*pi*(f0+res_freq)*(scat_range(1)))/10^(snr/10);
signal1(i)=signal1(i)+sqrt(scat_rcs(jscat))*exp(-2*j*pi*(f0+res_freq)*...
(2.*scat_range(jscat)/vc-2*(v1/vc)*((i-1)*PRI)))+noise*(randn(1)+j*randn(1));
signal2(i)=signal2(i)+sqrt(scat_rcs(jscat))*exp(-2*j*pi*(f0+(num_pulses-1)*freq_step-res_freq)*...
(2.*scat_range(jscat)/vc-2*(v1/vc)*((i+127)*PRI)))+noise*(randn(1)+j*randn(1));
phase1_con(i)=phase1_con(i)+exp(-j*2*pi*(f0+res_freq)*2*(vs/vc)*((i-1)*PRI));
phase2_con(i)=phase2_con(i)+exp(-j*2*pi*(f0+(num_pulses-1)*freq_step-res_freq)*2*(vs/vc)*((i+127)*PRI));
end
end
if(winid>=0)
window(1:num_pulses)=hamming(num_pulses);
else
window(1:num_pulses)=1;
end
signal1_com=signal1.*phase1_con;
signal1_com=signal1_com.*window;
signal1_ifft=(ifft(signal1_com,nfft));
signal1_db=db(signal1_ifft);
signal2_com=signal2.*phase2_con;
signal2_com=signal2_com.*window;
signal2_ifft=(ifft(signal2_com,nfft));
signal2_db=db(signal2_ifft);
[signal1max,max1i]=max(signal1_db);
[signal2max,max2i]=max(signal2_db);
mn=max1i+max2i;
minvs=(nfft-mn)*vc/(2*nfft*PRI*(2*f0+(nfft-1)*freq_step));