正弦波信源编码、译码，bch信道编码.m

以正弦波信号为例

%%%%%%%%%对正弦序列进行BPCM编码，即信源编码
predictor=[0,1];%设初值
partition=[-1:0.1:0.9];
codebook=[-1:0.1:1];
t=[0:pi/50:2*pi-pi/50];

x=2*sin(3*t);%取3个周期，幅度为2的正弦曲线
xinyuancode=dpcmenco(x,codebook,partition,predictor);%十进制信源编码
decodedx=dpcmdeco(xinyuancode,codebook,predictor);%十进制信源解码
%            pause             %按任意键可以查看下图
subplot(2,1,1);
plot(x-decodedx);%原信号与信源译码的差
title('原信号与信源译码的差')
%           pause              %按回车键可以查看下图
subplot(2,1,2);
plot(xinyuancode);%编码序列
plot(t,x,t,xinyuancode,'--',t,decodedx,'r');%原码与译码序列
title('信源编码、原信号、信源译码序列')
%%%%%%%%%%%%%%%%%%%%%%对信源编码使用BCH码进行信道编码
sizexinyuancode=size(dec2bin(xinyuancode));%求出信源编码的二进制位数
params=bchpoly(sizexinyuancode(1,2));%计算bch编码有效长度
n=params(1,1);
k=params(1,2);
xindaocode=encode(xinyuancode,n,k,'bch/decimal');%%%十进制编码