通信系统建模与仿真.docx
- 文档编号:11440435
- 上传时间:2023-03-01
- 格式:DOCX
- 页数:25
- 大小:1.66MB
通信系统建模与仿真.docx
《通信系统建模与仿真.docx》由会员分享,可在线阅读,更多相关《通信系统建模与仿真.docx(25页珍藏版)》请在冰豆网上搜索。
通信系统建模与仿真通信系统建模与仿真电子信息系统仿真课程设计电子信息系统仿真课程设计级电子信息工程专业班级题目FM调制解调系统设计与仿真姓名学号指导教师胡娟二0年月日内容摘要内容摘要频率调制(FM)通常应用通信系统中。
FM广泛应用于高保真音乐广播、电视伴音信号的传输、卫星通信和蜂窝电话系统等。
FM调制解调系统设计是对模拟通信系统主要原理和技术进行研究,理解FM系统调制解调的基本过程和相关知识,利用MATLABS成环境下的M文件,编写程序来实现FM调制与解调过程,并分别绘制出基带信号,载波信号,已调信号的时域波形;再进一步分别绘制出对已调信号叠加噪声后信号,非相干解调后信号和解调基带信号的时域波形;最后绘出FM基带信号通过上述信道和调制和解调系统后的误码率与信噪比的关系,并通过与理论结果波形对比来分析该仿真调制与解调系统的正确性及噪声对信号解调的影响。
在课程设计中,系统开发平台为WindowsXP使用工具软件为7.0。
在该平台运行程序完成了对FM调制和解调以及对叠加噪声后解调结果的观察。
通过该课程设计,达到了实现FM信号通过噪声信道,调制和解调系统的仿真目的。
了解FM调制解调系统的优点和缺点,对以后实际需要有很好的理论基础。
关键词关键词FM解调;调制;MKTLAB仿真;抗噪性一、一、MATLAB软件简介软件简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。
其特点是:
(1)可扩展性:
Matlab最重要的特点是易于扩展,它允许用户自行建立指定功能的M文件。
对于一个从事特定领域的工程师来说,不仅可利用Matlab所提供的函数及基本工具箱函数,还可方便地构造出专用的函数。
从而大大扩展了其应用范围。
当前支持Matlab的商用Toolbox(工具箱)有数百种之多。
而由个人开发的Toolbox则不可计数。
(2)易学易用性:
Matlab不需要用户有高深的数学知识和程序设计能力,不需要用户深刻了解算法及编程技巧。
(3)高效性:
Matlab语句功能十分强大,一条语句可完成十分复杂的任务。
如fft语句可完成对指定数据的快速傅里叶变换,这相当于上百条C语言语句的功能。
它大大加快了工程技术人员从事软件开发的效率。
据MathWorks公司声称,Matlab软件中所包含的Matlab源代码相当于70万行C代码。
二、理论分析二、理论分析2.1一般通信系统通信的目的是传输信息。
一般通信系统的作用就是将信息从信息源发送到一个或多个目的地。
对于任何一个通信系统,均可视为由发送端、信道和接收端三大部分组成。
图2.1通信系统一般模型2.2频率调制(FM)原理频率(FM)调制的名称源于m(t)与已调信号的频率呈线性关系FM调制就是将调制信号的变化映射到已调信号的频率大小。
加(f)SfM(f)FM图2.2FM调制原理图其中,m(t)为基带调制信号,设调制信号为m(t)二二Acos(2fnt)设正弦载波为c(t)二二C0S(2二二fct)可得到已调调频信号为Sfm(tpAccos2fct2kFMm(t)dt信号传输信道为高斯白噪声信道,其功率为二2。
2.3调制过程分析产生调频信号有直接调频法和间接调频法。
1、直接调频法aVCO吊応吊応0图2.3FM直接调频法模型2、间接调频法枳甘器0戦波1a盯宜畑竖Acsin2?
i、*f何硏帘怖匪调制图2.4FM窄带与宽带角调信号产生框图在调制时,调制信号的频率去控制载波的频率的变化,载波的瞬时频偏随调制信号m(t)成正比例变化,即dtdtKfmmt式中:
Kfm调频灵敏度(rad/s-v)相位偏移为X必K-mds(t)=Acos2得到调频信号为Sfm调制信号产生的M文件:
dt=O.OO1;t=0:
dt:
1.5;am=15;fm=20;mt=am*cos(2*pi*fm*t);fc=50;ct=cos(2*pi*fc*t);kf=10;int_mt
(1)=0;fori=1:
length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;endsfm=am*cos(2*pi*fc*t+2*pi*kf*int_mt);fJKfm“dJ%设定时间步长%产生时间向量%设定调制信号幅度J可更改%设定调制信号频率J可更改%生成调制信号%设定载波频率-可更改%生成载波%设定调频指数%对mt进行积分%调制,产生已调信号调制信号的时域图1100.5时间t11.5图2.5FM调制时域波形图图2.6FM调制频域波形图2.4FM解调模型的建立调频信号通常采用非相干解调方式来接收。
鉴频器图2.7FM非相干解调原理图其中,限幅器及带通滤波器的作用是让信号完全通过而最大限度地抑制噪声。
FM信号应该是一个等幅振荡信号,所以其幅度上的波动是外来的,可用限幅器将其抑制,而BPF抑制信号频带外的噪声。
由微分器和包络检波器级联的核心单元称为鉴频器。
微分器把调频信号变成调幅调频波,然后由包络检波器检出包络,最后通过低通滤波器取出调制信号。
此外,调频解调还有调频负反馈解调方案和利用锁相环作调频解调器。
2.5解调过程分析设输入调频信号为St=Sfmt二Acos2fctKfmmd调频到调幅的变换1微分微分器的作用是把调频信号变成调幅调频波。
微分器输出为Sd(t)=a2兀f广2兀Kfmm(t)Sin.岸f+2兀KfmtcM片J调频到调幅的变换2:
使用调谐电路的上升频率特性(线性区)平衡鉴频器:
展宽频率特性的线性范围图2.8平衡鉴频器及相应的频响图检波器输出为Sot二Kd2二fc2二KfmmJKd称为鉴频灵敏度(%z),是已调信号单位频偏对应的调制信号的幅度,经低通滤波器后加隔直流电容,隔除无用的直流,得mo(t)二二KdKfm(t)微分器通过程序实现,代码如下:
fori=1:
length(t)-1%接受信号通过微分器处理diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);%hilbert变换,求绝对值得到瞬时幅度(包络检波)通过M文件绘制出两种不同信噪比解调的输出波形如下:
2.6高斯白噪声信道特性设正弦波通过加性高斯白噪声信道后的信号为r(t)二Acos(ct二)n(t)其中,白噪声n(t)的取值的概率分布服从高斯分布。
MATLAB本身自带了标准高斯分布的内部函数randn。
randn函数产生的随机序列服从均值为口=0,方差二1的高斯分布。
正弦波通过加性高斯白噪声信道后的信号为r(t)二二Acosc
(2)nt(故其有用信号功率为2噪声功率为通过这个公式可以方便的设置高斯白噪声的方差。
我们选择了10db和30db两种不同信噪比以示区别,其时域图如图2.7和图2.8。
20:
1c;0鈕曲讪诃込也胡碟-201:
00.511.5时间t图2.10含小信噪比高斯白噪声已调信号的时域图图2.11含小信噪比高斯白噪声已调信号的时域图2.7调频系统的抗噪声性能分析图中带通滤波器的作用是抑制信号带宽以外的噪声。
n(t)是均值为零,单边功率谱密度为n0的高斯白噪声,经过带通滤波器后变为窄带高斯噪声ni(t)。
限幅器是为了消除接收信号在幅度上可能出现的畸变。
设调频信号为Sfmti=Acos2ftK怖mdSi=故其输入功率为输入噪声功率为M=nBFMA因此输入信噪比为Ni2Bfm在大信噪比条件下,信号和噪声的相互作用可以忽略,这时可以把信号和噪声分开来算,这里,我们可以得到解调器的输出信噪比So3A2K2m2(t)No82nof:
上式中,A为载波的振幅,Kf为调频器灵敏度,fm为调制信号m(t)的最高频率,n。
为噪声单边功率谱密度。
我们如若考虑m(t)为单一频率余弦波时的情况,可得到解调器的制度增益为考虑在宽带调频时,信号带宽为Bm(FM1B则可以得到则可以得到Gm3FMm(t)/m(t)max可以看出,大信噪比时宽带调频系统的信噪比增益是很高的,它与调P频指数的立方成正比。
可见,加大调频指数FM,可使调频系统的抗噪声性能迅速改善。
但加大调频指数FM,又会增大信号占用带宽Bfm,所以说FM调制的大信噪比是用带宽来换取的。
低信噪比下,解调输出信号与噪声相混合,无法区分,造成输出信噪比急剧恶化T门限效应。
cfNofPn(f)_2在调频解调中,解调输出有一抛物线噪声谱0A图2.13具有预加重和去加重滤波的调频系统f-W解调输出噪声功率在-WW频率范围内的噪声功率Pno这意味着在调频输出中,基带信号的高频分量比低频分量受到噪声的干扰更严重,为此,在实际系统中,采用预加重和去加重技术,以减小此抛物线噪声谱的影响,使系统总的输出信噪比得以改善。
设计一系统,使它对于基带信号低频段呈现为一对调频器及解调器,对于基带信号高频段呈现为一对调相及解调器。
即在发端将一高通滤波器与一调频器级联,该高通滤波器在低频有意恒定增益的频率特性,在高频有一近似于kf的微分器的频率特性。
收端的解调器是由普通鉴频器与一低通滤波器相级联组成。
图3.1程序流程图模型3.1仿真结果仿真结果图3.2原调制信号、载波信号和已调信号的曲线图图3.3原调制信号和已调信号在频域内的图图3.4原调制信号、无噪声条件下已调信号和解调信号的曲线图图3.5调制信号、小信噪比高斯白噪声下已调信号和解调信号的曲线图图3.6调制信号、大信噪比高斯白噪声下已调信号和解调信号的曲线图3.2MATLAB源代码源代码%FM调制解调系统.m%频率调制与解调的Matlab演示源程序%可以任意改原调制信号函数m(t)%*echooffcloseallclearallclc%*%*FM调制*%设定调制信号幅度%设定调制信号频率%生成调制信号%设定载波频率%设定调频指数%求信号m(t)的积分%调制,产生已调信号%调制信号dt=0.001;t=0:
dt:
2.0;ac=5;fm1=5;mt=ac*cos(2*pi*fm1*t);fc=50;ct=cos(2*pi*fc*t);kf=10;int_mt
(1)=0;fori=1:
length(t)-1int_mt(i+1)=int_mt(i)+mt(i)*dt;endsfm=ac*cos(2*pi*fc*t+2*pi*kf*int_mt);%*%*添加高斯白噪声*n=sqrt(db)*randn(size(t);nsfm=n+sfm;号(信号通%过信道传输)%*sn1=10;sn2=30;sn=0;db=acA2/(2*(10A(sn/10);差%设信躁比(小信噪比)%设信躁比(大信噪比)%设信躁比(无信噪比)%计算对应的高斯白躁声的方%生成高斯白躁声%生成含高斯白躁声的已调信%*FM解调*%接收信号通过微分器处理fori=1:
length(t)-1diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);(包络检波)%hilbert变换,求绝对值得到瞬时幅度zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;%*%*%*时域到频域转换ts=0.001;fs=1/ts;df=0.25;傅里叶变换%时,它表示FFT的最小频率间隔%*对调制信号m(t)求傅里叶变换*m=ac*cos(2*pi*fm1*t);fs=1/ts;ifnargin=2n1=0;elsen1=fs/df;endn2=length(m);n=2A(max(nextpow2(n1),nextpow2(n2);M=fft(m,n);m=m,zeros(1,n-n2);df1=fs/n;u求傅里变换M=M/fs;察f=0:
df1:
df1*(length(m)-1)-fs/2;*%抽样间隔%抽样频率%所需的频率分辨率,用在求%原调信号%以上程序是对调制后的信号%缩放,便于在频谱图上整体观%时间向量对应的频率向量%*对已调信号u求傅里变换*fs=1/ts;ifnargin=2n1=0;elsen1=fs/df;endn2=length(sfm);n=2A(max(nextpow2(n1),nextpow2(n2);U=fft(sfm,n);u=sfm,zeros(1,n-n2);df1=fs/n;变换%以上是对已调信号u求傅里%缩放U=U/fs;%*%*%*显示程序*disp(按任意键可以看到原调制信号、载波信号和已调信号的曲线)pause*figure
(1)*figure
(1)%*disp(按任意键可以看到原调制信号和已调信号在频域内的图形)pause%*figure
(2)*figure
(2)%fftshift:
将FFT中的DC分量移subplot(2,1,1)plot(f,abs(fftshift(M)到频谱中心xlabel(频率f)title(原调制信号的频谱图)subplot(2,1,2)plot(f,abs(fftshift(U)xlabel(频率f)title(已调信号的频谱图)%*disp(按任意键可以看到原调制信号、无噪声条件下已调信号和解调信号的曲线)pause%*figure(3)*figure(3)subplot(3,1,1);plot(t,mt);%绘制调制信号的时域图xlabel(时间t);title(调制信号的时域图);subplot(3,1,2);plot(t,sfm);%绘制已调信号的时域图xlabel(时间t);title(无噪声条件下已调信号的时域图);nsfm=sfm;fori=1:
length(t)-1%接受信号通过微分器处理diff_nsfm(i)=(nsfm(i+1)-nsfm(i)./dt;enddiff_nsfmn=abs(hilbert(diff_nsfm);%hilbert变换,求绝对值得到瞬时幅度(包络检波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn-zero;subplot(3,1,3);%绘制无噪声条件下解调信号的时域图plot(1:
length(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel(时间t);title(无噪声条件下解调信号的时域图);%*disp(按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线)pausefigure(4)diff_nsfmn1=abs(hilbert(diff_nsfm1);%hilbert变换,求绝对值得到瞬时幅度(包络检波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:
length(diff_nsfm),diff_nsfm);%绘制含小信噪比高斯白噪声已调信号的时域图xlabel(时间t);title(含小信噪比高斯白噪声已调信号的时域图);subplot(3,1,3);%绘制含小信噪比高斯白噪声解调信号的时域图plot(1:
length(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel(时间t);title(含小信噪比高斯白噪声解调信号的时域图);%*disp(按任意键可以看到原调制信号、大信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线)pause%*figure(5)*figure(5)subplot(3,1,1);plot(t,mt);xlabel(时间t);title(调制信号的时域图);db1=acA2/(2*(10A(sn2/10);躁声的方差n1=sqrt(db1)*randn(size(t);nsfm1=n1+sfm;号(信号通过信道传输)fori=1:
length(t)-1diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i)./dt;enddiff_nsfmn1=abs(hilbert(diff_nsfm1);瞬时幅度(包%络检波)zero=(max(diff_nsfmn)-min(diff_nsfmn)/2;diff_nsfmn1=diff_nsfmn1-zero;subplot(3,1,2);plot(1:
length(diff_nsfm1),diff_nsfm1);声已调信号%的时域图xlabel(时间t);title(含大信噪比高斯白噪声已调信号的时域图subplot(3,1,3);声解调信号%绘制调制信号的时域图%计算对应的大信噪比高斯白%生成高斯白躁声%生成含高斯白躁声的已调信%接受信号通过微分器处理%hilbert变换,求绝对值得到%绘制含大信噪比高斯白噪);%绘制含大信噪比高斯白噪%的时域图plot(1:
length(diff_nsfmn1)./1000,diff_nsfmn1./400,r);xlabel(时间t);title(含大信噪比高斯白噪声解调信号的时域图);四心得与体会五参考文献五参考文献1通信原理李晓峰等。
北京:
清华大学出版社。
2通信原理周炯磐等。
北京:
北京邮电大学出版社。
3MATLAB编程(第4版)StephenJ.Chapman著。
北京:
科学出版社。
4MATLAB在通信与电子中的应用徐明远著。
西安电子科技大学出版社。
5通信系统课程设计与实验教程雷菁主编。
科学出版社。
6MATLAB7.0在数字信号处理中的应用罗军辉编。
机械工业出版社。
指导教师评语:
课程设计成绩:
指导教师签名:
年月曰
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 通信 系统 建模 仿真