发明名称 一种改变批处理数据块块长的均衡系统及方法
摘要 本发明涉及一种改变批处理数据块块长的均衡系统及方法,所述均衡系统包括:批处理模块,预变换模块,频域输出信号模块,时域输出信号模块,批处理数据块块长变换模块和频域权系数变换模块;所述均衡方法是根据块长变化决定因子ΔL,更新第k+1次批处理数据块块长L<SUB>k+1</SUB>,并根据更新的批处理数据块块长L<SUB>k+1</SUB>,更新频域权系数W(k+1),所述更新频域权系数包括一次更新频域权系数W’(k+1)和二次更新频域权系数W(k+1)。本发明的优点在于:计算量小速度快,能应用于实时要求高的系统中;在矛盾的收敛速度和稳态误差之间找到了平衡点,既能够获得很快的收敛速度又能够获得较低的稳态误差;本方法实现了变块长在工程方面的应用。
申请公布号 CN101364840A 申请公布日期 2009.02.11
申请号 CN200810223644.1 申请日期 2008.09.28
申请人 中国科学院声学研究所 发明人 曾学文;杨群;肖琳;刘马飞;王劲林
分类号 H04B7/26(2006.01);H04Q7/22(2006.01) 主分类号 H04B7/26(2006.01)
代理机构 北京法思腾知识产权代理有限公司 代理人 杨小蓉
主权项 1、一种改变批处理数据块块长的均衡系统,所述均衡系统包括:一批处理模块,用于将接收到的时域输入信号分块处理,形成批处理数据块;一预变换模块,用于将第k批批处理数据块之前的、长度与该数据块块长相同的时域输入信号,与该数据块顺序组合,形成预变换数据块;一频域输出信号模块,用于对预变换数据块做快速傅立叶变换,通过计算输出频域信号;一时域输出信号模块,用于对输入的频域信号做快速反傅立叶变换,通过计算输出时域信号;其特征在于,还包括:一批处理数据块块长变换模块,用于根据块长变化决定因子ΔL,更新第k+1次批处理数据块块长L<sub>k+1</sub>;若块长变化决定因子ΔL的值大于设定的阈值范围的最大值,则L<sub>k+1</sub>=2L<sub>k</sub>,若块长变化决定因子ΔL的值小于设定的阈值范围的最小值,则<maths num="0001"><![CDATA[<math><mrow><msub><mi>L</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>L</mi><mi>k</mi></msub><mo>,</mo></mrow></math>]]></maths>否则L<sub>k+1</sub>=L<sub>k</sub>;其中,块长变化决定因子ΔL的阈值范围应该根据参数α、ρ<sub>1</sub>、ρ<sub>2</sub>和实际信道环境来做设定;一频域权系数变换模块,用于更新频域权系数W(k+1),所述更新频域权系数包括一次更新频域权系数W'(k+1)和二次更新频域权系数W(k+1);所述二次更新频域权系数根据更新的批处理数据块块长L<sub>k+1</sub>和一次更新频域权系数得到;所述块长变化决定因子ΔL,通过下述公式计算得到:如果e<sub>av,k</sub>=0、e<sub>av,k</sub>-e<sub>av,k-1</sub>≠0,则<maths num="0002"><![CDATA[<math><mrow><mi>&Delta;L</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mfrac><msub><mi>&rho;</mi><mn>2</mn></msub><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mfrac><mo>;</mo></mrow></math>]]></maths>如果e<sub>av,k</sub>-e<sub>av,k-1</sub>=0,e<sub>av,k</sub>=0,则L<sub>k+1</sub>=2L<sub>k</sub>;如果e<sub>av,k</sub>≠0,e<sub>av,k-1</sub>=0,则<maths num="0003"><![CDATA[<math><mrow><mi>&Delta;L</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mfrac><msub><mi>&rho;</mi><mn>2</mn></msub><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub></mfrac><mo>;</mo></mrow></math>]]></maths>当e<sub>av,k</sub>≠0,e<sub>av,k-1</sub>≠0时,如果e<sub>av,k</sub>&gt;=e<sub>av,k-1</sub>,则<maths num="0004"><![CDATA[<math><mrow><mi>&Delta;L</mi><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mfrac><msub><mi>&rho;</mi><mn>1</mn></msub><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub></mfrac><mo>-</mo><mi>&alpha;</mi><mfrac><mrow><msub><mi>&rho;</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mfrac><mo>;</mo></mrow></math>]]></maths>如果e<sub>av,k</sub>&lt;e<sub>av,k-1</sub>,则<maths num="0005"><![CDATA[<math><mrow><mi>&Delta;L</mi><mo>=</mo><mi>&alpha;</mi><mfrac><msub><mi>&rho;</mi><mn>1</mn></msub><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub></mfrac><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mi></mi><mfrac><msub><mi>&rho;</mi><mn>2</mn></msub><mrow><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>其中,参数α表示平衡系数,用于控制增长因子和缩短因子在块长变化决定因子ΔL中的比例,其取值范围是(0,1);参数ρ<sub>1</sub>表示增长步长,用于控制增长因子的步长,其取值范围是(0,1);参数ρ<sub>2</sub>表示缩短步长,用于控制缩短因子的步长,其取值范围是(0,1);所述参数α、ρ<sub>1</sub>、ρ<sub>2</sub>的值由外部给定;其中,e<sub>av,k-1</sub>、e<sub>av,k</sub>分别表示第k-1、k批次的平均绝对误差的值,e<sub>av,k</sub>通过下述公式计算得到:<maths num="0006"><![CDATA[<math><mrow><msub><mi>e</mi><mrow><mi>av</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>L</mi><mi>k</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>L</mi><mi>k</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover><mrow><mo>|</mo><msub><mi>e</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow></mrow></math>]]></maths>其中,e<sub>k</sub>(i)是e(k)=[e<sub>k</sub>(0),e<sub>k</sub>(1),...,e<sub>k</sub>(L<sub>k</sub>-1)]<sup>T</sup>式中的e<sub>k</sub>(0),e<sub>k</sub>(1),...,e<sub>k</sub>(L<sub>k</sub>-1);e(k)是误差信号,e(k)=d(k)—y(k);d(k)是时域输出信号值第k批的判决值,d(k)=[d<sub>k</sub>(0),d<sub>k</sub>(1),...,d<sub>k</sub>(L<sub>k</sub>-1)]<sup>T</sup>;所述一次更新频域权系数通过下述公式得到:<maths num="0007"><![CDATA[<math><mrow><msup><mi>W</mi><mo>,</mo></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><mfrac><msub><mi>&mu;</mi><mi>B</mi></msub><msub><mi>L</mi><mi>k</mi></msub></mfrac><mo>*</mo><mi>FFT</mi><mrow><mo>[</mo><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>]</mo></mrow></mrow></math>]]></maths>其中,μ<sub>B</sub>是块步长因子;Φ(k)=G*IFFT[X<sup>H</sup>(k)*E(k)];X<sup>H</sup>(k)是X(k)的复共轭;<maths num="0008"><![CDATA[<math><mrow><mi>G</mi><mo>=</mo><mrow><mfenced open='[' close=']' separators=','><mtable><mtr><mtd><msub><mi>I</mi><msub><mi>L</mi><mi>k</mi></msub></msub></mtd><mtd><msub><mi>O</mi><msub><mi>L</mi><mi>k</mi></msub></msub></mtd></mtr><mtr><mtd><msub><mi>O</mi><msub><mi>L</mi><mi>k</mi></msub></msub></mtd><mtd><msub><mi>O</mi><msub><mi>L</mi><mi>k</mi></msub></msub></mtd></mtr></mtable></mfenced></mrow><mo>,</mo></mrow></math>]]></maths><img file="A200810223644C00035.GIF" wi="64" he="53" />和<img file="A200810223644C00036.GIF" wi="47" he="54" />分别表示L<sub>k</sub>阶的零矩阵和单位矩阵;所述二次更新频域权系数通过下述公式得到:如果L<sub>k+1</sub>=L<sub>k</sub>,则W(k+1)=W’(k+1);如果L<sub>k+1</sub>=2L<sub>k</sub>,根据一次更新频域权系数W’(k+1),计算获得时域系数w’(k+1),表示为:w’(k+1)=IFFT(W’(k+1));利用在w’(k+1)末端内插L<sub>k</sub>个零,再进行快速傅立叶变换,获得批处理数据块块长为2L<sub>k</sub>的二次更新频域权系数W(k+1),表示为:<img file="A200810223644C00037.GIF" wi="616" he="112" />则,W(k+1)=FFT(w(k+1));如果<maths num="0009"><![CDATA[<math><mrow><msub><mi>L</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>L</mi><mi>k</mi></msub><mo>,</mo></mrow></math>]]></maths>根据一次更新频域权系数W’(k+1),抽取W<sub>2i</sub>’(k+1),i=0,…,L<sub>k</sub>-1,获得批处理数据块块长为<img file="A200810223644C00042.GIF" wi="81" he="103" />的二次更新频域权系数W(k+1)表示为:<maths num="0010"><![CDATA[<math><mrow><mi>W</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>[</mo><msubsup><mi>W</mi><mn>0</mn><mo>,</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msubsup><mi>W</mi><mn>2</mn><mo>,</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>W</mi><mrow><mn>2</mn><msub><mi>L</mi><mi>k</mi></msub><mo>-</mo><mn>2</mn></mrow><mo>,</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow><mi>T</mi></msup><mo>.</mo></mrow></math>]]></maths>
地址 100190北京市海淀区北四环西路21号中国科学院声学研究所