FM调制解调系统设计与仿真linWord文档下载推荐.docx
《FM调制解调系统设计与仿真linWord文档下载推荐.docx》由会员分享,可在线阅读,更多相关《FM调制解调系统设计与仿真linWord文档下载推荐.docx(26页珍藏版)》请在冰豆网上搜索。
2.3FM调制仿真结果(6
2.4FM解调模型的建立(6
2.5解调过程分析(7
2.6高斯白噪声信道特性(8
2.7信噪比分析(9
2.8调频系统的抗噪声性能分析(10
三.仿真实现(12
3.1MATLAB源代码(12
3.2MATLAB仿真结果及分析(12
四.心得体会(14
五.参考文献(14
摘要
调制在通信系统中有十分重要的作用。
通过调制,不仅可以进行频谱搬移,把调制信号的频谱搬移到所希望的位置上,从而将调制信号转换成适合于传播的已调信号,而且它对系统的传输有效性和传输的可靠性有着很大的影响,调制方式往往决定了一个通信系统的性能。
FM信号的调制属于频谱的非线性搬移,它的解调也有相干和非相干解调两种方式。
本课程设计使用的仿真软件为MATLAB,利用MATLAB集成环境下的M文件,编写程序来实现FM信号的仿真分析,并分别绘制出基带信号、载波信号、已调信号的时域波形和频域波形;
再进一步分别绘制出对已调信号叠加噪声后信号、同步解调前信号和解调后基带信号的时域波形;
最后绘出FM基带信号调制和解调系统后的输入输出信噪比的关系,并通过与理论结果波形对比来分析该仿真调制与解调过程的正确性及噪声对FM信号解调的影响。
在课程设计中,系统开发平台为Windows7,在该平台运行MATLAB程序完成了对FM调制和解调以及对叠加噪声后解调结果的观察和分析。
关键词:
FM调制解调MATLAB噪声
Abstract
Modulationhasaveryimportantroleinthecommunicationsystem.Modulation,canbeusedtonotonlymovethespectrum,thespectrumofthemodulatedsignalismovedtoadesiredpositionsoastoconvertthemodulatedsignalsuitablefortransmissionofmodulatedsignals,butalsoitsignificantlyaffectthevalidityofthetransmissionsystemandthetransmissionreliability.Oftendeterminestheperformanceofthemodulationschemeofacommunicationsystem.FMmodulatedsignalisNon-linearmovementofspectrum.italsohasthecoherentdemodulationandnon-coherentdemodulationoftwoways.
SimulationSoftwareThiscourseisdesignedforuseasMATLAB,theuseofM-filesMATLABintegratedenvironment,writeaprogramtoachievethesimulationofFMsignals,andweredrawnoutofthebasebandsignal,thecarriersignal,modulatedsignaldomainwaveformandfrequencydomainwaveforms;
weredrawnfurthermodulatedsignalsuperimposedonthenoisesignal,asynchronizationsignalanddemodulatesthefirsttimedomainwaveformofthedemodulatedbasebandsignal;
finaloutputSNRplotFMinputbasebandsignalafterthemodulationanddemodulationsystemofrelationsand,throughcomparisonwiththetheoreticalresultstoanalyzewaveformsandnoiseaffectthecorrectnessofthesimulationprocessofmodulationanddemodulationoftheFMsignaldemodulation.Inthecoursedesign,systemdevelopmentplatformforWindows7,theplatformrunningMATLABprogramcompletedandFMmodulationanddemodulationnoisesuperimposedonthedemodulationresultsafterobservationandanalysis.
Keywords:
FMModulationDemodulationMATLABNoise
一.通信系统简介
通信的目的是传输信息。
通信系统的作用就是将信息从信息源发送到一个或多个目的地。
对于任何一个通信系统,均可视为由发送端、信道和接收端三大部分组成(如图1所示。
图1通信系统一般模型
信息源(简称信源的作用是把各种信息转换成原始信号。
发送设备的作用产生适合传输的信号,即使发送信号的特性和信道特性相匹配,具有抗噪声的能力,并且具有足够的功率满足原距离传输的需求。
信息源和发送设备统称为发送端。
发送端将信息直接转换得到的较低频率的原始电信号称为基带信号。
通常基带信号不宜直接在信道中传输。
因此,在通信系统的发送端需将基带信号的频谱搬移(调制到适合信道传输的频率范围内进行传输。
这就是调制的过程。
信号通过信道传输后,具有将信号放大和反变换功能的接收端将已调制的信号搬移(解调到原来的频率范围,这就是解调的过程。
信号在信道中传输的过程总会受到噪声的干扰,通信系统中没有传输信号时也有噪声,噪声永远存在于通信系统中。
噪声对于信号的传输是有害的,它能使模拟信号失真。
在本仿真的过程中我们假设信道为高斯白噪声信道。
调制在通信系统中具有十分重要的作用。
一方面,通过调制可以把基带信号的频谱搬移到所希望的位置上去,从而将调制信号转换成适合于信道传输或便于信道多路复用的已调信号。
另一方面,通过调制可以提高信号通过信道传输时的抗干扰能力,同时,它还和传输效率有关。
具体地讲,不同的调制方式产生的已调信号的带宽不同,因此调制影响传输带宽的利用率。
可见,调制方式往往决定一个通信系统的性能。
在本仿真的过程中我们选择用调频调制方法进行调制。
信息源发送设备信道
接受设备信息源噪声源发送端
接收端信道
调制过程是一个频谱搬移的过程,它是将低频信号的频谱搬移到载频位置。
而解调是将位于载频的信号频谱再搬回来,并且不失真地恢复出原始基带信号。
在本仿真的过程中我们选择用相干解调方法进行解调。
角度调制及解调电路不同于频谱线性搬移电路。
它是用低频信号去调制高频振荡的相角,或是从已调波中解出调制信号所进行的频谱变换,这种变换不是线性变换,而是非线性变换。
因此,我们把角度调制及调角波的解调电路称为频谱非线性变换电路。
FM在通信系统中的使用非常广泛。
FM广泛应用于高保真音乐广播、电视伴音信号的传输、卫星通信和蜂窝电话系统等。
本次课设中使用功能强大的MATLAB仿真软件对FM信号进行仿真分析。
6max
(Kπ<
<
⎰dttmf二.FM调制解调系统设计
2.1MATBLAB简介
MATLAB的名称源自MatrixLaboratory,它的首创者是在数值线性代数领域颇有影响的CleveMoler博士,他也是生产经营MATLAB产品的美国Mathworks公司的创始人之一。
MATLAB是一种科学计算软件,专门以矩阵的形式处理数据。
MATLAB将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而使其被广泛地应用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作中。
MATLAB与其它计算机高级语言如C,C++等相比,MATLAB语言编程要简洁得多,编程语句更加接近数学描述,可读性好,其强大的图形功能和可视化数据处理能力也是其它高级语言望尘莫及的。
对于具有任何一门高级语言基础的学生来说,学习MATLAB十分容易。
MATLAB使得人们摆脱了常规计算机编程的繁琐,让人们能够将大部分精力投入到研究问题的数学建模上。
可以说,应用MATLAB这一数学计算和系统仿真的强大工具,可以使科学研究的效率得以成百倍的提高。
2.2FM调制模型的建立
从频率调制的相位与频率关系可以看出,调频信号可通过直接调频和间接调频两种方法得到,所谓间接调频就是先对调制信号积分再调相而得到。
同样,调相信号也可以通过直接调相和间接调相两种方法得到,间接调相就是先对调制信号进行微分再进行频率调制。
根据调制后已调信号的瞬时相位偏移的大小,可将角度调制分为宽带调制(宽带调频和宽带调相和窄带调制(窄带调频和窄带调相。
如果调频信号或调相信号的最大瞬时相位偏移保持在很小的范围内,一般小于30°
即满足条件:
时,则称为窄带调频或窄带调相。
当上述条件不满足时,就称为宽带调频或宽带调相。
图2FM间接调制模型
其中,m(t为基带调制信号,设调制信号为
m(t=cos(2*pi*fm*t
设正弦载波为
c(t=cos(2*pi*fc*t
在调制时,调制信号的频率去控制载波的频率的变化,载波的瞬时频偏随调制信号m(t成正比例变化,即:
dttmkttt
f⎰+=0
0((ωω
式中,fK为调频灵敏度。
这时相位偏移为
f((0
0⎰+=ωθ
则可得到调频信号为
(cos((u0
0FMdttmktUtt
fc⎰+=ω
FM的频谱的计算:
可以看出FM的频谱与
(fnmJ的值有关。
m
fmωω∆=
∑∞
1
204202cos(2(4cos(22cos(2(sincos(+=Ω+=+Ω+Ω+=Ωnfnffffft
nmJmJtmJtmJmJtm
12(53112sin((25sin(23sin(2sin(2sinsin(+=+Ω+=+Ω+Ω+Ω=Ωnfnfffft
nmJtmJtmJtmJtm
t
nmJAtVktAtcnfnfcFM(cos(sincos(U_
Ω+=⎥⎦
⎤
⎢⎣
⎡
ΩΩ+=∑∞
+∞=ω
ω
其信号带宽为Ω+=1(2fmB调制信号产生的部分M文件如下:
fc=30000;
%载波频率ct=cos(2*pi*fc*t;
%载波信号am=1;
%调制信号幅度fm=300;
%调制信号频率mt_sin=am*cos(2*pi*fm*t;
%正弦波mt_tri=sawtooth(2*pi*fm*t,0.5;
%三角波mf=1;
%调制指数,可调
int_mt_sin(1=0;
%正弦波积分fori=1:
length(t-1
int_mt_sin(i+1=int_mt_sin(i+mt_sin(i*dt;
end
int_mt_tri(1=0;
fori=1:
length(t-1%三角波积分int_mt_tri(i+1=int_mt_tri(i+mt_tri(i*dt;
sfm_sin=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_sin;
sfm_tri=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_tri;
2.3FM调制仿真结果
图3FM调制
2.4FM解调模型的建立
由前述公式知道输入调频信号为
相干载波为
乘法器的作用是把调频信号变成有多种频率的波的混合,乘法器输出为
[]
2cos1((2
2sin21(tdttmKttscfcpωω-*+-=⎰经低通滤波器后取出器低频分量为
dttmKtsfd⎰=
(2
(在经过微分器,即得出解调出的基带信号:
(0tmKtmf=
相干解调可以恢复出原来的基带信号,而且要求本地载波与调制载波同步,否则会使解调信号失真。
2.5解调过程分析
在这次的解调设计中,采用MATLAB自带的解调函数demod对已调信号进行解调。
该函数使用形式如下所示:
Z=demod(Y,Fc,Fs,Method
其中,Y为输入的调频波函数,Fc、Fs分别为载波频率和抽样频率,而且要满足Fc2Fs,Method为调制方式。
FM解调的M文件部分代码如下:
y_sin=demod(sfm_sin,fc,Fs,'
fm'
;
y_tri=demod(sfm_tri,fc,Fs,'
yyn_sin=demod(yn_sin,fc,Fs,'
yyn_tri=demod(yn_tri,fc,Fs,'
通过M文件绘制出大、小信噪比解调的输出波形如下:
图4小信噪比时的FM解调
图5大信噪比时的FM解调
2.6高斯白噪声信道特性
设正弦波通过加性高斯白噪声信道后的信号为
(cos((crtAtntωθ=++
其中,白噪声(nt的取值的概率分布服从高斯分布。
MATLAB本身自带了标准高斯分布的内部函数awgn。
awgn函数产生的随机序列服从
均值为0m=,方差2
1σ=的高斯分布。
正弦波通过加性高斯白噪声信道后的信号为
故其有用信号功率为
2
AS=
噪声功率为
Nσ
=
信噪比SN满足公式
1010log(SBN
则可得到公式
22
10
210
BAσ=
∙
我们可以通过这个公式方便的设置高斯白噪声的方差。
2.7信噪比分析
在本仿真过程中,选择了1,3,5三种不同调制指数以示区别,并绘制输入输出信噪比关系曲线,如下图6所示。
mf=1
mf=3
mf=5
图6正弦波输入输出信噪比
2.8调频系统的抗噪声性能分析
从前面的分析可知,调频信号的解调有相干解调和非相干解调两种。
相干解调仅适用于窄带调频信号,且需同步信号;
而非相干解调适用于窄带和宽带调频信号,而且不需同步信号,因而是FM系统的主要解调方式,所以这里仅仅讨论非相干解调系统的抗噪声性
能,其分析模型如图所示。
图7调频系统抗噪声性能分析模型
图中带通滤波器的作用是抑制信号带宽以外的噪声。
(nt是均值为零,单边功率谱密度为0n的高斯白噪声,经过带通滤波器后变为窄带高斯噪声(int。
限幅器是为了消除接收信号在幅度上可能出现的畸变。
设调频信号为
(cos((t
FMcfStAtKmdωττ-∞
=+⎰
故其输入功率为
iAS=
输入噪声功率为
ioFMNnB=
因此输入信噪比为
2iiFM
SANB=在大信噪比条件下,信号和噪声的相互作用可以忽略,这时可以把信号和噪声分开来算,这里,我们可以得到解调器的输出信噪比
《武汉理工大学MATLAB课程设计报告》
222
23
3(8o
f
o
om
AKmtSNnfπ=
上式中,A为载波的振幅,fK为调频器灵敏度,mf为调制信号(mt的最高频率,on为
噪声单边功率谱密度。
我们如若考虑(mt为单一频率余弦波时的情况,可得到解调器的制度增益为
232Ao
FM
iom
iSNGmSnfN==
考虑在宽带调频时,信号带宽为
2(12(FMfmmBmfff==+∆+
则可以得到
23(1FMffGmm=+
可以看出,大信噪比时宽带调频系统的信噪比增益是很高的,它与调频指数的立方成正比。
可见,加大调频指数fm,可使调频系统的抗噪声性能迅速改善。
三.仿真实现
本次设计是在MATLAB平台上,在前述FM通信系统的分析基础上,进行程序代码的编写,从而得到相关仿真波形,
由此进行分析讨论。
图8程序流程图
3.1MATLAB源代码
根据设计要求所编写的源代码已在在附页中给出,详情请看附页I。
3.2MATLAB仿真结果及分析
在MATLAB平台上运行程序得到各仿真图,由于篇幅有限,已在附页II给出。
这里针对部分仿真结果进行分析。
(1系统生成满足条件的信号波形,包括正弦波和三角波。
(2在未加入高斯白噪声时,调制系统输出理想的FM已调信号,已调信号的频率随着原信号的幅度大小而变。
解调系统也能理想地解调出原始输入信号,解调后的波形光滑,与
FM调制程序
添加高斯白噪声程序
FM解调程序
绘制频谱图程序
显示程序
初始化程序
原信号一致。
(3在加入较弱的高斯白噪声时,调制系统输出较为理想的FM已调信号,已调信号的频率随着原信号的幅度大小而变但并不是理想的。
解调系统也能正确地解调出原始输入信号,解调后的波形稍带毛刺,与原信号基本一致。
(4在加入较强的高斯白噪声时,调制系统输出不太理想的FM已调信号,已调信号的频率与原信号的幅度大小关系不明显。
解调后的波形有失真现象。
(5输入输出信噪比呈正向线性关系,输入信噪比越大,输出信噪比也越大。
总述:
FM调制解调系统设计成功。
结果表明:
信道噪声很大程度上影响着调制解调系统性能,实际运用中应尽量降低信道噪声;
为了少失真地得到原始信号,在提高系统性能的基础上,还应增大输入信噪比。
四.心得体会
课程设计是培养学生综合运用所学知识,发现,提出,分析和解决实际问题,锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程.回顾起此次通信原理课程设计,至今我仍感慨颇多,的确,从选题到定稿,从理论到实践,在整整一星期的日子里,可以说得是苦多于甜,但是可以学到很多很多的的东西,同时不仅可以巩固了以前所学过的知识,而且学到了很多在书本上所没有学到过的知识。
通过这次课程设计使我懂得了理论与实际相结合是很重要的,只有理论知识是远远不够的,只有把所学的理论知识与实践相结合起来,从理论中得出结论,才能真正为社会服务,从而提高自己的实际动手能力和独立思考的能力。
在设计的过程中遇到问题,可以说得是困难重重,这毕竟第一次做的,难免会遇到过各种各样的问题,同时在设计的过程中发现了自己的不足之处,对以前所学过的知识理解得不够深刻,掌握得不够牢固……通过这次课程设计之后,一定把以前所学过的知识重新温故。
这次课程设计终于顺利完成了,在设计中遇到了很多问题,最后在不懈的努力下终于迎刃而解。
感谢这次旅程,我又成长不少!
五.参考文献
[1]《通信原理(第六版》樊昌信等。
北京:
国防工业出版社。
[2]《MATLAB7.0在数字信号处理中的应用》罗军辉等。
机械工业出版社。
[3]《MATLAB程序设计教程》刘卫国等。
中国水利水电出版社。
附页I
仿真波形图
附页II
源代码
%**********本程序主要实现FM调制与解调系统设计*************%**********以正弦波和三角波进行设计**********************
%**********并考虑信道中高斯白噪声的影响******************
%************初始化***********************
dt=0.00001;
%积分步长
N=1000;
%采样点数
Fs=1/dt;
%抽样频率Fs=100000hz,满足抽样定理
n=0:
N-1;
t=n/Fs;
%***************信号设定***********************
%载波频率
ct=cos(2*pi*fc*t;
%载波信号
am=1;
%调制信号幅度
fm=300;
%调制信号频率
mt_sin=am*cos(2*pi*fm*t;
%正弦波
mt_tri=sawtooth(2*pi*fm*t,0.5;
%三角波
mf=3;
%调制指数
%******************调制模块*************************
%正弦波积分
end
length(t-1%三角波积分
int_mt_tri(i+1=int_mt_tri(i+mt_tri(i*dt;
sfm_tri=am*cos(2*pi*fc*t+mf*2*pi*fm*int_mt_tri;
%**********