LFM_huibo.rar_lfm回波_扩展目标_线性调频 雷达_雷达_雷达发射

function [Srt]=LFM_radar(T,B,Rmin,Rmax,R,RCS,S0,S1,fc,v) if nargin==0 S0=[0,0,0]; S1=[6000,5100,7000;6000,5000,6990;6000,5050,7000]; T=10e-6;%发射脉宽 B=30e6;%调频带宽 Rmin=10000; Rmax=13000;%距离门 %目标位置向量 RCS=[1 1 1 1 1 1];%目标相对有效反射面 fc=10e6; v=10; end R=Direction(S0,S1); Rmin=min(R); Rmax=max(R); R=R-Rmin;%目标位置向量，相对 C=3e8;%光速 fd=2*v/C*fc; K=B/T;%调频斜率 Rrec=Rmax-Rmin;%接收门宽/m Trec=2*Rrec/C;%接收门宽/s Fs=5*B;Ts=1/Fs;%计算机采样频率和周期 N0=ceil(T/Ts);%发射脉宽对应的采样 N=ceil(Trec/Ts);%接收窗对应的采样点数 %产生高采样率的chirp信号 t1=linspace(0,T,N0); Chirp=exp(1i*pi*K*t1.^2+1i*2*pi*t1*fc); figure(1); plot(t1*10^6,real(Chirp)); xlabel('\fontsize{9}时间（us）');title('\fontsize{9}线性调频信号'); ylabel('幅度') Chirp=exp(1i*pi*K*t1.^2+1i*2*pi*t1*(fc+fd)); M=length(R); Srt(1:N+N0)=0; %初始化接收信号 for i=1:1:M temp(1:N+N0)=0; tau=2*R(i)/C; K0=ceil(tau/Ts); %第i个回波的相对位置 temp(K0+1:K0+N0)=RCS(i)*Chirp(1:N0); a=rand(); a=0.3; Srt=Srt+temp*exp(1j*2*pi*a);%随机相位 将所有回波叠加 end %Srt figure(2); t=linspace(0,Ts*(N+N0),N+N0); plot(t*10^6,real(Srt)); xlabel('\fontsize{9}时间（us）');title('\fontsize{9}线性调频信号回波'); ylabel('幅度')