LFMCW雷达发射接收信号及测距matlab仿真.zip

freq_min=76.5e9; freq_max=77e9; t_cycle=200e-6; b_lfm=freq_max-freq_min; k_lfm=(freq_max-freq_min)/(t_cycle/2); sample_freq=2048/t_cycle; dt=1/(sample_freq); t_stop=400e-6; t_vec=0:dt:t_stop; len_t=length(t_vec); c=3e8; f_transt=(b_lfm/2)*sawtooth(2*pi*(1/t_cycle)*t_vec,0.5); f_trans=freq_min+f_transt+(b_lfm/2); %三角波 %plot(t_vec,f_trans) % 目标参数 r_obj=102; % 单位为m serie_delay=ceil((2*r_obj/(3e8))/(dt)); f_rece=zeros(1,len_t); f_rece(1:serie_delay)=freq_min; f_rece(serie_delay+1:len_t)=f_trans(1:len_t-serie_delay); % figure,plot(t_vec,f_trans),hold on % plot(t_vec,f_rece,'r'),hold off f_mix=abs(f_rece-f_trans); % f_mix(15:16) % figure,plot(t_vec,f_mix) sig_mix=zeros(1,len_t); for s=1:len_t sig_mix(s)=cos(2*pi*f_mix(s)*t_vec(s)); end % figure,plot(t_vec,sig_mix) temp_add=rand(1,len_t); temp_sig=temp_add-mean(temp_add); sig_noise=temp_sig/max(abs(temp_sig)); snr=0; % 单位为dB amp_noise=1/(10^(snr/20)); sig_real=sig_mix+amp_noise*sig_noise; temp_freq=fft(sig_real); len_f=fix(len_t/2); f_vec=(sample_freq/2)*linspace(0,1,len_f); % f_m=1/(2*dt); % 来自于matlab help的实例 % f_vec=linspace(0,f_m,len_f); r_vec=(c/2)*f_vec/k_lfm; % 传播时延与混频后的信号的最强频率，以及调频率的关系 sig_fft=2*abs(temp_freq(1:len_f))/len_t; %figure,plot(r_vec,sig_fft) % 以2048点fft为例，绘制每个fft输出的距离值。 len_fft=2048; len_f=fix(len_fft/2); f_vec=(sample_freq/2)*linspace(0,1,len_f); r_vec=((3e8)/2)*f_vec/k_lfm; % 传播时延与混频后的信号的最强频率，以及调频率的关系 velo_vec=(3e8)*f_vec/(4*(freq_min)); temp_num=fix(len_t/len_fft); range_vec=zeros(1,temp_num); % figure, for s=1:temp_num s % w_ham=hamming(len_fft);w_ham=w_ham.'; temp_freq=fft(sig_real(1:2048)); sig_fft=2*abs(temp_freq(1:len_f))/len_fft; plot(f_vec,sig_fft),axis tight; set(gca,'YTick',0:0.2:1); xlabel('频率'); ylabel('幅度'); temp_max=find(sig_fft==max(sig_fft)); range_vec(s)=r_vec(temp_max(1)); end time_vec1=1:temp_num; figure,scatter(time_vec1,range_vec,'.'); ylabel('距离/m');