通信原理课程设计DPCM编译码器设计及应用

(k)=eq(k)+Dr(k);РDr(k+1)=xr(k); %延迟器状态更新РendРsubplot(4,1,3);stairs(t,xr);axis([0 20*Ts,-2 2]);hold on; %解码输出Рsubplot(4,1,4);plot(t,x); %原信号Р4.数字信号编码解码设计РTs=0.002; %采样间隔Рt=0:Ts:20*Ts; %仿真时间序列Рx=[];Рfor i=1:21Р if rand()>0.5Р k=1;Р elseР k=-0;Р endР x=[x k];РendРdelta=0.4; %量化阶距РD(1+length(t))=0; %预测器初始状态Рfor k=1:length(t)Р e(k)=x(k)-D(k);Р e_q(k)=delta*(2*(e(k)>=0)-1);%量化器输出Р D(k+1)=e_q(k)+D(k); %延迟器状态更新Р codeout(k)=(e_q(k)>0); %编码输出РendРsubplot(3,1,1);plot(t,x,'-o');axis([0 20*Ts,0 2]);hold on;Рsubplot(3,1,2);stairs(t,codeout);axis([0 20*Ts,0 2]); %解码端РDr(1+length(t))=0; %解码端预测器初始状态Рfor k=1:length(t)Рeq(k)=delta*(2*codeout(k)-1); %解码Р xr(k)=eq(k)+Dr(k);РDr(k+1)=xr(k); %延迟器状态更新РendРsubplot(3,1,3);stairs(t,xr);hold on; %解码输出Рsubplot(3,1,3);plot(t,x); %原信号Р四.设计结果及其分析Р1.原始声音信号波形和频谱Р图2 模拟原始信号频谱