CN110046327B - 一种基于帕德逼近的通信误差函数逼近方法 - Google Patents
一种基于帕德逼近的通信误差函数逼近方法 Download PDFInfo
- Publication number
- CN110046327B CN110046327B CN201910299354.3A CN201910299354A CN110046327B CN 110046327 B CN110046327 B CN 110046327B CN 201910299354 A CN201910299354 A CN 201910299354A CN 110046327 B CN110046327 B CN 110046327B
- Authority
- CN
- China
- Prior art keywords
- approximation
- function
- specific
- communication error
- upper bound
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/17—Function evaluation by approximation methods, e.g. inter- or extrapolation, smoothing, least mean square method
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
- H04B1/06—Receivers
- H04B1/10—Means associated with receiver for limiting or suppressing noise or interference
- H04B1/1027—Means associated with receiver for limiting or suppressing noise or interference assessing signal quality or detecting noise/interference for the received signal
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Signal Processing (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Algebra (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种基于帕德逼近的通信误差函数逼近方法。现有逼近方法或逼近效果不够精确,或算式复杂度较高。本发明方法首先设定上界函数U(x),令x点第i阶导数U(i)(x)=Q(i)(x),得到含有n个未知数n个方程的方程组,解出的ai带入U(x),整理后得到UQ(x);然后设定下界函数L(x).建立方程组解出的bj带入L(x),整理后得到LQ(x);得到的UQ(x)、LQ(x)将Q(x)上下紧密包围,达到逼近效果。本发明通过对于Q(x)及其积分的高精度估算,有助于提高通信信号系统的性能分析精度。通信误差函数的逼近方法可以在各种通信系统计算机仿真软件中得到应用。相比较于常规的查表方法,更加精确,计算通信误码率更加准确。
Description
技术领域
本发明属于通信与电子信息技术领域,具体是属于计算机仿真领域,涉及一种基于帕德逼近的通信误差函数逼近方法。
背景技术
在通信领域,通信误差高斯Q函数在数字调制方案的符号错误概率(SEP)、加性高斯白噪声等的估计中起着十分关键的作用。高斯Q函数的逼近问题,在信号处理和通信原理中有广泛的应用。对于加性高斯白噪音、衰落信道,多个数字调制方案的误码率估值中应用了高斯Q函数,对符号差错概率的逼近扮演重要角色。
Q函数,即标准正态分布的右尾函数,又叫(标准正态分布的)互补累计分布函数,其表达式为:
由于它与正态分布的累积分布函数的关系,Q函数也可以用误差函数表示。Q(x)不是初等函数,通常需要用初等函数来逼近,以方便计算。为了更好地控制相应的逼近误差,通常需要求解Q(x)的上下界。目前已有一些结果,比如: 其中φ(x)为标准正态分布的密度函数。当前对于逼近误差有着越来越高的要求,改善Q(x)的逼近误差,对于提高通信系统的性能分析变的越发重要。
在近年一个关于一维高斯Q函数近似值的调查中明确指出,这个问题是当前普遍存在的。即使在近三年内,也有许多文献根据期望获取的精度提供了简单的闭合公式逼近Q函数。Q函数本身不是闭式解,逼近Q函数则为Q函数提供了一种简单的形式,提供了进一步的数学分析,进而方便了各种通信系统的性能分析。值得注意的是,高斯概率密度函数对语音信号、无线接收机中的信号、正交频分复用技术调制信号进行了表征,Q函数的逼近问题在很多应用领域都有出现,具有重要意义。下面为一些已有的Q函数逼近公式:
上述逼近方法或逼近效果不够精确,或算式复杂度较高。
在通信领域,在数字调制方案的符号错误概率(SEP)、加性高斯白噪声等的估计中高斯Q函数起着十分关键的作用。理论上,可用高斯Q函数改进的近似N点高斯Hermite求积法则来估算,其中增加N的值可以提高精度,但将相应地增加计算复杂性及其计算成本。
此外,各种数字调制方案的SEP表达式涉及到积分,其被积函数是Q函数和衰落概率密度函数(PDF)的乘积,表示为:其中代表衰落概率密度函数,m定义了衰落参数,范围从0.5到∞,γ是平均信噪比,Γ(·)表示Gamma函数,a和b是依赖于数字调制技术的实正常数,k为Q(x)的次数。
即若要精确计算数字调制方案的SEP表达式需要预先获取Q(x)的原函数,然而由于Q(x)不是初等函数,难以计算,因此必须采用Q函数的逼近函数来计算。
发明内容
本发明目的就是针对现有技术的不足,提供一种基于帕德逼近的通信误差函数逼近方法,以更好地方便各种通信系统的性能分析。
本发明方法步骤如下:
步骤(1).设定U(x)为Q(x)的上界函数,令ai为上界待定未知系数,满足U(0)=0.5=Q(0);n为根据上界具体逼近精度需要而定的参数,n为偶数,n∈[4,20];
步骤(2).令U(i)(x)=Q(i)(x),i=1,2,···,n,U(i)(x)、Q(i)(x)分别表示U(x)、Q(x)在x点第i阶导数;由此得到含有n个未知数n个方程的方程组:U(i)(x)=Q(i)(x),i=1,2,···,n;
步骤(3).求解步骤(2)的方程组,得到ai,i=1,2,···,n,解出的ai是与x相关的参数方程,将ai相应的函数式带入U(x),整理得到UQ(x),UQ(x)为带入具体的n值得到的上界函数;
步骤(4).设定L(x)为Q(x)的下界函数,令bj为下界待定未知系数,满足L(0)=0.5=Q(0);m为根据下界具体逼近精度需要而定的参数,m为奇数,m∈[5,21];
步骤(5).建立方程组j=1,2,···,m-1;其中L′(0)、Q′(0)为L(x)、Q(x)在0点的1阶导数,L(j)(x)、Q(j)(x)为L(x)、Q(x)在x点的j阶导数;
步骤(6).解步骤(5)的方程组,得到bj,j=1,2,···,m,解出的bj为与x相关的参数方程,将bj相应的函数式带入L(x),整理得到LQ(x),LQ(x)为带入具体的m值得到的下界函数;
步骤(7).得到UQ(x)、LQ(x)满足即UQ(x)和LQ(x)将Q(x)上下紧密包围,达到逼近效果。
本发明通过对于Q(x)及其积分的高精度估算,有助于提高通信信号系统的性能分析精度。
本发明通过对Q(x)的高精度快速估算,以及Q(x)的相关积分的高精度估算等两个内容的研究,获取Q(x)逼近精度更高的上下界,以及Q(x)相关积分高精度的估算,进而改善通信系统的性能及分析精度等。从理论上说,n、m的值可以取更大,以获得逼近Q函数效果更好的函数。
利用本发明方法取得的逼近函数来计算各种数字调制方案中的误码率,计算中用到的积分函数即高斯Q函数和衰落概率密度函数结合的函数。这个函数的形式如下:
将取得的逼近函数带入到上述公式计算误码率,取得了很好的效果。通信误差函数的逼近方法可以在各种通信系统计算机仿真软件中得到应用。相比较于常规的查表方法(采用其它逼近方法得到数值表),本发明方法逼近效果更加精确,计算通信误码率更加准确,并且计算复杂度不高,可以获得更加优良的结果。
附图说明
图1为用maple程序所作的Q(x)-UQ(x),Q(x)-LQ(x)的图像;
图2为Q(x)-LQ(x),Q(x)-UQ(x)和Q(x)-Qk(x),(i=1,2,3,4,5,6)的图像。
具体实施方式
以下结合附图对本发明进行进一步说明。
一种基于帕德逼近的通信误差函数逼近方法,具体步骤如下:
(1).设定U(x)为Q(x)的上界函数,n为根据上界具体逼近精度需要而定的参数,n为偶数,n∈[4,20],本实施例n=10。从理论上说,从理论上说,为了得到逼近精度更好的函数,n也可以取更大的偶数,n=10时取得的函数已取得较好的逼近效果。
因为为使U(x)满足U(0)=0.5=Q(0),令/>ai为上界待定未知系数。
(2).为使U(x)逼近Q(x)的效果尽可能好,令U(i)(x)=Q(i)(x),i=1,2,···,10,U(i)(x)、Q(i)(x)分别表示U(x)、Q(x)在x点第i阶导数;即在x点令U(x)与Q(x)的i阶导数分别相等,由此得到含有10个未知数10个方程的方程组:U(i)(x)=Q(i)(x),i=1,2,···,10。
(3).求解步骤(2)的方程组,得到ai,i=1,2,···,10,解出的ai是与x相关的参数方程,将ai相应的函数式带入U(x),整理得到UQ(x),UQ(x)为带入具体的n值得到的上界函数:。
该步骤可以通过计算机求解,求解软件包括maple,mathematica,matlab等。
(4).设定L(x)为Q(x)的下界函数,m为根据下界具体逼近精度需要而定的参数,m为奇数,m∈[5,21],本实施例m=11。从理论上说,从理论上说,为了得到逼近精度更好的函数,m也可以取更大的奇数,m=11时取得的函数已取得较好的逼近效果。
同理,令bj为下界待定未知系数,满足L(0)=0.5=Q(0)。
(5).建立方程组j=1,2,···,10;其中L′(0)、Q′(0)为L(x)、Q(x)在0点的1阶导数,L(j)(x)、Q(j)(x)为L(x)、Q(x)在x点的j阶导数。
(6).解步骤(5)的方程组,得到bj,j=1,2,···,11,解出的bj为与x相关的参数方程,将bj相应的函数式带入L(x),整理得到LQ(x),LQ(x)为带入具体的m值得到的下界函数。。
该步骤可以通过计算机求解,求解软件包括maple,mathematica,matlab等。
(7).得到UQ(x)、LQ(x)满足即UQ(x)和LQ(x)将Q(x)上下紧密包围,达到逼近效果。
如图1和图2所示,通过maple程序对比Q(x)-LQ(x)、Q(x)-UQ(x)和Q(x)-Qk(x),(k=1,2,3,4,5,6)的大小,图像可以得出LQ(x)和UQ(x)的逼近误差明显小于其他函数,则逼近效果明显优于其他函数,已取得很好的逼近效果。
图2中:
利用上述所得到的逼近函数来计算各种数字调制方案中的误码率,计算中用到的积分函数即高斯Q函数和衰落概率密度函数结合的函数。这个函数的形式如下:
通信误差函数的逼近方法可以在各种通信系统计算机仿真软件中得到应用。相比较于常规的查表方法(采用其它逼近方法得到数值表),本方法逼近效果更加精确,计算通信误码率更加准确,并且计算复杂度不高,可以获得更加优良的结果。
Claims (1)
1.通过逼近函数提高数字调制方案中的误码率准确度的方法,其特征在于该方法具体步骤如下:
步骤(1).设定U(x)为Q(x)的上界函数,令ai为上界待定未知系数,满足U(0)=0.5=Q(0);n为根据上界具体逼近精度需要而定的参数,n为偶数,n∈[4,20];
步骤(2).令U(i)(x)=Q(i)(x),i=1,2,…,n,U(i)(x)、Q(i)(x)分别表示U(x)、Q(x)在x点第i阶导数;由此得到含有n个未知数n个方程的方程组:U(i)(x)=Q(i)(x),i=1,2,…,n;
步骤(3).求解步骤(2)的方程组,得到ai,i=1,2,…,n,解出的ai是与x相关的参数方程,将ai相应的函数式带入U(x),整理得到UQ(x),UQ(x)为带入具体的n值得到的上界函数;
步骤(4).设定L(x)为Q(x)的下界函数,令bj为下界待定未知系数,满足L(0)=0.5=Q(0);m为根据下界具体逼近精度需要而定的参数,m为奇数,m∈[5,21];
步骤(5).建立方程组其中L′(0)、Q′(0)为L(x)、Q(x)在0点的1阶导数,L(j)(x)、Q(j)(x)为L(x)、Q(x)在x点的j阶导数;
步骤(6).解步骤(5)的方程组,得到bj,j=1,2,…,m,解出的bj为与x相关的参数方程,将bj相应的函数式带入L(x),整理得到LQ(x),LQ(x)为带入具体的m值得到的下界函数;
步骤(7).得到UQ(x)、LQ(x)满足即UQ(x)和LQ(x)将Q(x)上下紧密包围,达到逼近效果;
步骤(8).将取得的逼近函数确定各种数字调制方案中的误码率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910299354.3A CN110046327B (zh) | 2019-04-15 | 2019-04-15 | 一种基于帕德逼近的通信误差函数逼近方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910299354.3A CN110046327B (zh) | 2019-04-15 | 2019-04-15 | 一种基于帕德逼近的通信误差函数逼近方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110046327A CN110046327A (zh) | 2019-07-23 |
CN110046327B true CN110046327B (zh) | 2023-08-04 |
Family
ID=67277076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910299354.3A Active CN110046327B (zh) | 2019-04-15 | 2019-04-15 | 一种基于帕德逼近的通信误差函数逼近方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110046327B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110765506B (zh) * | 2019-09-30 | 2023-03-31 | 杭州电子科技大学上虞科学与工程研究院有限公司 | 实体模型的多分辨率等几何拓扑优化方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103618325A (zh) * | 2013-12-06 | 2014-03-05 | 武汉大学 | 计及电力系统运行方式变化的广域阻尼控制器设计方法 |
CN107612656A (zh) * | 2017-09-18 | 2018-01-19 | 北京航空航天大学 | 一种适用于极化码的高斯近似简化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8630409B2 (en) * | 2011-04-05 | 2014-01-14 | International Business Machines Corporation | Two-party private estimation of dataset similarity |
-
2019
- 2019-04-15 CN CN201910299354.3A patent/CN110046327B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103618325A (zh) * | 2013-12-06 | 2014-03-05 | 武汉大学 | 计及电力系统运行方式变化的广域阻尼控制器设计方法 |
CN107612656A (zh) * | 2017-09-18 | 2018-01-19 | 北京航空航天大学 | 一种适用于极化码的高斯近似简化方法 |
Non-Patent Citations (1)
Title |
---|
汪芳宗 等.基于帕德逼近的暂态稳定性快速数值计算方法.《电力系统保护与控制 》.2017,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN110046327A (zh) | 2019-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20140142931A1 (en) | Pitch detection method and apparatus | |
CN110046327B (zh) | 一种基于帕德逼近的通信误差函数逼近方法 | |
CN111860276A (zh) | 人体关键点检测方法、装置、网络设备及存储介质 | |
WO2022134995A1 (zh) | 一种时延调整方法、装置、存储介质及电子装置 | |
CN113780523A (zh) | 图像处理方法、装置、终端设备及存储介质 | |
US7800783B2 (en) | Auto white balancing by modified grey world and geometric distortion measurement | |
CN109635349B (zh) | 一种噪声增强最小化克拉美罗界的方法 | |
CN110095753B (zh) | 一种基于到达角度aoa测距的定位方法及装置 | |
CN106817756B (zh) | 一种上行通道功率的控制方法及装置 | |
CN110649913B (zh) | 用于晶振的频率调整方法、装置、存储介质及电子设备 | |
CN102346909B (zh) | 图像噪声强度处理方法 | |
CN115032667A (zh) | 高动态导航信号射频基带一体化方法、装置、设备及介质 | |
CN103514589A (zh) | 图像恢复中的滤波方法 | |
CN110209851B (zh) | 模型训练方法、装置、电子设备及存储介质 | |
CN113808011A (zh) | 一种基于特征融合的风格迁移方法、装置及其相关组件 | |
US20120084335A1 (en) | Method and apparatus of processing floating point number | |
CN113984208B (zh) | 一种光谱仪波长标定方法、系统、介质及电子终端 | |
CN112671679B (zh) | 应用于simo系统的信道估计方法、装置及存储介质 | |
WO2022061700A1 (en) | Method, apparatus, electronic device and readable storage medium for estimation of parameter of channel noise | |
CN111615714A (zh) | 一种rgb数据的色彩调整方法 | |
Ying et al. | Image sharpness metric based on algebraic multi-grid method | |
WO2017028615A1 (en) | System and method for representing complex numbers in fused floating point | |
CN114969631B (zh) | 基带芯片、信道估计方法、数据处理方法及设备 | |
CN114866171B (zh) | 信号统计方法、装置、计算机设备和存储介质 | |
CN112308199B (zh) | 数据块的处理方法、装置及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |