仲恺现代通信原理实验报告..docVIP

实验一 PCM脉冲编码调制 信息科学与技术学院学院（院、系） 网络工程 专业 132 班 现代通信系统 课 学号 姓名 实验日期 2015-4-26 教师评定 实验一：利用Matlab绘制带通信号x(t)=2sinc(20t)*cos[2π*100t+sinc(5t)],时间间隔为0.02s。 代码： 图像： ts=0.02; t=[-3:ts:3]; x=2*sinc(20*t).*cos(2*pi*100*t+sinc(5*t)); plot(t,x) 实验二：利用Matlab对模拟信源s=sint(0t2π)进行均匀量化，量化间隔为0.2s。 代码： 图像： t=[0:0.2:2*pi]; s=sin(t); partition=[-1:0.2:1]; codebook=[-5:1:5+1]; [index,quants]=quantiz(s,partition,codebook); subplot(2,1,1);plot(t,s); subplot(2,1,2);plot(t,quants) 实验三：编制一个函数实现均匀PCM量化编码，并计算量化噪声比(SQNR)。 代码： function [sqnr,a_quan,code]=upcm(a,n) %定义一个关于输入信号序列a和量化级数n的upcm函数 amax=max(abs(a)); %取变量amax等于序列a的绝对值 a_quan=a/amax; %对输入信号序列归一化，使信号幅度 取值范围为[-1,1] b_quan=a_quan; %令变量b_quan等于变量a_quan d=2/n; %取d=2/n为量化间隔 q=d.*[0:n-1]-(n-1)/2*d; %取q为每个量化区间对应的判决阈值 for i=1:n %对归一化后的输入信号序列进行量化 index=find((q(i)-d/2=a_quan)(a_quan=q(i)+d/2)); a_quan(index)=q(i)*ones(1,length(index)); b_quan(find(a_quan==q(i)))=(i-1).*ones(1,length(find(a_quan==q(i)))); end a_quan=a_quan*amax; %使量化后的归一化信号各点值变回 原来的值 nu=ceil(log2(n)); %设定给定量化级数所需比特数 code=zeros(length(a),nu); %取零矩阵，使其行数为序列a的长度，列数为量化所需比特数nu的矩阵 for i=1:length(a) %对输入信号序列量化后进行编码 for j=nu:-1:0 if(fix(b_quan(i)/(2^j))==1) code(i,nu-j)=1; b_quan(i)=b_quan(i)-2^j; end end end sqnr=20*log10(norm(a)./norm(a-a_quan)); %使公式计算量化噪声比(SQNR) 实验四：利用上题编制的函数，对正弦信号s=sint(0t2π)进行均匀PCM量化编码，并比较当量化级数分别为8、16时的量化噪声比的大小。 代码： 图像： t=[0: