t=[0:0.01:10];
a=sin(t);
k1=input('輸入第一個量化電平數:');
[aquana,codea]=u_pcm(a,k1);
k2=input('輸入第二個量化點平數:');
[aquanb,codeb]=u_pcm(a,k2);
plot(t,a,'-',t,aquana,' -.',t,aquanb,'-',t,zeros(1,length(t)))