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

1、%該程序完成16個脈沖信號的【 脈壓、動目標顯示/動目標檢測(MTI/MTD)】% =%程序中根據(jù)每個學生學號的末尾三位(依次為XYZ來決定仿真參數(shù),034%目標距離 為3000 8025 9000+(Y*10+Z)*200 8025，4個目標%目標速度 為50 0 (Y*10+X+Z)*6 100% close all ; %關(guān)閉所有圖形clear all ; %清除所有變量clc;%雷達參C=; %光速(m/s)RF=2; %雷達 射頻Lambda二二C/RF；雷達工作波長PulseNumber=16; %回波脈沖數(shù)Ban dWidth二；二；發(fā)射信號帶寬 帶寬B=1/T,T是脈沖寬度Ti

2、meWidth=; %發(fā)射信號時寬PRT=240e-6；%雷達發(fā)射脈沖重復(fù)周期(s),240us對應(yīng)1/2*240*300=36000米最大無模糊距離PRF=1/PRT;Fs=; %采樣頻率NoisePower=-12; %(dB);%噪聲功率(目標為OdB)% -%SampleNumber=fix(Fs*PRT); %計算 一個脈沖周期的采樣點數(shù)480；TotalNumber二二SampleNumber*PulseNumber%總的采樣點數(shù)480*16二二;BlindNumber=fix(Fs*TimeWidth);%計算一個脈沖周期的盲區(qū)-遮擋樣點數(shù)；%=目標參=%TargetNumber

3、=4; %目標個數(shù)SigPower(1:TargetNumber)=1 1 1 ;%目標功率,無量綱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); %把目標距離換算成 采樣點(距離門)fix函數(shù)向0靠攏取整TargetVelocity (1:TargetNumber)=50 0 204 100;%目標徑向

4、速度 單位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ù)目+1if rem(number,2)=0%em求余number=number+1;end%把numbe變?yōu)榕紨?shù)for i=-fix(number/2):fix(number/2)-1Chirp(i+fix (n umber/2)+1)=ex p(j*

5、( pi *(Ba ndWidth/TimeWidth)*(i/FsF2 );%exp(j*fi)*， 產(chǎn)生復(fù)數(shù)矩陣Chirp end coeff=conj(fliplr(Chirp);系數(shù)figure(1); %脈壓系數(shù)的實部plot(real(Chirp);axis(090 );title(脈壓系數(shù)實部 );%-產(chǎn)生目標回波串 -%-產(chǎn)生前3個目標的回波串 -%SignalAll=zeros(1,TotalNumber);%所有脈沖的信號,先填0for k=1:TargetNumber-1 %依次產(chǎn)生各個目標Sig nalTem p=zeros(1,Sa mp leNumber); %一個

6、PRTSignalTemp(DelayNumber(k)+1:DelayNumber(k)+number)=sqrt(SigPower( k)*Chirp; %一個脈沖的1個目標(未加多普勒速度)(DelayNumber(k)+1):(DelayNumber(k)+number)Signal=zeros(1,TotalNumber);for i=1:PulseNumber % 16個回波脈沖Signal(i-1)*SampleNumber+1:i*SampleNumber)=SignalTemp;個目標把16個Sig nalTe mp排在一起endFreqMove=exp(j*2*pi*Tar

7、getFd(k)*(0:TotalNumber-1)/Fs);多普勒速度*時間=目標的多普勒相移Signal=Signal.*FreqMove;SignalAll=SignalAll+Signal;end% %-fi=pi/3;SignalTemp=zeros(1,SampleNumber); %一個脈沖SignalTemp(DelayNumber(4)+1:DelayNumber(4)+number)=sqrt(SigPower(4)*exp(j*fi)*Chirp;%一個脈沖的1個目標(未加多普勒速度)Signal=zeros(1,TotalNumber);for i=1:PulseNum

8、berSignal(i-1)*SampleNumber+1:i*SampleNumber)=SignalTemp;end%把Chirp矩陣翻轉(zhuǎn)并把復(fù)數(shù)共軛，產(chǎn)生脈壓%每%目標的%加上多普勒速度后的16個脈沖1個目標%加上多普勒速度后的16個脈沖4個目標產(chǎn)生第4個目標的回波串產(chǎn)生系統(tǒng)噪聲信=%SystemNoise二二no rmmd(0,10八八(Noise Po wer/10),1,TotalNumber)+j* normrnd(0,10八八(Noise Power/10),1,TotalNumber);%勻值為0,標準差為10八八(NoisePower/10)的噪聲%=總的回波信=%Echo

9、=SignalAll+SystemNoise; % +SeaClutter+TerraClutter,加噪聲之后 的回波for i=1:PulseNumber%在接收機閉鎖期,接收的回波為0Echo(i-1)*SampleNumber+1:(i-1)*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 );FreqMove=exp(j

10、*2*pi*TargetFd(4)*(0:TotalNumber-1)/Fs);多普勒速度*時間=目標的多普勒相移Signal=Signal.*FreqMove;SignalAll=SignalAll+Signal;figure(2); subplot(2,1,1);plot(real(SignalAll), );gridon;zoom on;subplot(2,1,2);plot(imag(SignalAll);title( on;zoom on;%=%目標的r- );title( 目標信號的實部目標信號的虛部 );grid總回波信號的虛部,閉鎖期為%=時=域=脈壓=%pc_time0=co

11、nv(Echo,coeff);%pc_time0為Echo和coeff的卷積=頻=域=脈壓=%Echo_fft=fft(Echo,8192);%理應(yīng)進行TotalNumbe葉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é)果的幅度,有前暫態(tài)點 );subplot(2

12、,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-numbe葉1-TotalNumber)% =按=照脈沖號、距離門號重排數(shù)據(jù)=%for i=1:PulseNumber pc(i,1:SampleNumber)=pc_freq1(i-1)*SampleNumbe葉1:i*SampleN

13、umber); %每個PR為一行，每行480個采樣點 的數(shù)據(jù)endfigure(6);plot(abs(pc(1,:);title(頻域脈壓結(jié)果的幅度,沒有暫態(tài)點 );% =M(TI動目標顯示),對消靜止目標和低速目標-可抑 制雜波=%for i=1:PulseNumber-1%滑動對消，少了一個脈沖mti(i,:)=pc(i+1,:)-pc(i,:);endfigure(7);mesh(abs(mti);title( MTI result );pc_time1=pc_time0(number:TotalNumber+number-1);%去掉暫態(tài)點number-1個figure(4); %時

14、域脈壓結(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(態(tài)點);%pc_time1的模的曲線% =時域脈壓結(jié)果的幅度,無暫% =MTD目標檢測)，區(qū)分不同速度的目標，有測速作 用=% mtd=zeros(PulseNumber,SampleNumber);for i=1:SampleNumber buff(1:PulseNumber)=pc(1:PulseNumber,i);buff_fft=fft(b

15、uff);mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber);endfigure(8);mesh(abs(mtd);title(%=以下是為DSPS序提供回波數(shù)據(jù)、=% fo=fopen(F:my studyVisual_DSP_testtest_1,wt);%系數(shù)% for i=1:2*8192% fprintf(fo,%f,rn,coeff_fft_c(i);% end% fclose(fo);% fo=fopen(F:my studyVisual_DSP_testtest_1,wt);%16的% for i=1:2*TotalNumber% fprintf(fo,%f,rn,echo_c(i);% end% fclose(fo);MTD result );=虛=