CN1261475A - 借助卡尔曼滤波器源控制信道译码的方法和装置 - Google Patents
借助卡尔曼滤波器源控制信道译码的方法和装置 Download PDFInfo
- Publication number
- CN1261475A CN1261475A CN98806503A CN98806503A CN1261475A CN 1261475 A CN1261475 A CN 1261475A CN 98806503 A CN98806503 A CN 98806503A CN 98806503 A CN98806503 A CN 98806503A CN 1261475 A CN1261475 A CN 1261475A
- Authority
- CN
- China
- Prior art keywords
- reliability value
- value
- frame
- reliability
- binary digit
- 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
-
- 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/37—Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
- H03M13/39—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
- H03M13/41—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors
- H03M13/4138—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors soft-output Viterbi algorithm based decoding, i.e. Viterbi decoding with weighted decisions
- H03M13/4146—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors soft-output Viterbi algorithm based decoding, i.e. Viterbi decoding with weighted decisions soft-output Viterbi decoding according to Battail and Hagenauer in which the soft-output is determined using path metric differences along the maximum-likelihood path, i.e. "SOVA" 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/37—Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
- H03M13/3746—Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35 with iterative 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/37—Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
- H03M13/39—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
- H03M13/3961—Arrangements of methods for branch or transition metric calculation
-
- 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/37—Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
- H03M13/39—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
- H03M13/41—Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors
-
- 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/63—Joint error correction and other techniques
-
- 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/63—Joint error correction and other techniques
- H03M13/6312—Error control coding in combination with data compression
Landscapes
- Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Error Detection And Correction (AREA)
- Detection And Prevention Of Errors In Transmission (AREA)
- Complex Calculations (AREA)
Abstract
阐述用于处理通过接收机(16)接收数据的一种方法。借助于采用节拍增量的最大后验概率算法处理接收的数据。为了计算节拍增量,采用在使用卡尔曼滤波器的条件下在估算单元(20)中已求出的可靠性值(L(uk,l))。
Description
本发明涉及用于处理用接收器接收的数据的一种方法,在此方法上在相继的帧中接收经传输线段传输的数据。一个帧含有规定数目的二进制位。借助于采用节拍(Metrik)增量的最大后验概率算法或最大似然算法处理接收的数据,在此采用维特比算法。在此取决于对于此二进制位求出的一个现实的可靠性值来计算至少一个二进制位用的节拍增量,二进制位的值和/或二进制位的所谓的可靠性值从帧到帧相关。可靠性值是二进制位(1)的值具有一个规定的值的概率的一个尺度。
对于一种维特比算法例如从Joachim Hagenauer著的文章“Source-Coutrol led Channel Decoding(源控制信道译码)”在“IEEE通信学报”43卷,第9期,1995年9月,2449页至2457页中公开了一种这样的方法。在专利文献DE 42 24 214 C2中含有此文章的主要部分。尤其是文章的公式(21)说明新的节拍Mk(m)从旧的节拍Mk-1(m)和一个所属的节拍增量中的计算。在此m表示某个路径,而k表示一个帧中的某个二进制位。在计算节拍(Metrik)MK(m)时采用一个可靠性值L(uk)。在所述文章的2454和2455页上建议基于经验模型的所谓的Huk算法用于确定可靠性L(uk),在此模型上类似于在德国汽车责任保险上的所谓的点来计算可靠性值L(uk)。基于经验的特征可靠性值L(uk)的估算是不准确的。要么将节拍增量相加于已经计算的节拍上,要么相乘。
最大后验概率算法和最大似然算法是将经传输线段传输的数据序列与参考序列s相比较的和确定参考序列s的算法,此参考序列以最大的概率属于发明的数据序列。最大后验概率意味着,选择对其概率P(s|y)变成最大的参考序列s,在此y是接收的数据序列。最大似然性意味着,选择对其概率P(s|y)变为最大的参考序列s。传输路段例如是一个无线电信道,是一个传输线路或者甚至是从中读取数据的存储媒体。
本发明的任务在于提出用于求出可靠性值的一种简单的方法,此方法尽可能无差错和准确地计算现实可靠性值。
通过具有权利要求1特征的一种方法解决此任务。在归因于这些权利要求的从属权利要求中说明了有利的实施形式。
本发明根据的思路在于,用一种优化方法可以很准确地确定现实的可靠性值。可是困难在于,对于优化所采用的可靠性值基于伴有差错的观察之上。这些差错归因于在经传输线段传输数据时的干扰和维特比算法的伴有差错的结果。所以在按本发明的方法上,将一个帧用的基本上无差错的一个可靠值与为同一个帧求出的可靠性值的偏差的总和定义为目标函数。此偏差在此涉及在实际无差错可靠性值和求出的可靠性值之间的差值的数额。通过优化方法最小化此目标函数,使得为实施维特比算法求出的可靠性值很接近于无差错的可靠性值。
在按本发明的方法上尽管不知道无差错的可靠性值所以能够进行目标函数的最小化,因为在优化方法上采用一种数学模型,此模型在基于,在观察的可靠性值和无差错可靠性值之间建立一种联系。此外在模型中采用所谓的先验信息,意即已经在经传输线段传输数据之前已知的信息,例如像在相继帧中的关于数据相关性或可靠性值相关性的信息。在此意义上,数据源影响数据的其它处理。如果这种处理例如是一种译码的话,则人们也称之为源控制信道译码。
在无迭代的维特比算法上通常随着一个帧的延迟才存在基于观察的可靠性值。因此为了求出现实的可靠性值必须动用观察的可靠性值,从对于一个二进制位的一个已经处理的帧中已确定过此可靠性值。如果与此相反地进行带有迭代的维特比算法,则已经在第二个迭代步骤中存在对于瞬时处理帧的基于观察的一个可靠性值。在此情况下采用此观察的可靠性值用于求出现实的可靠性值。
相对于已知的Huk算法按本发明的方法提供较好的结果,因为用此方法可以处理强烈波动的源信号,在这些源信号情况下Huk算法是没有能力足够准确地估算可靠性值的。在按本发明的方法上,在进行维特比算法时采用的可靠性值接近于实际的可靠性值,使得在路径选择时降低差错的概率,因为所采用的可靠性值影响路径选择。例如在译码时正确的路径选择导致正确的译码的位串。
在本发明的一个进一步发展中按下式计算可靠性值:
在此ul是二进制位的优选可能为“+1”或“-1”的值。P(ul=+1)是所观察的二进制位1中的值是“+1”的概率。相应地P(ul=-1)是二进制位l中的值是“-1的概率。
在按本发明的一个其它特征的一种方法上,在优化时不用可靠性值直接计算,而是采用辅助可靠性值,下式优先适用于这些辅助可靠性值:m(ul)=P(ul=+1)-P(ul=-1)。
相对于可靠性值辅助可靠性值有这种优点,即这些辅助可靠性值对于在“0”和“1”之间的概率P有在“+1和“-1”之间的一个数值范围。与此相反地可靠性值在给定的函数范围中有从“-∞”到“+∞”的函数值,以致于给在有限的数值范围中的计算造成正如在机械计算机上存在的那样的困难。根据一种简单的数字关系进行从可靠性值向辅助可靠性值的,或从辅助可靠性值向可靠性值的换算,这种关系例如以表值形式存储到存储器中。
通过一种递归的优化可以减少优化时的计算工作量,在这种优化时从最后处理的帧用的可靠性值或辅助可靠性值中递归地计算现实的可靠性值或现实的辅助可靠性值。
为了优化,有利地采用在其上同样递归地进行优化的一种所谓的卡尔曼(Kalman)滤波器。例如在H.W.Sorenson著的文章“最小二乘方估算:从高斯到卡尔曼”中说明了一种这样的滤波器,在IEEE频谱,卷7,63至68页,1970年7月中阐述了此文章。基于卡尔曼滤波器的数学模型考虑这个事实,即只能间接地经基于伴随有差错的观察上的可靠性值求出实际的可靠性值。此外卡尔曼滤波器利用可靠性值和传输数据的统计性质,并且以足够的速度跟踪传输数据值的变化。
当在相继的帧的二进制位之间出现相关性时,则可有利地采用按本发明的方法。例如在一个帧中参数的较高值位用的二进制位上则是这种情况,按在移动无线电系统上的GSM标准传输此帧。在本发明方法的一个进一步发展中利用所述二进制位中的相关性用于这些二进制位的译码。
本发明此外涉及用于处理用接收机接收数据的和尤其是用于实施上述方法的一个装置。上述技术效应所以也适用于这种装置。
按本发明的一个其它特征,也通过具有权利要求15特征的一种方法解决所说明的任务。在此通过求平均值从至少两个观察的可靠性值中计算现实的可靠性值。这些措施基于这种认识,虽然每次观察的可靠值不与当时的实际可靠性值一致,可是通过求平均值相互消除观察时的差错。
在用求平均值的方法的一个进一步发展中仅针对一个规定数目的已经处理的帧来进行这种平均值计算,在此所规定的数目经越多次求平均值保持不变。通过这种措施,类似像通过一个窗口那样,仅优先考虑相继帧的序列的各自最后的区段。因此可以进行对变更着的条件的迅速适配。在递归地实施维特比算法时也将现实处理的帧计入窗口。
在用平均值计算的方法的一个另外的进一步发展上,采用例如由维特比算法处理的数据直接作为观察的可靠性值。所以这种措施导致现实可靠性值的一种很准确的计算,因为处理数据的平均值与辅助可靠性值一致,并且此平均值如已提及的那样对于实际无差错的可靠性值本身是一种良好的近似。
此外本发明涉及用于处理接收数据的和尤其是用于实施用按权利要求19的平均值计算的方法的一个装置。刚才所述的技术效应因此也适用于这个装置。
以下用附图阐述本发明的实施例。在其中:
图1展示具有数据传输用的主要功能单元的方框图,
图2展示两个由数据组成的应传输的帧,
图3展示卡尔曼滤波器中辅助可靠性值计算的一种第一方式用的示意图,
图4展示卡尔曼滤波器中辅助可靠性值计算的一种第二方式用的示意图,
图5展示来自GSM全价语音编码器的不同参数用的相关性系数的图表,和
图6展示一个GSM全价编码帧中的二进制位0至69用的相关性系数的图表。
图1展示具有数据传输用的主要功能单元的方框图。在此主要针对一个帧k的传输,在这方面为了清晰起见在图1中和在此图的以下阐述时仅在帧k的区分对于本发明是重要的情况下随同写上采用自然数的下标k。下面用图2阐述帧k的构造。
在发射机8中生成的一个源字符序列{q1′}由源字符q1′组成,这些源字符取决于要发射的信息,例如具有值“+1”和“-1”。下标1′对于在帧k中传输的源字符q1′从0进行至L′-1,在此L′是每个帧k的源字符q1′的数目。通过源编码器10,例如用GSM全价语音编码器压缩源字符序列{q1′}。在此由源编码的字符ul生成一个源编码序列{ul}。源编码的符号ul要到具有值“+1”,要么具有值“-1”。在一个帧之内下标l从0进行到L-1,在此L是一个帧中的源编码字符ul的数目。在此L′通常是大于L的。
然后在信道编码器12中防御信道干扰地编码源编码的序列{ul},在此例如采用一种卷积码。在此从码字x1中形成一种信道编码的序列{x1,n}。在码字x1之内通过下标n标记二进制位,下标在一个码字x1之内从零进行至N-1,在此N是一个码字x1中的二进制位的数目。码字x1的二进制位x1,n又具有要么数值“+1”,要么数值“-1”。在一个未表示的调制器中,进一步处理随后经传输线段14传输信道编码的序列{x1,n}。在传输时出现干扰,例如由衰落系数ak描述的衰落和由噪声系数No描述的噪声。
传输线段14位于发射机8和接收机16之间。接收机16必要时含有用于接收经传输线段14传输的信号的二个未表示的天线、一个扫描装置、用于解调信号的一个解调器和用于消除字符间干扰的校正器。出于简化原因,在图1中同样未表示这些装置。校正器输出接收序列{y1,n}的接收值y1,n。由于经传输线段14传输时的干扰,接收值y1,n拥有偏离“+1”和“-1”的,例如为“+0.2”或“-3.7”的值。
在信道译码器18中进一步处理接收值y1,n。在信道译码器18中重新取消由信道编码器12进行的卷积编码,在此应纠正传输差错。在卷积译码时例如采用已知的维特比算法。为了进行维特比算法也给信道译码器18输送一个信道状态信息Lc1,n。在信道译码器18中进行维特比算法时按以下的公式计算二进制位1用的路径m的一个所谓的节拍M1 (m): 在此M1-1 (m)是路径m的旧的节拍,意即对于二进制位1-1,
是属于路径m和二进制位1的码字的位,
是属于码字
的已译码的字符和L(uk,1)是求出的可靠性值,此可靠性值是对于源编码的字符ul等于“-1”的那个概率的一个尺度。
在下面用图3和4阐述其构造的一个计算单元20中求出可靠性值L(uk,1)。
隐含地确定信道状态信息Lc1,n,其办法是求出一个所谓的信道软件输出,此信道软件输出相当于由信道状态信息Lc1,n和各自接收值y1,n的乘积。对于此信道软件输出适用于关系式:LCy=L(x/y)-L(x) (2),在式中L(x/y)是一个可靠性值,它说明,当已接收了接收值y时,在码字的各自的二进制位中以何种概率出现值x,并且L(x)是一个可靠性值,它说明可以如何可靠地确定值x。Lc=4aEs/No适用于一种所谓的高斯信道/衰落信道,式中a为衰落系数,而Es/No为信噪比。
信道译码器18从源编码的或信道译码的接收字符
中生成一个接收的源编码的序列
。一个观察的可靠性值
属于每个源编码接收字符
,此可靠性值是信道译码器18的一个所谓的软输出(Soft-Ausgabe)。可靠性值L(uk-1,1)是对于用其已可以通过信道译码器18确定的源编码接收字符
的可靠性的一个尺度。为了区别于由计算单元20计算的可靠性值L(uk,1),以下将由信道译码器18生成的可靠性值称为基于观察的可靠性值或称为观察的可靠性值。通过信道译码器18以一个帧k的延迟才生成观察的可靠性值
。通过下标k-1延迟变得明显了。在以下阐述的本发明的实施例上,在计算单元20中采用观察的可靠性值
,以便求出现实处理帧k用的可靠性值L(uk,1)。
终止在一个加法单元28上的两个虚线箭头24和26应指明,在计算单元20中在可靠性值
之外,或者替代这些可靠性值,采用源编码的接收字符
和/或在译码时在源译码器22中生成的信息。计算单元20最后使得成为可能的是,接收的源字符序列
基本上与发射的源字符{q′1}一致,并且因此足够良好地纠正传输时的信道干扰的影响。
图2展示两个帧,即由源编码的字符ul组成的现实的帧k和先前的帧k-1。L个源编码的接收字符ul属于每个帧k,k-1,使得下标l从0进行至L-1。通过按图1的计算单元20才建立在帧k,k-1,k-2等等之间的一种联系。在此利用这个事实,即在相继帧k-1,k的某些源编码字符之间,例如字符uk,1和uk-1,1之间,或它们的可靠性值L(uk,1)和L(uk-1,1)之间出现一个明显的相关性。在本发明上通过一个数学模型描述此相关性,在模型的基础上求出可靠性值L(uk,1)。
如果对于属于接收字符
的状态来计算节拍增量M1 (m),则各自考虑按公式(1)的可靠性值L(uk-1,1)。在此时刻可靠性值
还是未知的。如果在此之后对于属于接收字符
的状态求出节拍M1 (m),计算单元20(请参阅图1)则提供从在此时刻已知的可靠性值
中计算出的一个可靠性值L(uk,1)。因此尽管延迟达一个帧,直到生成基于观察的可靠性值
时仍有帧k用的,具有足够精度的一个可靠性值L(uk,1)供支配。
图3展示计算卡尔曼滤波器中辅助可靠性值的一种第一方式的示意图。在首先详述了卡尔曼滤波器之后,在下面阐述辅助可靠性值。
以下采用一种书写方式,这种书写方式依据在阐述H.W.Sorenson著的文章“最小二乘方估算:从高斯到卡尔曼”中阐述卡尔曼滤波器时采用的书写方式,请比较IEEE频谱,卷7,63至68页,1970年7月。在方程(11′)和(10′)中采用的数学模型是多维的。为了阐述本发明的一个实施例却动用一维模型的较简单的情况。在此情况下矢量只有一个分量,并且因此是标量。产生以下的系统模型:xk=ρkxk-1+xk (3),式中xk为在时刻k的状态、xk-1为在时刻k-1的状态、ρk为xk和xk-1之间的相关系数,以及wk为在时刻k的系统干扰。采用此系统模型用于模拟所谓的一级马尔科夫过程。这样的系统模型用的,另外经常采用的一种关系是“一级自动回归的模型”。
所属的模型为:k=xk+xk (4),式中zk为在时刻k的测量值和vk为在时刻k的测量干扰。应与xk相乘的测量系数Hk具有值“1”,并且所以未随同写入公式(4)。干扰wk和vk通过所谓的白噪声形成,并且具有平均值零,此外wk的方差称为Qk,而vk的方差称为Rk。
根据公式(3)和(4)的模型,卡尔曼滤波器算法是与H.W.Sorenson一致地通过以下的公式定义的: Pk/k=(1-Kk)Pk/k-1 (9),式中
为状态xk用的估算,而ρk为差错
的协方差。在下标上第一下标说明各自的时刻k或k-1,在此时刻计算或已计算下标的量。由一斜杠分开的下标说明,用于计算各自已作下标的量所必要的值在何时刻供支配。如果第二下标对于理解是不必要的话,以下也采用在其上略去第二下标的书写方式。
方差Qk,Rk以及相关系数ρk优先是常数,并且在滤波开始时规定。在规定时有利地考虑事先已测量的测量值。
通过将公式(8)代入公式(5)产生
用的一种所谓的预测估算;
基于在处理时刻k+1存在的测量值zk+1的估算xk+1/k+1,通过将公式(5)代入公式(8)中产生:
通过数学的换算,在这些换算时在公式(7)和(9)中借助公式(6)代替协方差Pk/k-1,并且在这些换算时随后消除协方差Pk-1/k-1,产生计算放大系数K用的以下的公式:
从公式(12)中得出,当测量干扰vk与系统干扰wk相比是很小的和Rk因此大大小于Qk时,放大系数Kk大约取值1。如果放大系数Kk具有值1,则按公式(10)和(11)只通过现实的测量值zk或zk+1确定优化估算值
。如果与此相反地系统干扰wk与测量值vk相比是微小的话,放大系数Kk则趋于零。从方程(10)和(11)中得出结论,对于约为零的放大系数Kk所得出最佳估算值
仅与以前的估算值
有关。卡尔曼滤波器的这些性能保证无差错的实际状态xk的良好估算。
也可以如下地列出方程(11): 式中hj是与zj(j=1,…,k+1)无关的滤波器常数。估算值
因此是所有迄今存在的测量值zk的一种线性组合。从方程(10)、(11)和(13)中可以获知,在估算值
中既考虑测量值zk的统计相关性(一级统计)也考虑测量值zk的数值分布(零级统计)。
当由L(uk,1)替代
和由
替代zk时,可采用上面说明的公式直接用于可靠性值L(uk,1)的计算。可是出于下述原因合理的是,采用辅助可靠性值m(uk,1)和
代替可靠性值L(uk,1)和
,例如由可靠性值L(uk,1)与一个概率P的强烈非线性关系的线性化得出这些辅助可靠性值。
对于可靠性值通常适用以下的定义:
在此将ul看作为具有元素{+1,-1}的随机变量。P(ul=-1)=p是对于ul具有值“-1”的概率。与此相反地P(ul=+1)=1-p是对于ul具有值“+1”的概率。缩写“log”表示自然对数。可靠性值也称为软值或称为对数概率比。
按定义(14)可靠性值是[-∞,+∞]范围中的一个实数。如果用值L(ul)将ul分级为“+1或“-1,代数符号sign(L(ul))则说明所谓的硬抉择,而绝对值|L(ul)|则说明抉择的,也称为软值(Soft-Wert)的可靠性。
如果观察用于计算节拍MK (m)的,上面所说明的公式(1)的话,当可靠性值趋向于“+∞”或“-∞”时,则单独通过L(uk,1)确定此公式。正如从公式(14)中可获知的那样,这对于P趋向于零或P趋向于1是这种情况。此外在用处理器实现时,对数函数的转换常常是碍事的。在本发明的一个实施例上为了消除这些缺点,将在p和L(uk)之间的关系线性地近似:L(uk)≈
L(uk)≡K(1-2p)=Km(uk) (15).式中
L(uk)为近似的函数,K是一个常数,和m(uk)是已经提及的辅助可靠性值。常数K具有数值“z”,使得函数L(uk)的按概率p对于概率P=0.5的导数,等于函数L(uk)在位置p=0.5上按概率p的导数。
对于辅助可靠性值m(uk)适用以下的关系:m(uk)=1-2p=(-1)p+(+1)(1-p)=E{uk} (16).
从此公式中可以获知,辅助可靠性值m(uk)是uk的平均值。
从公式(14)和(16)中得出在辅助可靠性值和可靠性值之间的联系为: 或者为:
辅助可靠性值m(uk)具有像可靠性值L(uk)相似的意义,意即代数符号sign(m(uk))是所谓的硬抉择,而绝对值|m(uk)|是抉择的可靠性(所谓的软抉择)。按公式(15)或(17),例如借助于一个表格将由信道译码器18(请参阅图1)生成的观察的可靠性值
转换为一个辅助可靠性值m(uk-1,1)。然后同样地在采用一个表格的条件下,将由计算单元20(请参阅图1)求出的辅助可靠性值m(uk,1)转换为可靠性值L(uk,1)。按方程(15)通过辅助可靠性值m(uk)来近似可靠性值L(uk),在实践中同样导致良好的结果。
在用辅助可靠性值mk的按本发明的一个实施例上,在此未再随同写上与uk的关系,从方程(10)和(12)中对于卡尔曼滤波器得出以下的计算规范: 式中σ2 w,k相当于已经提及的方差Qk,而σ2 v,k相当于同样提及的方差Rk。
图3在一个方框图中展示按公式(19)的计算。用一个减法单元50计算在公式(19)中处于内括号中的表达式(m* k-1-mk-1/k-2)。在一个乘法单元52中用按公式(20)已计算的放大系数Kk-1秉此结果。在一个加法单元54中将辅助可靠性值mk-1/k-2叠加到相乘的结果上,以此确定了公式(19)中在外括号之内的表达式。随后在一个乘法单元56中将此表达式的计算结果与相关系数ρk相乘。相乘的结果是辅助可靠性值mk/k-1,如已提及的那样,在这个可靠性值由信道译码器18(请参阅图1)采用之前,将此辅助可靠性值mk/k-1转换为一个可靠性值。
延迟单元58从辅助可靠性值mk/k-1中通过延迟一个帧生成下一个处理步骤用的辅助可靠性值mk-1/k-2。图3中提出了一个乘法单元60,在此乘法单元中,在减法单元50中采用辅助可靠性mk-1/k-2之前,将此辅助可靠性值与测量系数Hk-1相乘。如上面所提及的那样,在实施例中测量系数Hk-1却具有值1。
在有关辅助可靠性值的一个另外的实施例中,在时刻k辅助可靠性值m* k已经是已知的,例如当进行迭代的译码时。在此情况下对于卡尔曼滤波器与公式(11)和(12)一致地适用以下两个公式:mk/k=ρkmk-1/k-1+Kk(m* k-ρkmk-1/k-1) (21)
图4展示公式(21)的计算。图4中的描述基本上符合图3中的描述,以致于不再次地阐述已经阐述过的功能单元,可是为了区别起见,所属的标记号获得一个位于上角的斜撇。与图3相比,在图4中下标k-1由下标k所代替。
作为初始条件在按图3或图4的实施例中例如选择:mo/-1=mo/o=0和k0=0。根据测量值来规定方差σ2 w,k和σ2 v,k。同样地根据例如借助于下述的图6求出的测量值来规定相关系数ρk。也可以在滤波过程期间按照规定的,必要时实现确定的函数,动态地计算方差σ2 w,k和σ2 v,k及相关系数ρk。
在图3和4中各自展示了帧k的二进制位1用的一个卡尔曼滤波器(请参阅图2)。通常有从帧到帧相关的多个二进制位1。在此情况下在计算单元20中对于这些二进制位1中的每一个采用按图3或按图4的一个卡尔曼滤波器,请参阅图1。
图5展示一个图表,在此图表上在横座标上展示了对于表示在纵座标上的参数号0至75的相关系数。具有号0至75的参数相当于由GSM全价语音编码器生成的参数,请参阅GSM建议06.10“欧洲数字蜂窝式电信系统;全价语音代码转换”,1995。GSM全价语音编码器例如是源编码器10,请参阅图1。
参数0至7是在LPC分析(线性预测编码)时生成的所谓的LAR系数(对数地区比例Logaritbmical area ratio)。这些系数的相关性是整个地大于0.3的。甚至是所谓的LTP(长期预测)用的相似度b的参数号9,26,43和60具有大于0.2的相关性。此外在每个帧中存在着来自所谓RPE分析(规则脉冲激发)的,从帧到帧仅少量变化的和具有大于0.7的相关系数的四个xMAX系数。所述参数的相关性也出现于在其中传输这些参数的较高值的二进制位中。
图6展示在一个GSM帧中的二进制位0至69用的相关系数的一个图表。许多二进制位的值具有一种可观的帧间相关性,意即在uk和uk-1之间的相关系数大于0.3。借助于在方程(16)中展示的平均值计算可以看出,在m(uk)和m(uk-1)之间的相关系数ρk必须更显著地大于uk和uk-1之间的相关系数,对于所述的值典型地为0.8至0.9,使得公式(3)的数学模型符合实际的情况,并且用卡尔曼滤波器生成的可靠性值L(uk,1)导致译码过程的改善。
在本发明的一个其它实施例中采用一个多维的模型代替在公式(3)和(4)中说明的一维的模型。在此情况下适用在已提到的Sorensen著的文章中说明的矢量书写方式。尤其当在帧间相关性之外也考虑帧内相关性时,则采用多维的模型。
也可以不用优化方法,用稍为有限的但是仍然对于许多用途是足够的精度求出可靠性值或辅助可靠性值的计算结果。在此利用上面在公式(16)中说明的联系,其办法是按下列公式求出辅助可靠性值m(uk): 式中N为帧的规定数目,经这些帧进行一种所谓的窗口化,意即只考虑最后的N帧的源编码的接收字符
通过这种措施也掌握应估算信号的迅速变化。
借助于辅助可靠性值的定义,即m=1-2p也可如下地表达公式(23):
Claims (19)
1.用于处理用接收器(16)接收的数据的方法,
在此方法上在相继的帧(K)中接收经传输线段(14)传输的数据,在此一个帧(K)含有规定数目的二进制位(1),
和在此方法上借助采用节拍增量的最大后验概率算法或最大似然算法处理接收的数据,
在此取决于对于现实处理的帧(K)求出的一个现实可靠性值(L(uk,1))地,来计算其值和/或其可靠性值(L(uk,1))从帧(k)到帧(k-1)相关的至少一个二进制位(1)用的节拍增量,此可靠性值是对于二进制位(1)的值有一个规定值的概率的一个尺度,其特征在于,采用至少一个观察的可靠性值
用于求出现实的可靠性值(L(uk,1)),从已经处理的一个帧(k-1)中或从二进制位(1)用的现实处理的帧(k)中确定此可靠性值,
和如此求出现实的可靠性值(L(uk,1)),使得对于多个帧(k)形成的,各自帧(k)用的无差错可靠性值的和对于同一个帧(k)求出的可靠性值(L(uk,1))的偏差的和变成最小。
2.按权利要求1的方法,其特征在于,按关系式: 求出可靠性值,式中L(ul)为可靠性值,ul表示帧之内的二进制位1的值,P(ul=+1)是对于二进制位1中的值有一个第一值的,优先有数值“+1”的概率,并且在此P(ul=-1)是对于二进制位1中的值有一个第二值的,优先有数值“-1的概率。
3.按权利要求1或2的方法,其特征在于,回归地,尤其在采用一个卡尔曼滤波器(20)的条件下从对于最后处理的帧(k-1)求出的可靠性值(L(uk,1))中计算现实的可靠性值(L(uk,1))。
4.按权利要求3的方法,其特征在于,按下列公式计算现实可靠性: 和/或按下列公式计算现实的可靠性值: 式中Lk为帧k用的可靠性值,ρk为规定的相关系数,L* k-1或L* k为帧k-1或k用的观察的可靠性值,和在此Kk或Kk-1是按下列公式回归地计算的放大系数: 式中σw,K-1为系统干扰的方差,而σ2 v,k-1为测量干扰的方差。
5.用于处理通过接收机(16)接收的数据的方法,
在此方法上在相继的帧(k)中接收经传输线段(14)传输的数据,在此一个帧(k)含有规定数目的二进制位(1),
和在此方法上借助于采用节拍增量的最大后验概率算法或最大似然算法处理接收的数据,
在此取决于对于现实处理的帧(k)求出的现实可靠性值(L(uk,1))地,来计算其值和/或可靠性值(L(uk,1))从帧(k)到帧(k-1)相关的至少一个二进制位(1)用的节拍增量,此现实的可靠性值是对于可以无差错地确定二进制位(1)的值的概率的一个尺度,
其特征在于,采用至少一个观察的可靠性值
用于求出现实的可靠性值(L(uk,1)),从一个已经处理的帧(k-1)中或从二进制位(1)用的现实处理的帧(k)中确定此观察的可靠性值,
如此求出现实的辅助可靠性值(m(uk,1)),使得对于多个帧(k)形成的,各自帧(k)用的无差错的辅助可靠性值和对于同一个帧(k)求出的辅助可靠性值(m(uk,1))的偏差的和变成最小,
和在采用换算函数或近似函数的条件下,从求出的辅助可靠性值(m(uk,1))中确定现实的可靠性值(L(uk,1))。
6.按权利要求5的方法,其特征在于,对于辅助可靠性值适用下列公式:m(ul)=P(ul=+1)-P(ul=-1),在此ul表示帧之内的二进制位1的值,P(ul=+1)是对于值ul有一个第一值,优先有数值“+1”的概率,而在此P(ul=-1)是对于值ul有一个第二值,优先有“-1的概率。
7.按权利要求5或6的方法,其特征在于,从对于最后处理的帧(k-1)求出的辅助可靠性值(m(uk-1,1))中回归地,尤其在采用一个卡尔曼滤波器(20)的条件下计算现实的辅助可靠性值(m(uk,1))。
8.按权利要求7的方法,其特征在于,按下列公式计算现实辅助可靠性值: 和/或按下列公式计算现实概率值mk:mk=ρkmk-1+Kk(m* k-ρkmk-1)式中mk为帧k用的辅助可靠性值,ρk为相关系数,m*k或m*k-1为观察的辅助可靠性值,和Kk或Kk-1为放大系数,优先按下列公式回归地确定此放大系数: 式中σ2 w,k-1为系统干扰的方差,而σ2 v,k-1为测量干扰的方差。
9.按权利要求8的方法,其特征在于,固定地规定相关系数ρk,系统干扰的方差σ2 w和测量干扰的方差σ2 v,在此在预先规定时优先考虑测量。
10.按以上权利要求之一的方法,其特征在于,在求取现实可靠性值时采用以下的数学模型:Xk=Φk,k-1Xk-1+Wk和Zk=HkXk+Vk式中Xk表示帧k的一个状态矢量,φk表示一个系统矩阵,Wk表示一个系统干扰矢量,Zk表示一个测量矢量,Hk表示一个测量矩阵,和Vk表示一个测量干扰矢量。
11.按以上权利要求之一的方法,其特征在于,采用此方法用于译码已编码地传输的数据。
12.按以上权利要求之一的方法,其特征在于,在优先按GSM标准工作的移动无线电系统中数据传输时采用此方法。
13.按权利要求12的方法,其特征在于,在二进制位(1)中优先传输在帧(k)中传输的参数的一个较高值的位。
14.用于处理接收数据,尤其是用于实施按以上权利要求之一的方法的装置(16),
具有用于接收经传输线段(14)传输的数据的一个接收单元,在由规定数目的二进制位组成的帧(k,k-1)中传输这些数据,
具有借助于采用节拍增量的最大后验概率算法或最大似然算法处理接收的数据的一个节拍增量单元(18),
和具有一个计算单元(20),此计算单元对于其值和/或可靠性值从帧(k)到帧(k-1)相关的至少一个二进制位(1),求出输出到节拍增量单元(18)上的一个现实可靠性值(L(uk,1)),
并且在于,计算单元(20)如此求取可靠性值(L(uk,1)),使得对于多个帧(k)形成的,一个帧的无差错可靠性值和对于同一个帧(k)求出的可靠性值(L(uk,1))的偏差的和变成最小,
或者在于,计算单元(20)如此求取可靠性值(L(uk,1)),使得对于多个帧(k)形成的,各自帧(k)用的基本上无差错的辅助可靠性值和对于同一个帧(k)求出的辅助可靠性值(m(uk,1))的偏差的和变成最小,在此借助于换算函数或近似函数从可靠性值(L(uk,1))中计算辅助可靠性值(m(uk,1))。
15.用于处理用接收机(16)接收的数据的方法,在此方法上在相继的帧(k)中接收经传输线段(14)传输的数据,在此一个帧(k)含有规定数目的二进制位(1),
和在此方法上借助于采用节拍增量的最大后验概率算法或最大似然算法处理接收的数据,
在此取决于对于现实处理的帧(k)求出的现实可靠性值(L(uk,1);m(uk,1))地,来计算其值和/或可靠性值从帧(k)到帧(k-1)相关的至少一个二进制位(1)用的节拍增量,此可靠性值是对于二进制位(1)有一个规定值的概率的一个尺度,
其特征在于,为了求取现实可靠性值(L(uk,1);m(uk,1))采用至少两个观察的可靠性值
或采用从二进制位(1)用的所属的帧(k,k-1)中确定的至少两个已经求出的可靠性值(L(uk,1),L(uk-1,1);m(uk,1),m(uk-1,1)),
和通过从观察的可靠性值
中的或从已经求出的可靠性值(L(uk,1),L(uk-1,1);m(uk,1),m(uk-1,1))中的平均值计算形成现实的可靠性值(L(uk,1),m(uk,1))。
16.按权利要求15的方法,其特征在于,对于一个恒定数目的已经处理的帧(k)形成平均值,在此优先采用最后处理的帧(k)的可靠性值。
17.按权利要求15或16的方法,其特征在于,可靠性值是辅助可靠性值,为了计算节拍增量用换算函数或近似函数将这些辅助可靠性值换算成可靠性值。
19.用于处理接收数据,尤其是用于实施按权利要求15至18之一的方法的装置(16),
具有用于接收经传输线段(14)传输的,在由规定数目的二进制位(1)组成的帧(k)中传输的数据,
具有借助于采用节拍增量的最大后验概率算法或最大似然算法处理接收的数据的一个节拍增量单元(18),
和具有一个计算单元(20),此计算单元对于其值和/或可靠性值从帧(k)到帧(k-1)相关的至少一个二进制位(1)求出输出到节拍增量单元(18)上的一个现实的可靠性值(L(uk,1))的一个计算单元(20),在此可靠性值(L(uk,1))是对于二进制位(1)的值有一个规定值的概率的一个尺度,
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE19726600 | 1997-06-23 | ||
DE19726600.2 | 1997-06-23 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1261475A true CN1261475A (zh) | 2000-07-26 |
CN1109409C CN1109409C (zh) | 2003-05-21 |
Family
ID=7833376
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN98806503A Expired - Fee Related CN1109409C (zh) | 1997-06-23 | 1998-06-04 | 借助卡尔曼滤波器源控制信道译码的方法和装置 |
Country Status (8)
Country | Link |
---|---|
US (1) | US6393075B1 (zh) |
EP (2) | EP1133063B1 (zh) |
JP (1) | JP2002505051A (zh) |
CN (1) | CN1109409C (zh) |
AU (1) | AU737227B2 (zh) |
DE (2) | DE59809365D1 (zh) |
ES (2) | ES2206364T3 (zh) |
WO (1) | WO1998059423A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104949093A (zh) * | 2014-03-26 | 2015-09-30 | 西门子公司 | 用于蒸汽发电厂的蒸汽产生器的状态观察器 |
CN104949094A (zh) * | 2014-03-26 | 2015-09-30 | 西门子公司 | 用于蒸汽发电厂的蒸汽产生器的多参量状态调节器 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6226318B1 (en) * | 1998-03-31 | 2001-05-01 | Zenith Electronics Corporation | Detection of code vectors in single frequency, multiple transmitter networks |
GB9824012D0 (en) * | 1998-11-03 | 1998-12-30 | Lucent Technologies Inc | Improved source controlled channel decoding |
EP1224783A2 (en) * | 2000-04-06 | 2002-07-24 | Lucent Technologies Inc. | Upstream data transmission |
US6798852B2 (en) | 2000-06-06 | 2004-09-28 | Ericsson Inc. | Methods and systems for extracting a joint probability from a map decision device and processing a signal using the joint probability information |
US6697441B1 (en) | 2000-06-06 | 2004-02-24 | Ericsson Inc. | Baseband processors and methods and systems for decoding a received signal having a transmitter or channel induced coupling between bits |
US6944206B1 (en) * | 2000-11-20 | 2005-09-13 | Ericsson Inc. | Rate one coding and decoding methods and systems |
US7139334B2 (en) * | 2001-06-21 | 2006-11-21 | Bartlett Alan M | Cooperative code-enhanced multi-user communications system |
US9331717B2 (en) * | 2009-02-27 | 2016-05-03 | Blackberry Limited | Forward error correction decoding avoidance based on predicted code block reliability |
CN104270119B (zh) * | 2014-09-22 | 2017-05-17 | 武汉理工大学 | 基于非线性未知随机偏差的两阶段容积卡尔曼滤波方法 |
US10958293B1 (en) * | 2020-03-02 | 2021-03-23 | GM Global Technology Operations LLC | System and method for near-lossless universal data compression using correlated data sequences |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4644564A (en) * | 1983-08-05 | 1987-02-17 | International Business Machines Corporation | Decoding the output signal of a partial-response class-IV communication or recording device channel |
EP0449327B1 (en) * | 1990-03-30 | 1998-07-15 | Nec Corporation | Noise-immune space diversity receiver |
US5432821A (en) * | 1992-12-02 | 1995-07-11 | University Of Southern California | System and method for estimating data sequences in digital transmissions |
US5673291A (en) * | 1994-09-14 | 1997-09-30 | Ericsson Inc. | Simultaneous demodulation and decoding of a digitally modulated radio signal using known symbols |
DE69526644T2 (de) * | 1994-10-31 | 2003-01-30 | Koninkl Philips Electronics Nv | Digitales übertragungs- und aufzeichungsystem mit einfacher fehlerkorrektur |
US6581179B1 (en) * | 1996-06-25 | 2003-06-17 | Ericsson Inc. | Methods for generating side information in the presence of time-selective fading |
US5953376A (en) * | 1996-09-26 | 1999-09-14 | Lucent Technologies Inc. | Probabilistic trellis coded modulation with PCM-derived constellations |
US6192503B1 (en) * | 1997-08-14 | 2001-02-20 | Ericsson Inc. | Communications system and methods employing selective recursive decording |
-
1998
- 1998-06-04 EP EP01108652A patent/EP1133063B1/de not_active Expired - Lifetime
- 1998-06-04 WO PCT/DE1998/001524 patent/WO1998059423A1/de active IP Right Grant
- 1998-06-04 CN CN98806503A patent/CN1109409C/zh not_active Expired - Fee Related
- 1998-06-04 DE DE59809365T patent/DE59809365D1/de not_active Expired - Fee Related
- 1998-06-04 ES ES01108652T patent/ES2206364T3/es not_active Expired - Lifetime
- 1998-06-04 DE DE59801798T patent/DE59801798D1/de not_active Expired - Fee Related
- 1998-06-04 EP EP98934858A patent/EP0992116B1/de not_active Expired - Lifetime
- 1998-06-04 ES ES98934858T patent/ES2166173T3/es not_active Expired - Lifetime
- 1998-06-04 AU AU84328/98A patent/AU737227B2/en not_active Ceased
- 1998-06-04 JP JP50353899A patent/JP2002505051A/ja active Pending
-
1999
- 1999-12-23 US US09/472,059 patent/US6393075B1/en not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104949093A (zh) * | 2014-03-26 | 2015-09-30 | 西门子公司 | 用于蒸汽发电厂的蒸汽产生器的状态观察器 |
CN104949094A (zh) * | 2014-03-26 | 2015-09-30 | 西门子公司 | 用于蒸汽发电厂的蒸汽产生器的多参量状态调节器 |
US10267512B2 (en) | 2014-03-26 | 2019-04-23 | Siemens Aktiengesellschaft | Multi-variable state closed-loop control for a steam generator of a thermal power plant |
Also Published As
Publication number | Publication date |
---|---|
CN1109409C (zh) | 2003-05-21 |
ES2206364T3 (es) | 2004-05-16 |
JP2002505051A (ja) | 2002-02-12 |
AU737227B2 (en) | 2001-08-16 |
EP0992116A1 (de) | 2000-04-12 |
ES2166173T3 (es) | 2002-04-01 |
DE59809365D1 (de) | 2003-09-25 |
AU8432898A (en) | 1999-01-04 |
WO1998059423A1 (de) | 1998-12-30 |
EP1133063B1 (de) | 2003-08-20 |
US6393075B1 (en) | 2002-05-21 |
EP1133063A1 (de) | 2001-09-12 |
EP0992116B1 (de) | 2001-10-17 |
DE59801798D1 (de) | 2001-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1109409C (zh) | 借助卡尔曼滤波器源控制信道译码的方法和装置 | |
CN1227816C (zh) | 有效的格子结构状态量度归一化 | |
CN1153399C (zh) | 用于提高当前语音帧质量的方法和无线系统 | |
US7978101B2 (en) | Encoder and decoder using arithmetic stage to compress code space that is not fully utilized | |
US7808406B2 (en) | Method and apparatus for realizing arithmetic coding/decoding | |
CN1421085A (zh) | 基于干扰消除及解码的组合软决策的方法和装置 | |
CN1198861A (zh) | 软判定卷积解码器中置信度及帧信号质量的检测 | |
CN110708079B (zh) | 一种极化码构造方法及装置 | |
CN1205759C (zh) | 用于计算回声测量值的方法和装置 | |
CN109525253B (zh) | 基于深度学习和集成方法的卷积码译码方法 | |
CN1901430A (zh) | 以循环冗余检查来进行盲目传输格式检测的系统与方法 | |
US7284185B2 (en) | Puncturing/depuncturing using compressed differential puncturing pattern | |
CN1983334A (zh) | 一种实现算术编解码的方法及装置 | |
CN1189264A (zh) | 降低了复杂度的信号传输系统 | |
CN1234898A (zh) | 具有改进语音编码器和解码器的发射机 | |
US6363341B1 (en) | Encoder for minimizing resulting effect of transmission errors | |
KR20040044589A (ko) | 다수결 논리를 이용한 rm 부호의 연판정 복호 방법 및그 장치 | |
RU2704722C2 (ru) | Перестановочный декодер с обратной связью | |
CA2054849C (en) | Speech parameter encoding method capable of transmitting a spectrum parameter at a reduced number of bits | |
EP3780607A1 (en) | Information maintaining encoding and decoding method and device | |
CN1338824A (zh) | 一种特博码解码的方法以及解码器 | |
CN1914670A (zh) | 预测编码方法 | |
CN112994842B (zh) | 一种软比特的量化方法、装置、解调器及计算机存储介质 | |
KR100590184B1 (ko) | 채널적응벡터양자화(covq)를 위한 코드북의 구현방법 | |
RU2672300C2 (ru) | Перестановочный декодер с памятью |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
ASS | Succession or assignment of patent right |
Owner name: NOKIA SIEMENS COMMUNICATION CO., LTD. Free format text: FORMER OWNER: SIEMENS AG Effective date: 20080404 |
|
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20080404 Address after: Munich, Germany Patentee after: Nokia Siemens Networks GmbH Address before: Munich, Germany Patentee before: Siemens AG |
|
C19 | Lapse of patent right due to non-payment of the annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |