蒙特卡罗法在数字调频信号解码方面的应用.pdf

蒙特卡罗法在数字调频信号解码方面的应用 马晓晖 北京邮电大学光通信中心，北京（100876 ） E-mail ：mxhui123@ 摘 要：在数字通信中，对于一个混有强线性调频干扰信号的数字调频信号进行解码，通过 应用蒙特卡罗法模拟仿真来估计该解码方法的误码率，并进行了讨论分析。 关键词：蒙特卡罗法，数字调频信号，误码率 中图分类号：TN92 1. 引言 蒙特卡罗方法，也叫蒙特卡罗分析，是一种使用随机抽样统计来估算数学函数的计算方 法。它需要一个良好的随机数源。这种方法往往包含一些误差，但是随着随机抽取样本数量 的增加，结果也会越来越精确。蒙特卡罗法可以说是把某个未知值取作某种概率分布或者概 率过程的未知参量，然后根据无作为标本抽出法对它进行统计、推定的一种方法。这个未知 值不一定是概率的值，也可以是推定的值[1] 。这种思想方法是一种在某种意义上可知的对策 (博弈)观点。随着方差减少法的开发和大型超高速计算机的出现，它不仅是仿真诸方法中最 重要的一种方法，而且也是系统工程诸方法中最有效的一种方法。 2. 提出问题 在数字通信中，对于数字调频信号，通常以频率250HZ 至4750HZ （每次加500HZ ）代 表0 至9 。现有一个五位码，例如数组[6，2 ，4 ，0，1]，调制成调频信号。接收到的信号中 混有 10 倍于原信号的强干扰，干扰信号为线性调频信号，其他噪声影响较小，可以忽略不 计，编程将信号解码，还原成原数组，应用蒙特卡罗法估计该方法的误码率。 3. 问题分析及实现方案 （1）生成实际接收信号：随机生成五位编码code，生成调频信号y1=cos （f/F*n ），F 为抽样频率。每个码元的长度有小范围随机抖动。录入强线性调频信号数据yint 。合成实际 信号y=y1+yint 。 （2 ）信号分析：将接收信号做时频分析，画出时频图，可见原信号与强干扰信号的频 域关系。并将信号按时间分段，加汉宁窗做傅立叶变换，得到五组频域数据，画成频域图。 由频域图可见强干扰信号对应一个高峰，原信号对应一个较矮的峰。 （3 ）信号处理：编程找出频谱图中信号的中心频率，并算出对应的码元PCODE ，即解 [2] 码过程 。 （4 ）误差估计：将解得的编码与原编码比较，计算误码率err 。 4. 数值分析及编程实现 （1）生成实际接收信号：编写函数[y，code]=tonegen （digits，yint ） 输入参数digits 为五位原编码，如果length(digits)5，则随机生成五位编码。输入参数 yint （取1111 个随机点）为干扰信号。默认值为强线性干扰信号。 函数返回值y 为生成的实际接收信号，code 为返回的原编码。 （2 ）分析接收信号：编写程序pic （digits ） -1- 调用tonegen 函数，返回实际接收信号和原编码。利用specgram 函数对实际信号做时频 分析，画出时频图，并将信号按时间分五段，加汉宁窗做傅立叶变换，得到五组频域数据， 画成五个频域图。 运行pic （0 ），得到PCODE=[0 ，3，8，0，1] （随机数） 时频图如图1 所示，五个频谱图如图2～图6 所示 图1 中桔色区域为能量集中区。可见强线性调频信号的阶梯状频谱，原编码信号的谱线 或与干扰信号谱线重叠（0，3 ）或远离干扰谱线（8，0，1）。 图1 时频图 图2 高峰对应干扰信号，矮峰对应码元0 图3 信号与干扰重叠，对应码元3 图4 矮峰对应码元8 图5 矮峰对应码元0