close all
clear all
echo off
i=1:500;
a=randn(1,500);%
mu=255;
[sqnr,a_quan,code]=mula_pcm(a,8,mu);%
f=a-a_quan;                 %
subplot(3,1,1);
%axis([0 500 -15 15]);
plot(i,f,'-');
title('非均勻PCM量化誤差,從上到下量化級數分別為8,16,64,')

[sqnr,a_quan,code]=mula_pcm(a,16,mu);%
f=a-a_quan;                 %
subplot(3,1,2);
%axis([0 500 -15 15]);
plot(i,f,'-');


[sqnr,a_quan,code]=mula_pcm(a,64,mu);%
f=a-a_quan;                 %
subplot(3,1,3);
%axis([0 500 -15 15]);
plot(i,f,'-');