怎么用matlab计算模拟信号,用MATLAB模拟实现数字信号的调制与频谱分析

1调制原理1.1双极性非归零电平码双极性非归零电平码(BipolarNRZ-L)是数字通信中用来构造指令码和数据的一种基本码型,用正电平对应数字1,等幅负电平对应数字“0”。采用BipolarNRZ-L码型。其中每一个数字“1”或“0”叫一个码元。这种码型抗干扰能力较强,有利于在实际信道中传输。在程序设计中,规定每个码元的抽样点数目为N_sample,N_sample是码元周期,为简便起见程序中所产生的基带信号的正负电平幅度都为“1”,并用an来表示。1.2调制方法由于数字基带信号是二进制的,因此调制方法都可以叫做二进制调制。1.2.12ASK调制2ASK调制是正弦载波的幅度随基带信号的变化而变化的数字调制,是最早出现的也是最简单的。抗噪声能力较差,实际应用并不广泛,但可以作为其他调制方式的研究基础。2ASK信号的时域表达式为:e2ASK=bnbAcos0t。式中bn为二进制数序列,当基带信号an为正电平时bn取“1”,当基带信号an为负电平时bn取“0”,由an转化为bn的硬件实现只需一个单向导通的二极管电路;A为调制信号的幅度;0为调制信号的频率。显然,当基带信号为正电平(+1)时,2ASK时域信号为一正弦波;当基带信号为负电平(-1)时,2ASK时域信号对应为“0”。正弦载波频率要比基带频率高的多得多,为了兼顾实际和视觉效果,在程序设计中,使一个码元长度对应于正弦函数的五个周期,故设计程序时频率设计为0=5*1/N_sample*2。1.2.22FSK调制2FSK是正弦载波的频率,随基带信号的变化而变化的数字调制,时域表达式为:an=1时,e2FSK=a*cos1t;an=1时e2FSK=a*cos2t。程序中设定基带信号为正电平时,对应于高频正弦载波,基带信号为负电平时对应低频正弦载波。同时考虑波形显示的实际性和简易性,设定“1”对应载波10个周期,“-1”对应载波5个周期,故表达式中参数设置为:1=10*1/N_sample*2*2,2=5;1/N_sample*2。1.2.32PSK调制2PSK是利用载波相位变化来传输信息的,通常用0或180来表示基带信号的两种状态,时域表达式为:an=1时,e2PSK=an*cos0t,an=-1时,e2PSK=an*cos(0t+),同样设定载波频率0=5*1/N_sample*2。2频谱分析2.1频谱计算MATLAB中的自带快速傅里叶变换函数fft为信号的频谱分析提供了方便,其调用格式为:Y=fft(X),其中X为时域信号序列,Y是X经过快速傅里叶变换后的序列。Y值相当于序列傅里叶变换在0-2上的N个抽样值,其中N为序列长度。在正常采样的情况下,序列Y的前半部与后半部分所含信息相同,分析时只需取出Y序列前半部分。2.2频谱分析数字基带信号,如计算机输出的二进制序列,或来自模拟信号经数字化处理后的PCM码组等。信号包含丰富的低频分量及直流分量,不适合长距离传输。为达到通信的目的,必须将低频的基带信号搬移到高频处,调制目的是实现信号的搬移,以适合信道中的长距离传输。3程序实现3.1双极性非归零电平码的可视化双极性非归零电平码(BipolarNRZ-L)可视化程序设计方法为:每一个码元显示由一定数目的“0”或“1”组成,如绘制[01]序列,设置密集度为5的话,则plot命令绘制图形参数为[-1-1-1-1-1]。functionA=B_NRZ_L(a,N_sample)%双极性非归零电平码的可视化%a为原码,N_sample为抽样点数n=length(a)Vb=find(a==0)Va(b)=-1Vp