通信原理软件实验.pdf

《通信原理软件实验.pdf》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

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

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

3、测4个及以上的连零来插入破坏符号（V）。倘若两个相邻破坏符号间的非零符号有偶数个时，还要插入B符号。2、数字带通调制仿真2DPSK数字调制系统的调制过程可以认为成差分双极性非归零矩形脉冲与高频载波的乘积。至于解调过程，我在本次实验中选取的为相位比较法。基本原理是将接收到的调制波形延迟一个码元周期，再与未延迟的调制波形相乘，分析该信号，从而解调码元。三、仿真方案2欢迎下载。精品文档1、基带码型仿真产生单双极键入占空比性归零码元产生相应码元数量的随机数将该数列采样，采样结果进行频输出到时域图表域处理，计算随上机码元的功率谱2、数字带通调制仿真调制波形输设定载波相乘观测波

4、形出产生基带码转换成差分双极性非归零矩形脉加入噪声元冲低通滤波相乘解调器输入（未加）延迟一个码抽样判别元周期得到解调码输出码元信元信号号四、实验结果1、基带码型仿真3欢迎下载。精品文档1）占空比为0.5占空比为0.754欢迎下载。精品文档占空比为1分析：由实验结果可知，随着占空比的增大，码元信号的直流分量会逐渐增多，频谱变窄。而且，相较于单极性，双极性的功率谱要更宽一些。2）（参数：HDB3码中V为1.5，B为2）5欢迎下载。精品文档分析：通过以上实验结果可见，AMI和HDB3码波形正确，而二者的功率谱在如此有限的码元数量内看不出来太多的区别。相比之下，HDB3的功

5、率谱成分要更复杂，更多。2、数字带通调制仿真上图为综合了调制和解调结果的总图。（仅取了一种信噪比）6欢迎下载。精品文档分析可得：调制完成后的2DPSK波形是正确的，叠加噪声后，送入解调器，得到的码元序列如图，可见解调结果还是比较正确和稳定的。上图为选取了100000个码元情况下的解调器信噪比和误码率之间的关系。欲要求误码率在104，则只需保证信噪比0.5dB即可。7欢迎下载。精品文档上图为选取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);forloop=1:1000a=(randn(1,M)>0);%产生单极性数据a1=sign(randn(1,M));%产生双极性数据tmp=zeros(L,M)

7、;tmp1=zeros(L,M);L1=L*D;%D是占空比8欢迎下载。精品文档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).^2/T;%样本信号的功率谱密度P1=abs(S1).^2/T;fin=fin*(1-1/loop)+P/loop;fin1=fin1*(1-1/loop)+P1/loop;endsubplot(2,2,1)plot(t,s);axis([-2,2