|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
fz=10e3%脉冲重复频率
7 T) @5 |+ E2 t% h, q4 otz=1/fz%脉冲重复周期
& k4 T9 {/ [0 i3 R7 Z+ T' Cfs=1e5%采样频率+ X Y% v; v+ @
ts=1/fs%采样周期
1 C( O, l2 Y3 f* wf_doppler=2.5e3%多普勒频率. ?4 N+ e% K a& j" j* d; u: D
N=tz/ts;%一个周期宽度; l1 [7 y. s6 K5 P
N=fix(N)%取接近且小于n的整数8 K3 S# O7 a- L" P% J0 c$ }
%%%产生雷达回波%%%%
3 o; o; b% I2 K8 aecho_mobj_pulse=[zeros(1,100),1,1,zeros(1,N-101)];/ |* |- }/ [3 ]) m
s_pc=repmat(echo_mobj_pulse,1,16);
& w4 [4 _ j5 N8 s. A% n4 Hn=1:16*N;4 _" r+ z+ ?4 U# N( E0 _8 c
s_doppler=cos(n*2*pi*f_doppler/fs);
7 B) g) h0 _$ fs_pc=s_pc.*s_doppler;9 ?" c- X% `" L6 z! W
s_noise=0.1*rand(1,N*16);: H: F: F" x+ Z$ Y. x5 a. t
s_pc=s_pc+s_noise;
/ a- }1 D+ e z0 ?! Z. bfigure,* G* y% c/ v p% F. Z3 y1 G
plot(0:ts16*N-1)*ts,s_pc),xlabel('t(单位:s)'),title('回波信号');, C; C( l! l% M9 b8 ^$ {
s_pc=reshape(s_pc,fix(N),16);# e/ e; N9 ^9 ~1 y( C/ E: O$ b6 A
|
|