通信原理软件实验

1、叮叮小文库雨流3汽2孝实验报告课程名称通信原理实验项目名称通信原理软件实验实验类型实验学时班级学号姓名指导教师实验室名称实验时间实验成绩实验原理（1分）实验步骤 （2分）实验结果 （2分）总成绩教师签字日期哈尔滨工程大学教务处制通信原理软件仿真实验一、实验题目1、基带码型仿真1）通过仿真观察占空比为 50%、75%以及100%的单、双极性归零码波形以 及其功率谱，分析不同占空比对仿真结果的影响。2）通过仿真产生一随机消息码序列，将其分别转换为 AMI码和HDB3码， 观察它们的波形及其功率谱密度。2、数字带通调制仿真设计一个采用2DPSK调制的数字通信系统：产生二进制随机数据，并仿 真其对应的

2、2DPSK调制波形，分析其频谱。所产生的调制波形加入不同信 噪比的白噪声，选取合适的接收方案，画出系统误码率曲线，并与理论误码 率进行对比。二、实验基本原理1、基带码型仿真想要产生不同占空比的单、双极性归零码波形，首先要确定码元序列，其次 要对码元序列进行采样并输出到图表上。确定码元序列的原理很简单，对于单极性码元，只需要产生一系列随机数并 判断随机数是否大于零即可。若大于等于零则码元为1,若小于零则码元为00对于双极性码元，可以通过判断随机数的极性来产生码元。AMI码为传号交替反转码，为1反转，否则归零。而HDB3码则是在AMI的基础上进行变换。通过检测 4个及以上的连零来 插入破坏符号（V

3、）。倘若两个相邻破坏符号间的非零符号有偶数个时，还要插 入B符号。2、数字带通调制仿真2DPSK数字调制系统的调制过程可以认为成差分双极性非归零矩形脉冲与 高频载波的乘积。至于解调过程，我在本次实验中选取的为相位比较法。基本原理是将接收到 的调制波形延迟一个码元周期，再与未延迟的调制波形相乘，分析该信号，从而18解调码元。三、仿真方案1、基带码型仿真四、实验结果1、基带码型仿真1)单极件功率谱双极性D.Ofi0.060040.0200.040.030 020.010-200-1000100200200-1C00100200刀段神.功率.诺0 080 060 040.020-200 -10001

4、00200占空比为0.5冷极性力率道0.2 -0.150.10.0So卜-200 -1000100200占空比为0.75m极性祖稠性功，酒极神功率酒-2C0 -100 Q 1Q0200JL.0.1S0.10.050-200 -1000100200占空比为1分析：由实验结果可知，随着占空比的增大，码元信号的直流分量会逐渐增 多，频谱变窄。而且，相较于单极性，双极性的功率谱要更宽一些。2）电极性（参数：HDB3码中V为1.5, B为2）分析：通过以上实验结果可见，AMI和HDB3码波形正确，而二者的功率谱 在如此有限的码元数量内看不出来太多的区别。 相比之下，HDB3的功率谱成分 要更复杂，更多。

5、2、数字带通调制仿真2DPSK调制波形噪声信噪比=48,53882F。QWWVaAWWVWVXAA2 -5-4-3-2-1012345c相位比较法解调出的码元识丹率付司=0.052 :Q .L .J L,_2 I| |5<-3-2X。1234s上图为综合了调制和解调结果的总图（仅取了一种信噪比）分析可得：调制完成后的2DPSK波形是正确的，叠加噪声后，送入解调器, 得到的码元序列如图，可见解调结果还是比较正确和稳定的。上图为选取了 100000个码元情况下的解调器信噪比和误码率之间的关系欲要求误码率在10 4,则只需保证信噪比 0.5dB即可上图为选取100000个码元情况下的误码率理论

6、曲线。可见仿真得到的曲线 要比理论曲线的误码率更小。原因可能是该处的信噪比较小，与书中理论公式所 要求的大信噪比条件不符，故而误差较大。五、程序源代码1、单双极性码元序列及可调的占空比：function output_args = exam_1( D )L=64; %每码元采样点数N=512; %采样点数M=N/L; %码元数Rs=2; %码元速率(个)/sTs=1/Rs;%比特间隔fs=L/Ts;%采样速率%外域采样点%产生单极性数据%产生双极性数据T=N/fs;%截短时间t=-(T/2):1/fs:(T/2-1/fs);fin=zeros(1,N);fin1=zeros(1,N);for

7、loop=1:1000a=(randn(1,M)>0);a1=sign(randn(1,M);tmp=zeros(L,M);tmp1=zeros(L,M);L1=L*D;%口是占空比tmp(1:L1,:)=ones(L1,1)*a;tmp1(1:L1,:)=ones(L1,1)*a1;s=tmp(:)'s1=tmp1(:)'f,S=T2F(t,s);f1,S1=T2F(t,s1);P=abs(S).A2/T;%羊本信号的功率谱密度P1=abs(S1).A2/T;fin=fin*(1-1/loop)+P/loop;fin1=fin1*(1-1/loop)+P1/loop;e

8、ndsubplot(2,2,1) plot(t,s);axis(-2,2,-2,2);title( '单极性);grid on;subplot(2,2,2)plot(t,s1);axis(-2,2,-2,2);title( '双极性);grid on;subplot(2,2,3)plot(f,fin)；title( '单极性功率谱')； subplot(2,2,4)plot(f1,fin1);title('双极性功率谱');end2、AMI和HDB3码的产生：function output_args = AMI_HDB3( D )L=64; %每

9、码元采样点数N=1280; %采样点数M=N/L; %码元数Rs=2; %码元速率(个)/sTs=1/Rs; %比特间隔fs=L/Ts;%采样速率T=N/fs;%截短时间t=-(T/2):1/fs:(T/2-1/fs);%外域采样点tmp=zeros(L,M);a=(randn(1,M)>0);%产生单极性数据L1=L*D; %口是占空比tmp(1:L1,:)=ones(L1,1)*a;s=tmp(:)'a_AMI=zeros(1,M);k=1;m=0;tmp_AMI=zeros(L,M);%生AMI码元数据for i=1:Mif (a(i)=1)if (k=1)a_AMI(i)

10、=1;m=i;k=k+1;elsea_AMI(i)=-a_AMI(m);m=i;k=k+1;end else a_AMI(i)=0; endend tmp_AMI(1:L1,:)=ones(L1,1)*a_AMI; s_AMI=tmp_AMI(:)'%产生HDB3码数据a_HDB3=a_AMI;V=1.5;B=2;m=0;k0=0; %连零计数f0=0; %V之间的非零计数v0=0; %V计数flag=0;%是否有B出现for i=1:Ma_HDB3(i尸a_HDB3(i)*(-1(flag);if (a_HDB3(i)=0) k0=k0+1;if (k0=4)if (m=0) a_H

11、DB3(i)=V;k0=0;v0=1;f0=0;elsea_HDB3(i)=a_HDB3(m)*V;v0=v0+1;k0=0;if (v0=1) f0=0;elseif (mod(f0,2)=0) f0=0;i1=m+1;a_HDB3(i1)=-a_HDB3(m)*B;i1=m+4;a_HDB3(i1)=sign(a_HDB3(m+1)*V; flag=flag+1;else f0=0;endend end end else k0=0;f0=f0+1; m=i;endendtmp_HDB3=zeros(L,M);tmp_HDB3(1:L1,:)=ones(L1,1)*a_HDB3; s_HDB

12、3=tmp_HDB3(:)'%求取功率谱f,S_AMI=T2F(t,s_AMI); f1,S_HDB3=T2F(t,s_HDB3);P=abs(S_AMI).A2/T; P1=abs(S_HDB3).A2/T;figure(1)subplot(3,1,1) plot(t,s);axis(-5,5,-2,2);title( '单极性);grid on;subplot(3,1,2) plot(t,s_AMI); axis(-5,5,-2,2);title( 'AMI 码')； grid on;subplot(3,1,3) plot(t,s_HDB3); axis(-

13、5,5,-2,2);title( 'HDB3 码)； grid on;figure(2) subplot(2,1,1) plot(f,P)；title( 'AMI 码功率谱');grid on;subplot(2,1,2) plot(f1,P1);title( 'HDB3 码功率谱');grid on; end3、2DPSK 调制：function wm,xzb = DPSK( D,M ) L=64; %每码元采样点数N=L*M; %采样点数 % M=N/L;Rs=2; %码元速率(个)/sTs=1/Rs;%比特间隔fs=L/Ts;%采样速率T=N/fs

14、;%截短时间t=-(T/2):1/fs:(T/2-1/fs);%外域采样点f0=2; zb=cos(2*pi*f0*t-pi/2);厅生载波数据tmp=zeros(L,M);a=(randn(1,M)>0);%产生单极性数据L1=L*1;tmp(1:L1,:)=ones(L1,1)*a; s=tmp(:)'a_AMI=zeros(1,M); tmp_AMI=zeros(L,M); a_AMI(1)=1;%产生差分双极性非归零码元数据for i=2:Mif (a(i)=1)a_AMI(i尸-a_AMI(i-1);elsea_AMI(i)=a_AMI(i-1); end end tm

15、p_AMI(1:L1,:)=ones(L1,1)*a_AMI; s_AMI=tmp_AMI(:)'s_jt=zeros(1,N);tz_DPSK=s_AMI.*zb;%产生调制波形noise=randn(1,N)*D;%D 为噪声比例P_noise=(sum(noise.A2)/T;%出噪声功率P_sig=(sum(tz_DPSK.A2)/T;%t 出信号功率tz_noise=tz_DPSK+noise;吸口入噪声xzb=P_sig/P_noise;for n=65:Ns_jt(n尸tz_noise(n-64);endjt_DPSK=s_jt;犷生解调波形(携带噪声)% 相位比较法解调

16、-fx=tz_noise.*jt_DPSK;jt_sig=zeros(1,N);%率调出的码元数据jt=zeros(1,M);%率调出的码元序列for loop=1:Nif (fx(loop)<0)jt_sig(loop)=1;endendk=1;cnt0=0;cnt1=0;m=1;for loop=1:Nif (jt_sig(loop)=0) cnt0=cnt0+1;k=k+1;if (k=65)k=1;if (cnt0>cnt1) jt(m)=0;m=m+1;cnt0=0;cnt1=0;else jt(m)=1; m=m+1; cnt0=0; cnt1=0;end end el

17、se cnt1=cnt1+1; k=k+1;if (k=65)k=1;if (cnt0>cnt1) jt(m)=0;m=m+1;cnt0=0;cnt1=0;else jt(m)=1; m=m+1; cnt0=0; cnt1=0;end end endendtmp_DPSK=zeros(L,M);tmp_DPSK(1:L1,:)=ones(L1,1)*jt;s_DPSK=tmp_DPSK(:)'%计算误码率x1=0;for x=1:Mif (jt(x)=a(x) x1=x1+1;end end wm=x1/M;figure(1)subplot(4,1,1) plot(t,s);ax

18、is(-T/2,T/2,-2,2);title('单极性码元');grid on;subplot(4,1,2)plot(t,tz_DPSK);axis(-T/2,T/2,-2,2);title( '2DPSK 调制波形');grid on;subplot(4,1,3) plot(t,tz_noise);title( '2DPSK 调制波形 + 噪声信噪比=',num2str(xzb);grid on;subplot(4,1,4) plot(t,s_DPSK); axis(-T/2,T/2,-2,2);title( '相位比较法解调出的码元误码率(dB尸，num2str(wm);grid on;end3、计算误码率图表：function output_args = cul( M )h=waitbar(0, 'Please wait.'); %&J建进度条P=19;x=zeros(1,P);y=zeros