%%?本函數(shù)實現(xiàn)國際通用的PCM量化A律13特性近似，x為輸入的序列A率自定義值，M為最大量化單位，yn為量化值，Cn為編碼值

function?[CnynQn]?=?APCM（xMmodel）
%%?輸入信號量化前歸一化處理
xt?=?x/max（abs（x））;?

s?=?sign（xt）;??%?符號函數(shù)，求得每一序列值的符號

Xt?=?abs（xt）;??%?求歸一化后序列的絕對值

Q?=M*Xt;?%?確定量化單位

Qn=zeros（length（x）2）;%Qn（:1）段數(shù)，Qn（:2）段內數(shù)

%%?均勻PCM編碼模式
a=M*[01/1281/641/321/161/81/41/21];
%段內最小量化間隔
for?ii?=?1:8;??an（ii）=?（a（ii+1）-a（ii））/16?;end

%%?選擇為PCM編碼或者對數(shù)PCM碼組的線性PCM的13位碼組
if?model==1
????
?????ypcm?=?zeros（length（x）8）;?
?????
?????for?i?=1:length（Q）
?????????
?????????for?j?=?1:8
?????????????
?????????????if?Q（i）>=a（j）?&&?Q（i）?????????????????
?????????????????Qn（i1）=j-1;?%?段數(shù)確定
?????????????????
?????????????????for?k?=0:15
?????????????????????
?????????????????????if?Q（i）>=a（j）+k*an（j）?&&?Q（i）??????????????????????????
?????????????????????????Qn（i2）?=?k;?%?段內數(shù)確定
?????????????????????????
??????????????????????????yn（i）?=?（（a（j）+（k+0.5）*an（j））*max（abs（x））*s（i））/M;
?????????????????????????
?????????????????????end
?????????????????end
?????????????end
?????????end
?????????
?????????if?Q（i）==a（9）
?????????????
?????????????Qn（i:）?=?[715]?;
?????????????
?????????????yn（i）?=?max（abs（x））*s（i）;
?????????????
?????????end
?????????
?????????%符號位的判斷
?????????if?s（i）>0
?????????????
?????????????ypcm（i1）=1;
?????????????
?????????elseif?s（i）<0
?????????????
?????????????ypcm（i1）=0;
?????????????
?????????end??
?????????ypcm（i2:4）?=?（numto2（Qn（i1）3））‘;
?????????
?????????ypcm（i5:8）?=??（numto2（Qn（i2）4））‘;
?????end
end???????

Cn?=?ypcm;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-07-04?09:55??音樂信號編碼matlab20180702\
?????文件????????1761??2018-07-03?21:53??音樂信號編碼matlab20180702\APCM.m
?????文件??????768044??2018-07-04?09:52??音樂信號編碼matlab20180702\Dquit1.wav
?????文件??????768044??2018-07-04?09:52??音樂信號編碼matlab20180702\Dquit2.wav
?????文件??????264155??2018-07-04?09:54??音樂信號編碼matlab20180702\SL_data.mat
?????文件?????3840044??2018-07-03?22:15??音樂信號編碼matlab20180702\SL_data.wav
?????文件?????????206??2018-07-02?22:16??音樂信號編碼matlab20180702\calsample.m
?????文件????????1098??2018-07-04?09:52??音樂信號編碼matlab20180702\cy_read01.m
?????文件?????????896??2018-07-03?19:01??音樂信號編碼matlab20180702\dec2binPN.m
?????文件?????????354??2018-07-03?21:42??音樂信號編碼matlab20180702\numto2.m
?????文件?????????275??2018-07-02?20:02??音樂信號編碼matlab20180702\wavread.m
?????文件??????????36??2018-07-02?22:51??音樂信號編碼matlab20180702\函數(shù)說明.txt
?????文件?????1536044??2018-07-04?09:51??音樂信號編碼matlab20180702\煙花易冷8s.wav