matlab--PCM编码.doc

PCM编码:clear all; close all; %建立原信号dt=0.002; %取时间间隔为0.01t=0:dt:10; %时域间隔dt为间隔从0到10画图fc=1 %xt里最大频率xt=sin(2*pi*fc*t)+cos(2*pi*fc*t); %xt方程 %采样：时间连续信号 变为 时间离散模拟信号fs=10; %抽样fs=2fc，每秒钟内的抽样点数目将等于或大于2fc个sdt=1/fs; %频域采样间隔0.1t1=0:sdt:10; %以sdt为间隔从0到10画图st=sin(2*pi*fc*t1)+cos(2*pi*fc*t1); %coswt=cos2pift,2pif=wfigure(1); subplot(311);plot(t,xt);title(原始信号); %条状图，连续图grid on %画背景subplot(312);stem(t1,st,.); %杆状图 ，离散图 title(抽样信号);grid on %画背景%量化编码一步完成n=length(st); %取st的长度为n，本题函数n为101 M=max(st); A=(st/M)*2048;%a1(极性码) a2a3a4（段落码）a5a6a7a8（段内电平码） code=zeros(i,8); %产生i*8的零矩阵%极性码a1 for i=1:n %if循环语句 if A(i)=0 code(i,1)=1; %代表正值 else code(i,1)=0; %代表负值 end % 段内码a2a3a4 if abs(A(i)=0&abs(A(i)16 code(i,2)=0;code(i,3)=0;code(i,4)=0;step=1;start=0;elseif 16=abs(A(i)&abs(A(i)32 code(i,2)=0;code(i,3)=0;code(i,4)=1;step=1;start=16;elseif 32=abs(A(i)&abs(A(i)64 code(i,2)=0;code(i,3)=1;code(i,4)=0;step=2;start=32;elseif 64=abs(A(i)&abs(A(i)128 code(i,2)=0;code(i,3)=1;code(i,4)=1;step=4;start=64;elseif 128=abs(A(i)&abs(A(i)256 code(i,2)=1;code(i,3)=0;code(i,4)=0;step=8;start=128;elseif 256=abs(A(i)&abs(A(i)512 code(i,2)=1;code(i,3)=0;code(i,4)=1;step=16;start=256;elseif 512=abs(A(i)&abs(A(i)1024 code(i,2)=1;code(i,3)=1;code(i,4)=0;step=32;start=512;elseif 1024=abs(A(i)&abs(A(i)=0 code(i,1)=1; else code(i,1)=0; end if abs(A(i)=0&abs(A(i)16 code(i,2)=0;code(i,3)=0;code(i,4)=0;step=1;start=0;elseif 16=abs(A(i)&abs(A(i)32 code(i,2)=0;code(i,3)=0;code(i,4)=1;step=1;start=16;elseif 32=abs(A(i)&abs(A(i)64 code(i,2)=0;code(i,3)=1;code(i,4)=0;step=2;start=32;elseif 64=abs(A(i)&abs(A(i)128 code(i,2)=0;code(i,3)=1;code(i,4)=1;step=4;start=64;elseif 128=abs(A(i)&abs(A(i)256 code(i,2)=1;code(i,3)=0;code(i,4)=0;step=8;start=128;elseif 256=abs(A(i)&abs(A(i)512 code(i,2)=1;code(i,3)=0;code(i,4)=1;step=16;start=256;elseif 512=abs(A(i)&abs(A(i)1024 code(i,2)=1;code(i,3)=1;code(i,4)=0;step=32;start=512;elseif 1024=abs(A(i)&abs(A(i)2048 code(i,2)=1;code(i,3)=1;code(i,4)=1;step=64;start=1024; endB=floor(abs(A(i)-start)/step); t=dec2bin(B,4)-48; code(i,