hbxh.m

有关雷达系统的仿真

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;