【正文】 t,y1)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Unmodulated pulse train39。)。 grid on。 axis([1e9,10e9 1 1]) subplot(2,2,2)。 plot(f,MY1)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLITUDE39。)。 grid on。 subplot(2,2,3)。 plot(f,20*log10(MY1))。xlabel(39。FREQUENCY39。)。ylabel(39。20LOG10=DB39。)。 grid on。 figure(3) subplot(2,2,1)。 plot(t,y2)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Receiver correlator output39。)。 grid on。 axis([1e9,10e9 1 1]) subplot(2,2,2)。 plot(f,MY2)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLITUDE39。)。 axis([0 7e9 0 .025])。 grid on。 subplot(2,2,3)。 plot(f,20*log10(MY2))。xlabel(39。FREQUENCY39。)。ylabel(39。20LOG10=DB39。)。 grid on。 24 0 5 10x 1 091 0 . 500 . 51T I M EAMPLITUDEU n m o d u l a t e d p u l s e t r a i n0 2 4 6x 1 01000 . 0 20 . 0 40 . 0 6FR E Q U E N C YAMPLITUDE0 2 4 6x 1 010 4 0 0 3 0 0 2 0 0 1 0 00FR E Q U E N C Y20LOG10=DB 圖 調(diào)制脈沖序列及其頻譜圖 0 5 10x 1 091 0 . 500 . 51T I M EAMPLITUDEM o d u l a t e d p u l s e t r a i n0 2 4 6x 1 01000 . 0 20 . 0 40 . 0 6FR E Q U E N C YAMPLITUDE0 2 4 6x 1 010 4 0 0 3 0 0 2 0 0 1 0 00FR E Q U E N C Y20LOG10=DB 圖 已調(diào)脈沖序列及其頻譜圖 25 0 5 10x 1 091 0 . 500 . 51T I M EAMPLITUDER e c e i v e r c o r r e l a t o r o u t p u t0 2 4 6x 1 0900 . 0 10 . 0 2FR E Q U E N C YAMPLITUDE0 2 4 6x 1 010 4 0 0 3 0 0 2 0 0 1 0 00FR E Q U E N C Y20LOG10=DB 圖 相關(guān)接收后的序列及其頻譜圖 在第三章中介紹了 PPW 的調(diào)制原理，本節(jié)將其進(jìn)行了仿真，由 圖 至圖 中可以看出，發(fā)送的信息比特是 [1 0 1 0],在調(diào)制序列中脈沖與脈沖之間的間隔是固定不變的，在時(shí)域圖上并沒(méi)有明顯 體現(xiàn)出信息，而已調(diào)序列中所有代表一個(gè)信息比特“ 1”的脈沖在每個(gè)幀中都位于相同的時(shí)隙，脈沖波形出現(xiàn)在每個(gè)時(shí)隙起始點(diǎn)后 ? 秒；而代表信息比特“ 0”的脈沖也是位于每個(gè)幀中相同的時(shí)隙，但是此時(shí)的脈沖波形卻出現(xiàn)在每個(gè)時(shí)隙起始點(diǎn)處，在時(shí)域圖上已經(jīng)較明顯的體現(xiàn)出了所帶信息。由于調(diào)制改變了脈沖與脈沖之間的間隔，信號(hào)能量更趨近平滑，信號(hào)更加隱蔽。但因數(shù)據(jù)調(diào)制使脈沖間隔變化不大，平滑改動(dòng)程度并不明顯，梳狀譜仍然存在，且發(fā)射的能量集中在脈沖重復(fù)頻率的整倍數(shù)處。功率譜密度的包絡(luò)的形 狀與基本脈沖波形（二階高斯微分波形）的傅氏變換相同。而對(duì)于采用相關(guān)接收的方式序列通過(guò)脈沖的長(zhǎng)短來(lái)體現(xiàn)信息，功率譜的也呈現(xiàn)梳 狀 。 BPSKUWB 仿真與分析 pw1=.2e9。 pw=pw1/。 Fs=100e9。 Fn=Fs/2。 t=1e9:1/Fs:30e9。 26 y =1*(1 4*pi.*((t)/pw).^2).* exp(2*pi.*((t)/pw).^2)。 noise=(1e50)*(randn(size(t)))。 figure(6) plot(t,noise)。 yp=1*y+ ... 1*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... 1*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... 1*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... 1*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)。 B=1。 yum=B*y+ ... B*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... B*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... B*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)+ ... B*(14*pi.*(()/pw).^2).*exp(2*pi.*(()/pw).^2)。 ym=yp+noise。 NFFYM=2.^(ceil(log(length(ym))/log(2)))。 FFTYM=fft(ym,NFFYM)。 NumUniquePts=ceil((NFFYM+1)/2)。 FFTYM=FFTYM(1:NumUniquePts)。 MYM=abs(FFTYM)。 MYM=MYM*2。 MYM(1)=MYM(1)/2。 MYM(length(MYM))=MYM(length(MYM))/2。 MYM=MYM/length(ym)。 f=(0:NumUniquePts1)*2*Fn/NFFYM。 NFFYUM=2.^(ceil(log(length(yum))/log(2)))。 FFTYUM=fft(yum,NFFYUM)。 NumUniquePts=ceil((NFFYUM+1)/2)。 FFTYUM=FFTYUM(1:NumUniquePts)。 27 MYUM=abs(FFTYUM)。 MYUM=MYUM*2。 MYUM(1)=MYUM(1)/2。 MYUM(length(MYUM))=MYUM(length(MYUM))/2。 MYUM=MYUM/length(yum)。 f=(0:NumUniquePts1)*2*Fn/NFFYUM。 NFFYC=2.^(ceil(log(length(yc))/log(2)))。 FFTYC=fft(yc,NFFYC)。 NumUniquePts=ceil((NFFYC+1)/2)。 FFTYC=FFTYC(1:NumUniquePts)。 MYC=abs(FFTYC)。 MYC=MYC*2。 MYC(1)=MYC(1)/2。 MYC(length(MYC))=MYC(length(MYC))/2。 MYC=MYC/length(yc)。 f=(0:NumUniquePts1)*2*Fn/NFFYC。 figure(1) subplot(2,2,1)。 plot(t,ym)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Modulated pulse train39。)。 grid on。 subplot(2,2,2)。 plot(f,MYM)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLITUDE39。)。 grid on。 subplot(2,2,3)。 plot(f,20*log10(MYM))。xlabel(39。FREQUENCY39。)。ylabel(39。20LOG10=DB39。)。 grid on。 figure(2) subplot(2,2,1)。 plot(t,yum)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Unmodulated pulse train39。)。 grid on。 28 axis([1e9,27e9 1 1]) subplot(2,2,2)。 plot(f,MYUM)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLITUDE39。)。 axis([0 10e9 0 .1])。 grid on。 subplot(2,2,3)。 plot(f,20*log10(MYUM))。xlabel(39。FREQUENCY39。)。ylabel(39。20LOG10=DB39。)。 grid on。 figure(3) subplot(2,2,1)。 plot(t,yc)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Receiver correlator outputno LPF39。)。 grid on。 subplot(2,2,2)。 plot(f,MYC)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLITUDE39。)。 grid on。 subplot(2,2,3)。 plot(f,20*log10(MYC))。xlabel(39。FREQUENCY39。)。ylabel(39。20LOG10=DB39。)。 grid on。 pt=.5 H=5。 L=0。 LEN=length(yc)。 for ii=1:LEN。 if yc(ii)=pt。 pv(ii)=H。 else。 pv(ii)=L。 end。 end 。 po=pv。 subplot(2,2,4)。 plot(t,po)。 29 axis([1e9 27e9 1 6]) title(39。Comparator output39。)。 xlabel(39。Frequency39。)。 ylabel(39。Voltage39。)。 grid on。 rc=.2e9。 ht=(1/rc).*exp(t/rc)。 ht=.2e9*ht。 ycfo=filter(yc,1,ht)/Fs。 yn=filter(noise,1,ht)/Fs。 NFFYCFO=2.^(ceil(log(length(ycfo))/log(2)))。 FFTYCFO=fft(ycfo,NFFYCFO)。 NumUniquePts=ceil((NFFYCFO+1)/2)。 FFTYCFO=FFTYCFO(1:NumUniquePts)。 MYCFO=abs(FFTYCFO)。 MYCFO=MYCFO*2。 MYCFO(1)=MYCFO(1)/2。 MYCFO(length(MYCFO))=MYCFO(length(MYCFO))/2。 MYCFO=MYCFO/length(ycfo)。 f=(0:NumUniquePts1)*2*Fn/NFFYCFO。 NFFYN=2.^(ceil(log(length(yn))/log(2)))。 FFTYN=fft(yn,NFFYN)。 NumUniquePts=ceil((NFFYN+1)/2)。 FFTYN=FFTYN(1:NumUniquePts)。 MYN=abs(FFTYN)。 MYN=MYN*2。 MYN(1)=MYN(1)/2。 MYN(length(MYN))=MYN(length(MYN))/2。 MYN=MYN/length(yn)。 f=(0:NumUniquePts1)*2*Fn/NFFYN。 30 figure(4) subplot(2,2,1)。 plot(t,ycfo)。xlabel(39。TIME39。)。ylabel(39。AMPLITUDE39。)。 title(39。Receiver filtered correlator output39。)。 grid on。 subplot(2,2,2)。 plot(f,MYCFO)。xlabel(39。FREQUENCY39。)。ylabel(39。AMPLI