(2021年整理)雷達系統(tǒng)仿真matlab代碼

1、編輯整理：尊敬的讀者朋友們：這里是精品文檔編輯中心，本文檔內(nèi)容是由我和我的同事精心編輯整理后發(fā)布的，發(fā)布之前我們對文中內(nèi)容進行仔細校對，但是難免會有疏漏的地方，但是任然希望（完整）雷達系統(tǒng)仿真matlab代碼）的內(nèi)容能夠給您的工作和學(xué)習(xí)帶來便利。同時也真誠的希望收到您的建議和反饋，這將是我們進步的源泉，前進的動力。本文可編輯可修改，如果覺得對您有幫助請收藏以便隨時查閱，最后祝您生活愉快業(yè)績進步，以下為（完整）雷達系統(tǒng)仿真matlab代碼的全部內(nèi)容。(完整)雷達系統(tǒng)仿真matlab代碼=%該程序完成16個脈沖信號的【脈壓、動目標(biāo)顯示/動目標(biāo)檢測（MTI/MTD）】=%程序中根據(jù)每個學(xué)生學(xué)號的末尾

2、三位（依次為XYZ）來決定仿真參數(shù)，034%目標(biāo)距離為300080259000+（Y10+Z）2008025,4個目標(biāo)%目標(biāo)速度為500（Y*10+X+Z）6100%=%closeall；%關(guān)閉所有圖形clearall;%清除所有變量clc；%雷達參數(shù)%C=3.0e8;%光速(m/s)RF=3.140e9/2;%雷達射頻1。57GHzLambda二C/RF；%雷達工作波長PulseNumber=16;%回波脈沖數(shù)BandWidth=2°0e6；%發(fā)射信號帶寬帶寬B=1/t,t是脈沖寬度TimeWidth=42。0e-6；%發(fā)射信號時寬PRT=240e6;%雷達發(fā)射脈沖重復(fù)周期(s),

3、240us對應(yīng)1/2大240*300=36000米最大無模糊距離PRF=1/PRT;Fs=2.0e6；%采樣頻率NoisePower=T2;%(dB);%噪聲功率(目標(biāo)為0dB)%-%SampleNumber二fix(Fs*PRT);%計算一個脈沖周期的采樣點數(shù)480;TotalNumber二SampleNumber大PulseNumber;%總的采樣點數(shù)480大16二；BlindNumber二fix(Fs*TimeWidth)；%計算脈沖周期的盲區(qū)一遮擋樣點數(shù)；%=%目標(biāo)參數(shù)%=%TargetNumber=4；%目標(biāo)個數(shù)SigPower(1:TargetNumber)=1110。25;%目標(biāo)

4、功率，無量綱TargetDistance(1:TargetNumber)二30008025158008025；%目標(biāo)距離，單位m距離參數(shù)為300080259000+(Y10+Z)2008025DelayNumber(1：TargetNumber)=fix(Fs*2TargetDistance(1:TargetNumber)/C);%把目標(biāo)距離換算成采樣點(距離門)fix函數(shù)向0靠攏取整TargetVelocity(1:TargetNumber)=500204100；%目標(biāo)徑向速度單位m/s速度參數(shù)為500(Y10+X+Z)6100TargetFd(1:TargetNumber)=2Target

5、Velocity(1:TargetNumber)/Lambda；計算目標(biāo)多卜勒頻移2v/入產(chǎn)生線性調(diào)頻信號%=number二fix(Fs*TimeWidth)；%回波的采樣點數(shù)二脈壓系數(shù)長度二暫態(tài)點數(shù)目+1ifrem(number,2)=0%rem求余number=number+1;end%把門umber變?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(flipl

6、r(Chirp)；%把Chirp矩陣翻轉(zhuǎn)并把復(fù)數(shù)共軛,產(chǎn)生脈壓系數(shù)figured);%脈壓系數(shù)的實部plot(real(Chirp)；axis(090-1.51.5);title('脈壓系數(shù)實部);%-產(chǎn)生目標(biāo)回波串-%-產(chǎn)生前3個目標(biāo)的回波串-%SignalAll二zeros(1,TotalNumber);%所有脈沖的信號，先填0fork=1:TargetNumber-1%依次產(chǎn)生各個目標(biāo)SignalTemp二zeros(1,SampleNumber)；%PRTSignalTemp(DelayNumber(k)+1:DelayNumber(k)+number)=sqrt(SigPow

7、er(k)*Chirp；%脈沖的1個目標(biāo)(未加多普勒速度)(DelayNumber(k)+1)：(DelayNumber(k)+number)Signal=zeros(1,TotalNumber);fori=1:PulseNumber%16個回波脈沖Signal(iT)大SampleNumber+1:i大SampleNumber)二SignalTemp;%每個目標(biāo)把16個SignalTemp排在一起endFreqMove二exp(j大2大pi大TargetFd(k)大(0:TotalNumber一1)/Fs);%目標(biāo)的多普勒速度大時間=目標(biāo)的多普勒相移Signal=Signal。*FreqMo

8、ve；%加上多普勒速度后的16個脈沖1個目標(biāo)SignalAll二SignalAll+Signal；%加上多普勒速度后的16個脈沖4個目標(biāo)end%一一一-一一一-一-一一一一產(chǎn)生第4個目標(biāo)的回波串一一一一-一-%fi=pi/3;SignalTemp二zeros(1,SampleNumber);%脈沖SignalTemp(DelayNumber(4)+1:DelayNumber(4)+number)=sqrt(SigPower(4)*exp(j*fi)*Chirp;%脈沖的1個目標(biāo)(未加多普勒速度)Signal=zeros(1,TotalNumber);fori=1：PulseNumberSign

9、al(i-1)SampleNumber+1：iSampleNumber)=SignalTemp；endFreqMove二exp(j*2*pi*TargetFd(4)大(0:TotalNumber一1)/Fs);%目標(biāo)的多普勒速度大時間二目標(biāo)的多普勒相移Signal=Signal。*FreqMove;SignalAll=SignalAll+Signal;figure(2)；subplot(2,1,1)；plot(real(SignalAll),r一);title('目標(biāo)信號的實部');gridon;zoomon；subplot(2,1,2);plot(imag(SignalAll

10、);title('目標(biāo)信號的虛部');gridon；zoomon;%=%號%產(chǎn)生系統(tǒng)噪聲信%SystemNoise=normrnd(0,10八(NoisePower/10),1,TotalNumber)+j*normrnd(0,10八(NoisePower/10),1,TotalNumber);%均值為0,標(biāo)準(zhǔn)差為10八(NoisePower/10)%的噪聲=%總的回波信號%=%Echo二SignalAll+SystemNoise；%+SeaClutter+TerraClutter,加噪聲之后的回波fori=1:PulseNumber%在接收機閉鎖期，接收的回波為0Echo(i

11、一1)*SampleNumber+1:(iT)大SampleNumber+number)=0;%發(fā)射時接收為0endfigure(3);%加噪聲之后的總回波信號subplot(2,1,1);plot(real(Echo),'r_');title('總回波信號的實部，閉鎖期為0');subplot(2,1,2);plot(imag(Echo)；title('總回波信號的虛部，閉鎖期為0');%=時域脈壓=%pc_time0=conv(Echo,coeff)；%pc_time0為Echo和coeff的卷積pc_time仁pc_time0(number

12、:TotalNumber+number1)；%去掉暫態(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(abs(pc_time1)；title('時域脈壓結(jié)果的幅度，無暫態(tài)點')；%pc_time1的模的曲線%=頻域脈壓=%Echo_fft=fft(Echo,8192);%理應(yīng)進行TotalNumber+number1點FFT,但為了提高運算速

13、度，進行了8192點的FFTcoeff_fft=fft(coeff,8192);（完整）雷達系統(tǒng)仿真matlab代碼pc_fft=Echo_fft.*coeff_fft;pc_freq0=ifft(pc_fft);figure(5)；subplot(2,1,1)；plot(abs(pc_freq0(1：TotalNumber+numberT)；title('頻域脈壓結(jié)果的幅度，有前暫態(tài)點');subplot(2,1,2);plot(abs(pc_time0(1:TotalNumber+number1)-pc_freqO(1:TotalNumber+number一1),'

14、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)大SampleNumber+1:i大SampleNumber)；%每個PRT為一行，每行480個采樣點的數(shù)據(jù)endfigure(6)；plot(abs(pc(1,：);title('頻域脈

15、壓結(jié)果的幅度，沒有暫態(tài)點');%=MTI(動目標(biāo)顯示)，對消靜止目標(biāo)和低速目標(biāo)-一可抑制雜波=%fori=1:PulseNumber-1%滑動對消，少了一個脈沖mti(i，：)=pc(i+1，：)-pc(i，：);endfigure(7)；mesh(abs(mti);title('MTIresult')；%=MTD(動目標(biāo)檢測)，區(qū)分不同速度的目標(biāo)，有測速作用=%mtd=zeros(PulseNumber，SampleNumber)；fori=1：SampleNumberbuff(1：PulseNumber)=pc(1：PulseNumber，i)；buff_fft=f

16、ft(buff);mtd(1：PulseNumber，i)=buff_fft(1:PulseNumber)；endfigure(8)；mesh(abs(mtd);title('MTDresult');%=虛實矩陣轉(zhuǎn)換=%coeff_fft_c=zeros(1，2*8192);fori=1：8192coeff_fft_c(2*i-1)=real(coeff_fft(i);coeff_fft_c(2i)=imag(coeff_fft(i);endecho_c=zeros(1，2*TotalNumber);fori=1:TotalNumberecho_c(2i-1)=real(Echo(i)；echo_c(2i)=imag(Echo(i);end%=二二二二二二二二二二二二二二以下是為dsp程序提供回波數(shù)據(jù)、脈壓系數(shù)=%(完整)雷達系統(tǒng)仿真matlab代碼%fo=fopen('F:mystudyVisual_DSP_testtest_1coeff_ff