huibo.rar_LFM雷达_lfm_lfm 回波_radar_雷达回波信号

%--------------线性调频信号载频f0=30MHz,脉冲重复周期30us,带宽B=4MHz,发射脉宽6us------% function [signal,M,N]=huibo(B,D) %( D=B*Tw=24) R=20000; %目标范围40公里 v=3000; %目标相对于径向运动的速度300m/s q=0; %初始相位 c=3e8; %光速 Vgain=2; %天线方向图调制，仿真近似处理为对脉冲的幅度的影响为常数% f0=30e6; %设置载波频率 L=c/f0; %雷达信号的波长10; fd=250; %多普勒频率250; Btar=4*pi*R/L; fs=40e6; Ts=1/fs; %仿真采样频率及周期 tao=100e-6 %发射脉宽 fr=1000;Tr=1/fr; %脉冲重复周期 B=5e6; %调频带宽4M k=B/tao; %调频斜率 M=floor(fs*Tr); %每个回波对应采样的数据个数点 m1=floor(2*fs*Tr/10); m2=m1+B/(k*Ts); %动目标开始点与结束点 m3=floor(5*fs*Tr/10); m4=m3+B/(k*Ts); %静目标开始点 f1=f0-B/2; N=10; % 回波数10个 s1=zeros(N,M); for n=1:N for m=1:M s1(n,m)=(u(m-m1)-u(m-m2)).*Vgain.*cos(2.*pi.*( f1.*(m-m1).*Ts + k/2.*(m-m1).*(m-m1).*Ts.*Ts) + 2*pi*fd.*n.*Tr -Btar +q); %----动目标存储矩阵 s2(n,m)=(u(m-m3)-u(m-m4)).*Vgain.*cos(2.*pi.*( f1.*(m-m3).*Ts + k/2.*(m-m3).*(m-m3).*Ts.*Ts) -Btar +q); %------静目标存储矩阵 end end s=s1+s2; [h,l]=size(s); ss=zeros(1,M*N); col=1; for i=1:h for j=1:l ss(1,col)=s(i,j); col=col+1; end end % % -----产生一组服从Gauss分布均值为0，方差为1的随机数,作为系统噪声------% % nois=0.1.*randn(1,length(ss)); % % % -----------------------地物杂波zb----------% % % cd=floor(length(ss)/2); %-------考虑到y1的长度与swerling的长度不一致，处理方法将swerling的杂波长度多取一位,叠加时再去除----% % cd=rand(1,length(ss)); % sigma=2; % z2=0.08*sqrt(2*log(1./cd))*sigma;%.*s; %产生瑞利分布信号1 % % z2=zabo(cd); %--------round函数:四舍五入------------% % for p=1:N % z2(1+(p-1).*M:(M-50)+(p-1).*M)=zeros(1,M-50); % end % % --------------线性调频信号js+系统噪声nois+地物杂波zb--------------------% % z=zeros(1,length(ss)) % for i=1:length(ss) % z(i)=ss(i)+nois(i)+z2(i); % end z=ss; signal=z; subplot(2,1,1); plot(z); %---------------------线性调频信号显示------------% title('目标回波的信号(含躁声)') zf1=z(1:M); Nf=M; zft=fft(zf1,Nf); n=0:Nf-1; subplot(2,1,2); q=(n./Nf).*(fs); plot(q,zft); title('一个周期内回波的频谱(含躁声)') %subplot(3,1,3); %plot(n,abs(zft));