資源描述:
《通信對(duì)抗原理大作業(yè)》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫(kù)。
1、通信對(duì)抗原理大作業(yè)題目:第五題、第八題學(xué)院:電子工程學(xué)院專業(yè):信息對(duì)抗技術(shù)學(xué)號(hào)姓名:題目五通信信號(hào)瞬時(shí)參數(shù)仿真。仿真模擬通信信號(hào)或者數(shù)字通信信號(hào)兩種以上,基于希爾伯特變化方法,分析不同通信信號(hào)的瞬時(shí)包絡(luò)、瞬時(shí)相位、瞬時(shí)頻率、瞬時(shí)譜。畫(huà)出信號(hào)的瞬時(shí)參數(shù)波形。根據(jù)題目要求在此選擇2ASK和2PSK兩種數(shù)字信號(hào),利用希爾伯特變化方法獲取信號(hào)的瞬時(shí)包絡(luò),瞬時(shí)相位,瞬時(shí)頻率,瞬時(shí)譜并給出仿真截圖。1、2ASK信號(hào)仿真對(duì)于二進(jìn)制情況2ASK信號(hào)可表示為:其中m(t)是單極性數(shù)字基帶信號(hào),取值為0和1,碼元寬度為,碼元速率是,是載波頻率。
2、m(t)是復(fù)包絡(luò)瞬時(shí)幅度a(t)和瞬時(shí)相位φ(t)分別為:2ASK信號(hào)瞬時(shí)參數(shù)仿真波形截圖此次仿真使用的信號(hào)參數(shù)為:N=150;fc=500;%載波頻率Fs=150;%采樣頻率t=0:1/fs:(N-1)/fs;fb=10;n_dB=150;在仿真過(guò)程中為了保證信號(hào)的真實(shí)性故添加了噪聲信號(hào)。仿真程序:closeall;clear;echoon;%---------------參數(shù)設(shè)置------------N=150;fc=500;%載波頻率fs=150;%采樣頻率t=0:1/fs:(N-1)/fs;fb=10;n_dB=1
3、50;echooff;%Modulatedsignals.n1=fs/fb;%ThenumberofsamplepersymbolL=N/n1;%ThenumberofsymbolsPN=randsrc(1,L);fork=0:L-1fori=1:n1a(n1*k+i)=PN(k+1);endendforj=1:Nifa(j)==1b(j)=1;elseb(j)=0;endendy=b.*cos(2*pi*fc*t);%signalof2ASKaskzerodb=y;Y=fft(askzerodb);pyy=Y.*conj(
4、Y);figure(1)subplot(211);plot(askzerodb);subplot(212);plot(pyy);[m,n]=size(askzerodb)ask_nos_db=awgn(askzerodb,25);sample_ask=reshape(ask_nos_db,n_dB,n/n_dB);ask_noise=sample_ask;%-----------------瞬時(shí)特征------------------askdata=askzerodb;data=(askdata);%第j行所有信號(hào)的希爾伯特
5、變化%求信號(hào)序列的希爾伯特變換hask=hilbert(data);%求信號(hào)序列的希爾伯特變換%求瞬時(shí)幅度%a_amplitude=sqrtm(y*y.+hask*hask);a_amplitude=abs(hask);%瞬時(shí)幅度f(wàn)igure(2)plot(b);title('瞬時(shí)幅度');phi_phase=angle(hask);%求相位修正值c_k(1)=0;%修正相位初始值fori=2:nifi+1>nbreak;endifphi_phase(i)-phi_phase(i-1)>pic_k(i)=c_k(i-1)-2
6、*pi;elseifphi_phase(i)-phi_phase(i+1)>pic_k(i)=c_k(i-1)+2*pi;elsec_k(i)=c_k(i-1);endendc_k(n)=phi_phase(n);c_k;%去卷疊相位phi_uw_phase=phi_phase+c_k;%去卷疊相位%去相位線性分量后真正的瞬時(shí)相位phi_NL_phase=phi_uw_phase-2*pi*(1:n)*fc/fs;%去相位線性分量后真正的瞬時(shí)相位figure(3)plot(phi_phase);title('瞬時(shí)相位');%
7、求瞬時(shí)頻率f_frequency=1/(2*pi)*diff(phi_uw_phase);%瞬時(shí)頻率figure(4)plot(f_frequency);title('瞬時(shí)頻率');信號(hào)截圖與各瞬時(shí)參量截圖如下:圖1:時(shí)域波形和頻域波形圖2:瞬時(shí)振幅波形圖3:瞬時(shí)相位波形圖4:瞬時(shí)相位波形從上面四張圖可以看出在添加了噪聲后以標(biāo)準(zhǔn)的2ASK信號(hào)出現(xiàn)了干擾,振幅是時(shí)變函數(shù),瞬時(shí)相位受噪聲影響不再是常數(shù)0,并且瞬時(shí)頻率也因噪聲影響在0附近波動(dòng)。2、2PSK信號(hào)瞬時(shí)參數(shù)仿真對(duì)于二進(jìn)制情況2PSK信號(hào)表示為:其中m(t)是雙極性數(shù)字基
8、帶信號(hào),其取值是-1和+1,碼元寬度為T(mén)b,碼元速率是Rb=1/Tb,fc是載波頻率,Dp頻率調(diào)制因子。其中瞬時(shí)幅度瞬時(shí)幅度a(t)和瞬時(shí)相位Φ(t)分別為:此次仿真使用的信號(hào)參數(shù)為:Fc=150;%載波頻率Fs=1200;%采樣頻率為了達(dá)到真實(shí)仿真效果在仿真過(guò)程中添加了高斯白噪聲。仿真程