数字通信系统的性能分析与仿真

淮海工学院

课程设计报告书

课程名称：综合专业课程设计

题目：数字通信系统的性能分析与仿真系（院）：电子工程系

学期：2011-2012-1

专业班级：通信081

姓名：宗淙

学号：510822114

一. 课程设计的目的

专业课程设计是本科教学中极为重要的实践性教学环节，它不但起着提高本科教学质量、水平和检验学生对课程内容掌握程度的作用，而且还起到从理论过渡到实践的桥梁作用。使学生得到一次工程技术实践方面的锻炼。为在下一阶段的教学环节(毕业设计)中,能够独立地进行工程设计工作，奠定一个良好的基础。

二.设计步骤

课程设计要求的主要步骤有：

1、根据设计课题的技术指标和给定条件，在教师指导下，能够独立而正确地进行方案论证，要求概念清楚、方案合理、方法正确、步骤完整。

2、要求了解有关通信电子工程中各种常用软件的使用方法以及常用的EDA手段。

3、通过专业课程设计中的方案选择和论证，了解、熟悉、并进一步掌握电子工程及其通信专业工程设计的特点、思路、以及具体的方法和步骤,以提高独立工程设计的能力。

4、通过专业课程设计中的设计计算、软件编制，产品图纸和系统线路的绘制与制作，提高计算机软、硬件设计，工程计算和工程制图的能力。

5、培养查阅相关参考资料和手册的能力。

6、通过专业课程设计，培养严谨的科学思维方法。

三.设计内容。

数字通信系统的性能分析与仿真。应该包括以下设计内容：

1.使用一种分组码或者卷积码进行信道纠错编码。

2.使用格雷码对数据进行映射。

3.使用MQAM进行调制，M是64，128，256，512。

4.选择合适的升余弦函数，使用升余弦对基带信号进行滤波。

5.在解调端，进行滤波，MQAM的解调，格雷码的逆映射，纠错解码。

6.改变信噪比，分析系统性能。

四.设计内容介绍：

MQAM是多进制正交幅度调制是在中、大容量数字微波通信系统中大量使用的一种载波控制方式。研究这种基本的数字调制信号的性能可以帮助学生理解数字通信的基本特点。

如信道纠错编码使用分组码，M必须选择数字64，128，256，512，以分析各种M 下的QAM系统性能。应用Matlab进行仿真，仿真采用蒙特卡罗模型。仿真的基本框图是：

信号源---信道编码---映射---调制---滤波—高斯信道---滤波---解调-----逆映射-----译码---信宿

信号源：函数信号发生器。

信道编码：可以选择卷积码，分组码。

调制：MQAM。M必须是数字64，128，256，512。在调制前使用格雷码进行映射。

信道：信号经过调制后。通过信道。信道可以选择高斯加性白噪声信道。设置不同的信道信噪比，对系统进行仿真，分析不同信噪比之下的系统性能。

解调：QAM。

译码：根据信道编码方式，选择对应的信道解码方式。

性能分析：信号经过调制，信道，解调过程。在接受端，将得到的数据于原来的信号源数据比较，得到在特定信噪比的误码率。改变系统信噪比，从而得到系统的误码率曲线图。

五设计过程：

1．仿真模型

（1）根据选择的调制信号方式，确定仿真流程及框图。

（2）确定仿真中所用的函数，模块的具体参数。

2．软件设计：

编写M文件，在程序中依次改变信噪比，得到在特定信噪比下的误码率。信噪比范围是-10dB---30dB之间，根据得到的误码率数据，给出误码率---信噪比的曲线图。3．仿真调试：

在完成设计草案后可以进行仿真，并且对所设计的系统进行修改，并同理论数据进行比较。完成以下内容：

1、选择一种在具体的信道编码，信道类型下的MQAM的数字调制方式。

2、结合具体数字调制方式，选择四种MQAM进行系统仿真。

3、绘出没有信道编码情况下的该数字调制系统下的误码率---信噪比曲线图。

4、绘出有信道编码情况下的该数字调制系统下的误码率---信噪比曲线图。

5、对有无信道编码情况下的系统误码率---信噪比曲线图进行比较，并分析。

5．1、分组信道编码、高斯白噪声信道下的16QAM数字调制方式

程序：

%产生二进制整数流（图1）

M=16; %调制系数

k=log2(M);

n=3e4; %数据流长度

nsamp=1;

data=randint(n,1); %产生随机二进制数据流

figure; %绘出前40个点的柱状图

stem(x(1:40),'filled');

title('随机二进制数');

xlabel('二进制的值');ylabel('比特数');

00.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1

随机二进制数

比特数二进

制的

值

图1 随机二进制数的柱状图

%定义滤波参数（图2）

nsamp=4;

filtorder=40;

delay=filtorder/(nsamp*2);

rolloff=0.25; %%滚降系数

rrcfilter=rcosine(1,nsamp,'fir/sqrt',rolloff,delay); %%升余弦平方根滤波

figure;impz(rrcfilter,1);

图2 初始化的升余弦滤波器的单位脉冲曲线图

%分组码纠错编码

genmat=[[1 0 1;0 1 1;0 1 0],eye(3)];

bianma=encode(x,6,3,'linear',genmat);

%准备调制，格雷码映射

mapping=[0 1 3 2 4 5 7 6 12 13 15 14 8 9 11 10].';

xsym=bi2de(reshape(bianma,k,length(bianma)/k).','left-msb'); xsym=mapping(xsym+1);

figure; （图3）

stem(xsym(1:10),'filled');

title('随机符号');

xlabel('符号的值');ylabel('整数值');

%进行16QAM调制（图4）

y=qammod(xsym,M);

%进行平方跟升余弦滤波

ytx=rcosflt(y,1,nsamp,'filter',rrcfilter); %yj表示传送的信号