【正文】 r(k)39。039。 %變?yōu)殡p極性 if SigCode_JiDai_In(j+k)==0 SigCode_JiDai_In(j+k)=0。 end k=k+1。 end j=j+8。 i=i+1。end%PCM抽樣周期Ts=8。%基帶信號(hào)碼元頻率 FdFd=8/Ts。%基帶時(shí)間軸t1=[0:(length(SigCode_JiDai_In)1)]./Fd。%濾波器采樣頻率 Fs_fFs_f=20*Fd。%滾降系數(shù) alpha、alpha=1。%濾波器群延遲 group_delay=3。%PCM編碼信號(hào)波形[Wave,tw]=RectangWaveofSqu3(SigCode_JiDai_In,Ts/8,1,1,Fs_f)。figure(2)plot(tw,Wave)grid onxlabel(39。時(shí)間(s)39。)。ylabel(39。幅度39。)。 title([39。PCM編碼信號(hào)波形39。])。axis([,])。Amax_D=Amax。%基帶信號(hào) SigCode_JiDai 變?yōu)?信號(hào)的PCM碼組 CodeofPCM_Out%串并轉(zhuǎn)換CodeofPCM_Out=cell(1,length(SigCode)/8)。j=0。for i=1:length(CodeofPCM_Out) temp=39。39。 for k=1:8 temp=strcat(temp,num2str(SigCode(j+k)))。 k=k+1。 end j=j+8。 CodeofPCM_Out(i)=cellstr(temp)。 i=i+1。end%信號(hào)譯碼 Sig_DSig_D=zeros(length(CodeofPCM_Out),1)。for i=1:length(Sig_D) Sig_D(i)=D_CodeofPCM(char(CodeofPCM_Out(i)))*Amax_D。 i=i+1。end%原始語音信號(hào)波形figure(1)subplot(2,1,1)plot(t,Sig_O), grid ontitle(39。原始語音信號(hào)波形39。)axis([1,1+200/Fs,1,1])xlabel(39。Time/S39。)。 ylabel(39。Amplitude/V39。)。zoom xon%語音信號(hào)譯碼輸出波形subplot(2,1,2)plot(t,Sig_D) ,grid ontitle(39。語音信號(hào)譯碼輸出波形39。)xlabel(39。Time/S39。)。 ylabel(39。Amplitude/V39。)。axis([1,1+200/Fs,1,1])zoom xon%語音信號(hào)編譯碼輸出波形對(duì)比figure(6)plot(t,Sig_O,39。r39。,t,Sig_D,39。b39。)grid onlegend(39。原始語音信號(hào)波形39。,39。譯碼輸出波形39。)title(39。語音信號(hào)編譯碼輸出波形對(duì)比39。)xlabel(39。Time/S39。)。 ylabel(39。Amplitude/V39。)。axis([1,1+200/Fs,1,1])zoom xonsound(Sig_D,Fs)。%% 遇到的問題d=0for i=1:length(Sig_D) d=d+(Sig_O(i)Sig_D(i))^2。 i=i+1。end%d_Sum_N=d。d_Sum=d。量化編碼C_CodeofPCMfunction CodeofPCM=C_CodeofPCM(A_Crent,Amax)% A率13折線編碼 %CodeofPCM=DianpingtoCodeofPCM(A_Crent,Amax)%入口參數(shù)%信號(hào)最大幅度 Amax （數(shù)值）%信號(hào)幅度 A_Crent （數(shù)值）%出口參數(shù)%信號(hào)的PCM碼組 CodeofPCM （字符串）%函數(shù)示例 CodeofPCM=DianpingtoCodeofPCM(635,2048)%返回 CodeofPCM 為 11100011（字符串）。%% 極性 JixingJixing=1。if A_Crent0 Jixing=0。endJixing=dec2bin(Jixing,1)。%% 段號(hào) DuanHaoDuanHao=0。%信號(hào)幅度歸一 A_Crentmaxto1 A_Crentmaxto1在 0~2048之間A_Crentmaxto1=floor(abs(A_Crent)/Amax*2048)。%段落起始 DuanluoQS%量化間隔 LinghuaJG LinghuaJG=0。%分段操作if A_Crentmaxto116 DuanHao=0。 DuanluoQS=0。 LinghuaJG=1。elseif A_Crentmaxto132 DuanHao=1。 DuanluoQS=16。 LinghuaJG=1。elseif A_Crentmaxto164 DuanHao=2。 DuanluoQS=32。 LinghuaJG=2。elseif A_Crentmaxto1128 DuanHao=3。 DuanluoQS=64。 LinghuaJG=4。elseif A_Crentmaxto1256 DuanHao=4。 DuanluoQS=128。 LinghuaJG=8。elseif A_Crentmaxto1512 DuanHao=5。 DuanluoQS=256。 LinghuaJG=16。elseif A_Crentmaxto11024 DuanHao=6。 DuanluoQS=512。 LinghuaJG=32。else DuanHao=7。 DuanluoQS=1024。 LinghuaJG=64。endDuanHao=dec2bin(DuanHao,3)。 %% 段內(nèi)碼 （自然碼） DuanNeima%異常處理if A_Crentmaxto1==2048 A_Crentmaxto1=A_Crentmaxto11。end DuanNeima1=floor((A_Crentmaxto1DuanluoQS)/LinghuaJG)。 DuanNeima=dec2bin(DuanNeima1,4)。 %% PCM碼組 CodeofPCM CodeofPCM=strcat(Jixing,DuanHao, DuanNeima)。end ——————————————————————————————————譯碼處理程序：function A_Crent=D_CodeofPCM(CodeofPCM)% A率13折線譯碼 % A_Crent=D_CodeofPCM(CodeofPCM,Amax)%入口參數(shù)%信號(hào)最大幅度 Amax （數(shù)值）%信號(hào)的PCM碼組 CodeofPCM （字符串）%出口參數(shù)%信號(hào)幅度 A_Crent （數(shù)值）%函數(shù)示例 A_Crent=D_CodeofPCM(39。1111001139。)%返回 A_Crent 為 1216%%字符串轉(zhuǎn)換為十進(jìn)制數(shù)處理CodeofPCM_int=str2num(CodeofPCM)。%% 極性判斷 JixingJixing=1。if floor(CodeofPCM_int/10^7)==0 Jixing=1。end%% 段號(hào)提取 DuanHao DuanHao=floor(mod(abs(CodeofPCM_int),10^7)/10^4)。DuanHao=num2str(DuanHao)。DuanHao=bin2dec(DuanHao)。%段起始電平提取 DuanluoQS%量化間隔提取 LinghuaJGswitch DuanHao case 0 DuanluoQS=0。 LinghuaJG=1。 case 1 DuanluoQS=16。 LinghuaJG=1。 case 2 DuanluoQS=32。 LinghuaJG=2。 case 3 DuanluoQS=64。 LinghuaJG=4。 case 4 DuanluoQS=128。 LinghuaJG=8。 case 5 DuanluoQS=256。 LinghuaJG=16。 case 6 DuanluoQS=512。 LinghuaJG=32。 case 7 DuanluoQS=1024。 LinghuaJG=64。end%% 段內(nèi)碼提取 （自然碼） DuanNeimaDuanNeima=mod(abs(CodeofPCM_int),10^4)。DuanNeima=num2str(DuanNeima)。DuanNeima=bin2dec(DuanNeima)。%信號(hào)幅度絕對(duì)值 absofA_Crent %absofA_Crent=DuanluoQS+LinghuaJG*DuanNeima。%信號(hào)幅度絕對(duì)值 absofA_Crent (7/12變換原理)absofA_Crent=DuanluoQS+LinghuaJG*DuanNeima+*LinghuaJG。%%信號(hào)幅度 A_CrentA_Crent=absofA_Crent*Jixing/2048。end