毕业论文数字信号载波传输的Matlab仿真.docx
《毕业论文数字信号载波传输的Matlab仿真.docx》由会员分享，可在线阅读，更多相关《毕业论文数字信号载波传输的Matlab仿真.docx（26页珍藏版）》请在冰点文库上搜索。
毕业论文数字信号载波传输的Matlab仿真
4.4.1FSK的产生21
附录28
摘要：
数字信号与信息科学的迅猛发展，计算机和通信系统的迅速发展的理论和实践，通信系统是将信息从信源发送到一个或者多个目的地，用来达到各种信息传输过程的一般技术体系。
数字信号调制与解调过程中起着重要的信息传递的作用。
为了结合计算机仿真，美国MathWork公司推出的MATLAB软件是一个最优秀的科技应用软件。
MATLAB完全可以将理论与实际应用结合在一起，以此激发的学生的学习兴趣，以此提升了学习的效果。
为了学习更方便的数字信号，美国公司将MATLAB特设了学习辅助软件作为学习数字信号基础，数字信号载波传输的分析及仿真完全充分的利用了MATLAB快速操作，文字，图形，声音及交互式人机界面等特点来进行的。
MATLAB是一种将强大功用结合为一体的科学计算机语言，它拥有了数值运算、制作图像、图像操作和系统仿真等，它工具箱和功能非常丰富，拥有强大的矩阵运算、图形可视化，因此提供了极大的方便在通信系统的调制和解调过程的分析中。
本论文大概介绍了数字信号系统的概念，接着介绍了调制和解调的方式，继而具体介绍了我们经常使用的几种调制和解调的方式和数字信号载波传输在MATLAB中的仿真。
数字信号课程也是我们专业的必修课程。
关键字：
数字信号载波传输，MATLAB，仿真
ABSTRACT：
Digitalsignaltheoryandpracticalapplicationwiththerapiddevelopmentofinformationscienceandcomputerandrapiddevelopment,thecommunicationsystemistosendinformationfromthesourcetooneormoredestinations,isusedtoimplementthefloorboardofthetechnologyofinformationtransmissionsystem.Intheprocessofthedigitalsignalmodulationanddemodulationplaysaimportantroleininformationtransfer.MathWorkcompaniesintheUnitedStateslaunchedinordertocombinetheresultofcomputersimulationofMATLABsoftwareisoneofthemostexcellenttechnologyapplicationsoftware.MATLABcanblendtheoryandpracticalapplication,toinspirethestudents'interestinlearning,toenhancethelearningeffect.
Inordertolearnmoreconvenientdigitalsignal,thecompanywillbeADhoclearningsoftwareMATLABasalearningbaseddigitalsignal,completelyfulluseofMATLABrapidoperation,text,graphics,soundandinteractivefeaturessuchasman-machineinterfacetoananalysisandasimulationofdigitalsignaltransmission.MATLABisapowerfulfunctionincombinationwithcomputerlanguagefortheintegrationofscience,ithasthenumericalcomputingandmakingimages,imagemanipulation,andsystemsimulation,etc.,anditstoolboxandfunctionisveryrich,hasastrongmatrixoperation,thegraphvisualization,thusprovidesagreatconvenienceinthecommunicationsystemofmodulationanddemodulationprocessanalysis.
Thispaperroughlyintroducestheconceptofdigitalsignalsystem,thenintroducessomemodesofthemodulationanddemodulation,andconcretelyintroducesthewayweoftenuseseveralkindsofmodulationanddemodulationanddigitalsignaltransmissionintheMATLABsimulation.Digitalsignalcourseisourprofessionalrequiredcourse.
KEYWORDS:
Digitalsignaltransmission,MATLAB,simulation
第一章前言
1.1课题背景
随着现代数字信号的快速发展，数字调制解调技术主要涉及在所有的工程技术领域中，数字信号处理技术的飞速成长，全球性信息性高速公路建设的浪潮也都正在进行当中。
今天的各种传输媒介，有线电视、卫星通信系统和微波中继等将会尽可能更多地选用了数字信号载波传输方式。
这些系统当中数字调制技术已经充分的在已经使用到了。
现在的人们都意识到：
人们的工作方式和生活方式随着数字信号的发展已经慢慢的改变，互联网和数字信号处理技术已经成为了影响上升的重要因素，以及成为了一个国家与地区、科学与文化发展的经济实力。
数字信号调制技术采取数字调制解调技术的加密和差错控制为了便于数字信号的整合，所以数字信号拥有了模拟信号不可替代的优势。
针对传统数字信号根据比较传统的意思就是信息的传输，在现代高度信息社会中，信息和数字信号俨然成为了当前我国社会的“命脉”。
通过社会广泛的传播与交流，信息作为这种需求性的资源，才能使其产生利用价值，增进社会成员之间的团队协作能力，推进社会生产力的急速发展，创造出了强大的经济效益。
随着数字通信的迅速发展，出现了许多新问题，对传输频带的限制和对传输质量的要求越来越高。
MATLAB作为一种跨平台的，简单和有效的数学语言的矩阵运算，为用户提供了非常直观和简单的程序开发环境。
数字信号载波的传输在仿真软件设计中选择采用了MATLAB作为仿真工具，使得结果更为直观清楚，对于现在学生的实践学习有很大的帮助。
1.2本课题主要内容
数字信号载波传输的调制方式有四种，分别是：
载波幅度调制（PAM/ASK），载波相位调制（PSK）、正交幅度调制（QAM）和载波频率调制（FSK）。
数字信号系统包含了很多内容，它们基本上可以用数字基带信号控制高频载波的三个参数里的振幅、频率和相位，从而获得数字带通载波信号。
然而美国的MathWorks公司研究了MATLAB，从而解决了这一问题。
本课题主要介绍数字信号载波调制的四种基本方式，并简单介绍数字通信仿真的过程。
第二章MATLAB介绍
2.1MATLAB简介
美国TheMathWorks公司出品的MATLAB是一款由具备商业价值的数学使用软件，它是MATrixLABoratory的缩写。
它已经发展称为一个非常灵活的计算系统，是用来解决数字信号传输中的各种技术问题。
它提供的数据分析，算法实现交互式开发环境下的应用于发展。
MATLAB具有广泛的功能定义库，适用于MATLAB编程语言。
MATLAB是把矩阵作为基本数据单元的一种计算机程序设计语言，它把具有高性能的数值运算和可视化这两种连接在一起，且提供了强大的内置函数和供数据分析，算法实现交互式开发环境下的应用与发展。
MATLAB能够方便地用于运算，也能够快速完成计算可视化，制作图像，建立和烘托场景，图象处理等分析处理工作。
现今已经有超过四十多个的MATLAB工具箱，其基本概括了数据获取，科学计算，控制系统设计与分析，数字信号处理，数字图像处理，财务分析、生物和遗传等专业范畴。
2.2特点
当前MATLAB作为世界顶级的数学信号载波传输的应用软件，它拥有着巨大的工程运算、算法钻研、工程制图、应用程序的建设和动态仿真等功能，从而为用户供给了十分明了的程序开发环境，MATLAB语言在MATLAB程序当中运行，且提供了一个极为普遍的预定义函数库，如此就使得技术工作变得简单高效。
MATLAB还拥有计算机强、很好的人机界面、编程效率极高、绘图功能强大和可扩展性强的特点。
丰富的功能应用工具箱给用户提供了大批简易适用的处理工具。
MATLAB的功能强大，也便于学生学习，且编程效率高。
2.3MATLAB界面
当MATLAB6.5.0程序启动的时候，MATLAB的桌面窗口出现。
MATLAB桌面结构如图2.3.1所示。
图2.1MATLAB桌面结构图
2.4MATLAB语言平台及基本语法
当MATLAB运行时，有许多不同类型的窗口出现在平台上，都有各自功能，即有一些用在接收命令，一些用在显示信息。
主要窗口有一个命令窗，工作空间，驱动平台，编辑和调试窗口，工作台窗口和数字编辑器，交互界面窗口，帮助导航浏览器图形窗口。
MATLAB拥有着高级的矩阵/阵列语言，它主要包括了控制语句，函数，数据结构，输入和输出及面向对象等特点。
MATLAB的语法结构：
M文件、程序流程控制、波形的绘制、数据的表示和运算、矩阵的表示和运算、数据的输入和输出、MATLAB程序编译步骤。
过程中需要继续不断运行和调试程序编辑文件。
2.5数据分析
MATLAB软件可以实现数据的分析和处理，是一种高性能的数值计算和可视化软件。
MATLAB是基本单元矩阵的数据操作，而不需要预先指定矩阵的维数（动态定维）。
MATLAB数值运算中有变量，数值，矩阵，阵列，符号运算、注释与标点。
MATLAB能够从其余的文献、应用程序和外部设备访问数据。
MATLAB中数值运算有以下几个相关概念：
1.变量。
变量命名原则：
（1）是以字母开头、后面可以跟字母、数字和下划线
（2）长度不超过63个字符
（3）字母区分大小写
2.数值。
MATLAB的数值选用了十进制数表示法，能够带小数点或符号，也判别复数。
MATLAB中的数值默认是双精度实数。
3.矩阵。
矩阵是MATLAB进行数据处理和运算的根本元素，建立矩阵方式有：
利用函数和M文件来建立。
4.数组。
在MATLAB中，数组是一个非常重要的概念，数组的要求有：
图形窗口的掌握、二维和三维绘图的使用等。
5.符号运算。
为用户提供了一个丰富的MATLAB内置函数。
有各种各样的MATLAB符号运算方式：
符号表达式的运算及复合、化简，符号微分方程求解等等。
第三章数字信号载波传输
3.1概述
1.数字信号的载波传输
数字信号的载波传输也被称为数字调制系统，数字基带信号的载波信息的复制处理过程。
在数字通信系统中，在数字通信系统中，数字测图需要在信道传输信号的输入序列，输出数字信号的信号映射后，适用于传输数字调制信号。
按照数字调制的方法能够分为基带传输和带通传输两种方式。
数字调制是用数字基带信号对载波的某些参量进行控制，使载波随基带信号变化；数字调制解调技术是在接收端已经搬到了基带信道铜频谱还原为基带信号的一种过程。
对数字信号传性的调制，可以方便于信号到数字信号的传输；信道复用的实现；信号变化占用的带宽；提高了系统的性能。
为了使传输信道中的各种元素，最好使用可组合成不同的调制，并能清楚地描述和表达它的数学模型。
数字调制技术是2ASK，2FSK，2PSK，QAM，数字通信的传输效率，可以实现且远比模拟通信，数字调制技术的种类也远远大于模拟通信，大大提高了用户按照现实情况使用需求来挑选系统配置的灵活性。
数字调制系统结构图如图3.1所示，它主要包括发送滤波器，调制器，信道，解调器和抽样判决器等部分组成。
图3.1数字调制系统结构图
2.为什进行频带调制？
首先，由于有限的频率资源有限，导致了我们不能在开放信道使用传输信息；基带传输损耗大、易误码，且方便使用不同类型的模拟信道资源传输数字信号；它不仅用于低速数据传输，但也为高速数据传输；在理论上来讲，无论是哪种线性载波传输系统都能够即是基带传输系统。
调制来实现多路复用，完成频率偏移和减少噪声的干扰等。
3.怎样进行频带调制？
数字调制的特点是信息传递的一些离散的载波信号。
在通信系统中，当一个正弦载波的振幅，频率和相位3参数，相应的也有3种基本类型的振幅调制，频率调制和相位调制。
频带调制的基本的四种方法有：
载波幅度调制（PAM）、载波相位调制（PSK）、正交幅度调制（QAM）及载波频率调制（FSK）。
第四章数字信号载波传输的几种方式
4.1载波幅度调制（ASK）
振幅键控是用于载波传输数字信息的变化，频率和初始相位稳定性。
载波幅度调制（ASK）在模拟信号中的振幅、频率和相位等作为常量，以振幅作为变量。
在2ASK中，在相应的二进制信息的变化，只有两个载流子速率0或1。
其信号表达式为（3-1）：
e0（t）=S（t）cos
ct式（4-1）
其中，e0（t）是一种脉冲，它的形状决定了传输信号的谱特性。
4.1.1ASK调制原理
ASK调制的调制信号（基带信号）过程中的的一些的转换规则来改变载波。
2ASK信号有着模拟调制法（相乘器法）和键控法这两种方式。
模拟调制法就是用基带信号与载波相乘，进而把基带信号调制到载波上进行传输。
键控法主要由s（t）来控制电路的开关继而来进行调制。
S（t）
e0（t）
cos
ct
图4.2模拟调制法（相乘器法）
4.1.2数字ASK相干解调方式
2ASK有相干解调法（同步检测法）和非相干解调法（包络检波法）这两种基本解调方法，它是将搬移到高频段的频谱再搬回到原点。
本次设计采用相干解调法。
相干解调就是同步检测法，同步检测器用于对载波抑制的双边带或单边带进行调制。
y1（t）y（t）s（t）v（t）s’（t）
定时
脉冲
图4.32ASK非相干解调流程图
输出
e2ASK（t）
cos
ct定时
脉冲
图4.42ASK相干解调法（同步检测法）
程序
clear
n=1:
8192;
m=1:
128;x（n）=randint（1,8192,2）;x=x（n）';
y（n）=zeros（1,8192）;z（m）=zeros（1,128）;
forn=1:
8192
form=1:
128
ifn==64*m-63
z（m）=x（n）;
ifm==ceil（n/64）
y（（（64*m-63）:
（64*m））'）=z（m）;
end
end
end
end
n=1:
8192;rm2=y（n）;
x2=rm2;%产生基带信号
n=1:
（2^13）;
x1=cos（n.*1e9*2*pi/4e9）;%载频1GHz
x=x1.*x2;%ASK频带调制
b=blackman（2^13）;%窗函数
X=b'.*x;%ASK频带调制加窗
x3=[ones（1,64）zeros（1,8128）];%独个基带信号码元
y1=X（1:
（2^13））;y4=x1.*x3;
Y1=fft（y1,（2^13））;magY1=abs（Y1（1:
1:
（2^12）+1））/（200）;%求FFT
Y4=fft（y4,（2^13））;magY4=abs（Y4（1:
1:
（2^12）+1））/（37）;%求FFT
k1=0:
（2^12）;w1=（2*pi/（2^13））*k1;
u=（2*w1/pi）*1e9;
figure
（1）
subplot（2,1,1）;
plot（u,magY1,'b',u,magY4,'r'）;grid
title（'ASKr'）;axis（[4e8,1.6e9,0,1.1]）
X2=b'.*x2;%基带信号加窗
y2=X2（1:
（2^13））;
Y2=fft（y2,（2^13））;magY2=abs（Y2（1:
1:
（2^12）+1））/（200）+eps;%求FFT
k1=0:
（2^12）;w1=（2*pi/（2^13））*k1;
u=（2*w1/pi）*1e9;
Y3=fft（x3,（2^13））;magY3=abs（Y3（1:
1:
（2^12）+1））/（35）+eps;
subplot（2,1,2）;semilogy（u,magY2,'b',u,magY3,'r'）;grid
title（'ASKr-modulation'）;axis（[0,1.2e9,3e-2,3]）
figure
（2）
subplot（2,1,1）;plot（n,x2）;title（'ASKr'）;
axis（[0,640,-0.2,1.2]）;grid
subplot（2,1,2）;plot（n,x）;
axis（[0,640,-1.2,1.2]）;grid
4.2基带移相键控（PSK）
基带移相键控（PSK）是利用载波的相位变化来传递数字信息，而振幅和频率保持稳定。
2PSK流程图:
S（t）
cos
ctn（t）cos
ct
图4.52PSK调制解调框图
4.2.1PSK调制原理
二进制移相键控（2PSK）信号的产生是在二进制数字调制中，当正弦载波的相位随着二进制数字基带信号离散的变化而变化的过程。
2PSK信号调制方式有模拟调制法和键控法这两种方法。
在PSK中，一般用初始相位0或者∏分别代表二进制1或者0。
载波相位变化的信号是以2PSK作为参考的基准，当基带信号为0时相位相对于初始相位为0°，当基带信号为1时相对于初始相位为180°。
键控法，是以载波的相位进行二进制信息的调制方式。
4.2.2数字PSK的相干解调法
PSK信号通常采用相干解调。
统一的参考相位解调载波PSK信号的参考相位和相干载波。
通过带通滤波的信号在相乘器中和本地载波相乘，而后用低通滤波器滤除掉了高频分量，再进行抽样判决。
判决器根据极性的基础上来进行判决的。
这是正抽样值判为1，负抽样值判为0。
图4.62PSK的相干解调原理图
4.2.3二进制差分相移键控（DPSK）
2DPSK流程图：
ctn（t）
ct
图4.72DPSK调制解调图
（1）DPSK调制原理
达到二进制差分相移键控的最通常的方式有：
首先对二进制数字基带信号进行差分编码，继而对变换出的差分码进行绝对调相便可。
（2）DPSK解调原理
2DPSK信号解调有相干解调方法和差分相干解调。
相干解调码转换法和相
干解调法的解调原理是：
首先对2DPSK信号进行了相干解调，恢复出相对应的码，之后通过码反变换器转换成为绝对码，以此来恢复出了发送出去的二进制数字信息。
2DPSK数据
信号输出
本地载波
图4.82DPSK差分相干解调原理图
acd
e
b
定时脉冲
图4.92DPSK相干解调原理图
由于DPSK系统的差分相干解调电路简单的多，因此2DPSK系统中大都采用差分相干解调。
PSK调制举例：
MATLAB实现如下（函数bpskdigital.m）
functionbpskdigital（s,f）
%本程序实现BPSK调制
%s-输出二进制序列，f-载波信号的频率
%调用举例：
（f0和f1必须是整数）bpskdigital（[10110],2）
t=0:
2*pi/99:
2*pi;%初始定义
cp=[];mod=[];mod1=[];bit=[];
forn=1:
length（s）;%调制过程
ifs（n）==0;
cp1=-ones（1,100）;
bit1=zeros（1,100）;
elses（n）==1;
cp1=ones（1,100）;
bit1=ones（1,100）;
end
c=sin（f*t）;
cp=[cpcp1];
mod=[modc];
bit=[bitbit1];
end
bpsk=cp.*mod;
subplot（2,1,1）;
plot（bit,'LineWidth',1.5）;gridon;%分别画出原信号，已调信号示意
ylabel（'BinarySignal'）;
axis（[0100*length（s）-2.52.5]）;
subplot（2,1,2）;
plot（bpsk,'LineWidth',1.5）;gridon;
ylabel（'BPSKmodulation'）;
axis（[0100*length（s）-2.52.5]）;
图4.10BPSK波形
DPSK举例程序：
其仿真MATLAB程序实现如下。
NRindB1=0:
2:
12;
SNRindB2=0:
0.1:
12;
fori=1:
length（SNRindB1）,
simu_err_prb（i）=snr2ps（SNRindB1（i））;
end;
fori=1:
length（SNRindB2）,
SNR=exp（SNRindB2（i）*log（10）/10）;
theo_err_prb（i）=2*Qfunct（sqrt（SNR））;
end;
semilogy（SNRindB1,smld_err_prb,'*'）;
hold
semilogy（SNRindB2,theo_err_prb）;
xlabel（'Eb/No（dB）'）;
ylabel（'误码率'）;
legend（'仿真符号误码率','理论符号误码率'）;
function[p]=snr2ps（snr_in_dB）
N=10000;
Es=1;
snr=10^（snr_in_dB/10）;
sgma=sqrt（Es/（4*snr））;
fori=1:
2*N,
temp=rand;
if（temp<0.5）,dsource（i）=0;
else
dsource（i）=1;
end;
end;
mapping=[0132];
M=4;
[diff_enc_output]=cm_dpske（E,M,mapping,dsource）;
fori=1:
N,[n
（1）n
（2）]=gausamp（sgma）;
r（i,:
）=diff_enc_output（i,:
）+n;
end;
numoferr=0;
prev_theta=0;
fori=1:
N,theta=angle（r（i,1）+j*r（i,2））;
if（（delta_theta
升级会员 