- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
数字通信理论实验报告
——〔实验4〕
实验四数字通信理论课程设计
4.1实验目的:
加深对AWGN信道下数字通信系统的理解。
掌握数字通信系统蒙特卡洛仿真的根本方法。
4.2实验内容:
在AWGN信道下,完成2FSK系统的误比特率性能仿真,绘制系统的BER曲线,并与理论计算的结果进行比照。具体包括如下内容:
编写程序生成随机的二元比特序列,该序列由{0,1}构成。
根据所选择的调制方式,将比特序列映射为星座图上的点。
将所生成的信号通过AWGN信道进行传输,编写程序实现随机的加性高斯白噪声过程,并完成对信号的加噪。
实现接收机的解调、检测与判决算法。要求使用相干接收机,最大似然检测。
在不同的比特信噪比〔〕的条件下统计系统的比特错误概率〔BER〕,画出BER随变化的曲线。的变化范围选为0~10dB。
在同一幅图中画出理论曲线,并将两者进行比照。
4.3实验结果及分析:
画出系统框图,说明仿真流程。
系统框图如下列图所示:
开始
开始
生成随机的二元比特序列
生成随机的二元比特序列
2FSK按照星座图映射
2FSK按照星座图映射
计算高斯白噪声的方差,产生噪声,叠加到调制信号上
计算高斯白噪声的方差,产生噪声,叠加到调制信号上
判决
判决
统计错误个数,计算
统计错误个数,计算BER
信噪比
信噪比10dB
N
画BER的随的变化曲线
画BER的随的变化曲线
结束
结束
仿真流程:编写程序生成随机的二元比特序列,按2FSK调制方式,将比特序列映射为星座图上的点,将所生成的信号通过AWGN信道进行传输,实现随机的加性高斯白噪声过程,并完成对信号的加噪,再使用相干接收机实现接收机的解调、最大似然检测与判决,在不同的比特信噪比〔〕的条件下统计系统的比特错误概率〔BER〕,画出BER随变化的曲线。的变化范围选为0~10dB。
给出仿真结果,提供理论的误码率结果,并进行必要的分析和讨论,
仿真结果如下列图所示:
分析:
如图中实线所示,为理论的信噪比与误码率的关系曲线;
图中所示为对信号进行2FSK调制及解调处理后得到的实际信噪比与误码率关系曲线;
有图可得,随着信噪比的增大,误码率呈下降趋势;
在图中观察可得,由于仿真码元数目的限制,误码率最小只到达了10^(-3)。
2FSK理论误码率推导
发送信号向量为,考虑到信道为高斯白噪声信道,噪声功率密度为。假设2FSK采用相干解调,那么接收信号向量为,,是相互独立的高斯变量,又因为采用最大似然检测,所以当接收端判决正确,那么判断错误的概率为其中所以其概率密度为
又Q函数为化简得2FSk的误码率为。
可见结果与理论值非常吻合。
程序源代码
clearall
snrdB_min=0;snrdB_max=10; %SNR(indB)limits
snrdB=snrdB_min:1:snrdB_max;
Nsymbols=10000;
snr=10.^(snrdB/10); %convertfromdB
h=waitbar(0,MonteCarloSimulation);
len_snr=length(snr);
errors=zeros(1,len_snr);
forj=1:len_snr %incrementSNR
waitbar(j/len_snr)
sigma=sqrt(1/(2*snr(j))); %noisestandarddeviation
error_count=0;
fork=1:Nsymbols %simulationloopbegins
d=round(rand(1)); %data
ifd==0
x_d=1; %directtransmitteroutput
x_q=0; %quadraturetransmitteroutput
else
x_d=0; %directtransmitteroutput
x_q=1; %quadraturetransmitteroutput
end
n_d=sigma*randn(1); %directnoisecomponent
n_q=sigma*randn(1); %quadraturenoisecomponent
您可能关注的文档
最近下载
- 多层住宅小区工程监理规划.docx VIP
- NY_T 3213-2023 植保无人驾驶航空器 质量评价技术规范.docx VIP
- 盘扣式悬挑脚手架施工方案施工计划.doc VIP
- 部编版小学一年级语文上册表格式教案全册.doc VIP
- 电线电缆抗张强度检测结果的测量不确定度评定.doc VIP
- (正式版)H-G-T 22820-2024 化工安全仪表系统工程设计规范.docx VIP
- 绘本故事PPT课件之我好害怕.pptx
- 广东省惠州市博罗县泰美镇城镇总体规划项目建议书.doc VIP
- 虚拟仪器及LabVIEW知到智慧树期末考试答案题库2025年山东科技大学.docx VIP
- 2025年全国Ⅰ卷英语听力评析和听力材料二次开发(课件).pptx
文档评论(0)