%--------------線性調(diào)頻信號載頻f0=30MHz,脈沖重復(fù)周期30us,帶寬B=4MHz,發(fā)射脈寬6us------%
function [signal,M,N]=huibo(B,D)     %( D=B*Tw=24)
R=20000;              %目標范圍40公里
v=3000;                %目標相對于徑向運動的速度300m/s
q=0;                  %初始相位
c=3e8;                %光速
Vgain=2;              %天線方向圖調(diào)制，仿真近似處理為對脈沖的幅度的影響為常數(shù)% 
f0=30e6;              %設(shè)置載波頻率   
L=c/f0;               %雷達信號的波長10;
fd=250;             %多普勒頻率250;
Btar=4*pi*R/L;          
fs=40e6; Ts=1/fs;    %仿真采樣頻率及周期
tao=100e-6           %發(fā)射脈寬
fr=1000;Tr=1/fr;             %脈沖重復(fù)周期
B=5e6;                %調(diào)頻帶寬4M   
k=B/tao;           %調(diào)頻斜率
M=floor(fs*Tr);       %每個回波對應(yīng)采樣的數(shù)據(jù)個數(shù)點
m1=floor(2*fs*Tr/10);  m2=m1+B/(k*Ts);   %動目標開始點與結(jié)束點
m3=floor(5*fs*Tr/10);  m4=m3+B/(k*Ts);   %靜目標開始點
f1=f0-B/2;
N=10;                 % 回波數(shù)10個
s1=zeros(N,M);
for n=1:N
  for m=1:M 
     s1(n,m)=(u(m-m1)-u(m-m2)).*Vgain.*cos(2.*pi.*( f1.*(m-m1).*Ts + k/2.*(m-m1).*(m-m1).*Ts.*Ts) + 2*pi*fd.*n.*Tr -Btar +q);  %----動目標存儲矩陣
     s2(n,m)=(u(m-m3)-u(m-m4)).*Vgain.*cos(2.*pi.*( f1.*(m-m3).*Ts + k/2.*(m-m3).*(m-m3).*Ts.*Ts) -Btar +q);    %------靜目標存儲矩陣
  end
end     
s=s1+s2;
[h,l]=size(s);
ss=zeros(1,M*N);
col=1;
for i=1:h
    for j=1:l
        ss(1,col)=s(i,j);
        col=col+1;
    end 
end
% % -----產(chǎn)生一組服從Gauss分布均值為0，方差為1的隨機數(shù),作為系統(tǒng)噪聲------%
% nois=0.1.*randn(1,length(ss)); 
% 
% % -----------------------地物雜波zb----------%
% % cd=floor(length(ss)/2);  %-------考慮到y(tǒng)1的長度與swerling的長度不一致，處理方法將swerling的雜波長度多取一位,疊加時再去除----%
% cd=rand(1,length(ss));
% sigma=2;
% z2=0.08*sqrt(2*log(1./cd))*sigma;%.*s;      %產(chǎn)生瑞利分布信號1
% % z2=zabo(cd); %--------round函數(shù):四舍五入------------%
% for p=1:N
%   z2(1+(p-1).*M:(M-50)+(p-1).*M)=zeros(1,M-50);
% end    
% % --------------線性調(diào)頻信號js+系統(tǒng)噪聲nois+地物雜波zb--------------------%
% z=zeros(1,length(ss))
% for i=1:length(ss)
%     z(i)=ss(i)+nois(i)+z2(i);
% end  
z=ss;
signal=z;
subplot(2,1,1);
plot(z);   %---------------------線性調(diào)頻信號顯示------------%
title('目標回波的信號(含躁聲)')
zf1=z(1:M);
Nf=M;
zft=fft(zf1,Nf);
n=0:Nf-1;
subplot(2,1,2);
q=(n./Nf).*(fs);
plot(q,zft);
title('一個周期內(nèi)回波的頻譜(含躁聲)')
%subplot(3,1,3);
%plot(n,abs(zft));