代码基于MATLAB的AM调制解调系统仿真报告.docx
- 文档编号:23447074
- 上传时间:2023-05-17
- 格式:DOCX
- 页数:23
- 大小:165.79KB
代码基于MATLAB的AM调制解调系统仿真报告.docx
《代码基于MATLAB的AM调制解调系统仿真报告.docx》由会员分享,可在线阅读,更多相关《代码基于MATLAB的AM调制解调系统仿真报告.docx(23页珍藏版)》请在冰豆网上搜索。
代码基于MATLAB的AM调制解调系统仿真报告
内部
基于MATLAB的AM调制解调系统
仿真报告
XXXX-XXXX-XXXX-XXXX
V1.0
天津市智能信号与图像处理重点实验室
2013年10月29日
修订历史记录
日期
版本
文档负责人
修改内容
2013-10-29
V1.0
刘亚洲
创建文档
编制
姓名
签字
日期
电话
审查
姓名
签字
日期
电话
审核
姓名
签字
日期
电话
批准
姓名
签字
日期
电话
文档评审负责人:
参加评审人员:
1引言
1.1设计目的
本报告依照传统模拟调制的规范,给出了AM调制解调的具体流程,重点研究了系统中各阶段信号时域和频谱波形以及频谱的搬移变化,为AM调制解调系统信号波形的进一步深入研究做了基础。
1.2术语定义
本文档使用以下关键术语和简略语。
英文缩写
英文全称
中文名称
AM
AmplitudeModulation
幅度调制
AWGN
AdditiveWhiteGaussianNoise
加性高斯白噪声
1.3参考资料
[1]通信原理(第六版)樊昌信曹丽娜编著国防工业出版社2007年1月
1.4文档组织
报告第二部分给出了AM调制解调的基本原理;第三部分给出了系统在MATLAB里面的程序调试及仿真;第四部分给出了各仿真模块输出时域和频域波形,并对比发射信号和接收信号的时域波形;第五部分对报告进行了总结。
2AM调制解调
图1AM调制解调系统框图
图1显示给出了用于AM调制解调的系统框图。
从图中可知,发送端信源信号
经AM调制器的调制后搭载高频载波发送出去,发送信道中经历加性高斯白噪声的干扰。
接收端信号经历AM解调器的解调输出,最终得到信宿信号
。
2.1AM调制
图2AM调制模型
图2显示给出了AM调制的原理模型。
从图中可知发送信号
和直流分量
叠加后乘以高频载波
后即可形成AM调制信号。
具体时域表波形为:
(1)
对应的频谱波形为:
(2)
2.2AM解调
图3AM解调模型
图3显示给出了AM解调的数学模型。
由上图可知,解调端信道输出信号
乘以跟发送端同频同相的高频载波
后,经低通滤波器提取低频分量,即可得到原始的基带调制信号。
具体理论推导如下:
送入解调器的AM的表达为
(3)
与同频同相的相干载波相
乘后得
(4)
经历低通滤波器滤除高频信号后得
(5)
再经过隔直流电容后
(6)
3基于MATLAB的AM仿真
本仿真利用MATLAB自带函数仿真了AM系统各阶段信号的时域和频域波形,并对比发射信号和解调恢复信号的波形。
3.1仿真基本参数
系统仿真前定义采样间隔ts,载波频率fc,噪声系数sigma,和时间长度t。
具体如下:
%----------------------------------------------------------------
ts=1.e-4;
t=-0.04:
ts:
0.04;
fc=500;
sigma=0.3;
3.2生成调制信号
仿真中采用时间长度[-0.020.02]高度为1的三角脉冲波形作为调制信号。
%----------------------------------------------------------------
%生成调制信号
m_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02);%调整信号m_sig由两个三角脉冲波形左右平移后叠加生成。
3.3AM调制器
在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.4相干解调器
相干解调依靠接收端提供跟发送端同频同相的高频载波
乘以接收信号,再经低通滤波器提取低频分量来恢复出原调制信号。
具体程序实现如下:
%----------------------------------------------------------------
%产生本地接收载波
s_carr=cos(2*pi*fc*t);
%同步解调
s_dem=s_am.*s_carr;
%定义lfft变量
Lfft=length(t);%取时间区间的长度并赋值给Lfft。
Lfft=2^ceil(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_dem);%理想低通滤波器filter滤除s_dem中的高频分量得恢复信号s_rec。
%绘制恢复信号频谱
S_rec=fftshift(fft(s_rec,Lfft)/(length(t)));
4仿真结果曲线
4.1发送信号波形和频谱
图4发送信号时域波形
图4显示给出了用于调制的发送信号时域波形。
图中横坐标和纵坐标分别对应表示时间和信号幅值。
从图中可以明显看出发送信号为三角波信号。
图5发送信号的频谱
图5显示给出了三角波发送信号对应的频谱。
图中横坐标表示频率,纵坐标表示频谱幅值。
从图中可以看出信号频谱主要集中在低频段,而且信号带宽较窄,大约在150Hz。
4.2载波信号波形和频谱
图6载波信号时域波形
图6显示给出了发送端调制基带信号所使用的载波波形。
图中横坐标表示载波持续时间,横坐标表示载波幅度。
从图中可知本次调制使用的载波是幅度为1,周期为
,频率为500Hz的余玄波。
图7载波信号频谱
图7显示给出了系统发射信号所需载波信号频谱。
图中横坐标表示频率,纵坐标表示信号频谱幅值。
从图中可以明显看出载波频谱的中心频率分别为+500Hz,-500Hz。
4.3AM信号波形和频谱
图8AM信号波形(SNR=19.72dB)
图8显示给出了系统中经过AM调制后信号的时域波形。
图中横坐标和纵坐标分别对应表示时间和信号幅值。
从图中可知,已调制信号的外包络仍然保持着跟调制信号相同的包络特性。
图9AM信号频谱波形(SNR=19.72dB)
图9显示给出了经历调制后信号所对应的频谱。
频谱图中横坐标表示频率,纵坐标表示频谱幅值。
从图中可知:
原发射信号经过调制后,频谱明显从原来的低频部分搬移到载波频率对应的高频部分。
但由于已调制信号中直流分量的作用,调制后信号频谱的幅值相比原调制信号频谱幅值并没有完全减半。
4.4相干解调波形和频谱
图10经过相干解调后的时域信号波形
图10显示给出了相干解调中信号乘以跟发送载波同频的本地余玄波后的时域波形。
图中横坐标为时间轴,纵坐标表示信号幅值。
从上述图中可知:
AM信号完整的输入波形转化为单一极性输出,同时信号的频率增加。
图11相干解调频谱
图11显示给出了信号通过相干解调乘以本地载波后对应频谱图。
图中横坐标表示频率变化范围,纵坐标表示频谱幅值。
由上图对比图9可知,乘以本地载波后的信号频谱将AM信号频谱又重新的搬移,而在高频段仍然保留频谱分量。
4.5恢复信号波形和频谱
图12经过低通滤波器后恢复信号时域波形(SNR=19.72dB)
图12显示给出了经过低通滤波器后恢复出的信号的时域波形。
图中横坐标表示信号长度时间坐标轴,纵坐标表示恢复信号的幅值。
从图中可看出,恢复的信号波形基本上跟发送端发送信号波形吻合。
但由于受到噪声的影响,信号的包络发生了抖动。
图13经历低通滤波器恢复信号频谱图
图13显示给出了相干解调信号经过低通滤波器后对应的频谱。
上述图像中横坐标代表频率变化,纵坐标表示频谱幅度变化。
对比原发射信号频谱可知,此恢复出的信号频谱基本跟发射信号频谱波形吻合并且带宽大约在150Hz,说明能较好的恢复出原调制信号。
5总结
本报告使用MATLAB软件仿真了AM调制解调系统中各阶段信号的时域和频谱波形,更直观的研究了信号在通过系统中各阶段的波形和频谱搬移变化情况。
6程序附录
%--------------------------------------------------------------------
%%AMModulation&Demodulation
%byyazooliu
%2013,10,28
%CopyrightbyN6_207lab
%--------------------------------------------------------------------
%thebasicparameterforAMModulation&Demodulationsystem
%--------------------------------------------------------------------
ts=1.e-4;
t=-0.04:
ts:
0.04;
fc=500;%frequencyofcarrier
sigma=0.2;%noiseindex
%--------------------------------------------------------------------
%Generatemessagesignalbyusingfunction’tripuls’
m_sig=tripuls(t+0.01,0.02)-tripuls(t-0.01,0.02);
%--------------------------------------------------------------------
%FFTofm_sig
Lfft=length(t);%thelengthoft
Lfft=2^ceil(log2(Lfft));%shiftthevariableof‘Lfft’
M_fre=fftshift(fft(m_sig,Lfft)/(length(t)));
freqm=(-Lfft/2:
Lfft/2-1)/(Lfft*ts);%frequencyrange
B_m=150;%BandwidthofthesignalisB_mHz
h=fir1(60,[B_m*ts]);%B_m*ts=0.015*5K=75Hz,thecut_offfreqofLPFfilter
%--------------------------------------------------------------------
%AmsignalgeneratedbyaddingacarriertoDSB-SC
s_carr=cos(2*pi*fc*t);%carriersignal
s_am=(1+m_sig).*cos(2*pi*fc*t)+sigma*randn(size(t));%modulatedsignal
%--------------------------------------------------------------------
%ffttos_am
%--------------------------------------------------------------------
Lfft=length(t);
Lfft=2^ceil(log2(Lfft)+1);
S_am=fftshift(fft(s_am,Lfft)/(length(t)));%fftto”s_am”
freqs=(-Lfft/2:
Lfft/2-1)/(Lfft*ts);%frequencyrange
%--------------------------------------------------------------------
%DisplaytheTime_domain&Frequency_domainwaveformofmessage_signal
%--------------------------------------------------------------------
figure
(1);
h1=plot(t,m_sig);
%title('MessageSignal');
xlabel('t(sec)');ylabel('m(t)');
axis([-0.0250.025-22]);
set(h1,'Linewidth',1.5);
figure
(2);
h2=plot(freqm,abs(M_fre));
%title('MessageSpectrum');
xlabel('f(Hz)');ylabel('M(f)');
axis([-20020000.5]);
set(h2,'Linewidth',1.5);
%--------------------------------------------------------------------
%DisplaytheTime_domaim&Frequency_domainwaveformOfCarrier
%--------------------------------------------------------------------
figure(3);
S_carr=fftshift(fft(s_carr,Lfft)/(length(t)));
h3=plot(freqs,abs(S_carr));
%title('CarrierSpectrum');
xlabel('f(Hz)');ylabel('{\itS}_{\itcarr}({\itf})');
axis([-70070000.5]);set(h3,'Linewidth',1.5);
figure(4);
h33=plot(t,s_carr);
%title('CarrierSignal');
xlabel('t(sec)');ylabel('{\its}_{carr}({\itt})');
axis([-0.0050.005-1.51.5]);set(h33,'Linewidth',1.5);
%--------------------------------------------------------------------
%DisplaytheTime_domain&Frequency_domainwaveofAMDodulatied_signal
%--------------------------------------------------------------------
figure(5);
h4=plot(t,s_am);xlabel('t(sec)');ylabel('{\its}_{\itAM}({\itt})');
%title('AMModulatedSignal');
axis([-0.0250.025-22]);set(h4,'Linewidth',1.5);
figure(6);
%S_amisgivenbefore
h5=plot(freqs,abs(S_am));xlabel('f(Hz)');ylabel('{\itS}_{\itAM}({\itf})');
%title('AMSpecturm');
axis([-70070000.5]);set(h5,'Linewidth',1.5);
%--------------------------------------------------------------------
%Demodulationbeginsbyusingthecoherent_demodulation
%--------------------------------------------------------------------
s_dem=s_am.*s_carr;
S_dem=fftshift(fft(s_dem,Lfft)/(length(t)));
%--------------------------------------------------------------------
%UsingtheidealLPFwithbandwidth150Hz
%--------------------------------------------------------------------
s_rec=filter(h,1,s_dem);
S_rec=fftshift(fft(s_rec,Lfft)/(length(t)));
%--------------------------------------------------------------------
%DisplaytheTime_domain&Frequency_domainwaveofCoherent_signal
%--------------------------------------------------------------------
figure(7);
h6=plot(t,s_dem);xlabel('t(sec)');ylabel('{\its}_{\itrect}({\itt})');
%title('CohernetSignal');
axis([-0.0250.025-22]);set(h6,'Linewidth',1.5);
figure(8);
h7=plot(freqs,abs(S_dem));xlabel('f(Hz)');ylabel('{\itS}_{\itrect}({\itf})');
%title('CoherentSpecturm');
axis([-1200120000.5]);set(h7,'Linewidth',1.5);
%--------------------------------------------------------------------
%DisplaytheTime_domain&Frequency_domainwaveofRcecovered_signal
%--------------------------------------------------------------------
figure(9);
h8=plot(t,s_rec);xlabel('t(sec)');ylabel('{\its}_{\itrec}({\itt})');
%title('RecoveredSignal');
axis([-0.0250.025-0.51]);set(h8,'Linewidth',1.5);
figure(10);
h9=plot(freqs,abs(S_rec));xlabel('f(Hz)');ylabel('{\itS}_{\itrec}({\itf})');
%title('RecoveredSpecturm');
axis([-20020000.5]);set(h9,'Linewidth',1.5);
%--------------------------------------------------------------------
%Displayalltime_domainwaveononefigure(10)
%--------------------------------------------------------------------
region1=[-0.0250.025-22];
figure(11);
subplot(221);
h_1=plot(t,m_sig);title('MessageSignal');xlabel('t(sec)');ylabel('m(t)');
set(h_1,'Linewidth',1.5);gridon;
axis(region1);
subplot(222);
h_2=plot(t,s_am);xlabel('t(sec)');ylabel('{\its}_{\itAM}({\itt})');
title('AMModulatedSignal');
set(h_2,'Linewidth',1.5);axis(region1);gridon;
subplot(223);
h_3=plot(t,s_dem);xlabel('t(sec)');ylabel('{\its}_{\itrect}({\itt})');
title('RectifiedSignal');
axis(region1);set(h_3,'Linewidth',1.5);gridon;
subplot(224);
h_4=plot(t,s_rec);xlabel('t(sec)');ylabel('{\its}_{\itrec}({\itt})');
title('RecoveredSignal');
axis([-0.0250.025-0.51]);set(h_4,'Linewidth',1.5);gridon;
%--------------------------------------------------------------------
%DisplayallFrequency_domainwaveononefigure(11)
%--------------------------------------------------------------------
region2=[-70070000.5];
figure(12);
subplot(221);
h_5=plot(freqm,abs(M_fre));title('MessageSpectrum');xlabel('
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 代码 基于 MATLAB AM 调制 解调 系统 仿真 报告
![提示](https://static.bdocx.com/images/bang_tan.gif)