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

实 验 报 告课程名称通信原理实验项目名称数字基带信号MATLAB软件仿真实验数字调制解调MATLAB软件仿真实验实验类型实验学时班级学号姓名指导教师实验室名称实验时间实验成绩实验原理（1分）实验步骤（2分）实验结果（2分）总成绩教师签字日期哈尔滨工程大学教务处 制实验一 基带码型仿真（1） 单、双极性归零码仿真1、 实验原理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占空比(Duty Ratio)在电信领域中有如下含义：例如：脉冲宽度1s，信号周期4s的脉冲序列占空比为0.25。在一段连续工作时间内脉冲占用的时间与总时间的比值。2、 仿真方案绘制图形产生双极性归零码并计算功率产生单极性归零码并计算功率定义占空比定义其他变量2.1程序流程图 绘制图形2.2参数设计分别设置占空比为0.5、0.75、1.0采样点数为2k，k取正整数每码元采样点数为64码元速率为2Mb/s3.实验程序global dt t df Nclose allk=14;Rt=0.5; %占空比N=2k; %采样点数L=64; %每码元的采样点数M=N/L; %码元数Rb=2; %码速率为2Mb/sTs=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); for jj=1:100a=round(rand(1,M); %产生M个取值0，1等概的随机码s=zeros(1,N); %产生一个N个元素的零序列for ii=1:Rt*Ts/dt s(ii+0:M-1*L)=a; %产生单极性归零码endQ=t2f(s); %傅氏变换P=Q.*conj(Q)/T; %P为单极性归零码的功率EP=(EP*(ii-1)+P)/ii; %累计平均aa=30+10*log10(EP+eps); %加eps以避免除以零end subplot(2,2,2)plot(f,aa,r)title(单极性归零码的功率谱)xlabel(f/MHZ)ylabel(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)grid for jj=1:100a=round(rand(1,M); %产生M个取值0，1等概的随机码a=1-2*a;s=zeros(1,N); %产生一个N个元素的零序列for ii=1:Rt*Ts/dt s(ii+0:M-1*L)=a; %产生双极性归零码endQ=t2f(s); %付氏变换P=Q.*conj(Q)/T; %P为双极性归零码的功率EP=(EP*(ii-1)+P)/ii; %累计平均aa=30+10*log10(EP+eps); %加eps以避免除以零end subplot(2,2,4)%set(2,position,10,50,750,350) %设定窗口位置及大小plot(f,aa,r)title(双极性归零码的功率)xlabel(f(MHZ)ylabel(Ps(f)(MHZ)axis(-15,+15,-50,50)gridsubplot(2,2,3)plot(t,s,g)title(双极性归零码)xlabel(t(ms)ylabel(s(t)(V)axis(-5,5,-1.5,1.5)grid子程序function X=t2f(x)global dt df N t f T%X=t2f(x)%X与x长度相同,并为2的整幂。%本函数需要一个全局变量dt(时域取样间隔)H=fft(x);X=H(N/2+1:N),H(1:N/2)*dt;end图一 占空比为0.5的波形图二 占空比为0.75的波形图三 占空比为1.0的波形（2） AMI码和HDB3码仿真一、实验原理1.1 AMI码1码通常称为传号，0码则叫空号，这是沿用了早期电报通信中的叫法。从形态上看，它已是三状态信号，所以AMI码是伪三进制码。（1）编码规则:消息代码中的0 传输码中的0消息代码中的1 传输码中的+1、-1交替例如:消息代码:1 0 1 0 1 0 0 0 1 0 1 1 1AMI码: +1 0 -1 0 +1 0 0 0 -1 0 +1 -1 +1（2）AMI码的特点:1 由AMI码确定的基带信号中正负脉冲交替，而0电位保持不变;所以由AMI码确定的基带信号无直流分量，且只有很小的低频分量;2 在接收端不易提取定时信号，由于它可能出现长的连0串;3 具有检错能力，如果在整个传输过程中，因传号极性交替规律受到破坏而出现误码时，在接收端很容易发现这种错误。（3）解码规则从收到的符号序列中将所有的-1变换成+1后，就可以得到原消息代码1.2 HDB3码三阶高密度双极性码（英语：High Density Bipolar of Order 3，简称：HDB3码）是一种适用于基带传输的编码方式，它是为了克服AMI码的缺点而出现的，具有能量分散，抗破坏性强等特点。三阶高密度双极性码用于所有层次的欧洲E-carrier系统，HDB3码将4个连续的“0”位元取代成“000V”或“B00V”。这个做法可以确保连续的violations are of differing polarity，即是相隔单数的一般B记号。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,符号的极性与前一非零符号的相反，并让后面的非零符号从符号开始再交替变化。例如:消息代码: 1 0 0 0 0 10 0 0 01 1 0 0 0 0 1 1AMI码: +1 0 0 0 0 -10 0 0 0+1 -1 0 0 0 0 +1 -1HDB3码：+1 0 0 0 +V -10 0 0 -V+1 -1 +B 0 0 +V -1 +1二、实验程序global dt df t f Nclose allN=214; %采样点数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);for ii=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; %标志前一个信号for ii=1:Mif a(ii)=1 sign1=0-sign1;ami(ii)=sign1;endendfor ii=1:M if b=3 %表示非0000if a(ii)=1 sign2=0-sign2;hdb(ii)=sign2; if c=0 %表示不是第一个0000 c=c+1; %用来计算相邻v之间的非0元素个数endelseif ii=M-3 & a(ii)+a(ii+1)+a(ii+2)+a(ii+3)=0 if mod(c,2)=1 %000Vhdb(ii+3)=sign2;else %B00V sign2=0-sign2;hdb(ii)=sign2;hdb(ii+3)=sign2;end c=0; b=0;endelseif b= Z %z即为阈值 Vs(ii)=H; elseVs(ii)=L; end end figure(4) subplot(311) plot(t,Vs) title(解调后差分信号) axis(0 10e-4 -0.2 1.2) grid on; c=randi(1,22,2); %产生解调后的差分码元 for f=0:19 c(f+1)=fix(Vs(f*200+50)+0.2) end d=randint(1,21,2); %定义差分译码后的码元 for l=1:21 %得到差分译码后的码元 if (l=1) if (s(1)=0) d(1)=0; else d(1)=1; end elseif (c(l)=c(l-1)d(l)=0; else d(l)=1; end end y=d(ceil(k*t+0.01); subplot(313); plot(t,y); axis(0 10e-4 -0.2 1.2); title(码反变换输出); %基带信号与解调后的信号对比 subplot(312) plot(t,m); axis(0 10e-4 -0.2 1.2); title(原始基带信号); %误码率与信噪比关系 figure(5) grid on initial_snr = 0; final_snr = 9; %snr信噪比 snr_step = 0.01; snr_in_dB = initial_snr:snr_step:final_snr; for i = 1:length(snr_in_dB) snr = 10(snr_in_dB(i)/10);%将信噪比单位dB转化一下 Pe(i) =erfc(sqrt(snr)*(1-1/2*erfc(sqrt(snr); %2dpSK相干解调 end semilogy(snr_in_dB,Pe