huibo.m

用MATLAB实现雷达回波信号的仿真

%--------------线性调频信号载频f0=30MHz,脉冲重复周期30us,带宽B=4MHz,发射脉宽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;              %天线方向图调制，仿真近似处理为对脉冲的幅度的影响为常数% 
f0=30e6;              %设置载波频率   
L=c/f0;               %雷达信号的波长10;
fd=250;             %多普勒频率250;
Btar=4*pi*R/L;          
fs=40e6; Ts=1/fs;    %仿真采样频率及周期
tao=100e-6           %发射脉宽
fr=1000;Tr=1/fr;             %脉冲重复周期
B=5e6;                %调频带宽4M   
k=B/tao;           %调频斜率
M=floor(fs*Tr);       %每个回波对应采样的数据个数点
m1=floor(2*fs*Tr/10);  m2=m1+B/(k*Ts);   %动目标开始点与结束点
m3=floor(5*fs*Tr/10);  m4=m3+B/(k*Ts);   %静目标开始点
f1=f0-B/2;
N=10;                 % 回波数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
% % -----产生一组服从Gauss分布均值为0，方差为1的随机数,作为系统噪声------%
% nois=0.1.*randn(1,length(ss)); 
% 
% % -----------------------地物杂波zb----------%
% % cd=floor(length(ss)/2);  %-------考虑到y1的长度与swerling的长度不一致，处理方法将swerling的杂波长度多取一位,叠加时再去除----%
% cd=rand(1,length(ss));
% sigma=2;
% z2=0.08*sqrt(2*log(1./cd))*sigma;%.*s;      %产生瑞利分布信号1
% % z2=zabo(cd); %--------round函数:四舍五入------------%
% for p=1:N
%   z2(1+(p-1).*M:(M-50)+(p-1).*M)=zeros(1,M-50);
% end    
% % --------------线性调频信号js+系统噪声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);   %---------------------线性调频信号显示------------%
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('一个周期内回波的频谱(含躁声)')
%subplot(3,1,3);
%plot(n,abs(zft));