【正文】 _2psk39。plot(t,e_psk)。figure(2)。)。title(39。subplot(414)。axis([0,5,2,2])。雙極性基帶信號(hào)st339。plot(t,st3)。st3=st1st2。)。title(39。subplot(412)。else江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 st2(k)=1。if st1(k)=1。for n=1:10 if a(n)st2=t。stem(a)。a=round(rand(1,i))。%碼元速率 B=2*fm。fc=。附錄1 2PSK調(diào)制解調(diào)程序及注釋?zhuān)?clear all close all i=10。橫向比較來(lái)看，對(duì)于同一種調(diào)制方式，當(dāng)信噪比相同時(shí)，采用相干解調(diào)方式的誤碼率低于非相干解調(diào)方式的誤碼率；縱向比較來(lái)看，對(duì)2PSK,2DPSK兩種調(diào)制方式若采用同一種解調(diào)方式江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告（相干解調(diào)或非相干解調(diào)），則2PSK的誤碼率最低，2DSPK的誤碼率次之。在信道高斯白噪聲的干擾下，各種二進(jìn)制數(shù)字調(diào)制系統(tǒng)的誤碼率取決于解調(diào)器輸入信噪比，而誤碼率表達(dá)式的形式則取決于解調(diào)方式。本地載波與發(fā)送端載波反向時(shí)，2PSK的解調(diào)波形與2DPSK完全相反、DPSK的誤碼率誤碼率是指接收的碼元數(shù)在傳輸總碼元數(shù)中所占的比例，即：誤碼率錯(cuò)誤碼元數(shù)傳輸總碼元數(shù) PSK相干解調(diào)：誤碼率（r/2）DPSK相干解調(diào)：誤碼率 DPSK差分非相干解調(diào)：誤碼率1/2e^(r)r為信噪比。2PSK信號(hào)載波恢復(fù)過(guò)程中，存在著180176。本次設(shè)計(jì)采用相干解調(diào)。在解調(diào)過(guò)程中，若相干載波產(chǎn)生180176。第二章 設(shè)計(jì)系統(tǒng)江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 兩種解調(diào)方式的原理框圖如圖18和圖19所示。絕對(duì)碼Dn相對(duì)碼BnCnS2dpsk(t)+延時(shí)Ts波形變換Coswc(t)圖17 2DPSK調(diào)制原理框圖 2DPSK解調(diào)原理2DPSK信號(hào)解調(diào)有相干解調(diào)方式和差分相干解調(diào)。實(shí)現(xiàn)二進(jìn)制差分相移鍵控的最常用的方法是：先對(duì)二進(jìn)制數(shù)字基帶信號(hào)進(jìn)行差分編碼，然后對(duì)變換出的差分碼進(jìn)行絕對(duì)調(diào)相即可。238。假設(shè)前后相鄰碼元的載波相位差為Dj，可定義一種數(shù)字信息與Dj之間的關(guān)系為：236。江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告圖16 2DPSK信號(hào)波形圖(a)絕對(duì)碼(b)相對(duì)碼10參考100011011(c) 2DPSK調(diào)制原理二進(jìn)制差分相移鍵控?，F(xiàn)假設(shè)用Φ表示本碼元初相與前一碼元初相之差，并規(guī)定：Φ＝0表示0碼，Φ＝π表示1碼。為了保證2PSK的優(yōu)點(diǎn)，又不會(huì)產(chǎn)生誤碼，將2PSK體制改進(jìn)為二進(jìn)制差分相移鍵控（2DPSK），及相對(duì)相移鍵控。傳輸系統(tǒng)中要保證信息的有效傳輸就必須要有較高的傳輸速率和很低的誤碼率。它不是利用載波相位的絕對(duì)數(shù)值傳送數(shù)字信息，而是用前后碼元的相對(duì)載波相位值傳送數(shù)字信息。這種現(xiàn)象稱(chēng)為2PSK方式的“倒π”現(xiàn)象或“反相工作”。由于2PSK信號(hào)的載波回復(fù)過(guò)程中存在著180176。即正抽樣值判為1，負(fù)抽樣值判為0。經(jīng)過(guò)帶通濾波的信號(hào)在相乘器中與本地載波相乘，然后用低通濾波器濾除高頻分量，在進(jìn)行抽樣判決。兩種方法原理圖分別如圖12和圖13所示。nang(tnTs)e2PSK=234。其時(shí)域233。通常用0176。鍵控法，是用載二進(jìn)制信息的調(diào)制和180176。和 二進(jìn)制數(shù)字基帶0，模擬調(diào)制法用波信號(hào)進(jìn)行調(diào)波的相位變化作為參考基準(zhǔn)的，當(dāng)基帶信號(hào)為0時(shí)相位相對(duì)于初始相位為0176。分別表示信號(hào)的 1 和 兩個(gè)反相的載制。2PSK信號(hào)調(diào)制有兩種方法，即模擬調(diào)制法和鍵控法。二進(jìn)制移相鍵控信號(hào)的典型時(shí)間波形如圖11。關(guān)鍵詞：數(shù)字調(diào)制與解調(diào)；MATLAB；2PSK；2DPSK；江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告第1章 基本工作原理 2PSK原理 2PSK基本原理二進(jìn)制移相鍵控，簡(jiǎn)記為2PSK或BPSK。進(jìn)一步學(xué)習(xí)了MATLAB編程軟件，將MATLAB與通信系統(tǒng)中數(shù)字調(diào)制知識(shí)聯(lián)系起來(lái)，為以后在通信領(lǐng)域?qū)W習(xí)和研究打下了基礎(chǔ)在計(jì)算機(jī)上，運(yùn)用MATLAB軟件來(lái)實(shí)現(xiàn)對(duì)數(shù)字信號(hào)調(diào)制技術(shù)的仿真。)title(‘HDB3碼功率譜密度’)第二篇：通信原理課程設(shè)計(jì)_(基于MATLAB的_2PSK_2DPSK仿真)江西農(nóng)業(yè)大學(xué)通信原理課程設(shè)計(jì)報(bào)告題 目 基于Matlab的相移鍵控仿真設(shè)計(jì)專(zhuān) 業(yè) 電子信息工程學(xué)生姓名 曾凡文學(xué) 號(hào) 20121206江西農(nóng)業(yè)大學(xué)課程設(shè)計(jì)報(bào)告 二 0 一五 年 六 月基于Matlab的2PSK,2DPSK仿真摘要：現(xiàn)代通信系統(tǒng)要求通信距離遠(yuǎn)、通信容量大、傳輸質(zhì)量好，作為其關(guān)鍵技術(shù)之一的調(diào)制技術(shù)一直是研究的一個(gè)重要方向。)ylabel(39。)grid on axis([5,5,50,50])xlabel(39。plot(f,aa,39。s(t)39。t(Tb)39。XTick39。b39。s(t)39。t(Tb)39。XTick39。m39。Ep=(Ep*(ii1)+P)/ii。HDB=t2f(hdb,dt)。end hdb=reshape(hdb,1,N)。for loop=1:L/2 init(loop,:)=a。hdb=zeros(L,M)。for ii=1:30 a=rand(1,M)。f=[Bs+df/2:df:Bs]。Bs=N*df/2。df=1/(N*dt)。Ts=1/Rb。M=N/L。end N=2^k。end end end %函數(shù)調(diào)用 clear all close all k=input(‘取樣點(diǎn)數(shù)=2^k，k=[4]’)。if(i4)if(changed_hdb3(i)*last_one==1)changed_hdb3(i3)=changed_hdb3(i)。changed_hdb3(i)=last_V。else count=count+1。last_one=changed_hdb3(i)。count=0。M=length(pcm)。*：function [changed_hdb3]=Hdb3_encoding(pcm)last_V=1。AMI(t)。AMI功率譜密度39。axis([5,+5,50,50])。subplot(2,1,2)。Ep=(Ep*(ii1)+P)/ii。end AMI=t2f(ami,dt)。for tt=1:Rt*Ts/dt。b(cc)=ones(size(cc))。cc=bb(1:2:length(bb))。n=length(find(a==1))。for ii=1:50 n=1。f=[Bs+df/2:df:Bs]。df=1/(N*dt)T=N*dt Bs=N*df/2。Rb=2。L=64。k=14。)。%采用title命令來(lái)實(shí)現(xiàn)標(biāo)記出各碼元對(duì)應(yīng)的二元信息title(39。plot(t,y)。m=min(y)。end end end y=[y,x(i)]。h=1。h=1。end else for j=1:grid y(j)=0。if(x(i)==1)for j=1:grid y(j)=1。t=0:1/grid:length(x)。、AMI碼及其功率譜的實(shí)現(xiàn) *：傅里葉變換 function X=t2f(x,dt)X=fftshift(fft(x))*dt。)%下面是mlm函數(shù)的調(diào)用 t=[1 1 0 1 0 0 1 0 1 0 0 1]。title(39。plot(x,miler)。x=x/pi。miler=miler./t。miler=(232*cos(x)22*cos(2*x)12*cos(3*x)+5*cos(4*x)+12*cos(5*x)+2*c 21 os(6*x)8*cos(7*x)+2*cos(8*x))./(17+8*cos(8*x))。x=1:1/fs:5。)。%采用title命令來(lái)實(shí)現(xiàn)標(biāo)記出各碼元對(duì)應(yīng)的二元信息title(39。plot(t,y)。m=min(y)。%所有時(shí)間與前一碼元后半時(shí)間值相反 end end end end y=[y,y(i*grid)]。%后半時(shí)間與前半時(shí)間相反 end else if(x(i1)==1)%反之，若前一信息為1，輸入為0 for j=1:grid y(grid*(i1)+j)=y(grid/2*(2*i3)+grid/4)。%碼元持續(xù)時(shí)間內(nèi)為0 end end for i=2:length(x)%開(kāi)始進(jìn)行密勒編碼if(x(i)==1)%若輸入信息為1 for j=1:grid/2 y(grid/2*(2*i2)+j)=y(grid/2*(2*i3)+grid/4)。%前半時(shí)間為0 y(grid/2*(2*i1)+j)=1。%定義時(shí)間序列 i=1。、密勒碼及其功率譜的實(shí)現(xiàn) function y=mlm(x)%該函數(shù)實(shí)現(xiàn)密勒編碼，輸入碼序列為x grid=100。%在命令窗口或M文件中調(diào)用函數(shù) t=[1 1 0 1 0 0 1 0 1 0 0 1]。數(shù)字雙相碼的功率譜39。axis([1,5,0,])。subplot(2,1,2)。mache=mache.*y。y(1)=1。19y=sin(pi*x/2)。fs=100。 1 1 0 1 0 0 1 0 1 0 0 1 39。axis([0,i,M+])。subplot(2,1,1)。M=max(y)。%用從0到1的電平跳變來(lái)表示‘0’ y(grid/2*(2*i1)+j)=1。%用從1到0的電平跳變來(lái)表示‘1’ y(grid/2*(2*i1)+j)=0。t=0:1/grid:length(x)。sjxglm(t)。)。title(39。plot(x,dan2)。dan2=dan2/4。y(1)=2。y=sin(pi*x)。fs=50。 1 1 0 1 0 0 1 0 1 0 0 1 39。axis([0,i,M+])。subplot(2,1,1)。M=max(y)。y(grid/2*(2*i1)+j)=0。y(grid/2*(2*i1)+j)=0。t=0:1/grid:length(x)。sjx(t)。)。title(39。plot(x,dan1)。dan1=dan1/4。y(1)=1。y=sin(pi*x)。fs=100。 1 1 0 1 0 0 1 0 1 0 0 1 39。axis([0,i,M+])。subplot(2,1,1)。M=max(y)。end else for j=1:grid y((i1)*grid+j)=1。t=0:1/grid:length(x)。djxglm(t)。)。title(39。plot(x,dan2)。dan2=dan2/4。y(1)=2。y=sin(pi*x)。fs=50。 1 1 0 1 0 0 1 150 1 0 0 1 39。axis([0,i,M+])。subplot(2,1,1)。M=max(y)。end else for j=1:grid/2 y(grid*(i1)+j)=0。for i=1:length(x)if(x(i)==1)for j=1:grid/2 y(grid/2*(2*i2)+j)=1。、單極性歸零碼及其功率譜實(shí)現(xiàn) function y=djxglm(x)grid=200。%函數(shù)調(diào)用t=[1 1 0 1 0 0 1 0 1 0 0 1]。單極性非歸零碼的功率譜39。axis([1,4,0,])。subplot(2,1,2)。dan1=dan1/4。14 y(1)=1。y=sin(pi*x)。fs=100。 1 1 0 1 0 0 1 0 1 0 0 1 39。axis([0,i,M+])。plot(t,y)。m=min(y)。end end end y=[y,x(i)]。for i=1:length(x)if(x(i)==1)for j=1:grid y((i1)*grid+j)=1。[h,l]=huffman(p)、單極性非歸零碼以及功率譜的實(shí)現(xiàn) function y=djx(x)%本函數(shù)實(shí)現(xiàn)將輸入的一段二進(jìn)制代碼編碼為相應(yīng)的單極性非歸零碼輸出 %輸入x為二進(jìn)制碼，輸出y為編好的碼%給出計(jì)算每一個(gè)碼元的點(diǎn)數(shù)，因?yàn)槲覀冎挥杏秒x散的點(diǎn)來(lái)得出連續(xù)的函數(shù)表示 grid=300。end l=sum(p.*ll)。end for i=1:n h(i,1:n)=c(1,n*(find(m(1,:)==i)1)+1:find(m(1,:)==i)*n)。for j=1:i1c(ni,(j+1)*n+1:(j+2)*n)=c(ni+1,n*(find(m(ni+1,:)==j+1)1)+1:n*find(m(ni+1,:)==j+1))。c(ni,2*n)=39。039。for i=2:n1c(ni,1:n1)=c(ni+1,n*(find(m(ni+1,:)==1))(n2):n*(find(m(ni+1,:)==1)))。c(n1,2*n)=39。end c(n1,n)=39。q=[q(1)+q(2),q(3:n),1]。for i=1:n1 [q,l]=sort(q)。q=p。Not a ,ponent do not add to 139。c1=hxhf1 %信道容量、哈弗曼編碼的實(shí)現(xiàn) p=[1/4,1/4,1/4,1/8,1/8]。hf1=hmessage(x,f1,4,6)。%利用函數(shù)dmessage來(lái)求信源的熵，利用函數(shù)hmessage來(lái)求平均互信息量 x=[,]。平均互信息量為：39。r=sum。end。%求平均互信息量 if t~=0 sum=sumt*log(f(i,j))/log(2)。r %直接在命令窗口中調(diào)用dmessage函數(shù)，也可以在其它函數(shù)中調(diào)用 dmessage([,],4)、離散信道容量的計(jì)算 %求互信息的函數(shù)文件如下 function r=hmessage(x,f,nx,my)%x為輸出的信源分布，f為轉(zhuǎn)移概率矩陣，nx為輸出的符號(hào)可選個(gè)數(shù) %my為輸出的符號(hào)個(gè)數(shù) sum=0。此離散信源的平均信息量為：39。for i=1:n r=rx(i)*log(x(i))/log(2)。%產(chǎn)生N（0，1）高斯源 px=1/sqrt(2*pi)*exp((x1).^2/2)。 HDB3碼及其功率譜 HDB3碼碼型和功率譜由此圖的功率譜與AMI碼的功率譜比較可知，HDB3碼與AMI碼是相似的，HDB3碼是AMI碼的變形。