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

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

2、idth=2.0e6;%發(fā)射信號帶寬帶寬B=1/ t, °是脈沖寬度TimeWidth=42.0e-6;%發(fā)射信號時寬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);%bf算一個脈沖周期的采樣點(diǎn)數(shù)480 ；TotalNumber=SampleNumber*PulseNumber;%總的采本點(diǎn)數(shù) 480*16=；BlindNumber=fix(Fs

3、*TimeWidth);%計算一個脈沖周期的盲區(qū)-遮擋樣點(diǎn)數(shù)；%= =% %目標(biāo)參數(shù)%= =% TargetNumber=4;%目標(biāo)個數(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 8025 DelayNumber (1:TargetNumber)=fix(Fs*2*TargetDistance(1:TargetNumber)/C); %把目標(biāo)距離換算成 采樣

4、點(diǎn)(距離門)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;目標(biāo)多卜勒頻移2v/人%十算%= =%產(chǎn)生線性調(diào)頻信號%= =%number=fix(Fs*TimeWidth);%回波的采樣點(diǎn)數(shù)=脈壓系數(shù)長度=暫態(tài)點(diǎn)數(shù)目+1if rem(number,2)=0%rem 求余number=number+1;end咐E num

5、ber變?yōu)榕紨?shù) for i=-fix(number/2):fix(number/2)-1Chirp(i+fix(number/2)+1)=exp(j*(pi*(BandWidth/TimeWidth)*(i/Fs)A2);%exp(j*fi)* ,產(chǎn)生復(fù)數(shù)矩陣Chirp endcoeff=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)生目標(biāo)回波串% 產(chǎn)生前3個目標(biāo)的回波串 %SignalAll=zeros(1,

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

7、)*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.*FreqMove;劾口上多普勒速度后的 16個脈沖1個目標(biāo)SignalAll=SignalAll+Signal;%;口上多普勒速度后的16個脈沖4個目標(biāo)end% % 產(chǎn)生第4個目標(biāo)的回波串%fi=pi/3;SignalTemp=zeros(1,SampleNumber);%

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

9、SignalAll=SignalAll+Signal;%目標(biāo)的多普勒速度figure(2);subplot(2,1,1);plot(real(SignalAll), on ;zoom on ;subplot(2,1,2);plot(imag(SignalAll);title( on ;'r-' );title( '目標(biāo)信號的實部);grid'目標(biāo)信號的虛部');grid on ;zoom%= =%產(chǎn)生系統(tǒng)噪聲信號%= =%SystemNoise=normrnd(0,10A(NoisePower/10),1,TotalNumber)+j*normrnd(0

10、,10A (NoisePower/10),1,TotalNumber);%勻值為0,標(biāo)準(zhǔn)差為10A(NoisePower/10)的噪聲%= =% %總的回波信號%= =% Echo=SignalAll+SystemNoise;% +SeaClutter+TerraClutter,加噪聲之后的回波for i=1:PulseNumber%接收機(jī)閉鎖期，接收的回波為0Echo(i-1)*SampleNumber+1:(i-1)*SampleNumber+number)=0;%發(fā)射時接收為0 end figure(3);%加噪聲之后的總回波信號subplot(2,1,1);plot(real(Echo

11、),'r-' );title( '總回波信號的實部，閉鎖期為 0');subplot(2,1,2);plot(imag(Echo);title('總回波信號的虛部，閉鎖期為 0');%= 時域脈壓 =% pc_time0=conv(Echo,coeff);%pc_time0 為 Echo 禾口 coeff 的卷積pc_time1=pc_time0(number:TotalNumber+number-1);%掉暫態(tài)點(diǎn) number-1 個figure(4);%寸域脈壓結(jié)果的幅度subplot(2,1,1);plot(abs(pc_time0),&#

12、39;r-' );title( '時域脈壓結(jié)果的幅度，有暫態(tài)點(diǎn)');%pc_time0 的模的曲線 subplot(2,1,2);plot(abs(pc_time1);title('時域脈壓結(jié)果的幅度，無暫態(tài)點(diǎn)');%pc_time1 的模的曲線 % =頻域脈壓=% Echo_fft=fft(Echo,8192);%3應(yīng)進(jìn)行 TotalNumber+number-1 點(diǎn) FFT,但為了提高運(yùn)算速度，進(jìn)行了 8192點(diǎn)的FFT coeff_fft=fft(coeff,8192); pc_fft=Echo_fft.*coeff_fft; pc_freq0=i

13、fft(pc_fft); figure(5); subplot(2,1,1);plot(abs(pc_freq0(1:TotalNumber+number-1);title('頻域脈壓結(jié)果的幅度，有前暫態(tài)點(diǎn))；subplot(2,1,2);plot(abs(pc_time0(1:TotalNumber+number-1)-pc_freq0(1:TotalNumber+number-1),'r' );title( '時域和頻域脈壓的差另1J');pc_freq1=pc_freq0(number:TotalNumber+number-1);%掉暫態(tài),點(diǎn) nu

14、mber-1個，后填充點(diǎn)若干(8192-number+1-TotalNumber) % =按照脈沖號、距離門號重排數(shù)據(jù)=% for i=1:PulseNumberpc(i,1:SampleNumber)=pc_freq1(i-1)*SampleNumber+1:i*SampleNumber);%5個 PRT為一行，每行 480 個采樣點(diǎn)的數(shù)據(jù)endfigure(6);plot(abs(pc(1,:);title('頻域脈壓結(jié)果的幅度，沒有暫態(tài)點(diǎn))；% =mti (動目標(biāo)顯示)，對消靜止目標(biāo)和低速目標(biāo)-可抑制雜波=%for i=1:PulseNumber-1%骨動對消，少了一個脈沖mti

15、(i,:)=pc(i+1,:)-pc(i,:);endfigure(7);mesh(abs(mti);title('MTI result' );% =mtd(動目標(biāo)檢測)，區(qū)分不同速度的目標(biāo)，有測速作用=% mtd=zeros(PulseNumber,SampleNumber); for i=1:SampleNumber buff(1:PulseNumber)=pc(1:PulseNumber,i); buff_fft=fft(buff); mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber); end figure(8);mesh(abs(

16、mtd);title('MTD result' );%= 虛實矩陣轉(zhuǎn)換 =% coeff_fft_c=zeros(1,2*8192); for i=1:8192 coeff_fft_c(2*i-1)=real(coeff_fft(i); coeff_fft_c(2*i)=imag(coeff_fft(i); end echo_c=zeros(1,2*TotalNumber); for i=1:TotalNumber echo_c(2*i-1)=real(Echo(i); echo_c(2*i)=imag(Echo(i); end %=以下是為DSP程序提供回波數(shù)據(jù)、脈壓系數(shù)=% % fo=fopen('F:my studyVisual_DSP_testtest_1coeff_fft_c.dat','wt');% 頻域脈壓系數(shù)