function [x,b2,f]=hbxh(n);
%
%假設發射n個脈沖,返回n個目標回波信號
%

%
%設置參數
%
fs=100*10^6;          %采樣頻率100M
Ts=1/fs;              %采樣周期
f1=28000000;          %
f2=32000000;          %發射信號頻率為30M，線性調頻帶寬為4M，范圍為28M~32M
fd=1000;              %多普勒頻移設為1000


Tr=0.00003;           %脈沖重復周期30微秒
mk=0.000005;          %脈寬為5微秒
t=0:Ts:Tr-Ts;
N=length(t);          %N為3000
f=0:fs/N:fs-fs/N;

n=5;                  %產生n個回波信號

a=zeros(1,N);         %在一個脈沖周期內進行采樣
b=zeros(n,N);         %線性調頻信號
c=zeros(n,N);
x=zeros(n,N);         %最終回波信號  
n1=zeros(n,N);        %靜目標雜波
n2=0.2*randn(n,3000);        %高斯白噪聲



a0=1;
a(a0:a0+mk*fs)=1;     % mk*fs＝500，前500次采樣值為1，其余為零


k=(f2-f1)/mk;
delay_n=2000;
for i=1:n
    b(i,:)=a.*exp(j*2*pi*((f1+k*t/2).*t+fd*i*Tr));   %線性調頻信號+多普勒頻移
    b1(i,:)=a.*exp(j*2*pi*((f1+k*t/2).*t));          %線性調頻信號
    n1(i,delay_n:delay_n+mk*fs)=b1(i,a0:a0+mk*fs);   %靜目標雜波                      
end
b2=a.*exp(j*2*pi*((f1+k*t/2).*t));

delay_s=1000;                                          %回波信號相對于發射信號的延遲
c(:,delay_s:delay_s+mk*fs)=b(:,a0:a0+mk*fs);             
x=c+n1+n2;        %最終回波信號      
% text=x;