下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
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% =%close all; %關(guān)閉所有圖形clear all; %清除所有變量clc;% =% 雷達(dá)參數(shù) % =%C=3.0e8; %光速(m/s)RF=3.140e9/2; %雷達(dá)射頻 1.57GHzLambda=C/RF;%雷達(dá)工作波長(zhǎng)PulseNumber=16; %回波脈沖數(shù)B
2、andWidth=2.0e6; %發(fā)射信號(hào)帶寬 帶寬B=1/,是脈沖寬度 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);%計(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)功率,無(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)距離換算成采樣點(diǎn)(距離門) fi
4、x函數(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)頻信號(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;end %把number變?yōu)榕紨?shù)for i
5、=-fix(number/2):fix(number/2)-1 Chirp(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,TotalNumber);%所
6、有脈沖的信號(hào),先填0for k=1:TargetNumber-1 % 依次產(chǎn)生各個(gè)目標(biāo) SignalTemp=zeros(1,SampleNumber);% 一個(gè)PRT SignalTemp(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è)回波脈沖 Signal(i-1)*SampleNu
7、mber+1:i*SampleNumber)=SignalTemp; %每個(gè)目標(biāo)把16個(gè)SignalTemp排在一起 end FreqMove=exp(j*2*pi*TargetFd(k)*(0:TotalNumber-1)/Fs);%目標(biāo)的多普勒速度*時(shí)間=目標(biāo)的多普勒相移 Signal=Signal.*FreqMove;%加上多普勒速度后的16個(gè)脈沖1個(gè)目標(biāo) SignalAll=SignalAll+Signal;%加上多普勒速度后的16個(gè)脈沖4個(gè)目標(biāo)end% %-產(chǎn)生第4個(gè)目標(biāo)的回波串-% fi=pi/3; SignalTemp=zeros(1,SampleNumber);% 一個(gè)脈沖 S
8、ignalTemp(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:PulseNumber Signal(i-1)*SampleNumber+1:i*SampleNumber)=SignalTemp; end FreqMove=exp(j*2*pi*TargetFd(4)*(0:TotalNumber-1)/Fs);%目標(biāo)的多普勒速度*時(shí)間=目標(biāo)的多普勒相移 Signal=Signal
9、.*FreqMove; SignalAll=SignalAll+Signal;figure(2);subplot(2,1,1);plot(real(SignalAll),r-);title(目標(biāo)信號(hào)的實(shí)部);grid on;zoom on;subplot(2,1,2);plot(imag(SignalAll);title(目標(biāo)信號(hào)的虛部);grid on;zoom on;%=% 產(chǎn)生系統(tǒng)噪聲信號(hào) %=%SystemNoise=normrnd(0,10(NoisePower/10),1,TotalNumber)+j*normrnd(0,10(NoisePower/10),1,TotalNumbe
10、r);%均值為0,標(biāo)準(zhǔn)差為10(NoisePower/10)的噪聲%=% 總的回波信號(hào) %=%Echo=SignalAll+SystemNoise;% +SeaClutter+TerraClutter,加噪聲之后的回波for i=1:PulseNumber %在接收機(jī)閉鎖期,接收的回波為0 Echo(i-1)*SampleNumber+1:(i-1)*SampleNumber+number)=0; %發(fā)射時(shí)接收為0endfigure(3);%加噪聲之后的總回波信號(hào)subplot(2,1,1);plot(real(Echo),r-);title(總回波信號(hào)的實(shí)部,閉鎖期為0);subplot(2
11、,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);%時(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(abs(pc_time1);title(時(shí)域
12、脈壓結(jié)果的幅度,無(wú)暫態(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é)果的幅度,有前暫態(tài)點(diǎn));subplot(2,1,2);plot(ab
13、s(pc_time0(1:TotalNumber+number-1)-pc_freq0(1:TotalNumber+number-1),r);title(時(shí)域和頻域脈壓的差別);pc_freq1=pc_freq0(number:TotalNumber+number-1);%去掉暫態(tài)點(diǎn) number-1個(gè),后填充點(diǎn)若干(8192-number+1-TotalNumber)% =按照脈沖號(hào)、距離門號(hào)重排數(shù)據(jù)=%for i=1:PulseNumber pc(i,1:SampleNumber)=pc_freq1(i-1)*SampleNumber+1:i*SampleNumber);%每個(gè)PRT為一行
14、,每行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:PulseNumber-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:SampleNumber buff(1
15、:PulseNumber)=pc(1:PulseNumber,i); buff_fft=fft(buff); mtd(1:PulseNumber,i)=buff_fft(1:PulseNumber);end figure(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);endecho_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_DS
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 花種銷售合同(2篇)
- 民辦學(xué)校勞務(wù)合同文本
- 2024年旅行中介項(xiàng)目建議書
- 2024年智能卡制作發(fā)行機(jī)項(xiàng)目合作計(jì)劃書
- 老舊小區(qū)水電安裝項(xiàng)目合同
- 2024版單位物業(yè)管理與服務(wù)合同(服務(wù)合同)
- 老年大學(xué)裝修家具合同
- 礦泉水運(yùn)輸外包合同樣本
- 食品飲料包機(jī)配送合同
- 生物科技傭金居間協(xié)議范本
- 宣傳片制作合同模板
- 蘑菇種植大棚承包合同范本
- 學(xué)校實(shí)驗(yàn)室安全應(yīng)急演練方案
- 消控室委托管理合同范本
- 小程序商城合作協(xié)議
- 防火道施工組織設(shè)計(jì)方案文件
- 中國(guó)銀行業(yè)務(wù)技能訓(xùn)練指南
- 人工智能概論P(yáng)PT全套完整教學(xué)課件
- 肺癌護(hù)理查房范文肺癌護(hù)理查房范文八篇
- 2021-2022學(xué)年廣東省深圳市南山外國(guó)語(yǔ)學(xué)校高二(上)期中英語(yǔ)試卷
- 前庭功能檢查技術(shù)的臨床應(yīng)用
評(píng)論
0/150
提交評(píng)論