基于MATLAB的AM调制解调系统仿真报告

1、基于基于 MATLABMATLAB 的的 AMAM 调制解调系统调制解调系统的仿真报告的仿真报告XXXX-XXXX-XXXX-XXXXXXXX-XXXX-XXXX-XXXX V1.0V1.0市智能信号与图像处理重点实验室市智能信号与图像处理重点实验室20132013 年年 1010 月月 2929 日日2 / 18修订历史记录修订历史记录日期版本文档负责人修改容2013-10-29V1.0V1.0创建文档编制签字日期审查签字日期审核签字日期批准批准3 / 18签字日期文档评审负责人：参加评审人员：目目 录录1 1 引言引言 5 51.1 设计目的 51.2 术语定义 51.3 参考资料 51.

2、4 文档组织 52 2 AMAM 调制解调调制解调 6 62.1 AM 调制 62.2 AM 解调 73 3 基于基于 MATLABMATLAB 的的 AMAM 仿真仿真 8 83.1 仿真基本参数 83.2 生成调制信号 83.3 AM 调制器 83.4 相干解调器 94 4 仿真结果曲线仿真结果曲线 10104.1 发送信号波形和频谱 104.2 载波信号波形和频谱 124.3 AM 信号波形和频谱 144.4 相干解调波形和频谱 164.5 恢复信号波形和频谱 185 5 总结总结 20206 6 程序附录程序附录 20204 / 181 1 引言引言1.11.1 设计目的设计目的本报告

3、依照传统模拟调制的规，给出了 AM 调制解调的具体流程，重点研究了系统中各阶段信号时域和频谱波形以与频谱的搬移变化，为 AM 调制解调系统信号波形的进一步深入研究做了基础。1.21.2 术语定义术语定义本文档使用以下关键术语和简略语。英文缩写英文缩写英文全称英文全称中文名称中文名称AMAmplitude Modulation幅度调制AWGNAdditive White Gaussian Noise加性高斯白噪声1.31.3 参考资料参考资料1通信原理（第六版）樊昌信 丽娜 编著 国防工业 2007 年 1 月21.41.4 文档组织文档组织报告第二部分给出了 AM 调制解调的基本原理；第三部分

4、给出了系统在 MATLAB 里面的程序调试与仿真；第四部分给出了各仿真模块输出时域和频域波形，并对比发射信号和接收信号的时域波形；第五部分对报告进行了总结。5 / 182 2 AMAM 调制解调调制解调信信源源信信号号信信宿宿信信号号AM调调制制AM解解调调信信道道加加性性噪噪声声图 1 AM 调制解调系统框图图 1 显示给出了用于 AM 调制解调的系统框图。从图中可知，发送端信源信号经( )m tAM 调制器的调制后搭载高频载波发送出去，发送信道中经历加性高斯白噪声的干扰。接收端信号经历 AM 解调器的解调输出，最终得到信宿信号。( )m t2.12.1 AMAM 调制调制+m(t)sAM(

5、t)Aocos(ct)图 2 AM 调制模型图 2 显示给出了 AM 调制的原理模型。从图中可知发送信号和直流分量叠加后( )m t0A乘以高频载波后即可形成 AM 调制信号。cos()ct具体时域表波形为： (1)00( )( ) cos( )cos( )( )cos( )AMcccstAm ttAtm tt对应的频谱波形为： (2)01( )()()()()22AMccccASfM ffM ffffff6 / 182.22.2 AMAM 解调解调低低通通滤滤波波器器sm(t)sd(t)cos(ct)图 3 AM 解调模型图 3 显示给出了 AM 解调的数学模型。由上图可知，解调端信道输出信

6、号乘以跟发( )mst送端同频同相的高频载波后，经低通滤波器提取低频分量，即可得到原始的基带cos()ct调制信号。具体理论推导如下：送入解调器的 AM 的表达为 (3)0( )( ) cos()mcstAm tt与同频同相的相干载波相乘后得( )=cos()cc tt (4)2000( )( ) cos ()11( )( ) cos(2)22pccstAm ttAm tAm tt经历低通滤波器滤除高频信号后得 (5)01( )( )2ds tAm t再经过隔直流电容后 (6)1( )( )2s tm t3 3 基于基于 MATLABMATLAB 的的 AMAM 仿真仿真本仿真利用 MATLA

7、B 自带函数仿真了 AM 系统各阶段信号的时域和频域波形，并对比发射信号和解调恢复信号的波形。7 / 183.13.1 仿真基本参数仿真基本参数系统仿真前定义采样间隔 ts,载波频率 fc,噪声系数 sigma,和时间长度 t。具体如下：%-ts=1.e-4;t=-0.04:ts:0.04 ;fc=500;sigma=0.3 ;3.23.2 生成调制信号生成调制信号仿真中采用时间长度-0.02 0.02高度为1的三角脉冲波形作为调制信号。%-% 生成调制信号m_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02)；%调整信号 m_sig 由两个三角脉冲波形

8、左右平移后叠加生成。3.33.3 AMAM 调制器调制器在 MATLAB 中使用如下的语句对发射信号进行调制：%-% 调制过程s_am=(1+m_sig).*cos(2*pi*fc*t)+sigma*randn(size(t)；%AM信号由调制信号m_sig叠加直流分量后乘以载波并混以噪声形成。cos(2*pi*fc*t)；% 载波信号，fc为载波频率。randn(size(t);%代表服从正态分布的噪声。3.43.4 相干解调器相干解调器相干解调依靠接收端提供跟发送端同频同相的高频载波乘以接收信号，再经cos()ct低通滤波器提取低频分量来恢复出原调制信号。具体程序实现如下：%-% 产生本地

9、接收载波8 / 18s_carr=cos(2*pi*fc*t)；% 同步解调s_dem=s_am.*s_carr；% 定义lfft变量Lfft=length(t)；%取时间区间的长度并赋值给Lfft。Lfft=2ceil(log2(Lfft)+1)；%将Lfft重新转化成为更大的偶数Lfft。ceil；%向正无穷大的方向取整% 绘制解调后信号频谱S_dem=fftshift(fft(s_dem,Lfft)/(length(t);% 生成低通滤波器h=fir1(60,B_m*ts)；%设计低通滤波器频域相应，其中滤波器阶数60,截至频率75Hz。% 低通滤波s_rec=filter(h,1,s_

10、dem)；%理想低通滤波器filter滤除s_dem中的高频分量得恢复信号s_rec。% 绘制恢复信号频谱S_rec=fftshift(fft(s_rec,Lfft)/(length(t)；9 / 184 4 仿真结果曲线仿真结果曲线4.14.1 发送信号波形和频谱发送信号波形和频谱-0.025-0.02 -0.015-0.01 -0.00500.0050.010.0150.020.025-2-1.5-1-0.500.511.52t(sec)m(t)图 4 发送信号时域波形图 4 显示给出了用于调制的发送信号时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图中可以明显看出发送信号为三

11、角波信号。10 / 18-150-100-5005010015000.050.10.150.20.250.30.350.40.450.5f(Hz)M(f)图 5 发送信号的频谱图 5 显示给出了三角波发送信号对应的频谱。图中横坐标表示频率，纵坐标表示频谱幅值。从图中可以看出信号频谱主要集中在低频段，而且信号带宽较窄，大约在 150Hz。11 / 184.24.2 载波信号波形和频谱载波信号波形和频谱-5-4-3-2-1012345x 10-3-1.5-1-0.500.511.5t(sec) scarr( t)图 6 载波信号时域波形图 6 显示给出了发送端调制基带信号所使用的载波波形。图中横坐

12、标表示载波持续时间，横坐标表示载波幅度。从图中可知本次调制使用的载波是幅度为 1，周期为，-32.0 10 s频率为 500Hz 的余玄波。12 / 18-600-400-200020040060000.050.10.150.20.250.30.350.40.450.5f(Hz) S carr( f)图 7 载波信号频谱图 7 显示给出了系统发射信号所需载波信号频谱。图中横坐标表示频率，纵坐标表示信号频谱幅值。从图中可以明显看出载波频谱的中心频率分别为+500Hz，-500Hz。13 / 184.34.3 AMAM 信号波形和频谱信号波形和频谱-0.025-0.02 -0.015-0.01 -

13、0.00500.0050.010.0150.020.025-2-1.5-1-0.500.511.52t(sec) s AM( t)图 8 AM 信号波形(SNR=19.72dB)图 8 显示给出了系统中经过 AM 调制后信号的时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图中可知，已调制信号的外包络仍然保持着跟调制信号一样的包络特性。14 / 18-600-400-200020040060000.050.10.150.20.250.30.350.40.450.5f(Hz) S AM( f)图 9 AM 信号频谱波形（SNR=19.72dB）图 9 显示给出了经历调制后信号所对应的频

14、谱。频谱图中横坐标表示频率，纵坐标表示频谱幅值。从图中可知：原发射信号经过调制后，频谱明显从原来的低频部分搬移到载波频率对应的高频部分。但由于已调制信号中直流分量的作用，调制后信号频谱的幅值相比原调制信号频谱幅值并没有完全减半。15 / 184.44.4 相干解调波形和频谱相干解调波形和频谱-0.025-0.02 -0.015-0.01 -0.00500.0050.010.0150.020.025-2-1.5-1-0.500.511.52t(sec) s rect( t)图 10 经过相干解调后的时域信号波形图 10 显示给出了相干解调号乘以跟发送载波同频的本地余玄波后的时域波形。图中横坐标为

15、时间轴，纵坐标表示信号幅值。从上述图中可知：AM 信号完整的输入波形转化为单一极性输出，同时信号的频率增加。16 / 18-1000-5000500100000.10.20.30.40.5f(Hz) S rect( f)图 11 相干解调频谱图 11 显示给出了信号通过相干解调乘以本地载波后对应频谱图。图中横坐标表示频率变化围，纵坐标表示频谱幅值。由上图对比图 9 可知，乘以本地载波后的信号频谱将 AM信号频谱又重新的搬移，而在高频段仍然保留频谱分量。17 / 184.54.5 恢复信号波形和频谱恢复信号波形和频谱-0.025-0.02 -0.015-0.01 -0.00500.0050.010.0150.020.025-0.500.51t(sec) s rec( t)图 12 经过低通滤波器后恢复信号时域波形（SNR=19.72dB）图 12 显示给出了经过低通滤波器后恢复出的信号的时域波形。图中横坐标表示信号长度时间坐标轴，纵坐标表示恢复信号的幅值。从图中可看出，恢复的信号波形基本上跟发送端发送信号波形吻合。但由于受到噪声的影响，信号的包络发生了抖动。18 / 18-200-150-100-5005010015020000.050.10.150.20.250.30