%>>>>>>>>>>PCM脈沖編碼調制（抽樣量化編碼）>>>>>>>>>>>>>>>>>>
%其中量化和編碼采用A率的近似十三折線法來實現
%模擬信號為Xa=sin（100*pi*t）
%抽樣頻率為fs=1/T=1000Hz

%>>>>>>>>>>>>>>>模擬信號時頻波形>>>>>>>>>>>>>>>>>>>>>>>>
A=1;
w0=100*pi;
dt=0.001;
t=0:dt:pi;
Xa=A*sin（w0*t）;
subplot（221）plot（tXa）
axis（[00.02-1.51.5]）title（‘a（1）模擬信號Xa=sin（100*pi*t）時域波形‘）
grid

%模擬信號頻譜
w=0:1000;
Xaw=fft（Xalength（w））;
subplot（222）;plot（wabs（Xaw））
axis（[01000600]）title（‘a（2）模擬信號Xa=sin（100*pi*t）幅頻特性‘）
grid

%>>>>>>>>>>>>>>>>>>>抽樣信號時頻波形>>>>>>>>>>>>>>>>>>>>>>>
T=0.001;n=0:21;
X=sin（w0*n*T）;
subplot（223）stem（X‘.‘）;
axis（[021-22]）title（‘b（1）抽樣信號X=sin（w0*n*T）‘）
grid

%抽樣信號頻譜
wn=0:1000;
Xwn=fft（Xlength（wn））;
subplot（224）;plot（wnabs（Xwn））
axis（[01000012]）title（‘a（2）模擬信號Xa=sin（100*pi*t）幅頻特性‘）
grid

%>>>>>>>>>>>>>>>>>>>>量化及編碼>>>>>>>>>>>>>>>>>>>>>>>>>>
c=zeros（1176）;??%量化碼組
d0=1/2048;??????%1/2048作為一個量化單位
X0=zeros（122）;?%22個抽樣點X0

for?n=1:22
?????X0（n）=fix（X（n）/d0）;??%抽樣值的量化單位數（取整）
?????%>>>>>>>>>>>>>>極性碼>>>>>>>>>>
????????if?X0（n）>0
????????????c（1+（n-1）*8）=1;
????????else
????????????c（1+（n-1）*8）=0;
????????end
????d=1024;
????q=2;???
????I=0;????%段落碼碼段
????J=0;????%段內碼碼段
????for?i=1:8???????????????%求段落碼碼段I
????????max=d*q;
????????min=d/（q^4）;
????????if?abs（X0（n））<=d*q?&&?abs（X0（n））>d
????????????I=8-i+1;
????????????for?j=1:16??????%求段內碼碼段J
????????????????if?abs（X0（n））<=max?&&?abs（X0（n））>max-min
?????????????????J=16-j+1;
????????????????end
????????????????max=max-min;
????????????end
????????end
????????d=d/q;
????end

????switch?I????????????%由段落碼碼段求段落碼
????????case?8
????????????c（2+（n-1）*8:4+（n-1）*8）=[111];
????????case?7
????????????c（2+（n-1）*8:4+（n-1）*8）=[110];
????????cas