pcm抽样量化编码程序

clear all; close all; %pcmcoding%设置输入信号 dt=0.1; %取时间间隔为0.1 t=0:dt:2; %时域间隔dt为间隔从0到2fc=1; %xt里最大频率xt=sin(2*pi*fc*t); %xt方程，T=1s %采样：时间连续信号变为时间离散模拟信号 fs=2.5; %抽样fs=2fc，每秒钟内的抽样点数目将等于或大于2fc个 sdt=1/fs; %频域采样间隔t1=0:sdt:2; %sdt为间隔从0到2st=sin(2*pi*fc*t1); %采样信号与输入信号对应 figure(1); subplot(311); plot(t,xt);title(输入信号); %绘出连续图形xlabel(t); %x轴显示时间t，y轴显示输入信号x(t)ylabel(x(t);grid on;axis(-inf,inf,-2,2); %inf:infinite,无穷大.-infxinf，-2y2 subplot(312); stem(t1,st,.); %绘出杆状离散图xlabel(t);ylabel(s(t); %y轴显示抽样信号s(t)axis(-inf,inf,-2,2); title(抽样信号); grid on %量化 n=length(st); %取st的长度为n，本题函数n为101 M=max(st); %M=1 c=(st/M)*2048 ; %编码 %c1(极性码)， c2c3c4(段落码），c5c6c7c8(段内电平码） code=zeros(n,8); %产生n*8的零矩阵%极性码c1 for n=1:n %if循环语句 if c(n)0 code(n,1)=1; %c为正值，c1=1 else code(n,1)=0; %c为负值，c1=0 end % 段落码c2c3c4 段落1-8，二进制表示即为0-7(000-111)if abs(c(n)=0&abs(c(n)16 code(n,2)=0;code(n,3)=0;code(n,4)=0;step=1;start=0; elseif 16=abs(c(n)&abs(c(n)32 code(n,2)=0;code(n,3)=0;code(n,4)=1;step=1;start=16; elseif 32=abs(c(n)&abs(c(n)64 code(n,2)=0;code(n,3)=1;code(n,4)=0;step=2;start=32; elseif 64=abs(c(n)&abs(c(n)128 code(n,2)=0;code(n,3)=1;code(n,4)=1;step=4;start=64; elseif 128=abs(c(n)&abs(c(n)256 code(n,2)=1;code(n,3)=0;code(n,4)=0;step=8;start=128; elseif 256=abs(c(n)&abs(c(n)512 code(n,2)=1;code(n,3)=0;code(n,4)=1;step=16;start=256; elseif 512=abs(c(n)&abs(c(n)1024 code(n,2)=1;code(n,3)=1;code(n,4)=0;step=32;start=512; elseif 1024=abs(c(n)&abs(c(n)=2048 code(n,2)=1;code(n,3)=1;code(n,4)=1;step=64;start=1024; end%段内码c5c6c7c8B=floor(abs(c(n)-start)/step); %对十进制段内码B进行floor取整，四舍五入 t=dec2bin(B,4)-48; %dec2bin定义,将B变为4位2进制码，-48改变格式 code(n,5:8)=t(1:4); %输出段内码 end code=reshape(code,1,8*n); %reshape代表从新塑形cod