1、%=%該程序完成 16 個脈沖信號的【脈壓、動目標顯示/ 動目標檢測（MTI/MTD ）】%=%程序中根據(jù)每個學(xué)生學(xué)號的末尾三位（依次為XYZ ）來決定仿真參數(shù)，034%目標距離 為 3000 8025 9000+(Y*10+Z)*200 8025，4個目標% 目標速度 為 50 0 (Y*10+X+Z)*6 100%=%closeall ; %關(guān)閉所有圖形clearall;%清除所有變量clc;%=%雷達參數(shù)%=%C=3.0e8;%光速 (m/s)RF=3.140e9/2;%雷達 射頻 1.57GHzLambda=C/RF; %雷達工作 波長PulseNumber=16;%回波脈沖數(shù)Band

2、Width=2.0e6;%發(fā)射信號帶寬帶寬 B=1/ ， 是脈沖寬度TimeWidth=42.0e-6;%發(fā)射信號時寬PRT=240e-6;% 雷達發(fā)射脈沖重復(fù)周期 (s),240us對應(yīng) 1/2*240*300=36000米最大無模糊距離PRF=1/PRT;Fs=2.0e6;%采樣頻率NoisePower=-12;%(dB);%噪聲功率 （目標為 0dB ）% -%SampleNumber=fix(Fs*PRT);%計算一個脈沖周期的采樣點數(shù)480 ；TotalNumber=SampleNumber*PulseNumber;%總的采樣點數(shù)480*16= ；BlindNumber=fix(Fs

3、*TimeWidth);%計算一個脈沖周期的盲區(qū)- 遮擋樣點數(shù) ；%=%目標參數(shù)%=%TargetNumber=4;%目標個數(shù)SigPower(1:TargetNumber)=1 1 1 0.25;%目標功率 , 無量綱TargetDistance(1:TargetNumber)=3000 8025 15800 8025;%目標距離 , 單位 m距離參數(shù)為 3000 8025 9000+(Y*10+Z)*200 8025DelayNumber(1:TargetNumber)=fix(Fs*2*TargetDistance(1:TargetNumber)/C);% 把目標距離換算成 采樣點 （距

4、離門） fix 函數(shù)向 0靠攏取整TargetVelocity (1:TargetNumber)=50 0 204 100; %目標徑向速度 單位 m/s 速度參數(shù)為 50 0 (Y*10+X+Z)*6 100TargetFd(1:TargetNumber)=2*TargetVelocity(1:TargetNumber)/Lambda;%計算目標多卜勒頻移2v/ %=%產(chǎn)生線性調(diào)頻信號%= =%number=fix(Fs*TimeWidth);%回波的采樣點數(shù)=脈壓系數(shù)長度=暫態(tài)點數(shù)目+1ifrem(number,2)=0%rem 求余number=number+1;end%把 number

5、 變?yōu)榕紨?shù)fori=-fix(number/2):fix(number/2)-1Chirp(i+fix(number/2)+1)=exp(j*(pi*(BandWidth/TimeWidth)*(i/Fs)2);%exp(j*fi)*，產(chǎn)生復(fù)數(shù)矩陣Chirpendcoeff=conj(fliplr(Chirp);%把 Chirp矩陣翻轉(zhuǎn)并把復(fù)數(shù)共軛，產(chǎn)生脈壓系數(shù)figure(1);%脈壓系數(shù)的實部plot(real(Chirp);axis(0 90 -1.5 1.5);title( 脈壓系數(shù)實部 );%-產(chǎn)生目標回波串-%-產(chǎn)生前 3個目標的回波串SignalAll=zeros(1,Total

6、Number);%所有脈沖的信號fork=1:TargetNumber-1% 依次產(chǎn)生各個目標SignalTemp=zeros(1,SampleNumber);% 一個 PRT-%, 先填0SignalTemp(DelayNumber(k)+1:DelayNumber(k)+number)=sqrt(SigPower(k)*Chirp;%一個脈沖的 1 個目標（未加多普勒速度）(DelayNumber(k)+1):(DelayNumber(k)+number)Signal=zeros(1,TotalNumber);fori=1:PulseNumber% 16個回波脈沖Signal(i-1)*S

7、ampleNumber+1:i*SampleNumber)=SignalTemp;%每個目標把16 個SignalTemp排在一起endFreqMove=exp(j*2*pi*TargetFd(k)*(0:TotalNumber-1)/Fs);%目標的多普勒速度*時間 = 目標的多普勒相移Signal=Signal.*FreqMove; SignalAll=SignalAll+Signal;%加上多普勒速度后的16 個脈沖 1個目標%加上多普勒速度后的16 個脈沖 4個目標end% %-產(chǎn)生第 4 個目標的回波串-%fi=pi/3;SignalTemp=zeros(1,SampleNumber

8、);% 一個脈沖SignalTemp(DelayNumber(4)+1:DelayNumber(4)+number)=sqrt(SigPower(4)*exp(j*fi)*Chirp;%一個脈沖的 1個目標（未加多普勒速度）Signal=zeros(1,TotalNumber);fori=1:PulseNumberSignal(i-1)*SampleNumber+1:i*SampleNumber)=SignalTemp;endFreqMove=exp(j*2*pi*TargetFd(4)*(0:TotalNumber-1)/Fs);%目標的多普勒速度*時間 = 目標的多普勒相移Signal=S

9、ignal.*FreqMove;SignalAll=SignalAll+Signal;figure(2);subplot(2,1,1);plot(real(SignalAll),r-);title( 目標信號的實部 );gridon ;zoomon ;subplot(2,1,2);plot(imag(SignalAll);title( 目標信號的虛部 );gridon ;zoomon ;%= =%產(chǎn)生系統(tǒng)噪聲信號%=%SystemNoise=normrnd(0,10(NoisePower/10),1,TotalNumber)+j*normrnd(0,10(NoisePower/10),1,To

10、talNumber);%均值為 0 ，標準差為 10(NoisePower/10) 的噪聲 %=%總的回波信號%= =%Echo=SignalAll+SystemNoise;% +SeaClutter+TerraClutter，加噪聲之后的回波fori=1:PulseNumber%在接收機閉鎖期, 接收的回波為0Echo(i-1)*SampleNumber+1:(i-1)*SampleNumber+number)=0;%發(fā)射時接收為0endfigure(3);%加噪聲之后的總回波信號subplot(2,1,1);plot(real(Echo), subplot(2,1,2);plot(imag

11、(Echo);title(r-);title( 總回波信號的實部, 閉鎖期為 總回波信號的虛部, 閉鎖期為 0 );0 );%=時域脈壓=%pc_time0=conv(Echo,coeff);%pc_time0pc_time1=pc_time0(number:TotalNumber+number-1);為 Echo和 coeff的卷積%去掉暫態(tài)點number-1個figure(4);%時域脈壓結(jié)果的幅度subplot(2,1,1);plot(abs(pc_time0),r-);title( 時域脈壓結(jié)果的幅度, 有暫態(tài)點 ); %pc_time0的模的曲線subplot(2,1,2);plot

12、(abs(pc_time1);title( 時域脈壓結(jié)果的幅度, 無暫態(tài)點 );%pc_time1的模的曲線% = =%頻域脈壓Echo_fft=fft(Echo,8192);%理應(yīng)進行 TotalNumber+number-1點 FFT, 但為了提高運算速度,進行了8192點的 FFTcoeff_fft=fft(coeff,8192);pc_fft=Echo_fft.*coeff_fft;pc_freq0=ifft(pc_fft);figure(5);subplot(2,1,1);plot(abs(pc_freq0(1:TotalNumber+number-1);title( 頻域脈壓結(jié)果的

13、幅度, 有前暫態(tài)點 );subplot(2,1,2);plot(abs(pc_time0(1:TotalNumber+number-1)-pc_freq0(1:TotalNumber+number-1),r);title(pc_freq1=pc_freq0(number:TotalNumber+number-1); 時域和頻域脈壓的差別 );%去掉暫態(tài)點number-1個 ,后填充點若干 (8192-number+1-TotalNumber)% =按照脈沖號、距離門號重排數(shù)據(jù)=%fori=1:PulseNumberpc(i,1:SampleNumber)=pc_freq1(i-1)*Sampl

14、eNumber+1:i*SampleNumber);%每個 PRT 為一行，每行480 個采樣點的數(shù)據(jù)endfigure(6);plot(abs(pc(1,:);title( 頻域脈壓結(jié)果的幅度, 沒有暫態(tài)點 );% =MTI（動目標顯示） , 對消靜止目標和低速目標 - 可抑制雜波=%fori=1:PulseNumber-1%滑動對消，少了一個脈沖mti(i,:)=pc(i+1,:)-pc(i,:);endfigure(7);mesh(abs(mti);title(MTI result);% =MTD（動目標檢測） , 區(qū)分不同速度的目標，有測速作用=%mtd=zeros(PulseNumb

15、er,SampleNumber); for i=1:SampleNumberbuff(1:PulseNumber)=pc(1:PulseNumber,i); buff_fft=fft(buff); mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber);endfigure(8);mesh(abs(mtd);title(MTD result);%=虛實矩陣轉(zhuǎn)換=%coeff_fft_c=zeros(1,2*8192);fori=1:8192coeff_fft_c(2*i-1)=real(coeff_fft(i);coeff_fft_c(2*i)=imag(coeff_fft(i);endecho_c=zeros(1,2*TotalNumber);fori=1:TotalNumberecho_c(2*i)=imag(Echo(i);end%=以下是為 DSP程序提供回波數(shù)據(jù)、脈壓系數(shù)=% fo=fopen(F:my studyVisua