基于MATLAB通信系统调制和解调实例分析.docVIP

基于MATLAB通信系统调制和解调实例分析 　　【摘 要】通信是通过某种媒体进行的信息传递，目的是传输信息，通信系统是用以完成信息传输过程的技术系统的总称，作用是将信息从信源发送到一个或多个目的地。调制与解调在信息的传输过程中占据着重要的地位，是不可或缺的，因此研究系统的调制和解调过程就极为重要。MATLAB是集数值计算、图形绘制、图像处理及系统仿真等强大功能于一体的科学计算语言，它强大的矩阵运算和图形可视化的功能以及丰富的工具箱，为通信系统的调制和解调过程的分析提供了极大的方便。 　　【关键词】通信系统 调制与解调 MATLAB 　　一、用MATLAB分析双边带幅度调制（DSM-AM） 　　我们可以利用 MATLAB 强大的符号运算功能来进行运算，再根据MATLAB的可视化结果进行分析。在DSB-AM中，已调信号的时域表示为：u（t）=m（t）c（t）=Acm（t） 式中，m（t）是消息信号，c（t）=Ac 为载波，fc是载波的频率（单位：HZ）， 是初始相位。为了讨论方便取初相=0（以下类似）。随u（t）作傅里叶变换，即可得到信号的频域表示：U（f）=Ac/2M（f-fc）+Ac/2M（f+fc）。传输带宽Bt是消息信号带宽W的两倍，即：Bt=2W。 　　某消息信号 　　用信号m（t）以DSB-AM方式调制载波c（t）=cos（2∏fct），所得到的已调制信号记为u（t）.设t。=0.15s，fc=250Hz。试比较消息信号与已调信号，并绘制它们的频谱。运用如下 MATLAB程序： 　　t0=0.15；%信号持续时间 　　ts=0.001；%采样时间间隔 　　Fc=250；%载波频率 　　FS=1/ts；%采样频率 　　df= 0.3；%频率分辨率 　　t=[0：ts：t0]；%时间矢量 　　m=[ones（1，t0/（3*ts））， -2*ones（1，t0/（3*ts）），zeros（1，t0/（3*ts）+1）]；%定义信号序列 　　C=cos（2*pi*Fc.*t）；%载波信号 　　u=m.*c；%调制信号 　　[M，m，df1]=fft_seq（m，ts，df）；%傅里叶变换 　　M=M/FS；[U，u，df1]=fft_seq（m，ts，df）； 　　U=U/FS；[C，c，df1]=fft_seq（c，ts，df）； 　　f=[0：df1：df1*（length（m）-1）]-Fs/2； %频率矢量 　　subplot（2，2，1）；plot（t，m（1：length（t）））；%未调制信号 　　title（‘未调制信号‘）； 　　subplot（2，2，2）；plot（t，u（1：length（t）））；%已调制信号 　　title（‘已调制信号‘）； 　　subplot（2，2，3）；plot（f，abs（fftshift（M）））； %未调制信号频谱 　　title（‘未调制信号频谱‘）；subplot（2，2，4）；plot（f，abs（fftshift（U）））； %已调制信号频谱 　　title（‘已调制信号频谱‘）； 　　傅里叶变换函数fft_seq，源代码如下： 　　function[M，m，df]=fft_seq（m，ts，df） 　　%[M，m，df]=fft_seq（m，ts，df） 　　%[M，m，df]=fft_seq（m，ts） 　　%M为输入序列m的傅里叶变换，ts为抽样间隔，输入df为频率分辨率 　　%输出序列m按要求的频率分辨率df进行补零后的序列 　　%输出df为最终的频率分辨率 　　fs=1/ts； 　　if nargin==2，n1=0； 　　else， n1=fs/df； 　　end 　　n2=length（m）； 　　n=2^（max（nextpow2（n1），nexpow2（n2）））； 　　M=fft（m，n）； 　　m=[m，zeros（n-n2）]； 　　df=fs/n； 　　运行后得到的信号和调制信号如图所示： 　　二、用 MATLAB 分析相干解调过程 　　例 　　用信号m（t）以DSB-AM方式调制载波c（t）=cos（2∏fct），所得到的已调制信号记为u（t）.设t0=0.15s，fc=250Hz。得到DSM-AM的调制信号，现试对该调制信号进行相干解调，并绘出消息信号的时域频域曲线。 　　MATLAB设计程序如下： 　　t0=0.15；%信号持续时间 　　ts=1/1500；%采样时间间隔 　　Fc=250；%载波频率 　　FS=1/ts；%采样频率 　　df= 0.3；%频率分辨率 　　t=[0：ts：t0]；%时间矢量 　　m=[on