利用Lloyd-Max，计算标准正态分布16层最佳量化.docx

题目一

 题目要求

利用 ，计算标准正态分布  层最佳量化。

 源代码

function main

clear all;

string1=input('Please input the pdf of the

signal: ');

xmin=input(' ');

disp('Please input the number of the quantization

levels:');

qsteps=input(' ');

lloyedmax(string2,xmax,xmin,qsteps)

function d=lloyedmax(pdf,xmax,xmin,qsteps)

g=inline(pdf);

step=abs((xmax-xmin))/qsteps;

for i=1:qsteps

xvalue(i)=xmin+(i-0.5)*step;

x(j)=dsum;

dold=dsum;

xdomin(1)=xmin;xdomin(qsteps+1)=xmax;

for i=2:qsteps

xdomin(i)=(xvalue(i-1)+xvalue(i))./2;

end

dsum=0;

for i=2:qsteps+1

xstemp=xdomin(i-1):abs(xdomin(i)-xdomin(i-1))/100:xdo

xstemp=xdomin(i-1):abs(xdomin(i)-xdomin(i-1))/100:xdo

min(i);

y1=xstemp.*g(xstemp);

y2=g(xstemp);

d1=trapz(xstemp,y1);

xlabel('regions (x) and outputlevels (*)')

end

outputlevel=xvalue;

regions=xdomin;

distortion=dsum;

end

j=j+1;

x(j)=dsum;

outputlevel=xvalue;

regions=xdomin;

for i=2:j

c(i-1)=x(i);

end

plot(1:1:j-1,c,'*')

xlabel('the compute numbers')

ylabel('everyone distortion')

 运行结果



 ! !"#$ 



%



%