[信息與通信]北郵通信原理軟件實(shí)驗(yàn)報(bào)告(編輯修改稿)
2025-04-19 00:44 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 Bs=fs/2。 //系統(tǒng)帶寬t=T/2+[0:N1]/fs。 //時(shí)域采樣點(diǎn)f=Bs+[0:N1]/T。 //頻域采樣點(diǎn)EP1=zeros(1,N)。EP2=zeros(1,N)。L0=input(39。您所需的占空比是？（輸入0~1范圍內(nèi)的數(shù)）：39。)for loop=1:1000 tmp1=zeros(L,M)。 tmp2=zeros(L,M)。 a=(rand(1,M))+0。 //單極性序列 b=sign()。 //雙極性序列 L1=L*L0。 //L0是占空比 //單極性碼 tmp1([1:L1],:)=ones(L1,1)*a。 s1=tmp1(:)39。 S1=t2f(s1,fs)。 //樣本信號(hào)的功率譜密度 P1=abs(S1).^2/T。 //隨機(jī)過(guò)程的功率譜是各個(gè)樣本功率譜的數(shù)學(xué)期望 EP1=EP1*(11/loop)+P1/loop。 //雙極性 tmp2([1:L1],:)=ones(L1,1)*b。 s2=tmp2(:)39。 S2=t2f(s2,fs)。 //樣本信號(hào)的功率譜密度 P2=abs(S2).^2/T。 //隨機(jī)過(guò)程的功率譜是各個(gè)樣本功率譜的數(shù)學(xué)期望 EP2=EP2*(11/loop)+P2/loop。end xset(window,1)plot(t,s1)set(gca(),grid,[1,1]) title(39。單極性碼時(shí)域圖39。)xlabel(39。t39。)ylabel(39。S(t)39。)mtlb_axis([3,3,])。xset(window,2) plot(f,10*log10(EP1+%eps))set(gca(),grid,[1,1]) title(39。單極性碼功率譜圖(dB)39。)xlabel(39。f39。)ylabel(39。功率39。) xset(window,3)plot(t,s2)set(gca(),grid,[1,1]) title(39。雙極性碼時(shí)域圖39。)xlabel(39。t39。)ylabel(39。S(t)39。)mtlb_axis([3,3,])。xset(window,4) plot(f,10*log10(EP2+%eps))set(gca(),grid,[1,1]) title(39。雙極性碼功率譜圖(dB)39。)xlabel(39。f39。)ylabel(39。功率39。) 波形圖：1）占空比為25%單極性歸零碼波形及其功率譜。2） 占空比為25%雙極性歸零碼波形及其功率譜。通過(guò)1）2）組對(duì)比發(fā)現(xiàn)，雙極性歸零碼的功率譜沒(méi)有離散的沖激，而單極性功率譜則有，這是因?yàn)閱螛O性信號(hào)的均值不為零。3） 占空比為50%單極性歸零碼波形及其功率譜。4） 占空比為50%雙極性歸零碼波形及其功率譜。5） 占空比為75%單極性歸零碼波形及其功率譜。6） 占空比為75%雙極性歸零碼波形及其功率譜。當(dāng)信號(hào)占空比增大，可以發(fā)現(xiàn)對(duì)應(yīng)的信號(hào)功率譜變得更加密集。7） 占空比為100%單極性歸零碼波形及其功率譜。8） 占空比為100%雙極性歸零碼波形及其功率譜。占空比為100%的單極性信號(hào)功率譜只在頻率為0處有沖激，這是因?yàn)槠溆嚯x散沖激剛好位于零點(diǎn)，所以看不出來(lái)。 實(shí)驗(yàn)心得：這次實(shí)驗(yàn)比較簡(jiǎn)單，雖然數(shù)據(jù)較多，但是原理相同。通過(guò)實(shí)驗(yàn)，我主要是對(duì)比了單雙極性信號(hào)波形以及信號(hào)功率譜之間的區(qū)別，比較了不同占空比的信號(hào)的功率譜之間的區(qū)別。實(shí)驗(yàn)四利用編程實(shí)現(xiàn)數(shù)字基帶傳輸系統(tǒng)（）1）畫出發(fā)送端輸入碼序列波形和功率譜、發(fā)送濾波器輸出波形和功率譜2）畫出接收端采樣判決后碼序列波形和功率譜、接收濾波器輸出波形和功率譜3）畫出接收濾波器輸出信號(hào)眼圖（在升余弦滾將系數(shù)分別為0、1的情況下）4）分別畫出升余弦滾將系數(shù)為0、1，采樣判決點(diǎn)在眼圖最大處的系統(tǒng)的實(shí)際誤碼曲線（Pe~s/n0曲線），并在同坐標(biāo)系中畫出理論誤碼曲線5）改變采樣點(diǎn)重復(fù)1）~4）。源代碼：clear allexec 。exec 。function eyes(result,L,Na) N=length(result)。 tt=0:1:Na*L。 set(gca(),auto_clear,off) for jj=1:Na*L:NNa*L plot(tt,result(jj:jj+Na*L))。 endset(gca(),auto_clear,on)endfunctionk=input(39。請(qǐng)輸入采樣點(diǎn)數(shù)2^k:（輸入k值）39。)。//輸入采樣點(diǎn)數(shù)aa=。N=2^k。L=8。M=N/L。Rs=2。 //碼率2MbpsTs=1/Rs。 //碼元間隔dt=Ts/L。 //時(shí)域采樣間隔fs=1/dt。df=1/(N*dt)。 //頻域采樣間隔T=N*dt。 //截?cái)鄷r(shí)間Bs=N*df/2。 //系統(tǒng)帶寬t=[T/2+dt/2:dt:T/2]。 //時(shí)域橫坐標(biāo)f=[Bs+df/2:df:Bs]。 //頻域橫坐標(biāo)alpha=input(39。所需要的滾降系數(shù)是（0~1范圍）：39。)。 // 設(shè)置滾降系數(shù)Hcos=zeros(1,N)。 // 升余弦濾波器算法i1=find(abs(f)(1alpha)/(2*Ts)amp。abs(f)=(1+alpha)/(2*Ts))。 Hcos(i1)=Ts/2*(1+cos(%pi*Ts/alpha*(abs(f(i1))(1alpha)/(2*Ts))))。 i2=find(abs(f)=(1alpha)/(2*Ts))。 Hcos(i2)=Ts。 //接受及發(fā)送的根升余弦濾波器GT=sqrt(Hcos)。 GR=GT。 // for loop1=1:20。 //用來(lái)改變信噪比的循環(huán)，方便之后畫信噪比誤碼率曲線 Eb_N0(loop1)=(loop11)。 // 信噪比 eb_n0(loop1)=10^(Eb_N0(loop1)/10)。 EP=1。 n0=EP/eb_n0(loop1)。 //信道的噪聲譜密度 sita=n0*Bs。 // 信道中噪聲功率 n_err=0。 // 誤碼數(shù) EP1=1。//發(fā)送序列功率譜 EP2=1。//zeros
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1