哈工程通信原理软件仿真实验报告.doc

试验讲演课程称号通讯道理数字基带旌旗灯号MATLAB软件仿真试验数字调制解调MATLAB软件仿真试验试验工程称号试验范例班级试验学时学号姓名指点老师试验时刻试验室称号试验道理〔1分〕试验步调〔2分〕试验后果〔2分〕总成果试验成果老师具名日期哈尔滨工程年夜学教务处制..试验一基带码型仿真〔一〕单、双极性归零码仿真一、试验道理1.1归零码归零码，是旌旗灯号电平在一个码元之内都要规复到零的编码方法，曼彻斯特编码两种编码方法它包含曼彻斯特编码跟差分1.2单、双极性归零码关于传输数字旌旗灯号来说，最常用的办法是用差别的电压电平来表现两个二进制数字，旌旗灯号由矩形脉冲构成即数字A)单极性不归零码，无电压表现〞0〞，恒定正电压表现〞1〞，每个码元时刻的两头点是采样时刻，裁决门限为半幅电平单极性归零码〔RZ〕等于以高电平跟零电中分不表现二进制码1跟0，并且在发送码1时高电平在全部码元时期T只继续一段时刻τ，其余时刻前往零电平．在单极性归零码中，τ/T称为占空比．单极性归零码的要紧长处是能够直截了当提取同步旌旗灯号，因而单极性归零码经常用作其余码型提取同步旌旗灯号时的过渡码型．也确实是说其余合适信道传输但不克不及直截了当提取同步信号的码型，可先变更为单极性归零码，而后再提取同步旌旗灯号B)双极性不归零码，〞1〞码跟〞0〞码都有电流，〞1〞为正电流，〞0〞为负电流，正跟负的幅度相称，裁决门限为零电平。

双极性归零码是二进制码0跟1分不对应于正跟负电平的波形的编码，在每个码之间都有空隙发生．这种码既存在双极性特征，又存在归零的特征．双极性归零码的特色是：接纳1比特的信息已接纳终了，而后预备下一比特信息信息．能够以为正负脉冲的前沿起了起动旌旗灯号的端依照接纳波形归于零电平就能够裁决的接纳，因而发送端不用按必定的周期发送感化，后沿起了停止旌旗灯号的感化．因而能够经常坚持准确的比特同步．即收发之间无需特不的准时，且各标记独破地形成起止方法，此方法也叫做自同步方法．因为这一特征，双极性归零码的使用非常普遍1.3功率谱密度求旌旗灯号的功率谱，功率谱=旌旗灯号的频率的相对平方/传输序列的继续时刻，求得的功率谱进展单元换算以dB值表现..1.4占空比(DutyRatio)在电信范畴中有如下含意：比方：脉冲宽度1，旌旗灯号周期4的脉冲序列占空比为0.25在一段延续任务时刻内脉冲占用的时刻与总时刻的比值二、仿真计划2.1顺序流程图发生单极性归零码并盘算功率界说占空比绘制图形界说其余变量发生双极性归零码并盘算功率绘制图形2.2参数计划分不设置占空比为0.5、0.75、1.0采样点数为2^k，k取正整数每码元采样点数为64码元速度为2Mb/s..3.试验顺序globaldttdfNcloseallk=14;Rt=0.5;N=2^k;L=64;%占空比%采样点数%每码元的采样点数%码元数M=N/L;Rb=2;%码速度为2Mb/s%码元距离Ts=1/Rb;dt=Ts/L;%时域采样距离%频域采样距离%截短时刻df=1/(N*dt);T=N*dt;Bs=N*df/2;%零碎带宽t=linspace(-T/2,T/2,N);%时域横坐标%频域横坐标f=linspace(-Bs,Bs,N);EP=zeros(1,N);forjj=1:100a=round(rand(1,M));s=zeros(1,N);%发生M个取值0，1等概的随机码%发生一个N个元素的零序列forii=1:Rt*Ts/dts(ii+[0:M-1]*L)=a;end%发生单极性归零码Q=t2f(s);%傅氏变更P=Q.*conj(Q)/T;EP=(EP*(ii-1)+P)/ii;aa=30+10*log10(EP+eps);end%P为单极性归零码的功率%累计均匀%加eps以防止除以零subplot(2,2,2)plot(f,aa,'r')title('单极性归零码的功率谱')xlabel(ylabel('f/MHZ')'Ps(f)/MHZ')axis([-15,+15,-50,50])gridsubplot(2,2,1)plot(t,s,'g')'单极性归零码')title(xlabel('t(ms)')..ylabel('s(t)(V)')axis([-5,5,-0.5,1.5])gridforjj=1:100a=round(rand(1,M));a=1-2*a;%发生M个取值0，1等概的随机码%发生一个N个元素的零序列%发生双极性归零码s=zeros(1,N);forii=1:Rt*Ts/dts(ii+[0:M-1]*L)=a;endQ=t2f(s);%付氏变更P=Q.*conj(Q)/T;EP=(EP*(ii-1)+P)/ii;%P为双极性归零码的功率%累计均匀aa=30+10*log10(EP+eps);end%加eps以防止除以零subplot(2,2,4)%set(2,'position',[10,50,750,350])%设定窗口地位及巨细plot(f,aa,'r')title('双极性归零码的功率')xlabel(ylabel('f(MHZ)')'Ps(f)(MHZ)')axis([-15,+15,-50,50])gridsubplot(2,2,3)plot(t,s,'g')'双极性归零码')title(xlabel(ylabel('t(ms)')'s(t)(V)')axis([-5,5,-1.5,1.5])grid子顺序functionX=t2f(x)globaldtdfNtfT%X=t2f(x)%X与x长度一样,并为2的整幂。

本函数需求一个全局变量H=fft(x);dt(时域取样距离)X=[H(N/2+1:N),H(1:N/2)]*dt;end..图一占空比为0.5的波形图二占空比为0.75的波形..图三占空比为1.0的波形〔二〕AMI码跟HDB3码仿真一、试验道理1.1AMI码1码平日称为传号，0码那么叫空号，这是相沿了晚期电报通讯中的叫法从形状上看，它已是三形状旌旗灯号，因而AMI码是伪三进制码〔1〕编码规那么:音讯代码中的0传输码中的0音讯代码中的1传输码中的+1、-1瓜代比方:音讯代码:1010100010111AMI码:+10-10+1000-10+1-1+1〔2〕AMI码的特色:1由AMI码断定的基带旌旗灯号中正负脉冲瓜代，而0电位坚持稳定;因而由AMI码断定的基带旌旗灯号无直流重量，且只要非常小的低频重量;2在接纳端不易提取准时旌旗灯号，因为它能够呈现长的连0串;..3存在检错才能，假如在全部传输进程中，因传号极性瓜代法则遭到毁坏而呈现误码时，在接纳端非常轻易发觉这种过错〔3〕解码规那么从收到的标记序列中将一切的-1变更成+1后，就能够失掉原音讯代码1.2HDB3码三阶高密度双极性码〔英语：HighDensityBipolarofOrder3，简称：HDB3码〕是一种实用于基带传输的编码方法，它是为了克制毁坏性强等特色。

AMI码的缺陷而呈现的，存在能量疏散，抗三阶高密度双极性码用于一切档次的欧洲E-carrier零碎，HDB3码将4个延续的“0位〞元代替成“000V〞或“B00V〞那个做法能够确保延续的是相隔双数的普通B暗号violationsareofdifferingpolarity，即1先将音讯代码变更成AMI码，假定AMI码中连0的个数小于4,如今的AMI码确实是HDB3码;2假定AMI码中连0的个数年夜于3,那么将每4个连0小段的第4个0变更成与前一个非0标记(+1或-1)同极性的标记,用表现(+1+,-1-);3为了不毁坏极性瓜代反转，当相邻V标记之间有偶数个非0标记时，再将该小段的第1个0变更成+B或-B,标记的极性与前一非零标记的相反，并让前面的非零标记从标记开始再瓜代变更比方:音讯代码:100001000011000011AMI码:+10000-10000+1-10000+1-1HDB3码：+1000+V-1000-V+1-1+B00+V-1+1二、试验顺序globalclosedtdftfNallN=2^14;L=64;%采样点数%每码元的采样点数%码元数M=N/L;Rb=2;..Ts=0.5;%码元宽度是0.5usdt=Ts/L;df=1/(N*dt);%MHzRT=0.5;%占空比T=N*dt;Bs=N*df/2;%截短时刻%零碎带宽t=[-T/2+dt/2:dt:T/2];f=[-Bs+df/2:df:Bs];%时域横坐标%频域横坐标EPAMI=zeros(size(f));EPHDB=zeros(size(f));forii=1:8ami=zeros(1,M);hdb=zeros(1,M);a=round(rand(1,M));b=3;%表现0000之间轮回个数c=-1;%记录相邻V之间的1元素个数sign1=-1;%标记前一个旌旗灯号%标记前一个旌旗灯号sign2=-1;forii=1:Mifa(ii)==1sign1=0-sign1;ami(ii)=sign1;endendforii=1:Mifb==3ifa(ii)==1%表现非0000sign2=0-sign2;hdb(ii)=sign2;ifc>=0%表现不是第一个0000c=c+1;%用来盘算相邻v之间的非0元素个数endelseifii<=M-3&a(ii)+a(ii+1)+a(ii+2)+a(ii+3)==0ifmod(c,2)==1%000Vhdb(ii+3)=sign2;else%B00Vsign2=0-sign2;hdb(ii)=sign2;hdb(ii+3)=sign2;endc=0;b=0;end..elseifb<3%对0000的轮回b=b+1;endendfori=[1:L]ami1(i+[0:M-1]*L)=ami;hdb1(i+[0:M-1]*L)=hdb;endAMI=t2f(ami1);PAMI=AMI.*conj(AMI)/T;HDB=t2f(hdb1);PHDB=HDB.*conj(HDB)/T;EPAMI=(EPAMI*(i。