Matlab雷达回波数据模拟.doc

clear,pactJ=sqrt(-1);closeall%Getrootfilenameforsavingresultsfile=input('Enterrootfilenamefordataandlistingfiles:','s');%formradarchirppulseT=10e-6;    %pulselength,secondsW=10e6;     %chirpbandwidth,Hzfs=12e6;    %chirpsamplingrate,Hz;oversamplebyalittlefprintf('\nPulselength=%gmicroseconds\n',T/1e-6)fprintf('Chirpbandwidth=%gMhz\n',W/1e6)fprintf('Samplingrate=%gMsamples/sec\n',fs/1e6)s=git_chirp(T,W,fs/W);% 120-by-1arrayplot((1e6/fs)*(0:length(s)-1),[real(s)imag(s)])title('RealandImaginaryPartsofChirpPulse')xlabel('time(usec)')ylabel('amplitude')gridNp=20;             %20pulsesjkl=0:(Np-1);      %pulseindexarray,慢时间采样的序列,注意第一个PRI标记为0是为了慢时间起始时刻从零开始PRF=;        %PRFinHzPRI=(1/PRF);       %PRIinsecT_0=PRI*jkl;       %relativestarttimesofpulses,insecg=ones(1,Np);      %gainsofpulsesT_out=[1240]*1e-6;%startandendtimesofrangewindowinsec,这个就是接收窗的时间宽度TrecT_ref=0;           %systemreferencetimeinusec,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;           %RFfrequencyinHz;10GHzisX-bandfprintf('\nWearesimulating%gpulsesatanRFof%gGHz',Np,fc/1e9)fprintf('\nandaPRFof%gkHz,givingaPRIof%gusec.',PRF/1e3,PRI/1e-6)fprintf('\nTherangewindowlimitsare%gto%gusec.\n',...   T_out(1)/1e-6,T_out(2)/1e-6)%ComputeunambiguousDopplerintervalinm/puteunambiguousrangeintervalinmetersvua=3e8*PRF/(2*fc);%第一盲速rmin=3e8*T_out(1)/2;rmax=3e8*T_out(2)/2;rua=3e8/2/PRF;fprintf('\nTheunambiguo