通信原理研究型学习报告 2.doc

通信 原理 设置目的： 检验学生对基带调制知识点的理解，掌握和运用程度，锻炼其分析问题，开展研究的能力。 内容要求： PCM，LPCM编译码仿真分析， 仿真实现A律13折线编译码，比较理论以及实际量化信噪比，信号源自定（包括语音信号的编译码）。 仿真结果： 假定抽样值Sa=150 （1）实现LPCM 正极性编码第一位为1，负极性第一位编码为0。理论信噪比采用SNR=3*M^2/Kcr^2 ，其中Kcr=x/(x/sqrt(2))；实际信噪比SN=S/N，其中S=sum(yk.^2,Sa+1)/(Sa+1) N=sum((yk-deyk).^2,Sa+1)/(Sa+1)，抽样点数越多，越接近理论量化信噪比。 对yt= cos(2000*t)进行抽样，抽样频率为fs=8000Hz，量化间隔delta=1/2^7。 信噪比：理论值SNR =43.9051dB ，实际值SN =44.4416dB .m程序如下： t=0:0.01:2*pi x=1 yt=x*cos(2000*t) M=2^7 delta=x/M fs=8000 Ts=1/fs Sa=150 ks=0:Sa yk=x*cos(ks*Ts*2000) code=zeros(1,8) CODE=zeros(Sa+1,8) for k=0:Sa if yk(k+1)=0 code(1)=1 i=2 a=floor(yk(k+1)/delta) if a1 code=[1 zeros(1,7)] elseif a==128 code=[1,ones(1,7)] elseif a1 while(a~=0) code(8+2-i)=mod(a,2) a=floor(a/2) i=i+1 end end elseif yk(k+1)0 code(1)=0 i=2 a=ceil(yk(k+1)/delta) if a-1 code=[0 zeros(1,7)] elseif a==-128 code=[0,ones(1,7)] elseif a-1 while(a~=0) code(8+2-i)=mod(a,2) a=ceil(a/2) i=i+1 end end end for i=1:8 CODE(k+1,i)=code(i) end code=zeros(1,8) end decode=zeros(1,8) deyk=double(zeros(1,Sa+1)) for k=0:Sa for i=1:8 decode(i)=CODE(k+1,i) end if decode(1)==1 deyk(k+1)=(decode(2)*2^6+decode(3)*2^5+decode(4)*2^4+decode(5)*2^3+decode(6)*2^2+decode(7)*2^1+decode(8)*2^0+0.5)*delta elseif decode(1)==0 deyk(k+1)=-(decode(2)*2^6+decode(3)*2^5+decode(4)*2^4+decode(5)*2^3+decode(6)*2^2+decode(7)*2^1+decode(8)*2^0+0.5)*delta end end plot(ks,yk,r,ks,deyk,b) figure plot(ks,yk-deyk) display(CODE) Kcr=x/(x/sqrt(2)) SNR=10*log10(3*M^2/Kcr^2) S=sum(yk.^2,Sa+1)/(Sa+1) N=sum((yk-deyk).^2,Sa+1)/(Sa+1) SN=10*log10(S/N) display(SNR) display(SN) 编译码图形 （红色是原码，蓝色的是译码之后的）差别较小，可以放大来看。 放大： 误差： （2）A率13折线编译码： 输入信号为yt= cos(2000*t)，正域值为极性值1，负域为0.理论信噪比采用分段求量化噪声，再求出总的量化噪声，进而得出量化信噪比。实际信噪比SN=vpa(Sq/N)其中，Sq=sum(yk.^2,Sa+1)/(Sa+1)，N=sum((deyk-yk).^2,Sa+1)/(Sa+1) 抽样频率f=8000Hz，最小量化间隔delta=1/