【導讀】脈沖編碼調制是將模擬信號變換成二進制信號的常用方法。于20世紀40年代，在。通信技術中就已經實現了這種編碼技術。增量調制可以看成是一種最簡單的DPCM，當DPCM. 系統(tǒng)中量化器的量化電平數取為2時，此DPCM系統(tǒng)就成為增量調制系統(tǒng)。是一個微分脈沖碼調制的簡化形式，它具有結構簡單，高效率的編碼，誤碼性能特點。本論文結合PCM與ΔM的抽樣、量化、編碼原理，利用MATLAB軟件編程和繪圖功能，完成了對脈沖編碼調制系統(tǒng)與增量調制系統(tǒng)（ΔM）的建模與仿真分析?；cA律13折線非均勻量化的量化性能及其差異。驗證了數字傳輸技術的優(yōu)越性。

　　

【正文】 27 title(39。抽樣正常 (fs=200Hz)時的信號波形 39。)。axis([,0,1])。 subplot(2,1,2) plot(f1,abs(fftshift(N1)))。 xlabel(39。頻率 39。)。ylabel(39。幅值 39。)。axis([500,500,0,])。 title(39。抽樣正常時的信號頻譜 39。)。axis([500,500,])。 t0=10。 %信號持續(xù)的時間 ts2=。 %不滿足抽樣條件的抽樣間隔 fs2=1/ts2。 t2=[t0/2:ts2:t0/2]。 %定義不滿足抽樣條件的時間序列 x2=sin(200*t2)。 m2=x2./(200*t2+eps)。 w2=t0/(2*ts2)+1。 m2(w2)=1。 %修正 t=0 時的信號值 m2=m2.*m2。 %定義信號 [M2,mn2,df2]=fft_seq(m2,ts2,df)。%對不滿足抽樣條件的信號進行傅立葉變換 M2=M2/fs2。N2=[M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2,M2]。 f2=[7*df2*length(mn2):df2:6*df2*length(mn2)df2]fs2/2。 figure(3) subplot(2,1,1)。 stem(t2,m2)。 xlabel(39。時間 39。)。ylabel(39。幅值 39。)。title(39。抽樣失真 (fs=100Hz)時的信號波形 39。)。 axis([,0,1])。subplot(2,1,2) plot(f2,abs(fftshift(N2)))。 xlabel(39。頻率 39。)。ylabel(39。幅值 39。)。axis([500,500,0,])。 title(39。抽樣失真時的信號頻譜 39。)。axis([500,500,])。 function [M,m,df]=fft_seq(m,ts,df) fs=1/ts。 if nargin==2 n1=0 else n1=fs/df end n2=length(m)。n=2^(max(nextpow2(n1),nextpow2(n2)))。 M=fft(m,n)。m=[m,zeros(1,nn2)]。df=fs/n 程序 2： function average() t=[0::4*pi]。 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 28 y=sin(t)。 w=jylh(y,1,64)。 subplot(2,1,1)。 plot(t,y)。 xlabel(39。時間 39。)。 ylabel(39。幅度 39。)。 axis([0,4*pi,])。 title(39。原始信號 39。)。 subplot(2,1,2)。 plot(t,w)。 xlabel(39。時間 39。)。 ylabel(39。幅度 39。)。 axis([0,4*pi,])。 title(39。均勻量化后的信號 39。)。 function h=jylh(f,V,L) n=length(f)。t=2*V/L。 p=zeros(1,L+1)。 for i=1:L+1,p(i)=V+(i1)*t。end for i=1:n if f(i)V,h(i)=V。end if f(i)=V,h(i)=V。end flag=0。 for j=2:L/2+1 if(flag==0) if(f(i)p(j)) h(i)=p(j1)。 flag=1。 end。 end。 end。 for j=L/2+2:L+1 if(flag==0) if(f(i)p(j)) h(i)=p(j)。 flag=1。 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 29 end end end end nq=V^2/(3*L^2)。 程序 3： function a_quantize() t=0::。 y=sin(8000*pi*t)。 figure subplot(2,1,1) plot(t,y) axis([0 ]) xlabel(39。時間 39。) ylabel(39。幅度 39。) title(39。原始信號 39。) z=a_pcm(y,)。 subplot(2,1,2) plot(t,z) axis([0 ]) xlabel(39。時間 39。) ylabel(39。幅度 39。) title(39。A律量化后的信號 39。) function y=a_pcm(x,a) t=1/a。 for i=1:length(x) if x(i)=0 if(x(i)=t) y(i)=(a*x(i))/(1+log(a))。 else y(i)=(1+log(a*x(i)))/(1+log(a))。 end else if(x(i)=t) 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 30 y(i)=(a*x(i))/(1+log(a))。 else y(i)=(1+log(a*x(i)))/(1+log(a))。 end end end 程序 4： function a_13code() t=0::。 y=sin(8000*pi*t) z=line13(y) c=pcmcode(z) function y=line13(x) x=x/max(x)。 z=sign(x)。 x=abs(x)。 for i=1:length(x) if((x(i)=0)amp。(x(i)1/64)) y(i)=16*x(i)。 else if((x(i)=1/64)amp。(x(i)1/32)) y(i)=8*x(i)+1/8。 else if((x(i)=1/32)amp。(x(i)1/16)) y(i)=4*x(i)+2/8。 else if((x(i)=1/16)amp。(x(i)1/8)) y(i)=2*x(i)+3/8。 else if((x(i)=1/8)amp。(x(i)1/4)) y(i)=x(i)+4/8。 else if((x(i)=1/4)amp。(x(i)1/2)) y(i)=1/2*x(i)+5/8。 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 31 else if((x(i)=1/2)amp。(x(i)=1)) y(i)=1/4*x(i)+6/8。 end end end end end end end end y=z.*y。 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 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 32 end 程序運行結果 ： a_13code y = Columns 1 through 7 0 Columns 8 through 11 z = Columns 1 through 7 0 Columns 8 through 11 c = 1 0 0 0 0 0 0 0 0 1 1 1 0 0 1 1 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 0 1 1 1 0 0 1 1 0 0 0 0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 1 1 1 0 0 0 0 0 0 0 程序 5： Ts=1e3。 t=0:Ts:20*Ts。 x=sin(2*pi*50*t)+*sin(2*pi*150*t)。 山東輕工業(yè)學院 2020 屆本科生畢業(yè)設計（論文） 33 delta=。 D(1+length(t))=0。 for k=1:length(t) e(k)=x(k)D(k)。 e_q(k)=delta*(2*(e(k)=0)1)。 D(k+1)=e_q(k)+D(k)。 codeout(k)=(e_q(k)0)。 end subplot(3,1,1)。plot(t,x,39。o39。)。axis([0 20*Ts,2 2])。hold on。 subplot(3,1,2)。stairs(t,codeout)。axis([0 20*Ts,2 2])。 Dr(1+length(t))=0。 for k=1:length(t) eq(k)=delta*(2*codeout(k)1)。 xr(k)=eq(k)+Dr(k)。 Dr(k+1)=xr(k)。 end subplot(3,1,3)。stairs(t,xr)。hold on。 subplot(3,1,3)。plot(t,x)。