CN111814107A - 一种高精度实现平方根倒数的计算系统及其计算方法 - Google Patents
一种高精度实现平方根倒数的计算系统及其计算方法 Download PDFInfo
- Publication number
- CN111814107A CN111814107A CN202010662685.1A CN202010662685A CN111814107A CN 111814107 A CN111814107 A CN 111814107A CN 202010662685 A CN202010662685 A CN 202010662685A CN 111814107 A CN111814107 A CN 111814107A
- Authority
- CN
- China
- Prior art keywords
- calculating
- unit
- data
- square root
- high precision
- 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.)
- Granted
Links
Images
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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/25—Integrating or interfacing systems involving database management systems
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种高精度实现平方根倒数的计算系统及其计算方法,在数字信号处理领域,比如在信号归一化等运算处常有运用。包括判零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元,所述零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元依次连接。本发明使用牛顿迭代与Tylor公式结合的方法计算平方根倒数。使用Tylor公式计算初值,给到牛顿迭代。Tylor公式,以48个数据的存储,及1个乘法3个加法的代价,相较于直接存储表的方法,节省了700个数据的存储资源。使用牛顿迭代计算,能够直接将初值的精度提高1倍。使得最大计算误差由2^(‑10)提升到2^(‑20),性能优越。
Description
技术领域
本发明涉及数字信号处理技术领域,具体为一种高精度实现平方根倒数的计算系统及其计算方法。
背景技术
在数字信号处理领域,常常需要求取平方根的倒数,比如在信号归一化时。但现有的方案常基于牛顿迭代方案,进行多次迭代。牛顿迭代时,需要选取合适的初值,初值一般由查表得到,为了减少运算量,则需要表格足够大才行,占用存储空间。且利用Tylor公式也可以求取平方根倒数,需要占用一定的存储资源及运算量。
发明内容
本发明的目的在于提供一种求取平方根倒数,资源消耗较少,时延较低,方便硬件实现的高精度实现平方根倒数的计算系统及其计算方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:一种高精度实现平方根倒数的计算系统,包括判零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元,所述零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元依次连接。
一种高精度实现平方根倒数的计算系统的计算方法,包括以下步骤:
步骤一:待求数据输入判零及符号判断单元,若为零或为负数,直接返回异常,否则,输入下一步;
步骤二:将步骤一输出的数据输入缩放模块,实现将数据缩放到[1~4)之间的功能,输出数据x,并记下缩放因子;
步骤三:将步骤二输出的数据x输入Tylor计算单元,计算得到y1,输出y1;
步骤四:将步骤三输出的结果y1和步骤二输出的x,给到牛顿迭代单元;进行一次牛顿迭代;
步骤五:将步骤四中进行牛顿迭代后输出的y1与步骤二输出的缩放因子输入,更新y1=y1*缩放因子,并输出结果。
进一步,步骤三的计算过程包括以下步骤:
(1)根据当前输入数据x,计算表格中的index;
(2)根据index计算得到x0,查得x1,x2,并更新x0=x0-x;
(3)根据(2)得到的x0,x1计算y0;
(4)根据(3)得到的y0及(2)得到的x2,计算得到y1,输出y1。
进一步,步骤四的计算过程包括以下步骤:将步骤三输出的结果y1,2输出的x,给到牛顿迭代单元,进行一次牛顿迭代:
(1)计算y2=y1*x;
(2)计算y3=y2*y1;
(3)更新y2=1-y3;
(4)更新y3=y2*y1/2;
(5)更新y1=y1+y3。
进一步,所述缩放模块内设置有辅助计算用的索引。
进一步,所述判零及符号判断单元基于BaseControl程序。
进一步,所述高精度实现平方根倒数的计算系统可基于LTE接收机使用。
进一步,所述高精度实现平方根倒数的计算系统还设置有数据库,所述数据库采用SybaseSQLAnywhere10,所述数据库连有数据后台服务程序,所述数据后台服务程序通过ODBC接口与数据库连接。
与现有技术相比,本发明的有益效果是:
1本发明使用牛顿迭代与Tylor公式结合的方法计算平方根倒数。
2使用Tylor公式计算初值,给到牛顿迭代。Tylor公式,以48个数据的存储,及1个乘法3个加法的代价,相较于直接存储表的方法,节省了700个数据的存储资源。
3使用牛顿迭代计算,能够直接将初值的精度提高1倍。使得最大计算误差由2^(-10)提升到2^(-20),性能优越。
附图说明
图1为本发明的系统整体框图;
图2为本发明的Tylor计算单元的整体框图;
图3为本发明的牛顿迭代单元的整体框图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下述实施例中的方法,如无特别说明,均为常规方法。
实施例:
对某一LTE接收机,设最大接收天线为2,最大接收层数为2,当实际接收天线为2,接收层数为2时,设其接收信号为对应的信道估计为则由y=H*x+n,n为2*1矩阵,需要求取x的ML解。将方程两边同时左乘2*2的矩阵则并同时乘以归一化参数则易由下三角矩阵特点,计算出x。此时,需要用到平方根的倒数求取模块。将a=|h11|2+|h01|2输入该装置,可以直接得到高精度的
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (8)
1.一种高精度实现平方根倒数的计算系统,其特征在于,包括判零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元,所述零及符号判断单元、缩放模块、Tylor计算单元和牛顿迭代单元依次连接。
2.一种高精度实现平方根倒数的计算系统的计算方法,其特征在于,包括以下步骤:
步骤一:待求数据输入判零及符号判断单元,若为零或为负数,直接返回异常,否则,输入下一步;
步骤二:将步骤一输出的数据输入缩放模块,实现将数据缩放到[1~4)之间的功能,输出数据x,并记下缩放因子;
步骤三:将步骤二输出的数据x输入Tylor计算单元,计算得到y1,输出y1;
步骤四:将步骤三输出的结果y1和步骤二输出的x,给到牛顿迭代单元;进行一次牛顿迭代;
步骤五:将步骤四中进行牛顿迭代后输出的y1与步骤二输出的缩放因子输入,更新y1=y1*缩放因子,并输出结果。
3.根据权利要求2所述的一种高精度实现平方根倒数的计算系统的计算方法,其特征在于,步骤三的计算过程包括以下步骤:
(1)根据当前输入数据x,计算表格中的index;
(2)根据index计算得到x0,查得x1,x2,并更新x0=x0-x;
(3)根据(2)得到的x0,x1计算y0;
(4)根据(3)得到的y0及(2)得到的x2,计算得到y1,输出y1。
4.根据权利要求2所述的一种高精度实现平方根倒数的计算系统的计算方法,其特征在于,步骤四的计算过程包括以下步骤:将步骤三输出的结果y1,2输出的x,给到牛顿迭代单元,进行一次牛顿迭代:
(1)计算y2=y1*x;
(2)计算y3=y2*y1;
(3)更新y2=1-y3;
(4)更新y3=y2*y1/2;
(5)更新y1=y1+y3。
5.根据权利要求1所述的一种高精度实现平方根倒数的计算系统,其特征在于,所述缩放模块内设置有辅助计算用的索引。
6.根据权利要求1所述的一种高精度实现平方根倒数的计算系统,其特征在于,所述判零及符号判断单元基于BaseControl程序。
7.根据权利要求1所述的一种高精度实现平方根倒数的计算系统,其特征在于,所述高精度实现平方根倒数的计算系统可基于LTE接收机使用。
8.根据权利要求1所述的一种高精度实现平方根倒数的计算系统,其特征在于,所述高精度实现平方根倒数的计算系统还设置有数据库,所述数据库采用SybaseSQLAnywhere10,所述数据库连有数据后台服务程序,所述数据后台服务程序通过ODBC接口与数据库连接。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662685.1A CN111814107B (zh) | 2020-07-10 | 2020-07-10 | 一种高精度实现平方根倒数的计算系统及其计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010662685.1A CN111814107B (zh) | 2020-07-10 | 2020-07-10 | 一种高精度实现平方根倒数的计算系统及其计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111814107A true CN111814107A (zh) | 2020-10-23 |
CN111814107B CN111814107B (zh) | 2021-03-12 |
Family
ID=72842673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010662685.1A Active CN111814107B (zh) | 2020-07-10 | 2020-07-10 | 一种高精度实现平方根倒数的计算系统及其计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111814107B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020116431A1 (en) * | 1999-07-30 | 2002-08-22 | Ho Ying-Wai | System and method for improving the accuracy of reciprocal operations performed by a floating-point unit |
CN1472636A (zh) * | 2002-07-29 | 2004-02-04 | 矽统科技股份有限公司 | 浮点数的对数运算方法和装置 |
CN1838024A (zh) * | 2005-09-09 | 2006-09-27 | 威盛电子股份有限公司 | 对数处理系统和方法 |
CN101865963A (zh) * | 2010-07-08 | 2010-10-20 | 山东电力研究院 | 一种新的继电保护装置的快速检测方法 |
CN102143103A (zh) * | 2011-05-19 | 2011-08-03 | 新邮通信设备有限公司 | 一种导频序列的生成方法 |
CN103975564A (zh) * | 2011-10-27 | 2014-08-06 | Lsi公司 | 具有拥有由用户定义的用于数字预失真(dpd)以及其它非线性应用的非线性函数的指令集的处理器 |
US8812575B2 (en) * | 2010-07-06 | 2014-08-19 | Silminds, Llc, Egypt | Decimal floating-point square-root unit using Newton-Raphson iterations |
CN104483539A (zh) * | 2015-01-08 | 2015-04-01 | 湖南大学 | 一种基于泰勒展开式的有功功率快速测量方法 |
CN104615404A (zh) * | 2015-02-15 | 2015-05-13 | 浪潮电子信息产业股份有限公司 | 一种基于查表操作的高速浮点除法部件装置 |
CN106227291A (zh) * | 2016-07-26 | 2016-12-14 | 中国科学院自动化研究所 | 基于分段式查表法的反正切函数的实现方法及实现装置 |
CN110688088A (zh) * | 2019-09-30 | 2020-01-14 | 南京大学 | 一种面向神经网络的通用非线性激活函数计算装置和方法 |
-
2020
- 2020-07-10 CN CN202010662685.1A patent/CN111814107B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020116431A1 (en) * | 1999-07-30 | 2002-08-22 | Ho Ying-Wai | System and method for improving the accuracy of reciprocal operations performed by a floating-point unit |
CN1472636A (zh) * | 2002-07-29 | 2004-02-04 | 矽统科技股份有限公司 | 浮点数的对数运算方法和装置 |
CN1838024A (zh) * | 2005-09-09 | 2006-09-27 | 威盛电子股份有限公司 | 对数处理系统和方法 |
US8812575B2 (en) * | 2010-07-06 | 2014-08-19 | Silminds, Llc, Egypt | Decimal floating-point square-root unit using Newton-Raphson iterations |
CN101865963A (zh) * | 2010-07-08 | 2010-10-20 | 山东电力研究院 | 一种新的继电保护装置的快速检测方法 |
CN102143103A (zh) * | 2011-05-19 | 2011-08-03 | 新邮通信设备有限公司 | 一种导频序列的生成方法 |
CN103975564A (zh) * | 2011-10-27 | 2014-08-06 | Lsi公司 | 具有拥有由用户定义的用于数字预失真(dpd)以及其它非线性应用的非线性函数的指令集的处理器 |
CN104483539A (zh) * | 2015-01-08 | 2015-04-01 | 湖南大学 | 一种基于泰勒展开式的有功功率快速测量方法 |
CN104615404A (zh) * | 2015-02-15 | 2015-05-13 | 浪潮电子信息产业股份有限公司 | 一种基于查表操作的高速浮点除法部件装置 |
CN106227291A (zh) * | 2016-07-26 | 2016-12-14 | 中国科学院自动化研究所 | 基于分段式查表法的反正切函数的实现方法及实现装置 |
CN110688088A (zh) * | 2019-09-30 | 2020-01-14 | 南京大学 | 一种面向神经网络的通用非线性激活函数计算装置和方法 |
Non-Patent Citations (1)
Title |
---|
刘苍: ""面向未来无线通信基带信号处理的QR分解关键技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111814107B (zh) | 2021-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112506935B (zh) | 数据处理方法、装置、电子设备、存储介质及程序产品 | |
CN108021537A (zh) | 一种基于硬件平台的softmax实现方式 | |
CN108196822A (zh) | 一种双精度浮点开方运算的方法及系统 | |
CN111813372B (zh) | 一种高精度低时延实现32位整数除法的方法及装置 | |
EP4195084A1 (en) | Method and device for adjusting model parameters, and storage medium and program product | |
CN110187866B (zh) | 一种基于双曲cordic的对数乘法计算系统及方法 | |
CN110222305B (zh) | 一种基于双曲cordic的对数函数计算系统及方法 | |
CN102566965B (zh) | 一种误差平坦的浮点数对数运算装置 | |
CN111814107B (zh) | 一种高精度实现平方根倒数的计算系统及其计算方法 | |
CN108959571B (zh) | Sql语句的运算方法、装置、终端设备及存储介质 | |
CN109614073A (zh) | 四象限反正切函数硬件实现电路 | |
CN109697234B (zh) | 实体的多属性信息查询方法、装置、服务器和介质 | |
CN112821392A (zh) | 一种基于全纯嵌入法的静态电压稳定边界计算方法 | |
CN111984226A (zh) | 一种基于双曲cordic的立方根求解装置及求解方法 | |
CN108319804B (zh) | 低资源调用的8192点基2 dit asic设计方法 | |
CN112100877A (zh) | 一种结构刚度高效拓扑优化方法及其应用 | |
CN111752532B (zh) | 一种高精度实现32位整数除法的方法、系统及装置 | |
CN115080650A (zh) | 一种数据类型转换方法、装置、设备及存储介质 | |
CN104734728A (zh) | 基于mwc系统的盲多带稀疏信号快速恢复算法 | |
CN107977587A (zh) | 一种平方根功能实现方法 | |
CN113515259B (zh) | 一种适用于浮点格式的复数的近似取模实现电路及方法 | |
CN113961168A (zh) | 数据处理方法、装置、电子设备及存储介质 | |
JP2023103419A (ja) | 演算方法、装置、チップ、電子機器及び記憶媒体 | |
CN117148269A (zh) | 一种基于多级维纳滤波的低计算复杂度直接位置估计方法 | |
CN101685384B (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 |