均匀与非均匀量化

1、实验二均匀、非均匀量化与编码实验目的1. 了解均匀与非均匀董化编码的原理与过程2. 掌握两种不同量化编码方式的特点并能加以比较3. 进一步熟悉mat I ab在通信方面的应用理论依据通原理论课中的标量董化原理，包括均匀量化和以PCM为例的非均匀量化（对数量化）。实验内容IS仿真理论及系统框图四. 实验数据1、均匀量化(曲线名称见各自曲线图)(1)n=8（2）非均匀量化发送端： 接收端非均匀董化（对数量化）原理框图U律压扩特性（本实验以U律为例）2、模块说明实脸程序设计依照量化编码原理可大体分为采样.量化.编码、作图及相关显示几部分。采样：本实验利用正弦函数的若千点（17个）。量化：根摇不同的输

2、入值.可以得到在不同量化级数下（量化级数n=8, 16, 64）的结 果，算法是程序的核心之一。编码：对董化结果进行二进制编码，本例进行自然编码，另一个算法核心。作图与显示：做出不同量化级数下量化曲线、量化误差曲线：显示量化信噪比、及部分 釆样点的编码结果。3、语句分析详细的分析见程序注释3468质姑愴号及月化由线<r» 1 11 1. 1 DSf化佑呢t1 0!> 18 42B30 15U"1J 0 0<. JI.I.Q 八Da06D.4D200609e10(2) n=16K<i及旨化出线9UnIu D I CIfDOOi>LLOdIIOC&

3、gt;ezzc=:1¾)S>I O I <>ILL，> I I I ()IlLl>oOOOOQV 9 br- rw D Cr 9 8Ktei及£化由线虽化值及城码显云<.IMIloIIIII<1110.<10 11)<>37 1207 >. d ( C U I I U?U1 U3 b I 1 1 t.Jhl-.八2、非均匀量化(U率对数量化)(1)n=8Bffei 号及fts-e琴W虧*tfgte及i码显示e盘化 1S¾ti=7.1<1 1 1<1 1 1<1 1 1(1 1 1I

4、<1 1 1<1 1 141 1 1 >D 1 11 D 0.(0 0»6E U n5O236(2)n=16s.½ 号及 gffi238曆化连盖曲妙B5SIW虽化這及城码显云(3)6 0-O . <1111<1111(1111<1111<1111 虧化恰呢1：.t. 1. .l.4.1.<I(>DlIl10 0 0弋(0 0八n=6423昱化误盖曲线量化值及i码显禾050.42 IOJ 2J 4 5 D O OODOO 马禦tw.<.1.IJIIIIl.QIIlll J<IlIII.rI1 Il-LlO<

5、;<S1Sbt-24 25B7I1110!0IllOIUCr>DlIlI1IOOOO01 a丨一"6 8-D-O5B9IQnum2、数据分析与说明：(1) n=8, 16, 64的均匀董化曲线对于均匀量化，董化级数越大，量化值与原始值约接近，董化误差越小，即量化效 果越好。(2) n=8, 16, 64的均匀董化信噪比对于均匀量化，量化级数越大，量化信噪比越大(3) n= 8, 16, 64的非均匀量化曲线对于非均匀量化，量化级数越大，量化值与原始值约接近，董化误差越小，即量化 效果越好。(4) n=8, 16, 64的非均匀量化信噪比对于非均匀量化，量化级数越大，量化信

6、噪比越大(5) 对于相同董化级的均匀与非均匀量化特性均匀量化的量化误差分布比较均匀，而非均匀量化的量化误差随信号幅度变化: 在大幅度处大，在小幅度处小。当n较小时，均匀量化的董化信噪比大于非均 匀量化，但是当n较大时，均匀量化的董化信噪比小于非均匀量化五、实验结果的分析与讨论通过对实验数据的分析，我们直观的发现了与通原理论课一样的结论。无论是均匀量化 还是非均匀量化其董化信噪比的是随着量化级数的增大而增大。而均匀量化与非均匀董 化相比，前者是先简单，后者相对复杂，但更适合小信号出现概率校大的信号，有实际 应用背景。六、实验总结与心得体会这次实验在利用老师给出的核心部分后相对简单，再加之对mat

7、Iab的掌握有了一定提 鬲完成得比上次顺利。其中对matIab编程的规范性有了进一步的强化，无论是注释 的规范性还是程序的灵活性都有了提鬲。七、思考题1、答：釆用奈圣斯特抽样频率为2fh,量化级数为n,则每码元比特数为Iog2,董化后的信号带宽为2wlog2n。2、答：非均匀PCM的量化信噪比之所以优于均匀PCM的董化信噪比是因为它具有较大的动态范国，结合实际语声信号的特点：大信号出现概率小.小信号出现槪率大，非均匀量化的压扩特性恰好适应这一条件。3、答：以n=8为例通过程序算出：祈者a 后者a均匀量化量化信噪比为非均匀量化量化信噪比为均匀量化量化信噪比为非均匀童化量化信噪比为虽然以上数据均为

8、不定值，但大小相对关系不会改变： 即：前者a均匀董化量化信噪比优于非均匀董化量化信噪比后者a非均匀董化量化信噪为优于均匀董化量化信噪比原因是前者a为服从正态分布，样值比较均匀，适合均匀量化而后者a的各个元素跳跃较大，分布不均，适合动态范国较大的非均匀量化附录：源程序：仁均匀量化N=2105L=26; M 二 N/L;%釆样点数a=sin(0:2*pi/M：2*pi);% 取样点n =input(,化级数=');%初始化，n为量化级数 if iSemPty(n),n=8;endamax=max (abs (a);a_quan=a/amax; % 归化 b_quan=a_quan;d=2/

9、n;%#化间隔q=d.*0:n-lJq=q- (n-1 )/2) *d; % 化电平%进行量化for i=1:na_quan(f ind(q(i)-d/2<=a_quan) & (aquan<=q(i) +d2) =.定位第 i 个董化 间隔码子q(i). *ones(1, Iength(f ind(q(i)-d2<=aquan)&(a_quan<=q(i) +d2) ;%赋值 为相应的量化电平b_quan (fi nd (a_quan=q (i) =(i-1)* OneS (1, Iength (f ind (a_qUan=q (i); enda_qu

10、an=a_quan*amax; %恢复原值(量化后)nu=ce i I (I og2 (n) ;%编码COde=ZerOS (Iength(a), nu);for i=1:Iength(a)for j=nu: T : 0%从离向低编码if (f iX (b_qUan (i)/(2 j) =1)COde(i, (nu-j)=1;b_quan (i)=b_quan(i)-2j;endendendSqnr=20* I Ogl 0 (norm (a) /norm (a-a_quan);%求量化信噪比t=(0:2*pi/N:2*pi);aq=sin(0:2*pi/N:2*pi)计算所有釆样点的量化值，为

11、作图服务 aq_quan=aq/amax; % 归一化aq=d.*0:n-1Jaq= (-1) /2) *daq j%,化电平for i=1:naq_quan (f i nd (q (i) -d/2<=aq_quan) & (aq-quan<=q (i) ÷d2) =. . %定位第 i 个 量化间隔码子q(i). *ones(1, Iength(f ind(q(i)-d/2<=aq_quan)&(aq_quan<=q(i)+d2) ;%賦 值为相应的量化电平endfigured)%作原始信号及量化曲线Set (1, ,position1, 1

12、0, 50, 300, 200)PIOt (t, s i n (t)gr id Onax i s (0,2*p i, *amax, *ama×)title('原始信号及量化曲线')Xlabel Ct')ylabel C原始信号及量化信号值text (2, 1, 'sin(t)')ho Id OnPlOt(匕 aq_quan,r,)figure (2) %做量化误差曲线(对所有采样点)Set(2, - position1, 350,50,300,200)PI Ot (t, aq_qUan-S i n (t)gr id Onax i s(0, 2*

13、p i, -2*amaxn, 2*amaxn)title('量化误差曲线')Xlabel ('t')ylabel C量化误差值')figure (3) %做前十个取样点的量化值曲线及表示量化信噪比(量化信噪比的计算依据一个周期内所有取样点，并非前十个)Set (IJPOSitiOr, 20, 60, 300, 200)Stem (a_quan)gr id Onax i s(0, -amax, amax)title(,量化值及编码显示JXlabel(,num,)ylabel ('量化值及编码显示')for k=1:MCODA=num2str

14、(COde(k,:);text (k, a_quan (k), CODA)text*amax,'量化信噪比=')text, *amax, num2str(Sqnr)end2、非均匀量化N=210jL=26; M 二 N/L;%采样点数aa=sin(0:2*pi/M:2*pi) 取样点 u=255;n=input(,化级数=');%初始化,n为量化级数if iSemPty(n),n=8;enda=mulaw (aa, U);amax=max (abs (a);a_quan=a/amax; % 归一化 b_quan=a_quan;d=2nj%化间隔q=d.*0:n-1Jq=

15、q- (nT) /2) *d; % 化电平%进行量化for i=1:na_quan (f i nd (q (i) -d/2<=a_quan) & (a_quan<=q (i) +d2) = . 定位第 i 个量化 间隔码子q(i). *ones (1, Iength (f i nd (q (i) -d/2<=a_quan) & (a_quan<=q (i) +d2) ; % 赋值 为相应的量化电平b_quan (f ind (a_quan=q (i) =(i-1) . * OneS (1, Iength (f ind (a_qUan=q (i); end

16、a_quan=a_quan*amax; %恢复原值(董化后)nu=ce i I (I og2 (n) 编码COde=ZerOS(Iength(a), nu);for i=1:Iength(a)for j=nu ：-1:0%从离向低编码if (f iX (b_quan (i)/(2 j) =1)COde(i, (nu-j)=1;b_quan (i)=b_quan(i)-2j;endendendaa_qUan=i nVmU I aw (a_qUan, U);Sqnr=20* I Ogl 0 (norm (aa)norm (aaaa_quan) ;%求董化信噪比t=(0:2*pi/N:2*pi);a

17、aq=s in (0:2*pi/N:2*pi) ;%计算所有采样点的量化值，为作图服务 aq=mulav(aaq, U);aq_qUan=aq amax; % 归一化aq=d.*0:n-1Jaq= (-1) /2) *d-aq j%,化电平for i=1:naq_quan (f i nd (q (i) -d/2<=aq_quan) & (aq-quan<=q (i) ÷d2) =. . %定位第 i 个 量化间隔码子q(i). *ones(1, Iength(f ind(q(i)-d/2<=aq_quan)&(aq_quan<=q(i)+d2)

18、;%赋值为相应的量化电平endaaq_quan=i nvmuI aw (aq_quan,U);figured)%作原始信号及量化曲线Set (1, -position1, 10, 50, 300, 200)PIOt (t, s i n (t)gr id Onax i s(0, 2*pi, *amax, *ama×)title('原始信号及量化曲线)Xlabel ('t')ylabel C原始信号及量化信号值')text (2, 1, ,sin(t)')ho Id OnPIOt (t, aaq_quan, , r,)figure (2) %做董化误差曲线(对所有釆样点)Set (2, -position-, 350,50,300,200)PIot(t, aaq_quansin(t)gr id Onax i s (0,2*p i,*amax, *amax)title('量化误差曲线)Xlabel Ct')ylabel C量化误差值')figure (3) %做前十个取样点的量化值曲线及表示