最优频带分配方案
最優頻帶分配問題
設兩個發送器分別具有P1P_1P1?,P2P_2P2?的功率,分別使用兩個不相交的頻帶帶寬 W1W_1W1?,W2W_2W2?,,其中 W1+W2=WW_1+W_2=WW1?+W2?=W(總帶寬),試計算最優頻帶分配方案,并畫出該信道的容量區域,進一步討論頻分多址(frequency-division multiplexing)系統的最優頻帶分配方案。
最優頻帶分配問題解決方案
利用單用戶的帶寬有限信道的容量公式,下面的碼率是可達的。
R1≤W1log?(1+P1NW1)                                                (1)R2≤W2log?(1+P2NW2)                                              (2)\begin{array}{l} {R_1} \le {W_1}\log (1 + \frac{{{P_1}}}{{N{W_1}}})\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(1)\\ {R_2} \le {W_2}\log (1 + \frac{{{P_2}}}{{N{W_2}}})\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;(2)\\ \end{array}R1?≤W1?log(1+NW1?P1??)(1)R2?≤W2?log(1+NW2?P2??)(2)?
滿足條件
W1+W2=W=constP1+P2=P=const\begin{array}{l} {W_1} + {W_2} = W = const\\ {P_1} + {P_2} = P = const\\ \end{array}W1?+W2?=W=constP1?+P2?=P=const?
則
R1+R2≤W1log?(1+P1NW1)+W2log?(1+P2NW2)                        (?){R_1}{\rm{ + }}{R_2} \le {W_1}\log (1 + \frac{{{P_1}}}{{N{W_1}}}) + {W_2}\log (1 + \frac{{{P_2}}}{{N{W_2}}})\;\;\;\;\;\;\;\;\;\;\;\;(*)R1?+R2?≤W1?log(1+NW1?P1??)+W2?log(1+NW2?P2??)(?)
下面證明上面(*)式右邊的碼率可達,并給出等號成立的條件。
首先從物理直觀上來看,由于R1R_1R1?和R2R_2R2?分別滿足(1)式和(2)式,這兩個不等式的右邊分別是它們的理論可達碼率,因此它們的加和一定小于各自的理論可達碼率的加和。
但是,這只是直觀上比較粗糙的理解,那么更嚴謹的問題在于: R1R_1R1?和R2R_2R2?的理論可達碼率是不是就是各自的理論可達碼率的加和呢?換句話說,它們的最大值能不能同時取到?如果可以,那么什么時候同時取到,等號成立的條件是什么?
畫出信道容量區域示意圖,大致的趨勢是實際的信道容量(曲線部分)小于等于理論最大值(折線部分),與橫坐標的交點表示信道的全部功率用來傳輸用戶1的信息( R2=0R_{2}=0R2?=0),與縱坐標的交點表示信道的全部功率用來傳輸用戶2的信息( R1=0R_1=0R1?=0),在這兩個點處顯然曲線部分與折線部分的取值相等。
命題轉化為證明:當 R1R_1R1?>0且 R2R_2R2?>0時,實際的信道容量(曲線部分)與理論最大值(折線部分)有且只有一個交點,且該交點對應的頻帶分配方案為
W1=P1P1+P2W,W2=P2P1+P2W                    (4){W_1} = \frac{{{P_1}}}{{{P_1} + {P_2}}}W,{W_2} = \frac{{{P_2}}}{{{P_1} + {P_2}}}W\;\;\;\;\;\;\;\;\;\;(4)W1?=P1?+P2?P1??W,W2?=P1?+P2?P2??W(4)
首先證明存在性,即:當 R1R_1R1?和R2R_2R2? 滿足(4)式時,(3)式等號成立。
R1max?=W1log?(1+P1NW1)=P1P1+P2Wlog?(1+P1NP1P1+P2W)=P1P1+P2Wlog?(1+P1+P2NW)R2max?=W2log?(1+P2NW2)=P2P1+P2Wlog?(1+P2NP2P1+P2W)=P2P1+P2Wlog?(1+P1+P2NW)(R1+R2)max?=R1max?+R2max?=Wlog?(1+P1+P2NW)\begin{array}{l} {R_{1\max }} = {W_1}\log (1 + \frac{{{P_1}}}{{N{W_1}}}) = \frac{{{P_1}}}{{{P_1} + {P_2}}}W\log (1 + \frac{{{P_1}}}{{N\frac{{{P_1}}}{{{P_1} + {P_2}}}W}}) = \frac{{{P_1}}}{{{P_1} + {P_2}}}W\log (1 + \frac{{{P_1} + {P_2}}}{{NW}})\\ {R_{2\max }} = {W_2}\log (1 + \frac{{{P_2}}}{{N{W_2}}}) = \frac{{{P_2}}}{{{P_1} + {P_2}}}W\log (1 + \frac{{{P_2}}}{{N\frac{{{P_2}}}{{{P_1} + {P_2}}}W}}) = \frac{{{P_2}}}{{{P_1} + {P_2}}}W\log (1 + \frac{{{P_1} + {P_2}}}{{NW}})\\ {({R_1} + {R_2})_{\max }} = {R_{1\max }} + {R_{2\max }} = W\log (1 + \frac{{{P_1} + {P_2}}}{{NW}})\\ \end{array}R1max?=W1?log(1+NW1?P1??)=P1?+P2?P1??Wlog(1+NP1?+P2?P1??WP1??)=P1?+P2?P1??Wlog(1+NWP1?+P2??)R2max?=W2?log(1+NW2?P2??)=P1?+P2?P2??Wlog(1+NP1?+P2?P2??WP2??)=P1?+P2?P2??Wlog(1+NWP1?+P2??)(R1?+R2?)max?=R1max?+R2max?=Wlog(1+NWP1?+P2??)?
再證明唯一性,即:實際的信道容量(曲線部分)與理論最大值(折線部分)有交點,而且只能有一個。
觀察碼率表達式
R1≤W1log?(1+P1NW1)R2≤W2log?(1+P2NW2)\begin{array}{l} {R_1} \le {W_1}\log (1 + \frac{{{P_1}}}{{N{W_1}}})\\ {R_2} \le {W_2}\log (1 + \frac{{{P_2}}}{{N{W_2}}})\\ \end{array}R1?≤W1?log(1+NW1?P1??)R2?≤W2?log(1+NW2?P2??)?
R1+R2R_1+R_2R1?+R2?的理論最大值可以寫成
R1+R2≤W1log?(1+P1NW1)+(W?W1)log?(1+P?P1N(W?W1)){R_1}{\rm{ + }}{R_2} \le {W_1}\log (1 + \frac{{{P_1}}}{{N{W_1}}}) + (W - {W_1})\log (1 + \frac{{P - {P_1}}}{{N(W - {W_1})}})R1?+R2?≤W1?log(1+NW1?P1??)+(W?W1?)log(1+N(W?W1?)P?P1??)
實際上,不等式右邊的理論可達碼率為凹函數,證明如下:
設PiN=ki>0,構造f(xi)=xilog?(1+kixi)f′(xi)=log?(1+kixi)+xi?(?kixi2)11+kixi=log?(1+kixi)?kixi+kif′′(xi)=11+kixi?(?kixi2)+ki(xi+ki)2=ki(1xi2+kixi+kixi+ki2?1xi2+kixi)<0∴f(x)為凹函數。\begin{array}{l} 設\frac{{{P_i}}}{N} = {k_i} > 0,構造f({x_i}) = {x_i}\log (1 + \frac{{{k_i}}}{{{x_i}}})\\ f'({x_i}) = \log (1 + \frac{{{k_i}}}{{{x_i}}}) + {x_i} \cdot ( - \frac{{{k_i}}}{{{x_i}^2}})\frac{1}{{1 + \frac{{{k_i}}}{{{x_i}}}}} = \log (1 + \frac{{{k_i}}}{{{x_i}}}) - \frac{{{k_i}}}{{{x_i} + {k_i}}}\\ f''({x_i}) = \frac{1}{{1 + \frac{{{k_i}}}{{{x_i}}}}} \cdot ( - \frac{{{k_i}}}{{{x_i}^2}}) + \frac{{{k_i}}}{{{{({x_i} + {k_i})}^2}}} = {k_i}(\frac{1}{{{x_i}^2 + {k_i}{x_i} + {k_i}{x_i} + {k_i}^2}} - \frac{1}{{{x_i}^2 + {k_i}{x_i}}}) < 0\\ \therefore f\left( x \right)為凹函數。 \end{array}設NPi??=ki?>0,構造f(xi?)=xi?log(1+xi?ki??)f′(xi?)=log(1+xi?ki??)+xi??(?xi?2ki??)1+xi?ki??1?=log(1+xi?ki??)?xi?+ki?ki??f′′(xi?)=1+xi?ki??1??(?xi?2ki??)+(xi?+ki?)2ki??=ki?(xi?2+ki?xi?+ki?xi?+ki?21??xi?2+ki?xi?1?)<0∴f(x)為凹函數。?
令g(x)=(x0?x)log(1+kix0?x),其中0<x<x0g′(x)=?log(1+kix0?x)+(x0?x)11+kix0?x?ki(x0?x)2              =?log(1+kix0?x)+kix0?x+kig′′(x)=?ki(x0?x)2?11+kix0?x+ki(x0?x+ki)2                  =ki[?1(x0?x)2+ki(x0?x)+1(x0?x+ki)2]                  =ki?(x0?x+ki)2+(x0?x)2+ki(x0?x)[(x0?x)2+ki(x0?x)](x0?x+ki)2                  =?ki2ki+(x0?x)[(x0?x)2+ki(x0?x)](x0?x+ki)2<0因此g(x)=(x0?x)log(1+kix0?x)仍為凹函數。\begin{array}{l} 令g(x) = ({x_0} - x)log(1 + \frac{{{k_i}}}{{{x_0} - x}}),其中0 < x < {x_0}\\ g'(x) = - log(1 + \frac{{{k_i}}}{{{x_0} - x}}) + ({x_0} - x)\frac{1}{{1 + \frac{{{k_i}}}{{{x_0} - x}}}} \cdot \frac{{{k_i}}}{{{{({x_0} - x)}^2}}}\\ \;\;\;\;\;\;\; = - log(1 + \frac{{{k_i}}}{{{x_0} - x}}) + \frac{{{k_i}}}{{{x_0} - x + {k_i}}}\\ g''(x) = - \frac{{{k_i}}}{{{{({x_0} - x)}^2}}} \cdot \frac{1}{{1 + \frac{{{k_i}}}{{{x_0} - x}}}} + \frac{{{k_i}}}{{{{({x_0} - x + {k_i})}^2}}}\\ \;\;\;\;\;\;\;\;\; = {k_i}[ - \frac{1}{{{{({x_0} - x)}^2} + {k_i}({x_0} - x)}} + \frac{1}{{{{({x_0} - x + {k_i})}^2}}}]\\ \;\;\;\;\;\;\;\;\; = {k_i}\frac{{ - {{({x_0} - x + {k_i})}^2} + {{({x_0} - x)}^2} + {k_i}({x_0} - x)}}{{[{{({x_0} - x)}^2} + {k_i}({x_0} - x)]{{({x_0} - x + {k_i})}^2}}}\\ \;\;\;\;\;\;\;\;\; = - {k_i}^2\frac{{{k_i} + ({x_0} - x)}}{{[{{({x_0} - x)}^2} + {k_i}({x_0} - x)]{{({x_0} - x + {k_i})}^2}}} < 0\\ 因此g(x) = ({x_0} - x)log(1 + \frac{{{k_i}}}{{{x_0} - x}})仍為凹函數。 \end{array}令g(x)=(x0??x)log(1+x0??xki??),其中0<x<x0?g′(x)=?log(1+x0??xki??)+(x0??x)1+x0??xki??1??(x0??x)2ki??=?log(1+x0??xki??)+x0??x+ki?ki??g′′(x)=?(x0??x)2ki???1+x0??xki??1?+(x0??x+ki?)2ki??=ki?[?(x0??x)2+ki?(x0??x)1?+(x0??x+ki?)21?]=ki?[(x0??x)2+ki?(x0??x)](x0??x+ki?)2?(x0??x+ki?)2+(x0??x)2+ki?(x0??x)?=?ki?2[(x0??x)2+ki?(x0??x)](x0??x+ki?)2ki?+(x0??x)?<0因此g(x)=(x0??x)log(1+x0??xki??)仍為凹函數。?
凹函數滿足
f[λx1+(1?λ)x2]≥λf(x1)+(1?λ)f(x2),?x1,x2∈I,0≤λ≤1f[\lambda {x_1} + (1 - \lambda ){x_2}] \ge \lambda f({x_1}) + (1 - \lambda )f({x_2}),\forall {x_1},{x_2} \in I,0 \le \lambda \le 1f[λx1?+(1?λ)x2?]≥λf(x1?)+(1?λ)f(x2?),?x1?,x2?∈I,0≤λ≤1
g[μx1+(1?μ)x2]≥μg(x1)+(1?μ)g(x2),?x1,x2∈I,0≤μ≤1g[\mu {x_1} + (1 - \mu ){x_2}] \ge \mu g({x_1}) + (1 - \mu )g({x_2}),\forall {x_1},{x_2} \in I,0 \le \mu \le 1g[μx1?+(1?μ)x2?]≥μg(x1?)+(1?μ)g(x2?),?x1?,x2?∈I,0≤μ≤1
由凹函數圖像特點可知實際的信道容量(曲線部分)與理論最大值(折線部分)有交點,而且只能有一個,唯一性得證。
由上述的推導可以得出結論:對于若干個電臺,只有當所有分配的帶寬與對應的功率成正比時,對應的頻帶分配方案是最優的。某子信道對應的功率越大,那么它分到的帶寬也越大,便于傳輸更多的信息。
對于FDMA系統,可以得出類似的結論。前面構造的 f(xi)f(x_i)f(xi?)可以推廣到n維,對應頻分多址的路數。
一般地,可以推廣為:對于m個具有功率為 P1,P2,P3,...,PmP_1,P_2,P_3,...,P_mP1?,P2?,P3?,...,Pm?的信源以及功率N的環境噪聲的高斯多接入系統,任何集合S對高斯公式平移為下列形式:
∑i∈SRi=S≤C(∑i∈SPiN)\sum\limits_{i \in S} {{R_i} = S} \le C(\frac{{\sum\limits_{i \in S} {{P_i}} }}{N})i∈S∑?Ri?=S≤C(Ni∈S∑?Pi??)
總結