《matlab在通信工程中的应用》由会员分享，可在线阅读，更多相关《matlab在通信工程中的应用（13页珍藏版）》请在金锄头文库上搜索。
1、安阳工学院专业课课程大作业课程名称（中文）MATLAB 在双边带调制与解调中的应用 成绩 姓名 夏娟班级 10 级通信工程学号 201002030026日期 2012 年 12 月 22 日1MATLAB 在双边带调制与解调中的应用【摘要摘要】MATLAB 是一个专门以矩阵的形式处理数据的软件，具有强大的矩阵计算和绘图功能，适用于科学计算、控制系统、信息处理等领域的分析、仿真和设计工作。双边带调制是一种通过抑制载波来实现高效率的调制方式，其解调方式为相干解调，是调制的逆过程，将已调信号中的调制信号恢复出来；调制与解调是通信系统中至关重要的方式。双边带调制与解调在 MATLAB 中的实现是一个在
2、现实中实现的基础仿真模型。【关关键键字字】MATLAB、双边带、调制、解调21.1. 调制与解调的概述调制与解调的概述所谓调制，就是把信号转换成适合在信道中传输的形式的一种过程，广义的调制分为基带调制和带通调制（也称为载波调制）。在无线通信中和其他大多数场合，调制一词均指载波调制。载波调制就是用调制信号去控制载波的参数的过程，使载波的某一个或某几个参数按照调制信号的规律而变换。调制信号是指来自信源的消息信号（基带信号），这些信号可以是数字的也可以是模拟的。未受调制的周期性振荡信号称为载波，它可以是正弦信号，也可以是非正弦波。载波调制后称为已调信号，它含有调制信号的全部特征。所谓解调（检波），是
3、调制的逆过程，作用是将已调信号从中的调制信号恢复出来。信号要进行调制的原因有三：第一，在无线传输中，信号是以电磁波的形式通过天线辐射到空间的，为了获得较高的辐射频率，天线的尺寸必须与发射信号波长相比拟。而基带信号包含较低的频率分量的波长较长，使得在现实生活中不易实现。第二，把多个基带信号分别搬至到不同的载频处，以实现信道的多路复用，提高信道利用率。第三，扩展信号宽带，提高系统抗干扰，抗衰落能力，还可实现传输带宽与信噪比之间的互换。调制的方式有很多种，最常用和最重要的模拟调制是用正弦波作为载波的幅度调制和角度调制，常见的幅度调制为 AM、DSB、SSB、VSB 等；而解调方式有检波和同步检波（相
4、干解调）。2.幅度调制与解调的设计思路幅度调制与解调的设计思路2.1 调制与解调的基本模型调制与解调的基本模型调制在通信过程中起着极其重要的作用：无线电通信是通过空间辐射方式传输信号的，调制过程可以将信号的频谱搬移到容易以电磁波形式辐射的较高频范围；此外，调制过程可以将不同的信号通过频谱搬移托付至不同频率的载波上，实现多路复用，不至于互相干扰。振幅调制是一种实用很广的连续波调制方式。调幅信号 X(t)主要有调制信号和载波信号组成。调幅器原理如图 1 所示。3其中载波信号 C(t)用于搭载有用信号，其频率较高。幅度调制信号 g(t)含有有用信息，频率较低。运用 MATLAB 信号 g(t)处理工
5、具箱的有关函数可以对信号进行调制。对于信号 x(t)，通信系统就可以有效而可靠的传输了。在接收端，分析已调信号的频谱，进而对它进行解调，以恢复原调制信号。解调器原理如图 2 所示。对于调制解调的过程以及其中所包含的对于信号的频谱分析均可以通过MATLAB 的相关函数实现。2.2 频谱及功率谱密度分析频谱及功率谱密度分析当调制信号 f(t)为确定信号时，已调信号的频谱为.抑制载波的双边带调幅虽然节省了载波功率，ccSDSB=1/2F-+1/2F( +) 但已制的频带宽度仍为调制信号的两倍，与常规双边带调幅时相同。 通信中，调制信号通常是平稳随机过程。其功率谱密度与自相关函数之间是一对傅氏变换关系
6、。这样就可以先找到信号的自相关函数，然后通过氏傅变换来实现信号的功率谱密度。3.3.设计流程设计流程 3.1 利用利用 matlab 绘制已知信号绘制已知信号 f(t)根据 f(t) 表达示 f(t)=sin(t)，0t2*pi。由于函数是正弦函数，故利用时间t 与 f(t)的关系，再利用 plot 函数实现图的画法，并且图标为幅度。3.23.2 利用利用 matlabmatlab 绘制已知信号绘制已知信号 f(t)f(t)的频谱的频谱根据 f(t)的表达示，通过求傅立叶变换来实现信号的频谱，具体可以取40000 个点来实现。并且运用算法 yw=2*/40000*abs(fftshift(yk
7、)，fw=-425000:24999/50000*fs。这样再利用 plot 函数实现图的画法，并且对所画的图做标识，如标题，幅度。3.33.3 利用利用 matlabmatlab 绘制载波信号绘制载波信号由给定的载波为，fc=100hz，的出余弦信号的画法，这样再利用cos2cf tplot 函数实现图的画法，并且对所画的图做标识，如标题，幅度，时间。3.4 利用利用 matlab 绘制已调信号绘制已调信号由调制信号知：抑制载波双边带调幅的调制过程实际上就是调制信号与载波的相乘运算。故此时将上述两个信号相乘，就可以得出已调信号 y4, y4=sin (t).*cos(2*.*fc.*t).
8、这样再利用 plot 函数实现图的画法，并且对所画的图做标识，如标题，幅度，时间。3.5 利用利用 matlab 绘制已调信号的频谱绘制已调信号的频谱根据已调信号的表达示，提高求傅立叶变换来实现信号的频谱，具体可以取 4000 个点来实现。并且运用算法 yw=2* /4000*abs(fftshift(yk)，fw=-2500:2499/5000*fs。这样再利用 subplot 函数实现子图的画法，并且对所画的图做标识，如标题，幅度。3.6 利用利用 matlab 绘制绘制 DSB-SC 调制信号的功率谱密度调制信号的功率谱密度通信中，调制信号通常是平稳随机过程。其功率谱密度与自相关函数之间
9、是一对付氏变换关系。此时先求调制信号的自相关函数，利用命令c,lags=xcorr(y4,20)以及 plot(lags/fs,c)就可以实现调制信号的自相关函数，此时将自相关函数求付氏变换。利用 SDSBp=fft(c,5000； fw=-2500:2499/5000*fs；yw=2*/4000*abs(fftshift(SDSBp)即可实现，此时用 figure 和 subplot 可以在另一页画出自相关函数波形和功率谱密度波形。3.7 利用利用 matlab 绘制相干解调后的信号波形绘制相干解调后的信号波形由抑制载波双边带调幅的解调过程实际上实际是将已调信号乘上一个同频同相的载波。即 y
10、=sin (t).*cos(2* *fc*t).*cos(2* *fc*t)，此时解调图形如图 10。再用一个低通滤波器就可以恢复原始的调制信号，这种调制方法称为相干解调。主要程序语句为n,Wn=ellipord(Wp,Ws,Rp,Rs)；b,a=ellip(n,Rp,Rs,Wn)；这样可以实现求取阶数 n 和传递函数的分子分母 b,a；Wp=0.5/50; Ws=1/50; Xl=5*filter(b,a,y)。通过这样可以使滤波后的波形失真更小。此时可得相干解调5后的信号波形，具体波形如图 11。4.4.源代码及仿真结果源代码及仿真结果4.1 调制信号的波形及仿真结果（如图调制信号的波形及
11、仿真结果（如图 3）t=0:0.001:2*pi;m=sin(t);plot(t,m)title(调制信号)ylabel(幅度)xlabel(时间：S)grid on图 3：已知信号波形4.2 调制信号的频谱及仿真结果（如图调制信号的频谱及仿真结果（如图 4）t=0:0.001:2*pi;fs=300;m=sin(t);yk=fft(m,50000); %对信号做傅里叶变换yw=2*pi/40000*abs(fftshift(yk) ; %频谱搬移fw=-25000:24999/50000*fs ; plot(fw,yw)title(调制信号的频谱)xlabel(频率：HZ)ylabel(幅度
12、)grid6xlim(-30,30)图 4：已知信号的频谱4.3 载波信号的波形及仿真结果（如图载波信号的波形及仿真结果（如图 5）t=0:0.001:2*pi;fc=100;y=cos(2*pi*fc*t); %载波信号plot(t,y)title(载波信号)xlabel(时间：s)ylabel(幅度)gridxlim(0,0.1)图 5 ：载波信号4.4 已调信号的波形及仿真结果（如图已调信号的波形及仿真结果（如图 6）t=0:0.001:2*pi;m=sin(t);7fc=100;y=cos(2*pi*fc*t); %载波信号u=m.*y;plot(t,u,b)title(已调信号)xl
13、abel(时间：s)ylabel(幅度)gridxlim(0,2*pi)图 6：已调信号4.5 已调信号的频谱及仿真结果（如图已调信号的频谱及仿真结果（如图 7）t=0:0.001:2*pi;m=sin(t);fs=200;fc=100;y=cos(2*pi*fc*t); %载波信号u=m.*y;yk=fft(u,50000); %对信号做傅立叶变换yw=2*pi/40000*abs(fftshift(yk) ; %频谱搬移fw=-25000:24999/50000*fs ; plot(fw,yw)title(已调信号的频谱)xlabel(频率：hz)8ylabel(幅度)gridxlim(-
14、30,30)图 7：已调信号波形的频谱4.6 已调信号的功率谱密度及仿真结果（如图已调信号的功率谱密度及仿真结果（如图 8）t=0:0.001:2*pi;m=sin(t);fs=200;fc=100;y=cos(2*pi*fc*t); %载波信号u=m.*y;c,lags=xcorr(u,200) %DSB信号自相关函数subplot(211) plot(lags/fs,c)title(DSB信号自相关函数)xlabel(t)ylabel(Rxx(t)grid on;SDSBp=fft(c,5000) %DSB功率谱fw=-2500:2499/5000*fsyw=2*pi/4000*abs(f
15、ftshift(SDSBp)%频谱搬移subplot(212)9plot(fw,yw)title(DSB信号功率谱)xlabel(w)ylabel(Rxx(t)grid图8 调制信号自相关函数波形和功率谱密度波形4.7 解调信号波形及仿真结果（如图解调信号波形及仿真结果（如图 9）t=0:0.001:2*pi;m=sin(t);fs=200;fc=100;y=cos(2*pi*fc*t); %载波信号v=m.*y.*y;plot(t,v)title(解调信号)xlabel(时间：s)ylabel(幅度)grid on10图9 解调信号波形4.8 经过滤波后的波形及仿真结果（如图经过滤波后的波形及仿真结果（如图 11）t=0:0.001:2*pi;m=sin(t);%fs=200;fc=100;y=cos(2*pi*fc*t); %载波信号v=m.*y.*y;Rp=0.5; %滤波后的f(t)信号Rs=20; %信号衰减幅度Wp=0.5/50; %通带截止频率Ws=1/50; %阻带截止频率,100为载波频率的一半n,Wn=ellipord(Wp,Ws,Rp,Rs); %阶数nb,a=ellip(n,Rp,Rs,Wn); %传递函数分子分母b,aXl=5*filter(b,a,v);%f
