【正文】 1 1 1 1 0 0 1 0 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 19 1 1 1 1 1 1 1 0 1 1 1 1 1 0 0 1 1 1 1 1 0 0 1 0 1 1 1 0 0 0 1 1 1 0 0 0 0 1 0 1 PCM 解碼結(jié)果 dn = Columns 1 through 11 0 0 Columns 12 through 21 圖 48 到達(dá)信宿波形 對(duì)比原始信號(hào)正弦波型可知，到達(dá)信宿的波形在小信號(hào)部分失真小而在大信號(hào)部分失真明顯，這符合 A 律 13 折線(xiàn)的編碼規(guī)則，同時(shí)信源編碼和信道編碼在加入加性高斯白噪聲的情況下，譯碼信號(hào)基本正確，實(shí)現(xiàn)了整個(gè)通信系統(tǒng)的工作流程，建模測(cè)試成功。 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 20 5 分析與小結(jié) 本設(shè)計(jì)結(jié)合 PCM 的抽樣、量化、編碼原理把模擬信號(hào)轉(zhuǎn)換數(shù)字信號(hào)、漢明碼編譯碼原理實(shí)現(xiàn)信道編譯碼、 HDB3 編譯碼原理用于實(shí)現(xiàn)信源編譯碼、 PSK 調(diào)制解調(diào)原理用于實(shí)現(xiàn)數(shù)字調(diào)制與解調(diào)、 AWGN 原理用于信道和噪聲源的近似，利用 MATLAB 軟件編程，完成了整個(gè)通信系統(tǒng)編碼仿真分析 。 本次課程設(shè)計(jì)主要使用了 MATLAB 軟件，使我對(duì) MATLAB 的運(yùn)用更加熟練。對(duì)循環(huán)碼的設(shè)計(jì)也使我投入到對(duì)課本知識(shí)的再一次探討中去，從不同的角度講所學(xué)應(yīng)用到實(shí)踐中去。加深了我對(duì)這部分知識(shí)的理解，鍛煉了對(duì)于資料的收集整理運(yùn)用的能力。對(duì) MATLAB 自帶幫助文檔的運(yùn)用更是提高了我對(duì)英語(yǔ)學(xué)習(xí)重要性的認(rèn)知。 設(shè)計(jì)初，查閱資料，構(gòu)建基本模型，完成程序主要流程并成功運(yùn)行下來(lái)耗時(shí)一天，這個(gè)過(guò)程并不漫長(zhǎng)，在學(xué)校基于 MATLAB 的 DSP 實(shí)驗(yàn)課上，我積累了一定的基礎(chǔ)認(rèn)知，可以初步簡(jiǎn)單運(yùn)用該軟件，為這一環(huán)節(jié)打下了基礎(chǔ)。然后接下來(lái)日 子將重點(diǎn)放在各個(gè)模塊的精細(xì)改進(jìn)，各種程序之間的優(yōu)劣，編程效率也納入自己的考慮之中，這一階段我忙于對(duì)各個(gè)部分功能以及網(wǎng)上資料的測(cè)試，真正印證了“矩陣實(shí)驗(yàn)室”的意義。此階段耗時(shí)最長(zhǎng)，收獲也最多，改變了平時(shí)學(xué)習(xí)中浮于表面，毛毛躁躁的習(xí)慣。 撰寫(xiě)報(bào)告過(guò)程是自己對(duì)這次設(shè)計(jì)的一個(gè)完整回顧，又從一個(gè)新的角度整理實(shí)驗(yàn)數(shù)據(jù)與資料。將林林總總的實(shí)驗(yàn)程序、數(shù)據(jù)的整合過(guò)程中鍛煉了自己的綜合素質(zhì)。對(duì)于問(wèn)題的考慮還不夠全面，諸如沒(méi)有備份程序，修改后程序已經(jīng)面目全非也使自己的心情極為煩躁。克服這些問(wèn)題其實(shí)也可以學(xué)到其他方面的知識(shí)，終將使 自己的能力得到全方位的提高。 誠(chéng)然，自己的設(shè)計(jì)還有許許多多不完善的地方，各個(gè)方面也不一定面面俱到，但總體來(lái)講，過(guò)程中的收獲大于結(jié)果，更是培養(yǎng)了自己對(duì)于學(xué)習(xí)的興趣。 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 21 參考文獻(xiàn) [1] 樊昌信 曹麗娜，《通信原理第六版》，國(guó)防工業(yè)出版社， 2021 [2] 周開(kāi)利，《鄧春暉主編 MATLAB 基礎(chǔ)及其應(yīng)用教程》，北京大學(xué)出版社， 2021 [3] 董振海 ,《 精通 MATLAB 7 編程與數(shù)據(jù)庫(kù)應(yīng)用 》， 電子工業(yè)出版社 ,2021 [4] 陳懷琛，《 MATLAB 在電子信息課程中的應(yīng)用（第二版）》，電子工業(yè)出版社， 2021 [5] 張德豐 主編，《 MATLAB 通信工程仿真》， 機(jī)械工業(yè)出版社 ， 2021 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 22 附錄 一、主程序： close all。 clear all。 figure(1)。set(gcf,39。color39。,39。white39。) t0=0::。 x=sin(8000*pi*t0)。 subplot(2,1,1) plot(t0,x) axis([0 1 1])。 title(39。原始信號(hào)波形 39。)。 t=0::。 y=sin(8000*pi*t)。 subplot(2,1,2) stem(t,y) axis([0 1 1])。 title(39。采樣 39。)。 for i=1:21 gn(i)=y(i)。 end gn z=line13(y)。 xn=pcmcode(z)。 K=4。N=7。 msg=reshape(xn,[],4)。 yn=encode(msg,N,K)。 yn1=reshape(yn,1,[]) 。 zn=hdb3(yn1)。 zn1=abs(zn)。 pn=psk(zn1)。 SNR=100。 an = awgn(pn,SNR)。 subplot(2,1,2) plot(abs(fft(an))) axis([0 100*length(zn) 0 400])。 title(39。加噪后信號(hào)頻譜 39。)。 bn=depsk(an,zn1)。 mn=dehdb3(zn)。 mn1=reshape(mn,[],7)。 kn=decode(mn1,N,K)。 kn1=reshape(kn,[],8)。 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 23 dn=depcm(kn1)。 out=deline13(dn)。 figure(9)。set(gcf,39。color39。,39。white39。) plot(out)。 m=0。 for i=1:21 e(i)=dn(i)gn(i)。 m=m+e(i)。 end jun=m/21。 s=0。 for i=1:21 s=(e(i)jun)^2+s。 end 二、相關(guān)函數(shù)： %13 折線(xiàn)量化 function y=line13(x) x=x/max(x)。 z=sign(x)。 x=abs(x)。 for i=1:length(x) if((x(i)=0)amp。amp。(x(i)1/64)) y(i)=16*x(i)。 else if((x(i)=1/64)amp。amp。(x(i)1/32)) y(i)=8*x(i)+1/8。 else if((x(i)=1/32)amp。amp。(x(i)1/16)) y(i)=4*x(i)+2/8。 else if((x(i)=1/16)amp。amp。(x(i)1/8)) y(i)=2*x(i)+3/8。 else if((x(i)=1/8)amp。amp。(x(i)1/4)) y(i)=x(i)+4/8。 else if((x(i)=1/4)amp。amp。(x(i)1/2)) y(i)=1/2*x(i)+5/8。 else if((x(i)=1/2)amp。amp。(x(i)=1)) y(i)=1/4*x(i)+6/8。 end 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 24 end end end end end end end y=z.*y。 %PCM 編碼 程序 function f=pcmcode(y) f=zeros(length(y),8)。 z=sign(y)。 y=y.*128。 y=fix(y)。 y=abs(y)。 for i=1:length(y) if(y(i)==128) y(i)=。 end end for i=1:length(y) for j=6:1:0 f(i,8j)=fix(y(i)/2^j)。 y(i)=mod(y(i),(2^j))。 end end for i=1:length(y)。 if(z(i)==1) f(i,1)=0。 else f(i,1)=1。 end end %HDB3 編碼程序 function y=hdb3(x) n=length(x)。 last_V=1。 last_one=1。 y=zeros(size(x))。 count=0。 for i=1:n if x(i)==1 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 25 y(i)=last_one。 last_one=y(i)。 count=0。 else count=count+1。 if count==4 count=0。 y(i)=last_V。 last_V=y(i)。 if y(i)*last_one==1 y(i3)=y(i)。 end last_one=y(i)。 end end end figure(2)。set(gcf,39。color39。,39。white39。) subplot(2,1,1)。 a=x。 i=0:n1。 stairs(i,a)。 axis([0,n,])。 title(39。原碼型 39。)。 xlabel(39。x39。)。 ylabel(39。y39。)。 grid on subplot(2,1,2)。 a=y。 i=0:n1。 stairs(i,a)。 axis([0,n,2,2])。 title(39。HDB3碼型 39。)。 xlabel(39。x39。)。 ylabel(39。y39。)。 grid on %PSK 調(diào)制程序 function p=psk(g) cp=[]。 mod1=[]。 f=2*2*pi。 t=0:2*pi/99:2*pi。 for n=1:length(g) if g(n)==0 武漢理工大學(xué)《 學(xué)科基礎(chǔ)課群綜合訓(xùn)練 》 報(bào)告 26 A=1.*ones(1,100)。 elseif g(n)==1 A=ones(1,100)。 end cp=[cp A]。 c=cos(f*t)。 mod1=[mod1 c]。 end figure(3)。set(gcf,39。color39。,39。white39。) subplot(2,1,1) plot(cp)。 grid on axis([0 100*length(g) 2 2])。 title(39。雙極型信號(hào)序列 39。) p=cp.*mod1。 subplot(2,1,2) plot(p) grid on axis([0 2021 2 2]) title(39。psk調(diào)制信號(hào) 39。) figure(4)。set(gcf,39。color39。,39。white39。) subplot(2,1,1) plot(abs(fft