通讯系统仿真综合实训报告.docx
《通讯系统仿真综合实训报告.docx》由会员分享，可在线阅读，更多相关《通讯系统仿真综合实训报告.docx（15页珍藏版）》请在冰豆网上搜索。
通讯系统仿真综合实训报告
天津电子信息职业技术学院
综合实训报告
课题名称通讯系统仿真综合实训报告
姓名温智
学号11
班级通信S09-5
专业通信技术
所在系电子技术系
指导教师张澄
完成日期2011.11.08
目录
摘要………………………………………
一实训目的1
二实训原理1
三实训正文2
实例一.2
实例二.3
实例三.4
自选题：
6
四.心得体会14
五.参考文献15
摘要
通过MATLAB仿真，具体培养对仿真软件的应用，通过仿真锻炼加深对知识的应用及作用。
一实训目的
通过课程设计进一步理解扩展频谱通信的基本概念及其系统模型；重点是伪随机编码的基本原理，m序列、Gold序列的性质及特点；扩展频谱信号的相关解扩、基带解调与载波同步，跳频信号的解跳和解调等等。
要求学生在课程设计中建立基本的扩频系统模型，仿真计算出伪随机编码的相关特性，通过扩频调制的解扩仿真系统的抗干扰性能。
二实训原理
Gold序列是m序列的组合码，由优选对的两个m序列逐位模2加得到，当改变其中的一个m序列的相位（向后移位）时，可以得到一新的Gold序列。
Gold序列虽然是由m序列模2加得到的，但它已不是m序列，不过它具有与m序列优选对类似的自相关性和互相关特性，而且构造简单，比m序列所产生的序列码组多的多，因而获得广泛的应用。
周期均为N=2*n-1的m序列优选对{an}和{bn}，{an}与后移τ位的{bn+τ}（τ=0，1，，…，N-1）逐位模2加所得的序列{an+bn+τ}即得到Gold序列，改变序列移位值τ，可以得到不同的Gold序列。
而m序列的优选对是指在m序列集中，互相关函数绝对值的最大值|RXY（τ）|max最接近或者达到相关下限（最小值）的一对m序列。
三实训正文
实例一.判断特征多项式F（x）=x9+x6+x4+x3+1是否可生成m序列，并建模验证。
F（x）对应的系数二进制表示为，相应的十进制数是601。
测试模型如下：
测试结果为：
实例二.计算特征多项式为F（x）=x9+x6+x4+x3+1的m序列的自相关系数。
对于周期N的序列，其自相关系数是偶函数，即ρ（-j）=ρ（j），而且也是以N为周期的周期函数。
周期为N的m序列自相关系数理论值为
ρ（j）=
其中k为整数。
本例中m序列的周期为N=29-1=511，首先计算出一个周期的m
序列，然后再根据自相关系数的定义进行计算，计算中应注意将二进制输出的m序列转换为取值{±1}的双极性序列，然后再求相关函数。
程序如下：
reg=ones（1,9）;%寄存器初始状态：
全1，寄存器级数为9
coeff=[1,0,0,1,0,1,1,0,0,1];%抽头系数a0a1...ar，取决于特征多项式
N=2^length（reg）-1;%周期
fork=1:
N%计算一个周期的m序列输出
a_n=mod（sum（reg.*coeff（1:
length（coeff）-1））,2）;%反馈系数
reg=[reg（2:
length（reg））,a_n];%寄存器移位,反馈
out（k）=reg
（1）;%寄存器最低位输出
end
out=2*out-1;%转换为双极性序列
forj=0:
N-1
rho（j+1）=sum（out.*[out（1+j:
N）,out（1:
j）]）/N;
end
j=-N+1:
N-1;
rho=[fliplr（rho（2:
N））,rho];
plot（j,rho）;axis（[-1010-0.11.2]）;
title（'洪松-实例2'）
实例三.计算r=6本原多项式（八进制表示）103和147对应的两个m序列的互相关函数序列。
八进制数103和147转换为二进制分别是：
和。
对应m序列的特征多项式以向量形式表示为
[1,0,0,0,0,1,1]和[1,1,0,0,1,1,1]
编写程序如下：
clear;
reg=ones（1,6）;%寄存器初始状态：
全1，寄存器级数为6
coeff=[1,0,0,0,0,1,1];%抽头系数cr...c1c0，取决于特征多项式
N=2^length（reg）-1;%周期
fork=1:
N%计算一个周期的m序列输出
a_n=mod（sum（reg.*coeff（1:
length（coeff）-1））,2）;%反馈
reg=[reg（2:
length（reg））,a_n];%寄存器移位,反馈
out1（k）=2*reg
（1）-1;%寄存器最低位输出,转换为双极性序列
end
reg=ones（1,6）;
coeff=[1,1,0,0,1,1,1];%抽头系数
fork=1:
N%计算一个周期的m序列输出
a_n=smod（um（reg.*coeff（1:
length（coeff）-1））,2）;%反馈
reg=[reg（2:
length（reg））,a_n];%寄存器移位,反馈
out2（k）=2*reg
（1）-1;%寄存器最低位输出,转换为双极性序列
end
%得出两个双极性电平的m序列
forj=0:
N-1
R（j+1）=sum（out1.*[out2（1+j:
N）,out2（1:
j）]）;%计算相关函数
end
j=-N+1:
N-1;%相关函数自变量
R=[fliplr（R（2:
N））,R];%利用相关函数的偶函数特性计算j为负值的情况
plot（j,R）;axis（[-NN-2020]）;xlabel（‘j’）;ylabel（‘R（j）’）;%作图
max（abs（R））%计算相关函数绝对值的最大值
title（'洪松-实例3'）
自选题：
基于Matlab的CDMA通信系统仿真
设计原理：
CDMA是码分多址的英文缩写（CodeDivisionMultipleAccess），它是在数字技术的分支--扩频通信技术基础上发展起来的一种崭新而成熟的无线通信技术。
CDMA技术的原理是基于扩频技术，即将需传送的具有一定信号带宽信息数据，用一个带宽远大于信号带宽的高速伪随机码进行调制，使原数据信号的带宽被扩展，再经载波调制并发送出去。
接收端使用完全相同的伪随机码，与接收的带宽信号作相关处理，把宽带信号换成原信息数据的窄带信号即解扩，以实现信息通信。
1．CDMA是扩频通信的一种，他具有扩频通信的以下特点：
（1）抗干扰能力强。
这是扩频通信的基本特点，是所有通信方式无法比拟的。
（2）宽带传输，抗衰落能力强。
（3）由于采用宽带传输，在信道中传输的有用信号的功率比干扰信号的功率低得多，因此信号好象隐蔽在噪声中；即功率谱密度比较低，有利于信号隐蔽。
（4）利用扩频码的相关性来获取用户的信息，抗截获的能力强等
。
所有CDMA类型都使用扩频过程增益来允许接受者部分衰减非期望信号。
若具有期望扩频码的信号和定时则被接受，如果信号有不同的扩频码（或者相同扩频码但是不同的时间偏移）将被过程增益认为随机噪声衰减掉。
这项操作的方法是给每一个站点分配一个扩频码或者芯片序列.这些芯片序列被表示成由＋1和－1组成的序列。
每个芯片序列和本身点乘得到＋1,（和补码点乘得到－1）,反之点乘不同的芯片序列得到0。
这种特性叫做正交性.这种序列叫做Walsh码可以从一个二进Walsh矩阵导出。
当多个终端发送多个片码时，信号就会在空中叠加。
例如芯片序列是（-1,-1,-1,-1）和（+1,-1,+1,-1），叠加后变成（0,-2,0,-2）。
接受方只要计算发送信号到空中的终端目点值。
例如（-1,-1,-1,-1）.（0,-2,0,-2）=+1。
就可以正确接收所需信息，而把其他信号当作噪声过滤掉。
CDMA通信系统原理框图
根据上述原理搭建仿真模块如下：
本实验中采用PNSequenceGenerator来产生m系列，PNSequenceGenerator的参数有Generatorpolynomial（特征多项式）：
[11001]，对应上图中a0=1,C1=1,C2=0,C3=0,C4=1;Initialstate（初始状态）:
[0100],对应a1=0,a2=1,a3=0,a4=0.周期为15，Sampletime（采样时间）即码元宽度为2e^-5,可以得到一个周期的系列为1011。
上图中三个BernoulliRandomBinaryGenerator表示三个不同的通信用户发射各自的通信信息，码元宽度为6x10^-6s,PNSequenceGenerator产生用于直接扩频的正交码组，它产生的是m系列，实验中m系列的周期为15，码元宽度为2x10^-5s,基带信号码元宽度是m系列码元宽度的30倍，正好是两个m系列的周期。
延迟4个码元及8个码元的两个码组与原始的码组构成三个正交的码组，它们分别对三个用户的信号进行直接扩频。
扩频的操作是将转换为二进制双极性信号的基带信号和用于扩频的码组直接相乘。
扩频后的信号在sum中与GaussianNoiseGenerator产生的高斯白噪声混合，这表现了码分多址通信的特点：
在同一时间、同一频段利用正交的码组承载不同通信用户的信息传输。
其余模块的功能及参数设置如下
1.
模块名称：
AWGNChannel
参数名称
参数值
功能说明：
加性高斯白噪声信道用高斯分布的噪声信号叠加在通过它的信号上，使通过该信道的信号产生与噪声均值相应的偏移，并且围绕平均值做随机波动。
Initialseed（初始化种子）
1567
SNR（dB）（信噪比）
Signaltonoiseratio（SNR）
SNR（dB）（信噪比）
8
Inputsignalpower（watts）（输入信号功率）
1
2.模块名称：
BernoulliRandomBinaryGenerator
参数名称
参数值
功能说明：
伯努利二进制随机数产生器用于产生二进制的数字系列，系列中0和1出现的概率服从伯努利分布
Probabilityofazero（0出现的概率）
0.5
0.5
0.5
Initialseed（初始化种子）
12345
54321
12345
Sampletime（采样时间）
6e-4
3.模块名称：
PNSequenceGenerator
参数名称
参数值
功能说明：
PN码序列产生器产生具有伪随机特性的二进制序列，用于扩频。
Generatorpolynomial（特征多项式）
[11001]
Initialstate（初始状态）
[0100]
Shift（ormask）（移动）
0
Sampletime（采样时间）
2e-5
4.模块名称：
GaussianNoiseGenerator
参数名称
参数值
功能说明：
加性高斯白噪声产生器产生具有高斯分布的噪声信号
Mean（均值）
0
Variance（方差）
0.1
Initialseed（初始化种子）
12345
Sampletime（采样时间）
2e-5
5.模块名称：
Relay
参数设置
参数值
功能说明：
扩频之前用于将单极性信号转换为双极性信号
上三个
下三个
Switchonpoint（大于等于时，开关打开）
0.9
0.9
Switchoffpoint（小于等于时，开关关闭）
0.1
0.1
Outputwhenon（开关打开时输出值）
1
1
Outputwhenoff（开关关闭时输出值）
-1
-1
6.模块名称：
Relay
参数设置
参数值
功能说明：
接在低通滤波器后用于对信号的抽样判决
Switchonpoint（大于等于时，开关打开）
0.51
Switchoffpoint（小于等于时，开关关闭）
0.49
Outputwhenon（开关打开时输出值）
1
Outputwhenoff（开关关闭时输出值）
-1
7.模块名称：
DigitalFilter
参数名称
参数值
功能说明：
数字滤波器实现滤波功能，解扩后经该滤波器得到一个携带信息的中频信号
Numerator（分子系数）
[0.00040.00170.00250.00170.0004]
Denominator（分母系数）
[1.000-3.18113.8623-2.11300.4385]
Initialconditions（初始条件）
0
其他模块：
模块名称：
UnitDelay
位置：
Simulink/Discrete
模块名称：
Display
位置：
Simulink/Sinks
模块名称：
DotProduct
位置：
Simulink/MathOperations
模块名称：
Sum
位置：
Simulink/MathOperations
模块名称：
Scope
位置：
Simulink/Sinks
模块名称：
ErrorRateCalculation
位置：
CommunicationsBlockset/CommSinks
运行仿真模块绘出信号的变化几个关键阶段的波形：
1.伯努利二进制信号波形
2.伯努利二进制信号双极性转换后波形
3.扩频后信号波形
4.三路信号与噪声信号叠加后波形
5.经过加性高斯白噪声信道后的信号波形
6.解扩后的信号波形
7.经过数字滤波器滤波后的波形
8.抽样判决后的波形
以上波形图可看出信号的接收端收到的信号大致和原始的信号（由发送端发送的信号）的波形一致，有些许失真。
由误码率计算知道误差是非常小的，在允许的范围。
误码率曲线的绘制过程：
（用m文件绘制误码率曲线）
程序如下：
snR1=-20:
1:
20;%信噪比范围
forn=1:
length（snR1）
SNR=snR1（n）;
sim（'cdma'）
s2（n）=[mean（s）]‘;%求s的均值
A1（n）=s2（n）+eps;%计算误码率
%eps:
从1.0到下一个最大浮点数的距离，常用来作为浮点计算相对误差使用
s3（n）=[mean（s1）]';%求s1的均值
A2（n）=s3（n）+eps;%计算误码率
s3（n）=[mean（s2）]';%求s2的均值
A3（n）=s3（n）+eps;%计算误码率
EN（n）=[snR1（n）]';
end
figure
（1）
semilogy（EN,（A1））;grid;%绘制误码率曲线
title（'信噪比与误码率关系曲线1'）;
figure
（2）;
semilogy（EN,（A2））;grid;%绘制误码率曲线
title（'信噪比与误码率关系曲线2'）;
figure（3）;
semilogy（EN,（A3））;grid;%绘制误码率曲线
title（'信噪比与误码率关系曲线3'）
经过无数次的调试，仿真系统终于基本达到的理想的要求，误码率很小，满足了通信的要求。
通过本实验可以详细具体地实现CDMA通信的过程，通过跟踪各个阶段的信号波形，对于加深CDMA通信原理的理解相当形象简明，再通过各个模块参数的设定，可以了解在实际通信过程中对信号产生影响的各种因数及影响的程度。
由此，利用这种方法对于我们加深学习其他通信技术是十分有裨益的。
四.心得体会
通过这段时间的实训，对MATLAB仿真软件有了更深的了解。
对此次实训，颇有体会。
本次实习主要任务是进一步理解扩展频谱通信的基本概念及其系统模型；重点是伪随机编码的基本原理，m序列、Gold序列的性质及特点，对此专门找了《扩展频谱通信》进行了学习，从原理到实现，一步一步从零开始，期间多次经过了同学们的帮助答疑解惑，使得这次实训能顺利完成。
经过理论与实践相结合，进一步加深我们的理论知识。
要想在短暂的实训时间内，尽可能能多的学一些东西，这就需要我们认真听取老师的讲解，按计划完成相应的任务。
这次实训让我学到的东西太多，使我受益非浅。
实训，虽然辛苦了点，但能让我学到不同的东西，我心里还是高兴的。
人非生而知之，要学得知识，一靠学习，二靠实践。
没有实践，学习就是无源之水，无本之木。
以上就是我在成都的进行实训的心得和感受。
不到半年的时间就将步入社会的我们，面临是就业的压力，我想我们更应该把握住最后的一段时间，充实、完善自我，争取做一名好的学生！
五.参考文献
[1]徐明远邵玉斌.MATLAB仿真在通信与电子工程中的应用[M].西安：
西安电子科技大学出版社，2005.6
[2]樊昌信曹丽娜.通信原理（第六版）[M].北京：
国防工业出版社，2007.3
[3]刘焕淋向劲松代少升.扩展频谱通信[M].北京：
北京邮电大学出版社，2008
[4]李建东郭梯云邬国扬.移动通信（第四版）[M].西安：
西安电子科技大学出版社，2006.12
关键词：
MATLAB仿真、AWGNChannel、BernoulliRandomBinaryGenerator、
PNSequenceGenerator、GaussianNoiseGenerator、Relay
升级会员 