[go: up one dir, main page]

CN116527207B - 一种北斗三号通信基带自适应迭代次数的turbo译码方法 - Google Patents

一种北斗三号通信基带自适应迭代次数的turbo译码方法 Download PDF

Info

Publication number
CN116527207B
CN116527207B CN202310810437.0A CN202310810437A CN116527207B CN 116527207 B CN116527207 B CN 116527207B CN 202310810437 A CN202310810437 A CN 202310810437A CN 116527207 B CN116527207 B CN 116527207B
Authority
CN
China
Prior art keywords
alpha
beta
delta
max
value
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
Application number
CN202310810437.0A
Other languages
English (en)
Other versions
CN116527207A (zh
Inventor
林仁杰
张勇鹏
杨华炜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujian Fuda Beidou Communication Technology Co ltd
Original Assignee
Fujian Fuda Beidou Communication Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Fujian Fuda Beidou Communication Technology Co ltd filed Critical Fujian Fuda Beidou Communication Technology Co ltd
Priority to CN202310810437.0A priority Critical patent/CN116527207B/zh
Publication of CN116527207A publication Critical patent/CN116527207A/zh
Application granted granted Critical
Publication of CN116527207B publication Critical patent/CN116527207B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/29—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes combining two or more codes or code structures, e.g. product codes, generalised product codes, concatenated codes, inner and outer codes
    • H03M13/2957—Turbo codes and decoding
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00—Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/35—Unequal or adaptive error protection, e.g. by providing a different level of protection according to significance of source information or by adapting the coding according to the change of transmission channel characteristics
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0408—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas using two or more beams, i.e. beam diversity
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0047—Decoding adapted to other signal detection operation
    • H04L1/005—Iterative decoding, including iteration between signal detection and decoding operation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00—Arrangements for detecting or preventing errors in the information received
    • H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045—Arrangements at the receiver end
    • H04L1/0054—Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Error Detection And Correction (AREA)

Abstract

本发明属于通信基带技术领域，特别涉及一种北斗三号通信基带自适应迭代次数的turbo译码方法，turbo译码方法由以下单元实现包含原始待译码码元信息、主码码元获取单元、分量码码元获取单元、当前分支度量值、前向分支度量值Alpha计算器等实现。本发明采用自适应迭代次数的turbo译码方法，实现接收机收到不同个数波束的情况下，译码器根据系统时钟主频率自动调整迭代次数，实现波束个数多时译码器能实现多波束通信，波束个数少时译码器可以增加迭代次数，实现性能最优化，且系统主频率和捕获波束个数计算最优迭代次数并用于控制译码器迭代次数，实现自适应迭代次数的turbo译码方法，大大平衡了系统资源和功耗之间的矛盾。

Description

一种北斗三号通信基带自适应迭代次数的turbo译码方法
技术领域
本发明属于通信基带技术领域，特别涉及一种北斗三号通信基带自适应迭代次数的turbo译码方法。
背景技术
随着卫星通信产业的不断发展，卫星通信核心模组中的基带技术，已经成为产业的关键核心技术，全世界都在为卫星无线通信领域基带技术找寻新的突破，尤其是北斗通信基带译码器算法，其设计性能、资源利用率、设计功耗等直接影响了整个北斗行业终端产品的性能；另外北斗通信基带译码器算法中，单个译码器能解调译码的信息数量，取决于系统时钟的主频率和译码迭代的次数，现有技术中turbo译码方法还不能够实现接收机收到不同个数波束的情况下，译码器根据系统时钟主频率自动调整迭代次数，实现波束个数多时译码器能实现多波束通信，波束个数少时译码器可以增加迭代次数，实现性能最优化的情况。
发明内容
（一）要解决的技术问题
为了克服上述现有技术中的不足，本发明提供一种北斗三号通信基带自适应迭代次数的turbo译码方法，能够实现接收机收到不同个数波束的情况下，译码器根据系统时钟主频率自动调整迭代次数，实现波束个数多时译码器能实现多波束通信，波束个数少时译码器可以增加迭代次数，实现性能最优化的情况。
（二）技术方案
本发明通过如下技术方案实现：本发明提出了一种北斗三号通信基带自适应迭代次数的turbo译码方法，turbo译码方法由以下单元实现包含原始待译码码元信息、主码码元获取单元、分量码码元获取单元、当前分支度量值、前向分支度量值Alpha计算器、后向分支度量值Beta计算器、Delta存储单元、Alpha存储单元、Beta存储单元、最大似然比Max_likehood计算器、外信息Le计算器、3GPP交织器、3GPP解交织器、最优递推迭代次数niter计算器实现，其所述turbo译码方法具体步骤如下：
1st)：根据接收的信息确定本帧为8Kbps信息帧、16Kbps信息帧、24Kbps信息帧，分别记为kind = 8Kbps ，kind = 16Kbps， kind = 24Kbps；
2nd)：根据系统的主频率f_sys和当前总接收波束个数cnt_b计算每个波束可以占用译码器的时钟周期个数为cyc_clk =f_sys/8/cnt_b；根据卫星导航电文为每个波束每125ms一帧循环的特点，通过系统主频率f和接收波束个数cnt计算在当前系统主频率下，可以分配给解调每一帧的时钟周期数cyc_clk = f/8/cnt；
例如系统主频率130MHz，则125ms共可以提供(130*10^6)/8个时钟周期，若此时接收波束为3个，则这个波束可以占用译码器的时钟周期个数为:
cyc_clk =(130*10^6)/8/3；
3rd)：查询时钟周期表，根据kind值查询相应的列，找到满足时钟周期表值小于cyc_clk的最大迭代次数，迭代索引niter值作为最优的迭代索引niter值；
4th)：主码码元获取待译码的码元信息和外信息La；
5th)：计算当前分支度量值Delta，并存储在Delta存储单元；
6th)：读取Delta存储单元的值计算前向分支度量值Alpha值，存储在Alpha存储单元；
7th)：读取Delta存储单元的值计算后向分支度量值Beta值，存储在Beta存储单元；
8th)：同时读取Delta存储单元、Alpha存储单元、Beta存储单元值计算最大似然比Max_likelyhood；
9th)：根据Max_likelyhood与输入的外信息La共同计算新的外信息Le；
10th)：Le经过3GPP标准交织器进行交织生成新的La作为副码或者分量码的外信息输入；副码也经过同样的过程输出新的外信息Le；
11th)：判断是否达到最优的迭代次数niter，是就结束，否则将Le经过3GPP解交织器进行交织后作为新一轮主码的外信息La输入，继续新的一轮循环，跳转步骤4th）。
该turbo译码方法进一步地包括有，所述步骤3rd）时钟周期表中的建立是从主码码元获取单元达到最优迭代次数过程中，译码器依次计算分别接收1帧8Kbps、16Kbps或24Kbps电文信息时，且经过不同迭代次数时，译码器所需要的时钟周期，并根据该信息建立时钟周期表，建立结果如图3所示。
该turbo译码方法进一步地包括有，所述步骤4th）第一次迭代外信息La默认为0。
该turbo译码方法进一步地包括有，所述步骤5th）中的分支度量值Delta为计最优递推迭代次数niter计算器计算，计算方式是根据turbo编码器建时钟计数值和实际符号值对应关系表评估输入码元符号的方式来实现；例如：第一个时钟周期计数值为000，对应实际符号值为000，首先假设当前获取的3个码元信息分别为X，A，B，则Delta = X+A+B；第二个时钟周期计数值为001，对应实际符号值为111，首先假设当前获取的3个码元信息分别为X，A，B，则Delta = -X-A-B；符号值0表示正，1表示负，建立结果如图4所示。
该turbo译码方法进一步地包括有，所述步骤6th）前向分支度量值Alpha值为前向分支度量值Alpha计算器计算，且采用迭代法来实现；具体步骤如下：
首先初始化第1列的8个Alpha值：
Alpha(1,1) = 0；
Alpha(1,2) = -Max_num；
Alpha(1,3) = -Max_num；
Alpha(1,4) = -Max_num；
Alpha(1,5) = -Max_num；
Alpha(1,6) = -Max_num；
Alpha(1,7) = -Max_num；
Alpha(1,8) = -Max_num；
当k＞1时，依次计算2~n组Alpha值，方法如下：
Alpha(k,1) = Max(Delta (k-1,1)+Alpha(k-1,1), Delta (k-1,2)+Alpha(k-1,2))；
Alpha(k,2) = Max(Delta (k-1,4)+Alpha(k-1,3), Delta (k-1,3)+Alpha(k-1,4))；
Alpha(k,3) = Max(Delta (k-1,5)+Alpha(k-1,5), Delta (k-1,6)+Alpha(k-1,8))；
Alpha(k,4) = Max(Delta (k-1,8)+Alpha(k-1,7), Delta (k-1,7)+Alpha(k-1,6))；
Alpha(k,5) = Max(Delta (k-1,2)+Alpha(k-1,3), Delta (k-1,1)+Alpha(k-1,2))；
Alpha(k,6) = Max(Delta (k-1,3)+Alpha(k-1,1), Delta (k-1,4)+Alpha(k-1,4))；
Alpha(k,7) = Max(Delta (k-1,6)+Alpha(k-1,5), Delta (k-1,5)+Alpha(k-1,6))；
Alpha(k,8) = Max(Delta (k-1,7)+Alpha(k-1,7), Delta (k-1,8)+Alpha(k-1,8))；
其中，Max_num为Alpha对应二进制位数下的最大值，例如10位对应Max_num=1023，16位对应Max_num=65535；n为本帧待译码信息电文总bit数。
该turbo译码方法进一步地包括有，所述步骤7th）后向分支度量值Beta值为后向分支度量值Beta计算器计算，且采用迭代法来实现；具体步骤如下；
首先初始化最后1列的8个Beta值：
Beta(n,1) = 0；
Beta(n,2) = -Max_num；
Beta(n,3) = -Max_num；
Beta(n,4) = -Max_num；
Beta(n,5) = -Max_num；
Beta(n,6) = -Max_num；
Beta(n,7) = -Max_num；
Beta(n,8) = -Max_num；
当k<n时，依次计算1~n-1组Beta值，方法如下：
Beta(k,1) = Max(Delta (k+1,1) + Beta(k+1,1), Delta (k+1,2) + Beta(k+1,5))；
Beta(k,2) = Max(Delta (k+1,1) + Beta(k+1,6), Delta (k+1,2) + Beta(k+1,1))；
Beta(k,3) = Max(Delta (k+1,3) + Beta(k+1,5), Delta (k+1,4) + Beta(k+1,2))；
Beta(k,4) = Max(Delta (k+1,3) + Beta(k+1,2), Delta (k+1,4) + Beta(k+1,6))；
Beta(k,5) = Max(Delta (k+1,5) + Beta(k+1,3), Delta (k+1,6) + Beta(k+1,4))；
Beta(k,6) = Max(Delta (k+1,5) + Beta(k+1,7), Delta (k+1,6) + Beta(k+1,3))；
Beta(k,7) = Max(Delta (k+1,7) + Beta(k+1,8), Delta (k+1,8) + Beta(k+1,7))；
Beta(k,8) = Max(Delta (k+1,7) + Beta(k+1,4), Delta (k+1,8) + Beta(k+1,8))；
其中，Max_num为Beta对应二进制位数下的最大值，例如10位对应Max_num=1023，16位对应Max_num=65535；n为本帧待译码信息电文总bit数。
该turbo译码方法进一步地包括有，所述步骤9th）最大似然比Max_likelyhood采用最大似然比Max_likelyhood计算器计算和对数运算计算，最大似然比Max_likelyhood计算器由对数运算器和线性计算器构成，其对数运算计算步骤如下：
获取当前信息bit位数且读地址，通过对数运算器分别对Delta、Alpha、Beta值做以10为底的对数运算，得到lg（Delta）、lg（Alpha）、lg（Beta），再然后利用线性组合器获取得到lg（Delta）+ lg（Alpha）+ lg（Beta）其值等于最大似然比Max_likelyhood±1%的误差值范围。
（三）有益效果
本发明相对于现有技术，具有以下有益效果：
本发明采用自适应迭代次数的turbo译码方法，实现接收机收到不同个数波束的情况下，译码器根据系统时钟主频率自动调整迭代次数，实现波束个数多时译码器能实现多波束通信，波束个数少时译码器可以增加迭代次数，实现性能最优化，且系统主频率和捕获波束个数计算最优迭代次数并用于控制译码器迭代次数，实现自适应迭代次数的turbo译码方法，大大平衡了系统资源和功耗之间的矛盾。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明的turbo译码方法流程图。
图2为本发明Max_likelyhood似然比的计算流程图。
图3为本发明时钟周期表表格图。
图4为本发明turbo编码器建时钟计数值和实际符号值对应关系表评估输入码元符号的表格图。
具体实施方式
本技术方案中：
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参照图1所示，本发明提出了一种北斗三号通信基带自适应迭代次数的turbo译码方法，turbo译码方法由以下单元实现包含原始待译码码元信息、主码码元获取单元、分量码码元获取单元、当前分支度量值、前向分支度量值Alpha计算器、后向分支度量值Beta计算器、Delta存储单元、Alpha存储单元、Beta存储单元、最大似然比Max_likehood计算器、外信息Le计算器、3GPP交织器、3GPP解交织器、最优递推迭代次数niter计算器实现，其所述turbo译码方法具体步骤如下：
1st)：根据接收的信息确定本帧为8Kbps信息帧、16Kbps信息帧、24Kbps信息帧，分别记为kind = 8Kbps ，kind = 16Kbps， kind = 24Kbps；
2nd)：根据系统的主频率f_sys和当前总接收波束个数cnt_b计算每个波束可以占用译码器的时钟周期个数为cyc_clk =f_sys/8/cnt_b；根据卫星导航电文为每个波束每125ms一帧循环的特点，通过系统主频率f和接收波束个数cnt计算在当前系统主频率下，可以分配给解调每一帧的时钟周期数cyc_clk = f/8/cnt；
例如系统主频率130MHz，则125ms共可以提供(130*10^6)/8个时钟周期，若此时接收波束为3个，则这个波束可以占用译码器的时钟周期个数为:
cyc_clk =(130*10^6)/8/3；
3rd)：查询时钟周期表，根据kind值查询相应的列，找到满足时钟周期表值小于cyc_clk的最大迭代次数，迭代索引niter值作为最优的迭代索引niter值；
4th)：主码码元获取待译码的码元信息和外信息La；
5th)：计算当前分支度量值Delta，并存储在Delta存储单元；
6th)：读取Delta存储单元的值计算前向分支度量值Alpha值，存储在Alpha存储单元；
7th)：读取Delta存储单元的值计算后向分支度量值Beta值，存储在Beta存储单元；
8th)：同时读取Delta存储单元、Alpha存储单元、Beta存储单元值计算最大似然比Max_likelyhood；
9th)：根据Max_likelyhood与输入的外信息La共同计算新的外信息Le；
10th)： Le经过3GPP标准交织器进行交织生成新的La作为副码或者分量码的外信息输入；副码也经过同样的过程输出新的外信息Le；
11th)： 判断是否达到最优的迭代次数niter，是就结束，否则将Le经过3GPP解交织器进行交织后作为新一轮主码的外信息La输入，继续新的一轮循环，跳转步骤4th）。
该turbo译码方法进一步地包括有，所述步骤3rd）时钟周期表中的建立是从主码码元获取单元达到最优迭代次数过程中，译码器依次计算分别接收1帧8Kbps、16Kbps或24Kbps电文信息时，且经过不同迭代次数时，译码器所需要的时钟周期，并根据该信息建立时钟周期表，建立结果如图3所示。
该turbo译码方法进一步地包括有，所述步骤4th）第一次迭代外信息La默认为0。
该turbo译码方法进一步地包括有，所述步骤5th）中的分支度量值Delta为计最优递推迭代次数niter计算器计算，计算方式是根据turbo编码器建时钟计数值和实际符号值对应关系表评估输入码元符号的方式来实现；例如：第一个时钟周期计数值为000，对应实际符号值为000，首先假设当前获取的3个码元信息分别为X，A，B，则Delta = X+A+B；第二个时钟周期计数值为001，对应实际符号值为111，首先假设当前获取的3个码元信息分别为X，A，B，则Delta = -X-A-B；符号值0表示正，1表示负，建立结果如图4所示。
该turbo译码方法进一步地包括有，所述步骤6th）前向分支度量值Alpha值为前向分支度量值Alpha计算器计算，且采用迭代法来实现；具体步骤如下：
首先初始化第1列的8个Alpha值：
Alpha(1,1) = 0；
Alpha(1,2) = -Max_num；
Alpha(1,3) = -Max_num；
Alpha(1,4) = -Max_num；
Alpha(1,5) = -Max_num；
Alpha(1,6) = -Max_num；
Alpha(1,7) = -Max_num；
Alpha(1,8) = -Max_num；
当k＞1时，依次计算2~n组Alpha值，方法如下：
Alpha(k,1) = Max(Delta (k-1,1)+Alpha(k-1,1), Delta (k-1,2)+Alpha(k-1,2))；
Alpha(k,2) = Max(Delta (k-1,4)+Alpha(k-1,3), Delta (k-1,3)+Alpha(k-1,4))；
Alpha(k,3) = Max(Delta (k-1,5)+Alpha(k-1,5), Delta (k-1,6)+Alpha(k-1,8))；
Alpha(k,4) = Max(Delta (k-1,8)+Alpha(k-1,7), Delta (k-1,7)+Alpha(k-1,6))；
Alpha(k,5) = Max(Delta (k-1,2)+Alpha(k-1,3), Delta (k-1,1)+Alpha(k-1,2))；
Alpha(k,6) = Max(Delta (k-1,3)+Alpha(k-1,1), Delta (k-1,4)+Alpha(k-1,4))；
Alpha(k,7) = Max(Delta (k-1,6)+Alpha(k-1,5), Delta (k-1,5)+Alpha(k-1,6))；
Alpha(k,8) = Max(Delta (k-1,7)+Alpha(k-1,7), Delta (k-1,8)+Alpha(k-1,8))；
其中，Max_num为Alpha对应二进制位数下的最大值，例如10位对应Max_num=1023，16位对应Max_num=65535；n为本帧待译码信息电文总bit数。
该turbo译码方法进一步地包括有，所述步骤7th）后向分支度量值Beta值为后向分支度量值Beta计算器计算，且采用迭代法来实现；具体步骤如下；
首先初始化最后1列的8个Beta值：
Beta(n,1) = 0；
Beta(n,2) = -Max_num；
Beta(n,3) = -Max_num；
Beta(n,4) = -Max_num；
Beta(n,5) = -Max_num；
Beta(n,6) = -Max_num；
Beta(n,7) = -Max_num；
Beta(n,8) = -Max_num；
当k<n时，依次计算1~n-1组Beta值，方法如下：
Beta(k,1) = Max(Delta (k+1,1) + Beta(k+1,1), Delta (k+1,2) + Beta(k+1,5))；
Beta(k,2) = Max(Delta (k+1,1) + Beta(k+1,6), Delta (k+1,2) + Beta(k+1,1))；
Beta(k,3) = Max(Delta (k+1,3) + Beta(k+1,5), Delta (k+1,4) + Beta(k+1,2))；
Beta(k,4) = Max(Delta (k+1,3) + Beta(k+1,2), Delta (k+1,4) + Beta(k+1,6))；
Beta(k,5) = Max(Delta (k+1,5) + Beta(k+1,3), Delta (k+1,6) + Beta(k+1,4))；
Beta(k,6) = Max(Delta (k+1,5) + Beta(k+1,7), Delta (k+1,6) + Beta(k+1,3))；
Beta(k,7) = Max(Delta (k+1,7) + Beta(k+1,8), Delta (k+1,8) + Beta(k+1,7))；
Beta(k,8) = Max(Delta (k+1,7) + Beta(k+1,4), Delta (k+1,8) + Beta(k+1,8))；
其中，Max_num为Beta对应二进制位数下的最大值，例如10位对应Max_num=1023，16位对应Max_num=65535；n为本帧待译码信息电文总bit数。
该turbo译码方法进一步地包括有，所述步骤9th）最大似然比Max_likelyhood采用最大似然比Max_likelyhood计算器计算和对数运算计算，最大似然比Max_likelyhood计算器由对数运算器和线性计算器构成，其对数运算计算步骤如下（图2所示）：
获取当前信息bit位数且读地址，通过对数运算器分别对Delta、Alpha、Beta值做以10为底的对数运算，得到lg（Delta）、lg（Alpha）、lg（Beta），再然后利用线性组合器获取得到lg（Delta）+ lg（Alpha）+ lg（Beta）其值等于最大似然比Max_likelyhood±1%的误差值范围；此处对Max_likelyhood做近似运算牺牲约0.2dB信噪比，但是大幅度降低了运算所需要的资源。
本发明turbo译码方法由以下单元实现包含有原始待译码码元信息、主码码元获取单元、分量码码元获取单元、当前分支度量值、前向分支度量值Alpha计算器、后向分支度量值Beta计算器、Delta存储单元、Alpha存储单元、Beta存储单元、最大似然比Max_likehood计算器、外信息Le计算器、3GPP交织器、3GPP解交织器、最优递推迭代次数niter计算器实现；且时钟周期表的建立已主码码元获取单元达到最优迭代次数过程中，译码器依次计算分别接收1帧8Kbps、16Kbps或24Kbps电文信息时，且经过不同迭代次数时，译码器所需要的时钟周期，并根据该信息建立的，且根据卫星导航电文为每个波束每125ms一帧循环的特点，通过系统主频率f和接收波束个数cnt计算在当前系统主频率下，可以分配给解调每一帧的时钟周期数cyc_clk = f/8/cnt；并根据当前系统参数计算的cyc_clk，匹配时钟周期表中该码率下满足cyc_clk<时钟周期的最大迭代次数niter值，作为次译码最优迭代次数；从而实现北斗卫星接收机在实际工作中，根据当前接收卫星波束的个数，实时调整译码器的最优迭代次数，从而达到一个译码解决所有波束译码且性能最优化的问题，且根据系统主频率和捕获波束个数计算最优迭代次数并用于控制译码器迭代次数，实现自适应迭代次数的turbo译码方法，大大平衡了系统资源和功耗之间的矛盾。
本发明中的分支度量值Delta计算器的计算方法采用建查找表评估输入码元符号的方式来实现、前向分支度量值ALpha计算器和后向分支度量值Beta计算器的计算均采用迭代法来实现，降低了计算所需的时钟周期数。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述turbo译码方法具体步骤如下：
1st)：根据接收的信息确定本帧为8Kbps信息帧、16Kbps信息帧、24Kbps信息帧，分别记为kind = 8Kbps ，kind = 16Kbps， kind = 24Kbps；
2nd)： 根据系统的主频率f_sys和当前总接收波束个数cnt_b计算每个波束可以占用译码器的时钟周期个数为cyc_clk =f_sys/8/cnt_b；
3rd)：查询时钟周期表，根据kind值查询相应的列，找到满足时钟周期表值小于cyc_clk的最大迭代次数，迭代索引niter值作为最优的迭代索引niter值；
4th)：主码码元获取待译码的码元信息和外信息La；
5th)：计算当前分支度量值Delta，并存储在Delta存储单元；
6th)：读取Delta存储单元的值计算前向分支度量值Alpha值，存储在Alpha存储单元；
7th)：读取Delta存储单元的值计算后向分支度量值Beta值，存储在Beta存储单元；
8th)：同时读取Delta存储单元、Alpha存储单元、Beta存储单元值计算最大似然比Max_likelyhood；
9th)：根据Max_likelyhood与输入的外信息La共同计算新的外信息Le；
10th)：Le经过3GPP标准交织器进行交织生成新的La作为副码或者分量码的外信息输入；副码也经过同样的过程输出新的外信息Le；
11th)：判断是否达到最优的迭代次数niter，是就结束，否则将Le经过3GPP解交织器进行交织后作为新一轮主码的外信息La输入，继续新的一轮循环，跳转步骤4th）；
所述步骤5th）中的分支度量值Delta为最优递推迭代次数niter计算器计算，计算方式是根据turbo编码器建立时钟计数值和实际符号值对应关系表评估输入码元符号的方式来实现。
2.根据权利要求1所述的一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述步骤3rd）时钟周期表中的建立是从主码码元获取单元达到最优迭代次数过程中，译码器依次计算分别接收1帧8Kbps、16Kbps或24Kbps电文信息时，且经过不同迭代次数时，译码器所需要的时钟周期，并根据该信息建立时钟周期表。
3.根据权利要求1所述的一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述步骤4th）第一次迭代外信息La默认为0。
4.根据权利要求1所述的一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述步骤6th）前向分支度量值Alpha值为前向分支度量值Alpha计算器计算，且采用迭代法来实现；具体步骤如下：
首先初始化第1列的8个Alpha值：
Alpha(1,1) = 0；
Alpha(1,2) = -Max_num；
Alpha(1,3) = -Max_num；
Alpha(1,4) = -Max_num；
Alpha(1,5) = -Max_num；
Alpha(1,6) = -Max_num；
Alpha(1,7) = -Max_num；
Alpha(1,8) = -Max_num；
当k＞1时，依次计算2~n组Alpha值，方法如下：
Alpha(k,1) = Max(Delta (k-1,1)+Alpha(k-1,1), Delta (k-1,2)+Alpha(k-1,2))；
Alpha(k,2) = Max(Delta (k-1,4)+Alpha(k-1,3), Delta (k-1,3)+Alpha(k-1,4))；
Alpha(k,3) = Max(Delta (k-1,5)+Alpha(k-1,5), Delta (k-1,6)+Alpha(k-1,8))；
Alpha(k,4) = Max(Delta (k-1,8)+Alpha(k-1,7), Delta (k-1,7)+Alpha(k-1,6))；
Alpha(k,5) = Max(Delta (k-1,2)+Alpha(k-1,3), Delta (k-1,1)+Alpha(k-1,2))；
Alpha(k,6) = Max(Delta (k-1,3)+Alpha(k-1,1), Delta (k-1,4)+Alpha(k-1,4))；
Alpha(k,7) = Max(Delta (k-1,6)+Alpha(k-1,5), Delta (k-1,5)+Alpha(k-1,6))；
Alpha(k,8) = Max(Delta (k-1,7)+Alpha(k-1,7), Delta (k-1,8)+Alpha(k-1,8))；
其中，Max_num为Alpha对应二进制位数下的最大值， n为本帧待译码信息电文总bit数。
5.根据权利要求1所述的一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述步骤7th）后向分支度量值Beta值为后向分支度量值Beta计算器计算，且采用迭代法来实现；具体步骤如下；
首先初始化最后1列的8个Beta值：
Beta(n,1) = 0；
Beta(n,2) = -Max_num；
Beta(n,3) = -Max_num；
Beta(n,4) = -Max_num；
Beta(n,5) = -Max_num；
Beta(n,6) = -Max_num；
Beta(n,7) = -Max_num；
Beta(n,8) = -Max_num；
当k<n时，依次计算1~n-1组Beta值，方法如下：
Beta(k,1) = Max(Delta (k+1,1) + Beta(k+1,1), Delta (k+1,2) + Beta(k+1,5))；
Beta(k,2) = Max(Delta (k+1,1) + Beta(k+1,6), Delta (k+1,2) + Beta(k+1,1))；
Beta(k,3) = Max(Delta (k+1,3) + Beta(k+1,5), Delta (k+1,4) + Beta(k+1,2))；
Beta(k,4) = Max(Delta (k+1,3) + Beta(k+1,2), Delta (k+1,4) + Beta(k+1,6))；
Beta(k,5) = Max(Delta (k+1,5) + Beta(k+1,3), Delta (k+1,6) + Beta(k+1,4))；
Beta(k,6) = Max(Delta (k+1,5) + Beta(k+1,7), Delta (k+1,6) + Beta(k+1,3))；
Beta(k,7) = Max(Delta (k+1,7) + Beta(k+1,8), Delta (k+1,8) + Beta(k+1,7))；
Beta(k,8) = Max(Delta (k+1,7) + Beta(k+1,4), Delta (k+1,8) + Beta(k+1,8))；
其中，Max_num为Beta对应二进制位数下的最大值， n为本帧待译码信息电文总bit数。
6.根据权利要求1所述的一种北斗三号通信基带自适应迭代次数的turbo译码方法，其特征在于：所述步骤9th）最大似然比Max_likelyhood采用最大似然比Max_likelyhood计算器计算和对数运算计算，最大似然比Max_likelyhood计算器由对数运算器和线性计算器构成，其对数运算计算步骤如下：
获取当前信息bit位数且读地址，通过对数运算器分别对Delta、Alpha、Beta值做以10为底的对数运算，得到lg（Delta）、lg（Alpha）、lg（Beta），再然后利用线性组合器获取得到lg（Delta）+ lg（Alpha）+ lg（Beta）其值等于最大似然比Max_likelyhood±1%的误差值范围。
CN202310810437.0A 2023-07-04 2023-07-04 一种北斗三号通信基带自适应迭代次数的turbo译码方法 Active CN116527207B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310810437.0A CN116527207B (zh) 2023-07-04 2023-07-04 一种北斗三号通信基带自适应迭代次数的turbo译码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310810437.0A CN116527207B (zh) 2023-07-04 2023-07-04 一种北斗三号通信基带自适应迭代次数的turbo译码方法

Publications (2)

Publication Number Publication Date
CN116527207A CN116527207A (zh) 2023-08-01
CN116527207B true CN116527207B (zh) 2023-08-29

Family

ID=87403290

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310810437.0A Active CN116527207B (zh) 2023-07-04 2023-07-04 一种北斗三号通信基带自适应迭代次数的turbo译码方法

Country Status (1)

Country Link
CN (1) CN116527207B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101674094A (zh) * 2009-10-16 2010-03-17 中山大学 一种Turbo码信道自适应实现方法
CN102148681A (zh) * 2011-05-19 2011-08-10 北京邮电大学 连续相位调制信号实现迭代定时同步的系统与方法
CN111262656A (zh) * 2020-01-16 2020-06-09 哈尔滨工业大学 低轨星地链路下高谱效稀疏码多址接入多用户迭代检测译码方法
CN115883298A (zh) * 2022-11-11 2023-03-31 中国船舶重工集团公司第七一五研究所 一种基于Haar分布域编码分集的水声通信方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7346833B2 (en) * 2002-11-05 2008-03-18 Analog Devices, Inc. Reduced complexity turbo decoding scheme
JP5479580B2 (ja) * 2009-06-18 2014-04-23 ゼットティーイー コーポレーション Ｌｔｅにおける並列ｔｕｒｂｏデコーディングの方法及び装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101674094A (zh) * 2009-10-16 2010-03-17 中山大学 一种Turbo码信道自适应实现方法
CN102148681A (zh) * 2011-05-19 2011-08-10 北京邮电大学 连续相位调制信号实现迭代定时同步的系统与方法
CN111262656A (zh) * 2020-01-16 2020-06-09 哈尔滨工业大学 低轨星地链路下高谱效稀疏码多址接入多用户迭代检测译码方法
CN115883298A (zh) * 2022-11-11 2023-03-31 中国船舶重工集团公司第七一五研究所 一种基于Haar分布域编码分集的水声通信方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种应用于卫星导航信号的Turbo译码方案及仿真;薛志芹等;《第十三届中国卫星导航年会》;第18-22页 *

Also Published As

Publication number Publication date
CN116527207A (zh) 2023-08-01

Similar Documents

Publication Publication Date Title
CN1937473A (zh) 随机存取多向cdma2000 turbo编码交织器
KR102547476B1 (ko) 경로 메트릭 값 기반의 디코딩 프로세스 제어 방법, 연산 장치 및 모바일 장치
US20090094505A1 (en) High speed turbo codes decoder for 3g using pipelined siso log-map decoders architecture
US6868518B2 (en) Look-up table addressing scheme
CN116527207B (zh) 一种北斗三号通信基带自适应迭代次数的turbo译码方法
US20030018942A1 (en) Common memory device and controlling method thereof
US20030091129A1 (en) Look-up table index value generation in a turbo decoder
US7340670B2 (en) Decoding apparatus, communication apparatus and decoding method
US20200028526A1 (en) Convolutional code decoder and convolutional code decoding method
US6886127B2 (en) Implementation of a turbo decoder
CN1988391A (zh) 运算电路
US7716554B2 (en) System and method for blind transport format detection with cyclic redundancy check
US20030023919A1 (en) Stop iteration criterion for turbo decoding
US6954507B2 (en) Method for estimating the bit error rate in a radio receiver and corresponding radio receiver
CN110768748B (zh) 回旋码解码器及回旋码解码方法
Chaikalis et al. Reconfigurable turbo decoding for 3G applications
CN100527647C (zh) 传输格式组合指示器解码电路及其解码方法
US6754503B1 (en) Method for adaptively selecting a soft symbol for a subsequent operation a communication device
CN119729619B (zh) 用于无线通信的方法、装置、芯片及终端
CN112564859A (zh) 一种软比特量化方法、装置、终端及存储介质
CN214707679U (zh) 一种基于MAP算法的卷积Turbo码译码器
CN112865814B (zh) 一种卷积码的Viterbi译码方法及译码器
CN110768747B (zh) 回旋码解码器及回旋码解码方法
US9136880B2 (en) Method for stopping iteration in an iterative turbo decoder and an iterative turbo decoder
CN116566403A (zh) 一种北斗三号通信基带可变信息率的turbo译码方法

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