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

1、=%該程序完成16個(gè)脈沖信號(hào)的【脈壓、動(dòng)目標(biāo)顯示/動(dòng)目標(biāo)檢測(cè)(MTI/MTD )】% =% 程序中根據(jù)每個(gè)學(xué)生學(xué)號(hào)的末尾三位(依次為 XYZ)來(lái)決定仿真參數(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 ; 聯(lián)閉所有圖形clearall ; %青除所有變量clc;% =%雷達(dá)參數(shù)% =%C=3.0e8;%t速(m/s)RF=3.140e9/2;%M 達(dá)射頻 1.57GHzLambda=C/RF; %雷達(dá)工作波長(zhǎng)PulseNumber=16;刎波脈沖數(shù)BandW

2、idth=2.0e6;%發(fā)射信號(hào)帶寬帶寬B=1/ t, °是脈沖寬度TimeWidth=42.0e-6;%發(fā)射信號(hào)時(shí)寬PRT=240e-6;%雷達(dá)發(fā)射脈沖重復(fù)周期(s),240us對(duì)應(yīng)1/2*240*300=36000米最大無(wú)模糊距離PRF=1/PRT;Fs=2.0e6;%采樣頻率NoisePower=-12; %(dB);% 噪聲功率(目標(biāo)為 0dB)%SampleNumber=fix(Fs*PRT);%bf算一個(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)功率，無(wú)量綱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;%bf 算目標(biāo)多卜勒頻移2v/人 %= =% %產(chǎn)生線性調(diào)頻信號(hào)%= =% number=fix(Fs*TimeWidth);%回波的采樣點(diǎn)數(shù)=脈壓系數(shù)長(zhǎng)度=暫態(tài)點(diǎn)數(shù)目+1if rem(number,2)=0%rem 求余number=number+1; en

5、d咐E number變?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ù)的實(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=ze

6、ros(1,TotalNumber);湎有脈沖的信號(hào)，先填 0for k=1:TargetNumber-1%依次產(chǎn)生各個(gè)目標(biāo)SignalTemp=zeros(1,SampleNumber);% 一個(gè) PRTSignalTemp(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);for i=1:PulseNumber % 16 個(gè)回波脈沖Sig

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

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

9、=Signal.*FreqMove; SignalAll=SignalAll+Signal; figure(2); subplot(2,1,1);plot(real(SignalAll),'r-' );title( '目標(biāo)信號(hào)的實(shí)部');gridon ;zoom on ; subplot(2,1,2);plot(imag(SignalAll);title('目標(biāo)信號(hào)的虛部');grid on ;zoomon ; %= =% %產(chǎn)生系統(tǒng)噪聲信號(hào)%= =% SystemNoise=normrnd(0,10A(NoisePower/10),1,Tota

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

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

12、plot(abs(pc_time0),'r-' );title( '時(shí)域脈壓結(jié)果的幅度，有暫態(tài)點(diǎn)');%pc_time0的模的曲線subplot(2,1,2);plot(abs(pc_time1);title('時(shí)域脈壓結(jié)果的幅度，無(wú)暫態(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.*coef

13、f_fft; pc_freq0=ifft(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( '時(shí)域和頻域脈壓的差另1J');pc_freq1=pc_freq0(number:TotalNumber+

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

15、er-1%骨動(dòng)對(duì)消，少了一個(gè)脈沖mti(i,:)=pc(i+1,:)-pc(i,:);endfigure(7);mesh(abs(mti);title('MTI result' );% =MTD(動(dòng)目標(biāo)檢測(cè))，區(qū)分不同速度的目標(biāo)，有測(cè)速作用=% mtd=zeros(PulseNumber,SampleNumber);for i=1:SampleNumberbuff(1:PulseNumber)=pc(1:PulseNumber,i);buff_fft=fft(buff);mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber);endfigure

16、(8);mesh(abs(mtd);title('MTD result' );%=虛實(shí)矩陣轉(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');% 頻域脈