【正文】 %把矩陣變成 1 行 800 列的矩陣 errors2=reshape(errors2,1。%發(fā)現(xiàn)錯誤讓其值為 1 errors2(find(c33~=c7))=1。 errors2=zeros(100,8)。)。%hamming 信道譯碼 200*4 c33=decode(rx2,15,8,39。hamming/binary39。rx2(find(rx20))=0。rx1(find(rx10))=0。%加噪聲 rx2=tx2+sigma*randn(100,15)。 bit1(q)=biterrors/(100*8)。%發(fā)現(xiàn)錯誤讓 error 為 1 errors=reshape(errors,1,800)。%把 m3 變成 8 行 100 列的矩陣 errors=zeros(100,8)。 通信原理課程 設(shè)計 第 33 頁 共 33 頁 m3=reshape(m3,8,100)。)。 m3=de2bi(qq2,2,39。 qq2(find(qq2==3))=2。(qq20)))=1。%判決，解調(diào) qq2(find(qq2=2))=3。%加噪聲 qq2(find((qq2=0)amp。 sigma=sqrt(r1)。 r1=10.^(dB(q)/10)。biterrors1=0。%調(diào)制 errorbit=0。tx2=c3。%(15,8)循環(huán)碼編碼 tx1=c2。cyclic/binary39。)。%c1=200*4,前兩行對應(yīng) c 第一行 c2=encode(c1,7 ,4,39。 c1=reshape(c1,4,200)。 c1=c.39。 m2(find(m2==0))=1。leftmsb39。m1=m1.39。%調(diào)制 m1=m.39。else 通信原理課程 設(shè)計 第 32 頁 共 33 頁 u(5:8)=[1 1 1 1] end c(i,1:8)=u(1:8)。 p=dec2bin(t,4)48。st=1024。u(4)=1。y=2048) u(2)=1。st=512。u(4)=0。y1024) u(2)=1。st=256。u(4)=1。y512) u(2)=1。st=128。u(4)=0。y256) u(2)=1。st=64。u(4)=1。y128) u(2)=0。st=32。u(4)=0。y64) u(2)=0。st=16。u(4)=1。y32) u(2)=0。st=0。u(4)=0。y16)%段落碼判斷 u(2)=0。 end y=abs(y)。 u=[0 0 0 0 0 0 0 0 ] if(y0)%極值碼判斷 u(1)=1。%歸一化 s2=s1./(1/2048)。 s=A*sin(w*n*T)。fs=8*10^3。%語音信號頻率 2K，幅值 t=0:: 通信原理課程 設(shè)計 第 31 頁 共 33 頁 y=A*sin(w*t)。f=2020。 out(i)=ss*(st+dt)/4096*v。 st=slot(tmp)。 for i=1:n/8 ss=2*in(1,i)1。step(7)=64。 step(5)=16。step(3)=4。 step(1)=2。slot(7)=1024。 slot(5)=256。slot(3)=64。 slot(1)=0。 6. 非均勻量化 PCM 譯碼 function[out]=pcm_decode(in,v) n=length(in)。 end end out=reshape(out39。 t=dec2bin(tmp,4)48。 通信原理課程 設(shè)計 第 30 頁 共 33 頁 end if abs(x(i))=4096 out(i,2:8)=[1 1 1 1 1 1 1]。step=128。out(i,3)=1。st=2048。out(i,4)=1。abs(x(i))4096 out(i,2)=1。 elseif abs(x(i))=2048amp。step=64。out(i,3)=1。amp。st=512。out(i,4)=1。abs(x(i))1024 out(i,2)=1。 elseif abs(x(i))=512amp。step=16。out(i,3)=0。amp。st=128。out(i,4)=1。abs(x(i))256 out(i,2)=0。 elseif abs(x(i))=128amp。step=4。out(i,3)=1。amp。st=32。out(i,4)=1。abs(x(i))64 out(i,2)=0。 elseif abs(x(i))=32amp。step=2。out(i,3)=0。amp。 else out(i,1)=0。 end end 通信原理課程 設(shè)計 第 29 頁 共 33 頁 5. 非均勻量化 PCM 編碼 function[out]=pcm_encode(x) n=length(x)。if(z(i)==1) f(i,1)=0。 y(i)=mod(y(i),(2^j))。for i=1:length(y) if(y(i)==128) y(i)=。y=fix(y)。z=sign(y)。 end end end end end end end end y=z.*y。amp。(x(i)1/2)) y(i)=1/2*x(i)+5/8。 else if((x(i)=1/4)amp。amp。(x(i)1/8)) y(i)=2*x(i)+3/8。 else if((x(i)=1/16)amp。amp。(x(i)1/32)) 通信原理課程 設(shè)計 第 28 頁 共 33 頁 y(i)=8*x(i)+1/8。 else if((x(i)=1/64)amp。amp。 x=abs(x)。 y=sin(8000*pi*t) z=line13(y) c=pcmcode(z) function y=line13(x) x=x/max(x)。 else y(i)=(1+log(a*x(i)))/(1+log(a))。 else y(i)=(1+log(a*x(i)))/(1+log(a))。) function y=a_pcm(x,a) t=1/a。) title(39。) ylabel(39。subplot(2,1,2) plot(t,z) 通信原理課程 設(shè)計 第 27 頁 共 33 頁 axis([0 ]) xlabel(39。原始信號 39。原始信號 39。幅度 39。時間 39。 y=sin(8000*pi*t)。 end end end end nq=V^2/(3*L^2)。 end end end for j=L/2+2:L+1 if(flag==0) if(f(i)p(j)) h(i)=p(j)。for j=2:L/2+1 if(flag==0) if(f(i)p(j)) h(i)=p(j1)。end if f(i)=V,h(i)=V。 for i=1:L+1,p(i)=V+(i1)*t。t=2*V/L。均勻量化后的信號 39。 axis([0,4*pi,])。幅度 39。)。 xlabel(39。 subplot(2,1,2)。原始信號 39。 axis([0,4*pi,])。幅度 39。)。 xlabel(39。subplot(2,1,1)。y=sin(t)。df=fs/n。 M=fft(m,n)。end n2=length(m)。if nargin==2 n1=0。axis([500,500,])。抽樣失真時的信號頻譜 39。axis([500,500,0,])。幅值 39。)。 xlabel(39。 axis([,0,1])。抽樣失真時的信號頻域波形 (fs=100Hz)39。)。ylabel(39。時間 39。%stem(t2,m2)。bo39。 figure(3) 通信原理課程 設(shè)計 第 25 頁 共 33 頁 subplot(2,1,1)。N2=[M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2]。 %定義信號 [M2,mn2,df2]=fft_seq(m2,ts2,df)。 m2(w2)=1。 m2=x2./(200*t2+eps)。 t2=[t0/2:ts2:t0/2]。 %信號持續(xù)的時間 ts2=。axis([500,500,])。抽樣正常時的信號頻譜 39。axis([500,500,0,])。幅值 39。)。 xlabel(39。axis([,0,1])。抽樣正常時的信號頻域波形 (fs=200Hz)39。)。ylabel(39。時間 39。%stem(t1,m1)。bo39。 figure(2) subplot(2,1,1)。N1=[M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1,M1]。 %定義信號 [M1,mn1,df1]=fft_seq(m1,ts1,df)。 m1(w1)=1。 m1=x1./(200*t1+eps)。 t1=[t0/2:ts1:t0/2]。 %信號持續(xù)的時間 ts1=。)。title(39。)。ylabel(39。頻率 39。 plot(f,abs(fftshift(M)))。 axis([,0,])。原始信號的頻域圖 (fh=200/2piHz)39。)。ylabel(39。時間 39。 %legend(t,m)。 %定義頻率序列 figure(1) subplot(2,1,1)。 %傅立葉變換 M=M/fs。 %修正 t=0 點的信號值 m=m.*m。 w=t0/(2*ts)+1。 %定義頻率分辨率 x=sin(200*t)。 t=[t0/2:ts:t0/2]。 %定義時間長度 ts=。其次，我要感謝幫助過我的同學(xué)，他們也為我解決了不少我不太明白的設(shè)計上的難題。首先我要感謝我的老師在課程設(shè)計上給予我的指導(dǎo)、提供給我的支持和幫助，這是我能順利完成這次設(shè)計的主要原因，更重要的是老師幫我解決了許多技術(shù)上的難題，讓我能把系統(tǒng)做得更加完善。 通信原理課程 設(shè)計 第 21 頁 共 33 頁 致 謝 本次課設(shè)從開始到結(jié)束經(jīng)歷了三周時間，在這三周的時間里對程序進行了無數(shù)次的修改，而且還有好多次問過老師，剛開始做時一點思緒都沒有，老師講解后有了大體的操作輪廓，然后翻閱了好多資料后確定了這個方案并完成了它。掌握了對其進行信號分析的基本方法。 選擇這個課題總的來說還是稍微有點難度的，剛開始編程時一點思路都沒有，最后翻閱各種資料，方案，最后才的到了這樣的結(jié)果。 PCM 編碼采用 A 律十三折線編碼與譯碼 ，不僅簡便，且可減少量化誤差。 均勻量化輸出波形圖清晰地顯示處均勻量化的特征，每個量階都是均勻分布的，每個間隔都是相等的。 由于在基帶傳輸時誤碼率降到 以下是不難的，所以此時通常用式 (1)來估算 PCM 系統(tǒng)的性能。 （ 3） PCM 系統(tǒng)接收端輸出信號的總信噪比 ? ? NeNNeeqq ppMNN tmENS 222220241 2241)( ?????? 在接收端輸入大信噪比的情況下，誤碼率 ep 將極小，于是 124 2 ??Nep ，所以總信噪比近似為NqNS 20 2? 與只 考慮 量化噪 聲情況下的系統(tǒng)輸出信噪比是相同的。 eN ：信道加性噪聲的平均功率 。 )(tne ：由信道加性噪聲引起的輸出噪聲成分。在圖 1中的 PCM 系統(tǒng)的低通濾波器的輸出信號為 )()()()(? tntntmtm eq ??? 其 中 )(tm ：接收端輸出的信號成分 。 PCM 通信系統(tǒng)的主要參數(shù)設(shè)置 信噪比（ SNR）范圍： dB=[25:5:25]； 信源模塊：我們設(shè)定 原始信號的 幅值 A 為 ，頻率 f 為 2KHZ； 信源編碼 /譯碼：取樣時 間間隔 t： ，范圍 T 是 ；抽樣頻率 fs=8K，抽樣數(shù)量n=100，進行非均勻量化 。 通信原理課程 設(shè)計 第 18 頁 共 33 頁 4. 系統(tǒng) 性能分析 PCM 通信系統(tǒng)的性能指標 誤碼率：錯誤接收的碼元數(shù)在傳送總碼元數(shù)中所占的比例，誤碼率是碼元在傳輸系統(tǒng)中被傳錯的概率。 兩種信道編碼方式在此傳輸系統(tǒng)中，性能相當。 通信