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