您好,欢迎访问三七文档
当前位置:首页 > 电子/通信 > 综合/其它 > 基于MATLAB模拟调制系统的仿真
11线性模拟调制1.1模拟调制原理模拟调制是指用来自信源的基带模拟信号去调制某个载波,而载波是一个确知的周期性波形。模拟调制可分为线性调制和非线性调制,本文主要研究线性调制。线性调制的原理模型如图1.1所示。设c(t)=Acos2tfo,调制信号为m(t),已调信号为s(t)。图1.1线性调制的远离模型调制信号m(t)和载波在乘法器中相乘的结果为:tAtmtswocos)()(',然后通过一个传输函数为H(f)的带通滤波器,得出已调信号为。从图1.1中可得已调信号的时域和频域表达式为:(1-1)式(1-1)中,M(f)为调制信号m(t)的频谱。由于调制信号m(t)和乘法器输出信号之间是线性关系,所以成为线性调制。带通滤波器H(f)可以有不同的设计,从而得到不同的调制种类。1.2双边带调制DSB的基本原理在幅度调制的一般模型中,若假设滤波器为全通网络,调制信号m(t)中无直流分量,则输出的已调信号就是无载波分量的双边带调制信号,或称抑制载波双边带(DSB)调制信号,简称双边带(DSB)信号。设正弦型载波c(t)=Acos(t),式中:A为载波幅度,为载波角频率。根据调制定义,幅度调制信号(已调信号)一般可表示为:(t)=Am(t)cos(t)(1-2))()]()([21)()(*]cos)([)(fHffMffMfsthttmtsooowm(t)H(t)Aostwocs(t))('ts2其中,m(t)为基带调制信号。设调制信号m(t)的频谱为M(),则由公式2-2不难得到已调信号(t)的频谱:)]()([2)(ccmMMAsωωωωω(1-3)由以上表示式可见,在波形上,幅度已调信号随基带信号的规律呈正比地变化;在频谱结构上,它的频谱完全是基带信号频谱在频域内的简单搬移。标准振幅就是常规双边带调制,简称调幅(AM)。假设调制信号m(t)的平均值为0,将其叠加一个直流偏量后与载波相乘,即可形成调幅信号。其时域表达式为:)cos())(()(0ttmtcAMAs(1-4)式中:为外加的直流分量;m(t)可以是确知信号,也可以是随机信号。若为确知信号,则AM信号的频谱为:(1-5)AM信号的频谱由载频分量、上边带、下边带三部分组成。AM信号的总功率包括载波功率和边带功率两部分。只有边带功率才与调制信号有关,也就是说,载波分量并不携带信息。因此,AM信号的功率利用率比较低。AM调制器模型如下图所示。图1.2AM调制器模型AM信号的时域和频域表达式分别为(1-6)(1-7)式中,Ao为外加的直流分量;m(t)可以是确知信号也可以是随机信号,但通常认为其平均值为0,即0)(tm—。由频谱可以看出,AM信号的频谱由载波分量、上边带、下边带三部分组成。上边带的频谱结构与原调制信号的频谱结构相同,下边带是上边带的镜像。因此,AM信号是带有载波分量的双边带信号,他的带宽是基带信号带宽的2倍,即)(cos)()(cos)(cos)]([)(twctmtwcAtwctmAotsoAM)]()([21)]()([)(wcwMwcwMwcwwcwAotsAM)]()([21)]()([)(0ccccmMMAsfHfBHAM23(1-8)AM调制典型波形和频谱如图1-1所示:图1.3AM调制典型波形和频谱如果在AM调制模型中将直流A去掉,即可得到一种高调制效率的调制方式—抑制载波双边带信号(DSB—SC),简称双边带信号。其时域表达式为)cos()()(ttmtcDSBs(1-9)式中,假设的平均值为0。DSB的频谱与AM的谱相近,只是没有了在±ω处的函数δ,即(1-10)与AM信号比较,因为不存在载波分量,DSB信号的调制效率是100,即全部效率都用于信息传输。但由于DSB信号的包络不再与调制信号的变化规律一致,因而不能采用简单的包络检波来恢复调制信号。DSB信号借条是需采用相干解调,也称同步检波(比包络检波器复杂得多)。)]()([21)(ccmMMs4其典型波形和频谱如图1-2所示:图1.4DSB调制典型波形和频谱1.3DSB解调原理与抗噪性能解调是调制的逆过程,其作用是从接收的已调信号中恢复原基带信号(即调制信号)。解调的方法可分为两类:相干解调和非相干解调(包络检波)。相干解调,也称同步检波,为了无失真地恢复原基带信号,接收端必须提供一个与接收的已调载波严格同步(同频同相)的本地载波(称为相干载波),它与接受的已调信号相乘后,经低通滤波器取出低频分量,即可得到原始的基带调制信号。包络检波器就是直接从已调波的幅度中提取原调制信号,通常由半波或全波整流器和低通滤波器组成。由于DSB信号的包络不再与调制信号的变化规律一致,因而不能采用简单的包络检波来恢复调制信号。DSB信号解调时需采用相干解调。DSB相干解调性能分析模型如图2.5所示:图1.5DSB相干解调性能分析模型设解调器输入信号为:(t)=m(t)cos(t)(1-11)cos0tOttOm(t)sDSB(t)OtO-ccM()OH-HSDSB()O-cc载波反相点2H带通滤波器sm(t)sm(t)n(t)ni(t)mo(t)no(t)低通滤波器cosct+5与相干载波cos(t)相乘后,得(1-12)经低通滤波器后,输出信号为:(1-13)因此,解调器输出端的有用信号功率为:(1-14)解调DSB信号时,接收机中的带通滤波器的中心频率与调制频率相同,此解调器输入端的窄带噪声:)sin()()cos()()(tttttcsccinnn(1-15)它与相干载波cos(t)相乘后,得:(1-16)经低通滤波器后,解调器最终输出噪声为:)(21)(0ttnnc(1-17)故输出噪声功率为:(1-18)式中,B=2,为DSB的带通滤波器的带宽,为噪声单边功率谱密度。解调器输入信号平均功率为:)(212tmsi(1-19)可得解调器的输入信噪比:BtmnNSii02)(21(1-20)解调器的输出信噪比:BtmnNS0200)((1-21)因此制度增益为:NSNSGiiDSB00(1-22)也就是说,DSB信号的解调器使信噪比改善一倍。)2cos()(21)cos()(21)(cos)(2ttmttmttmccc)(21)(0tmtm)(41)(2200ttmms)]2sin()()2cos()([21)(21)cos()(tttttttcsccccinnnnBtntnnNNic022004141)(41)(62DSB双边带信号的模拟与仿真2.1DSB调制解调分析的MATLAB实现信号DSB调制采用MATLAB函数modulate实现,其函数格式为:Y=MODULATE(X,Fc,Fs,METHOD,OPT)X为基带调制信号,Fc为载波频率,Fs为抽样频率,METHOD为调制方式选择,DSB调制时为’am’,OPT在DSB调制时可不选,Fs需满足Fs2*Fc+BW,BW为调制信号带宽。DSB信号解调采用MATLAB函数demod实现,其函数使用格式为:X=DEMOD(Y,Fc,Fs,METHOD,OPT)Y为DSB已调信号,Fc为载波频率,Fs为抽样频率,METHOD为解调方式选择,DSB解调时为’am’,OPT在DSB调制时可不选。观察信号频谱需对信号进行傅里叶变换,采用MATLAB函数fft实现,其函数常使用格式为:Y=FFT(X,N),X为时域函数,N为傅里叶变换点数选择,一般取值2'。频域变换后,对频域函数取模,格式:Y1=ABS(Y),再进行频率转换,转换方法:f=(0:length(Y)-1)’*Fs/length(Y)分析解调器的抗噪性能时,在输入端加入高斯白噪声,采用MATLAB函数awgn实现,其函数使用格式为:Y=AWGN(X,SNR),加高斯白噪声于X中,SNR为信噪比,单位为dB,其值在假设X的功率为0dBM的情况下确定。信号的信噪比为信号中有用的信号功率与噪声功率的比值,根据信号功率定义,采用MATLAB函数var实现,其函数常使用格式为:Y=VAR(X),返回向量的方差,则信噪比为:SNR=VAR(X1)/VAR(X2)。绘制曲线采用MATLAB函数plot实现,其函数常使用格式:PLOT(X,Y),X为横轴变量,Y为纵轴变量,坐标范围限定AXIS([x1x2y1y2]),轴线说明XLABEL(‘‘)和YLABEL(‘‘)。程序设计流程图见附录。2.2仿真程序用频率300HZ正弦波调制频率30KHZ的正弦波,采用同步解调,观察调制信号、已调信号、解调信号的波形、频谱以及解调器输入输出信噪比的关系。编程如下:Fs=100000;%抽样频率Fc=30000;%载波频率N=1000;%FFT长度n=0:N-1;t=n/Fs;%截止时间和步长x=sin(2*pi*300*t);%基带调制信号y=modulate(x,Fc,Fs,'am');%抑制双边带振幅调制7yn=awgn(y,4);%加入高斯白噪声yn1=awgn(y,10);yn2=awgn(y,15);yn3=awgn(y,20);yn4=awgn(y,25);y1=demod(y,Fc,Fs,'am');%无噪声已调信号解调yyn=demod(yn,30000,Fs,'am');%加噪声已调信号解调yyn1=demod(yn1,30000,Fs,'am');yyn2=demod(yn2,30000,Fs,'am');yyn3=demod(yn3,30000,Fs,'am');yyn4=demod(yn4,30000,Fs,'am');dy1=yn-y;%高斯白噪声snr1=var(y)/var(dy1);%输入信噪比dy2=yyn-y1;%解调后噪声snr2=var(y1)/var(dy2);%输出信噪比dy11=yn1-y;snr11=var(y)/var(dy11);dy21=yyn1-y1;snr21=var(y1)/var(dy21);dy12=yn2-y;snr12=var(y)/var(dy12);dy22=yyn2-y1;snr22=var(y1)/var(dy22);dy13=yn3-y;snr13=var(y)/var(dy13);dy23=yyn3-y1;snr23=var(y1)/var(dy23);dy14=yn4-y;snr14=var(y)/var(dy14);dy24=yyn4-y1;snr24=var(y1)/var(dy24);in=[snr1,snr11,snr12,snr13,snr14];8out=[snr2,snr21,snr22,snr23,snr24];ff1=fft(x,N);%傅里叶变换mag1=abs(ff1);%取模f1=(0:length(ff1)-1)'*Fs/length(ff1);%频率转换ff2=fft(y,N);mag2=abs(ff2);f2=(0:length(ff2)-1)'*Fs/length(ff2);ff3=fft(y1,N);mag3=abs(ff3);f3=(0:length(ff3)-1)'*Fs/length(ff3);figure(1);subplot(221)%绘制曲线plot(t,x)xlabel('调制信号波形')subplot(222)plot(f1,mag1)axis([0100001000])xlabel('调制信号频谱')subplot(223)plot(t,y)xlabel('已调信号波形')subplot(224)plot(f2,mag2)axis([0400000500])xlabel('已调信号频谱')figure(2);subplot(311)plot(t,yyn)xlabel('加噪声解调信号波形')subplot(313)plot(f3,mag3)axis([010000600]
本文标题:基于MATLAB模拟调制系统的仿真
链接地址:https://www.777doc.com/doc-7308886 .html