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

上传人:桔**** 文档编号:492890843 上传时间:2023-08-28 格式:DOC 页数:12 大小:205.50KB
返回 下载 相关 举报
基于MATLAB的AM调制解调系统仿真报告_第1页
第1页 / 共12页
基于MATLAB的AM调制解调系统仿真报告_第2页
第2页 / 共12页
基于MATLAB的AM调制解调系统仿真报告_第3页
第3页 / 共12页
基于MATLAB的AM调制解调系统仿真报告_第4页
第4页 / 共12页
基于MATLAB的AM调制解调系统仿真报告_第5页
第5页 / 共12页
点击查看更多>>
资源描述

《基于MATLAB的AM调制解调系统仿真报告》由会员分享,可在线阅读,更多相关《基于MATLAB的AM调制解调系统仿真报告(12页珍藏版)》请在金锄头文库上搜索。

1、-内部基于MATLAB的AM调制解调系统仿真报告*-*-*-* V1.0*市智能信号与图像处理重点实验室2013年10月29日修订历史记录日期版本文档负责人修改容2013-10-29V1.0亚洲创立文档编制签字日期审查签字日期审核签字日期批准签字日期文档评审负责人:参加评审人员:目 录1引言51.1设计目的51.2术语定义51.3参考资料51.4文档组织52 AM调制解调62.1 AM调制62.2 AM解调73基于MATLAB的AM仿真83.1仿真根本参数83.2生成调制信号83.3 AM调制器83.4相干解调器94 仿真结果曲线104.1发送信号波形和频谱104.2载波信号波形和频谱124.

2、3 AM信号波形和频谱144.4相干解调波形和频谱164.5恢复信号波形和频谱185总结206程序附录20. z.-1引言1.1设计目的本报告依照传统模拟调制的规,给出了AM调制解调的具体流程,重点研究了系统中各阶段信号时域和频谱波形以及频谱的搬移变化,为AM调制解调系统信号波形的进一步深入研究做了根底。1.2术语定义本文档使用以下关键术语和简单语。英文缩写英文全称中文名称AMAmplitude Modulation幅度调制AWGNAdditive White Gaussian Noise加性高斯白噪声1.3参考资料1通信原理第六版樊昌信 丽娜 编著 国防工业 2007年1月1.4文档组织报告

3、第二局部给出了AM调制解调的根本原理;第三局部给出了系统在MATLAB里面的程序调试及仿真;第四局部给出了各仿真模块输出时域和频域波形,并比照发射信号和接收信号的时域波形;第五局部对报告进展了总结。2AM调制解调图1 AM调制解调系统框图图1显示给出了用于AM调制解调的系统框图。从图中可知,发送端信源信号经AM调制器的调制后搭载高频载波发送出去,发送信道中经历加性高斯白噪声的干扰。接收端信号经历AM解调器的解调输出,最终得到信宿信号。2.1AM调制图2 AM调制模型图2显示给出了AM调制的原理模型。从图中可知发送信号和直流分量叠加后乘以高频载波后即可形成AM调制信号。具体时域表波形为: (1)

4、对应的频谱波形为: (2)2.2AM解调图3 AM解调模型图3显示给出了AM解调的数学模型。由上图可知,解调端信道输出信号乘以跟发送端同频同相的高频载波后,经低通滤波器提取低频分量,即可得到原始的基带调制信号。具体理论推导如下:送入解调器的AM的表达为(3)与同频同相的相干载波相乘后得 (4)经历低通滤波器滤除高频信号后得 (5)再经过隔直流电容后 (6)3基于MATLAB的AM仿真本仿真利用MATLAB自带函数仿真了AM系统各阶段信号的时域和频域波形,并比照发射信号和解调恢复信号的波形。3.1仿真根本参数系统仿真前定义采样间隔ts,载波频率fc,噪声系数sigma,和时间长度t。具体如下:%

5、-ts=1.e-4;t=-0.04:ts:0.04;fc=500;sigma=0.3;3.2生成调制信号仿真中采用时间长度-0.02 0.02高度为1的三角脉冲波形作为调制信号。%-%生成调制信号m_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02);%调整信号m_sig由两个三角脉冲波形左右平移后叠加生成。3.3 AM调制器在MATLAB中使用如下的语句对发射信号进展调制:%-% 调制过程s_am=(1+m_sig).*cos(2*pi*fc*t)+sigma*randn(size(t);%AM信号由调制信号m_sig叠加直流分量后乘以载波并混以噪声形

6、成。cos(2*pi*fc*t);% 载波信号,fc为载波频率。randn(size(t);%代表服从正态分布的噪声。3.4相干解调器相干解调依靠接收端提供跟发送端同频同相的高频载波乘以接收信号,再经低通滤波器提取低频分量来恢复出原调制信号。具体程序实现如下:%-%产生本地接收载波s_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;%向正无穷大的方向取整%绘制解调后

7、信号频谱S_dem=fftshift(fft(s_dem,Lfft)/(length(t);% 生成低通滤波器h=fir1(60,B_m*ts);%设计低通滤波器频域相应,其中滤波器阶数60,截至频率75Hz。% 低通滤波s_rec=filter(h,1,s_dem);%理想低通滤波器filter滤除s_dem中的高频分量得恢复信号s_rec。% 绘制恢复信号频谱S_rec=fftshift(fft(s_rec,Lfft)/(length(t);4仿真结果曲线4.1发送信号波形和频谱图4发送信号时域波形图4显示给出了用于调制的发送信号时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图

8、中可以明显看出发送信号为三角波信号。图5发送信号的频谱图5显示给出了三角波发送信号对应的频谱。图中横坐标表示频率,纵坐标表示频谱幅值。从图中可以看出信号频谱主要集中在低频段,而且信号带宽较窄,大约在150Hz。4.2载波信号波形和频谱图6 载波信号时域波形图6显示给出了发送端调制基带信号所使用的载波波形。图中横坐标表示载波持续时间,横坐标表示载波幅度。从图中可知本次调制使用的载波是幅度为1,周期为,频率为500Hz的余玄波。图7 载波信号频谱图7显示给出了系统发射信号所需载波信号频谱。图中横坐标表示频率,纵坐标表示信号频谱幅值。从图中可以明显看出载波频谱的中心频率分别为+500Hz,-500H

9、z。4.3 AM信号波形和频谱图8 AM信号波形(SNR=19.72dB)图8显示给出了系统中经过AM调制后信号的时域波形。图中横坐标和纵坐标分别对应表示时间和信号幅值。从图中可知,已调制信号的外包络仍然保持着跟调制信号一样的包络特性。图9 AM信号频谱波形SNR=19.72dB图9显示给出了经历调制后信号所对应的频谱。频谱图中横坐标表示频率,纵坐标表示频谱幅值。从图中可知:原发射信号经过调制后,频谱明显从原来的低频局部搬移到载波频率对应的高频局部。但由于已调制信号中直流分量的作用,调制后信号频谱的幅值相比原调制信号频谱幅值并没有完全减半。4.4相干解调波形和频谱图10经过相干解调后的时域信号

10、波形图10显示给出了相干解调号乘以跟发送载波同频的本地余玄波后的时域波形。图中横坐标为时间轴,纵坐标表示信号幅值。从上述图中可知:AM信号完整的输入波形转化为单一极性输出,同时信号的频率增加。图11相干解调频谱图11显示给出了信号通过相干解调乘以本地载波后对应频谱图。图中横坐标表示频率变化围,纵坐标表示频谱幅值。由上图比照图9可知,乘以本地载波后的信号频谱将AM信号频谱又重新的搬移,而在高频段仍然保存频谱分量。4.5恢复信号波形和频谱图12 经过低通滤波器后恢复信号时域波形SNR=19.72dB图12显示给出了经过低通滤波器后恢复出的信号的时域波形。图中横坐标表示信号长度时间坐标轴,纵坐标表示

11、恢复信号的幅值。从图中可看出,恢复的信号波形根本上跟发送端发送信号波形吻合。但由于受到噪声的影响,信号的包络发生了抖动。图13 经历低通滤波器恢复信号频谱图图13显示给出了相干解调信号经过低通滤波器后对应的频谱。上述图像中横坐标代表频率变化,纵坐标表示频谱幅度变化。比照原发射信号频谱可知,此恢复出的信号频谱根本跟发射信号频谱波形吻合并且带宽大约在150Hz,说明能较好的恢复出原调制信号。5总结本报告使用MATLAB软件仿真了AM调制解调系统中各阶段信号的时域和频谱波形,更直观的研究了信号在通过系统中各阶段的波形和频谱搬移变化情况。6程序附录%-% AM Modulation &Demodula

12、tion% by yazooliu% 2013,10,28% Copyright by N6_207lab%-% the basic parameter for AM Modulation &Demodulation system%-ts=1.e-4;t=-0.04:ts:0.04;fc=500; %frequency of carriersigma=0.2;%noise inde*%-% Generate message signal by using functiontripulsm_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02);%-% FFT of m_sigLfft=length(t);%the length of tLfft=2ceil(log2(Lfft);% shift the variable of LfftM_fre=fftshift(fft(m_sig,Lfft)/(length(t); freqm=(-Lfft/2:Lfft/2-1)/(Lfft*ts);%frequency rangeB_m=1

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 建筑/环境 > 施工组织

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号