[通信工程電子信息工程信息工程論文]超寬帶調制技術研究[論文]-資料下載頁

　　

【正文】 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 圖 調制脈沖序列及其頻譜圖 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 圖 已調脈沖序列及其頻譜圖 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 圖 相關接收后的序列及其頻譜圖 在第三章中介紹了 PPW 的調制原理，本節將其進行了仿真，由 圖 至圖 中可以看出，發送的信息比特是 [1 0 1 0],在調制序列中脈沖與脈沖之間的間隔是固定不變的，在時域圖上并沒有明顯 體現出信息，而已調序列中所有代表一個信息比特“ 1”的脈沖在每個幀中都位于相同的時隙，脈沖波形出現在每個時隙起始點后 ? 秒；而代表信息比特“ 0”的脈沖也是位于每個幀中相同的時隙，但是此時的脈沖波形卻出現在每個時隙起始點處，在時域圖上已經較明顯的體現出了所帶信息。由于調制改變了脈沖與脈沖之間的間隔，信號能量更趨近平滑，信號更加隱蔽。但因數據調制使脈沖間隔變化不大，平滑改動程度并不明顯，梳狀譜仍然存在，且發射的能量集中在脈沖重復頻率的整倍數處。功率譜密度的包絡的形 狀與基本脈沖波形（二階高斯微分波形）的傅氏變換相同。而對于采用相關接收的方式序列通過脈沖的長短來體現信息，功率譜的也呈現梳 狀 。 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