close all
clear all

mu=255;
m=randn(1,500);
n=10;%量化級數的和數目
[f,maximum]=mulaw(m,mu);%產生歸一化的壓縮曲線f
[sqnr,y_q,code]=u_pcm(f,n);%對歸一化的壓縮曲線進行均勻量化
m_quan=invmulaw(y_q,mu);%恢復量化后的歸一化序列
m_quan=maximum*m_quan;%恢復量化后的序列
sqnr=20*log10(norm(m)/norm(m-m_quan))%求量化信噪比