Matlab雷达回波数据模拟..docVIP

clear, hold offformat compactJ = sqrt(-1);close all% Get root file name for saving resultsfile=input(Enter root file name for data and listing files: ,s); % form radar chirp pulseT = 10e-6;???? % pulse length, secondsW = 10e6;????? % chirp bandwidth, Hzfs = 12e6;???? % chirp sampling rate, Hz; oversample by a littlefprintf(\nPulse length = %g microseconds\n,T/1e-6)fprintf(Chirp bandwidth = %g Mhz\n,W/1e6)fprintf(Sampling rate = %g Msamples/sec\n,fs/1e6)s = git_chirp(T,W,fs/W); %? 120-by-1 arrayplot((1e6/fs)*(0:length(s)-1),[real(s) imag(s)])title(Real and Imaginary Parts of Chirp Pulse)xlabel(time (usec))ylabel(amplitude)gridNp = 20;????????????? % 20 pulsesjkl = 0:(Np-1);?????? % pulse index array, 慢时间采样的序列，注意第一个PRI标记为0是为了慢时间起始时刻从零开始PRF = 10.0e3;???????? % PRF in HzPRI = (1/PRF);??????? % PRI in secT_0 = PRI*jkl;??????? % relative start times of pulses, in secg = ones(1,Np);?????? % gains of pulsesT_out = [12 40]*1e-6; % start and end times of range window in sec， 这个就是接收窗的时间宽度TrecT_ref = 0;??????????? % system reference time in usec， T_ref = 0指T_0=0时，r_at_T_0 = ri ；当T_0 ~= 0时，r_at_T_0 = ri - vi*T_0(j)fc = 10e9;??????????? % RF frequency in Hz; 10 GHz is X-bandfprintf(\nWe are simulating %g pulses at an RF of %g GHz,Np,fc/1e9)fprintf(\nand a PRF of %g kHz, giving a PRI of %g usec.,PRF/1e3,PRI/1e-6)fprintf(\nThe range window limits are %g to %g usec.\n, ...??? T_out(1)/1e-6,T_out(2)/1e-6)% Compute unambiguous Doppler interval in m/sec% Compute unambiguous range interval in metersvua = 3e8*PRF/(2*fc); %第一盲速rmin = 3e8*T_out(1)/2;rmax = 3e8*T_out(2)/2;rua = 3e8/2/PRF;fprintf(\nThe unambiguous velocity interval is %g m/s.,vua)fprintf(\nThe range window starts at %g km.,rmin/1e3)fprintf(\nThe range window ends at %g km.,rmax/1e3)fprintf(\nThe unambiguous range interval is %g km.\n\n,rua/1e3)% Define number of targets, then range, SNR, and% radial velocity of each.? The SNR will be the actual SNR of the target in% the final data; it will no