matlab radar simulink system matlab雷达系统仿真程序

matlab雷达系统仿真程序，包含回波信号模拟产生、脉压、积累；可很大帮助初学者认识雷达系统。

close all
clear all
clc

%%产生发射信号%%%
code=[1 1 1 1 1 -1 -1 1 1 -1 1 -1 1];
tao=10e-6;
fc=28e6;
f0=30e6;
fs=100e6;
ts=1/fs;
B=4e6;
t_tao=0:1/fs:tao-1/fs;
N=length（t_tao）;
k=B/fs*2*pi/max（t_tao）;
n=length（code）;
pha=0;

s=zeros（1n*N）;
for i=1:n
    if code（i）==1
        pha=pi;
    else pha=0;
    end
    s（1（i-1）*N+1:i*N）=cos（2*pi*fc*t_tao+k*cumsum（t_tao）+pha）;
end
t=0:1/fs:13*tao-1/fs;
% figure
% subplot（211）
% plot（ts）;
% title（‘13位巴克码+线性调频‘）
% s_fft_result=abs（fft（s（（1:N））））;
% subplot（212）
% plot（（0:fs/N:fs/2-fs/N）abs（s_fft_result（1:N/2）））;
% title（‘码内信号频谱‘）


%%%产生脉压系数%%%
%---zjjt---%
N=tao/ts;
n=0:N-1;
s1=s（1:N）;
local_oscillator_i=cos（n*f0/fs*2*pi）;
local_oscillator_q=sin（n*f0/fs*2*pi）;
fbb_i=local_oscillator_i.*s1;
fbb_q=local_oscillator_q.*s1;
window=chebwin（5140）;
[ba]=fir1（502*B/fswindow）;
fbb_i=[fbb_izeros（125）];
fbb_q=[fbb_qzeros（125）];
fbb_i=filter（bafbb_i）;
fbb_q=filter（bafbb_q）;
fbb_i=fbb_i（26:end）;
fbb_q=fbb_q（26:end）;
fbb=fbb_i+j*fbb_q;

%---产生线性调频脉压匹配系数---%
M=131072;
D=B*tao;
match_filter_1=ts*fliplr（conj（fbb））*sqrt（D）*2/tao;
match_filter_1_fft=fft（match_filter_1M）;
% figure
% subplot（211）
% plot（real（match_filter_1_fft））;
% title（‘脉压系数实部‘）
% subplot（212）
% plot（imag（match_filter_1_fft））;
% title（‘脉压系数虚部‘）

N=length（s）;
n=0:N-1;
local_oscillator_i=cos（n*f0/fs*2*pi）;
local_oscillator_q=sin（n*f0/fs*2*pi）;
fbb_i=local_oscillator_i.*s;
fbb_q=local_oscillator_q.*s;
window=chebwin（5140）;
[ba]=fir1（500.5window）;
fbb_i=[fbb_izeros（125）];
fbb_q=[fbb_qzeros（125）];
fbb_i=filter（bafbb_i）;
fbb_q=filter（bafbb_q）;
fbb_i=fbb_i（26:end）;
fbb_q=fbb_q（26:end）;
signal=fbb_i+j*fbb_q;
clear fbb_i;
clear fbb_q;

signal_fft=fft（signalM）;
pc_result_fft=signal_fft.*match_filter_1_fft;
pc_result=ifft（pc_result_fftM）;
% figure
% plot（（0:ts:length（signal）*ts-ts）pc_result（1:length（signal）））;
% title（‘回波脉冲压缩处理结果‘）

clear local_oscillator_i;
clear local_oscillator_q;
t=tao*length（code）;
match_filter_2=2*ts*fliplr（conj（pc_result））*2/t;
match_filter_2_fft=fft（match_filter_2M）;
% figure
% subplot（211）
% plot（real（match_filter_2_fft））;
% title（‘脉冲压缩系数实部‘）
% subplot（212）
% plot（imag（match_filter_2_fft））;
% title（‘脉冲压缩系数虚部‘）
clear fbb；
clear match_filter_1;
clear match_filter_2;
clear signal;
clear signal_fft;
clear pc_result;
clear pc_result_fft;

%%%产生雷达回波%%%
f_frame=1e3;
T_frame=1/f_frame;
N_echo_frame=18;
f_doppler=3.5e3;
t_frame=0:ts:T_frame-ts;
t_mobj=200e-6;
echo_mobj_pulse=[zeros（1t_mobj/ts）szeros（1（T_frame-t_mobj）/ts-length（s））];
echo_mobj=repmat（echo_mobj_pulse1N_echo_frame）;
t_doppler=0:ts:N_echo_frame*T_frame-ts;
s_doppler=cos（2*pi*f_doppler*t_doppler）;
s_echo_mobj=echo_mobj.*s_doppler;
t_fobj=450e-6;
echo_fobj_pulse=[zeros（1t_fobj/ts）szeros（1（T_frame-t_fobj）/ts-length（s））];
s_echo_fobj=repmat（echo_fobj_pulse1N_echo_frame）;

t_clutter=700e-6;
t_clutter_pulse=39e-6;
sigm

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       6738  2014-04-18 19:28  radarsystem.m

----------- ---------  ---------- -----  ----

                 6738                    1