雷达系统仿真设计报告
pxx2=abs(fft(xdata(401:800),Ns).^2)/Ns; pxx3=abs(fft(xdata(801:1200),Ns).^2)/Ns; pxx4=abs(fft(xdata(1201:1600),Ns).^2)/Ns; pxx5=abs(fft(xdata(1601:2000),Ns).^2)/Ns; pxx=10*log10((pxx1+pxx2+pxx3+pxx4+pxx5)/5); f=(0:length(pxx)-1)/length(pxx);
subplot(313),plot(f,pxx);title('功率谱估计');ylabel('幅度(取对数)');xlabel('f');grid on;
第 17 页 共 26 页
雷达系统仿真设计报告
设计报告三实现代码
close all; clear all; clc;
%%%%%%产生雷达发射信号%%%%%%%%%%%%%%%%%% f_frame=1e3; %雷达发射信号重频
T_frame=1/f_frame; %雷达发射脉冲重复周期 code=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1]; 位巴克码 tao=10e-6; %脉冲宽度10us fc=28e6; %调频信号起始频率 f0=30e6; %中频频率 fs=50e6; %采样频率50MHz 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(1,n*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
%%%%%%%产生脉冲压缩系数%%%%%%%%%%%%%%%%%%%
第 18 页 共 26 页
雷达系统仿真设计报告
%对单个脉宽内线性LMF信号进行正交解调 n=0:N-1; s1=s(1:N);
local_oscillator_i=cos(n*f0/fs*2*pi); %i路本振信号 local_oscillator_q=sin(n*f0/fs*2*pi); %q路本振信号 fbb_i=local_oscillator_i.*s1; %i路解调 fbb_q=local_oscillator_q.*s1; %q路解调 window=chebwin(51,40); [b,a]=fir1(50,2*B/fs,window); fbb_i=[fbb_i,zeros(1,25)]; fbb_q=[fbb_q,zeros(1,25)]; fbb_i=filter(b,a,fbb_i); fbb_q=filter(b,a,fbb_q);
fbb_i=fbb_i(26:end); %截取有效信息 fbb_q=fbb_q(26:end); fbb=fbb_i+j*fbb_q;
%%产生理想线性调频脉冲压缩匹配系数%%%%%%%%%%%%%% M=65536;% 在一个周期内采样了50000个点，所以利用fft做65536点fft（2^16) D=B*tao;
match_filter_1=ts*fliplr(conj(fbb))*sqrt(D)*2/tao;
match_filter_1_fft=fft(match_filter_1,M);%第一次脉冲压缩匹配处理系数 %%%%对混合调制信号进行正交解%%%%%%%%%%%%%%% N=length(s); n=0:N-1;
local_oscillator_i=cos(n*f0/fs*2*pi); %i路本振信号 local_oscillator_q=sin(n*f0/fs*2*pi); %q路本振信号 fbb_i=local_oscillator_i.*s; %i路解调 fbb_q=local_oscillator_q.*s; %q路解调 window=chebwin(51,40); [b,a]=fir1(50,0.5,window);
第 19 页 共 26 页
雷达系统仿真设计报告
fbb_i=[fbb_i,zeros(1,25)]; fbb_q=[fbb_q,zeros(1,25)]; fbb_i=filter(b,a,fbb_i); fbb_q=filter(b,a,fbb_q);
fbb_i=fbb_i(26:end); %截取有效信息 fbb_q=fbb_q(26:end); signal=fbb_i+j*fbb_q;
signal_fft=fft(signal,M);
pc_result_fft=signal_fft.*match_filter_1_fft; pc_result=ifft(pc_result_fft,M); t=tao*length(code);
match_filter_2=2*ts*fliplr(conj(pc_result))*2/t;
match_filter_2_fft=fft(match_filter_2,M);%第二次脉冲压缩匹配处理系数
clear local_oscillator_i; clear local_oscillator_q; clear fbb;clear fbb_q;clear fbb_i; clear match_filter_1; clear match_filter_2; clear signal; clear signal_fft; clear pc_result; clear pc_result_fft;
%%%产生雷达回波信号：动目标反射回波+静目标反射回波+瑞利杂波+噪声 N_echo_frame=18; %本程序设计为对18个脉冲进行采样
%%动目标反射回波信号%%%%%%%%%%%%%%%%%%%%
第 20 页 共 26 页