雷達(dá)系統(tǒng)仿真matlab代碼

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

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

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

4、離門） fix 函數(shù)向 0靠攏取整TargetVelocity (1:TargetNumber)=50 0 204 100; %目標(biāo)徑向速度 單位 m/s 速度參數(shù)為 50 0 (Y*10+X+Z)*6 100TargetFd(1:TargetNumber)=2*TargetVelocity(1:TargetNumber)/Lambda;%計(jì)算目標(biāo)多卜勒頻移2v/ %=%產(chǎn)生線性調(diào)頻信號%= =%number=fix(Fs*TimeWidth);%回波的采樣點(diǎn)數(shù)=脈壓系數(shù)長度=暫態(tài)點(diǎn)數(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ù)的實(shí)部plot(real(Chirp);axis(0 90 -1.5 1.5);title( 脈壓系數(shù)實(shí)部 );%-產(chǎn)生目標(biāo)回波串-%-產(chǎn)生前 3個(gè)目標(biāo)的回波串SignalAll=zeros(1,Total

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

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

8、);% 一個(gè)脈沖SignalTemp(DelayNumber(4)+1:DelayNumber(4)+number)=sqrt(SigPower(4)*exp(j*fi)*Chirp;%一個(gè)脈沖的 1個(gè)目標(biāo)（未加多普勒速度）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);%目標(biāo)的多普勒速度*時(shí)間 = 目標(biāo)的多普勒相移Signal=S

9、ignal.*FreqMove;SignalAll=SignalAll+Signal;figure(2);subplot(2,1,1);plot(real(SignalAll),r-);title( 目標(biāo)信號的實(shí)部 );gridon ;zoomon ;subplot(2,1,2);plot(imag(SignalAll);title( 目標(biāo)信號的虛部 );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 ，標(biāo)準(zhǔn)差為 10(NoisePower/10) 的噪聲 %=%總的回波信號%= =%Echo=SignalAll+SystemNoise;% +SeaClutter+TerraClutter，加噪聲之后的回波fori=1:PulseNumber%在接收機(jī)閉鎖期, 接收的回波為0Echo(i-1)*SampleNumber+1:(i-1)*SampleNumber+number)=0;%發(fā)射時(shí)接收為0endfigure(3);%加噪聲之后的總回波信號subplot(2,1,1);plot(real(Echo), subplot(2,1,2);plot(imag

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

12、(abs(pc_time1);title( 時(shí)域脈壓結(jié)果的幅度, 無暫態(tài)點(diǎn) );%pc_time1的模的曲線% = =%頻域脈壓Echo_fft=fft(Echo,8192);%理應(yīng)進(jìn)行 TotalNumber+number-1點(diǎn) FFT, 但為了提高運(yùn)算速度,進(jìn)行了8192點(diǎn)的 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)點(diǎn) );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); 時(shí)域和頻域脈壓的差別 );%去掉暫態(tài)點(diǎn)number-1個(gè) ,后填充點(diǎn)若干 (8192-number+1-TotalNumber)% =按照脈沖號、距離門號重排數(shù)據(jù)=%fori=1:PulseNumberpc(i,1:SampleNumber)=pc_freq1(i-1)*Sampl

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