主权项 |
基于误码率优化的OFDM系统资源分配算法,其特征在于:基于共轭梯度提出最小误码率算法,把共轭性与最速下降法相结合,利用已知点处的梯度构造一组共轭方向,并沿着此方向进行搜索,求出目标函数极小值;首先构造拉格朗日函数:<img file="FDA0000997280320000011.GIF" wi="1011" he="278" />其中,λ<sub>k</sub>、μ<sub>k</sub>表示Lagrange乘子,<img file="FDA0000997280320000012.GIF" wi="798" he="175" />原问题的对偶函数是:D(λ,μ)=min L(λ,μ),对偶问题为:max D(λ,μ)s.t.λ≥0,μ≥0对偶函数为:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>D</mi><mrow><mo>(</mo><mi>P</mi><mo>,</mo><mi>λ</mi><mo>,</mo><mi>μ</mi><mo>)</mo></mrow><mo>=</mo><mi>min</mi><mo>[</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>Q</mi><mrow><mo>(</mo><msub><mover><mi>c</mi><mo>~</mo></mover><mi>k</mi></msub><mo>(</mo><mi>p</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>λ</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>μ</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>C</mi><mo>-</mo><mi>C</mi><mo>(</mo><mi>P</mi><mo>)</mo><mo>)</mo></mrow><mo>]</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000997280320000013.GIF" wi="830" he="271" /></maths>用共轭梯度优化方法来求解非线性的最优化问题;通过惩罚函数的方法将目标函数转化为:I(p;ε<sub>n</sub>)=P<sub>E</sub>(p)+ε<sub>n</sub>B(p),B(p)=‑log(1‑p<sup>T</sup>p)约束条件:<maths num="0002"><math><![CDATA[<mrow><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>P</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>≤</mo><msub><mi>P</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>≥</mo><mn>0</mn><mo>,</mo><mo>∀</mo><mi>n</mi><mo>,</mo><mi>k</mi></mrow>]]></math><img file="FDA0000997280320000014.GIF" wi="757" he="127" /></maths><maths num="0003"><math><![CDATA[<mrow><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mfrac><msub><mi>ρ</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mi>K</mi></mfrac><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>KP</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><msup><msub><mi>h</mi><mrow><mi>n</mi><mo>,</mo><mi>k</mi></mrow></msub><mn>2</mn></msup></mrow><mrow><msub><mi>BN</mi><mn>0</mn></msub></mrow></mfrac><mo>)</mo></mrow><mo>≤</mo><mi>C</mi><mrow><mo>(</mo><mi>P</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000997280320000015.GIF" wi="870" he="175" /></maths>其中下标n为迭代的次数;B(p)为惩罚函数,取ε<sub>1</sub>>0,c>1,令ε<sub>n+1</sub>=ε<sub>n</sub>/c;当n→∞的时候,ε<sub>n</sub>→0,I(p;ε<sub>n</sub>)→P<sub>E</sub>(p)。 |