Chow比特功率分配算法及MATLAB实现

原理介绍

Chow算法是根据各个子信道的信道容量来分配比特的。它的优化准则是在维持目标误比特率的前提下，使系统的余量最大。该算法通过迭代过程，逐步分配比特，同时使系统的余量逐步增大，直到所有的比特都分配完毕。

比特分配

（1）计算各个子信道的信噪比，$SNR\left(i\right),\forall i$， 假设所有子信道上的信号能量都是归一化的，$\epsilon \left(i\right)=1,\forall i$
（2）令${\gamma }_{m\arg in}=0\text{dB}$ ，迭代次数为$Itercount=0$ ，已使用的子信道数为 $Usedcarriers=N$，其中 $N$为可用的子信道数的最大数目
（3）依次计算$b(i)$ ，$\hat{b}(i)$ ，$diff(i)$ 和 $Usedcarriers$
$$b(i)={\log }_2\left(1+\frac{SNR(i)}{\Gamma +{\gamma }_{m\arg in}(dB)}\right)$$
$$\hat{b}(i)=round(b(i))$$
$$diff(i)=b(i)-\hat{b}(i)$$
若 $\hat{b}(i)\text{=}0$， $Usedcarriers\text{=}N-1$，$\Gamma$ 为 $SNR$间距，定义为系统达到传输容量Ｃ时的理想$SNR$ 和实际传输容量 $R$ 的 $SNR$之比。 $\Gamma$是$BER$ 和信道编码方式的函数，假设采用QAM调制，则
$$\Gamma \text{=}\frac{-\ln \left(5BE{R}_{t\arg et}\right)}{1.5}$$
（4）计算 $R\text{=}sum\left(\left\{\hat{b}\left(i\right)\right\}\right)$，若$R=0$ ,则信道太差，无法使用
（5）计算 ${\gamma }_{margin}$
$${\gamma }_{margin}={\gamma }_{margin}+10{\log }_{10}\left(2^{\frac{R-R_T}{USedCarriers}}\right)$$
（6）迭代次数 $Itercount=Itercount+1$
（7）若 $R\ne R_T$且$Itercount<Maxcount$ ，令$Usedcarriers=N$ ，转步骤（3），否则转步骤（8）
（8）若$R>R_T$ ，找到最小的$diff\left(i\right)$ ，相应地 $\hat{b}(i)\text{=}\hat{b}(i)\text{-}1$，$diff\left(i\right)=diff\left(i\right)+1$ ，重复此步骤直至 $R\text{=}{R}_T$
（9）若$R<R_T$ ，找到最大的 $diff\left(i\right)$，相应地$\hat{b}(i)\text{=}\hat{b}(i)\text{+}1$ ，$diff\left(i\right)=diff\left(i\right)\text{-}1$ ，重复此步骤直至 $R\text{=}{R}_T$
（10）调整每个子信道上的发射功率，使得相应于 $b(i)$
$$P_e\left(i\right)=P_{e,t\arg et},\forall i$$
（11）调整发射总功率。对所有已使用的子信道乘以相同的比例因子，
使总的信号功率 $P=P_T$

仿真参数设置

参数名称	参数值
平均信噪比	20
误比特率	$10^{-4}$
子载波数目	32
平均发射功率	1
比特总数上限	128

实验结果如下

从上图中可以看出，通过逐步调整分配比特，从而使系统最终的比特数达到最大。同时，也可以看出，Chow算法分配的结果有一定的优势，它是将信道状况较差的子载波舍弃，也就是说系统最终传输的所用的子载波一定有信息的传输，从仿真的结果可能证实这一点。
代码如下：
主函数

clear ;
close all;
clc;
N_subc=32; 
BER=1e-4; 
gap=-log(5*BER)/1.5; 
P_av=1; 
Pt=P_av*N_subc; 
SNR_av=20; 
noise=P_av./10.^(SNR_av./10); 
Rt=128; 
subcar_gains=random('rayleigh',1,1,N_subc); 
SNR=(subcar_gains.^2)./(noise*gap);  
[bit_alloc,power_alloc,Iterate_count]=Chow(SNR,N_subc,gap,Rt); 
power_alloc=Pt.*(power_alloc./sum(power_alloc));
figure(1); 
subplot(2,1,1); 
plot(subcar_gains,'-r'); 
legend('信道增益'); 
hold on; 
stem(bit_alloc); 
title('Chow比特功率分配'); 
ylabel('比特分配结果'); 
xlabel('子载波序号'); 
subplot(2,1,2); 
stem(power_alloc); 
ylabel('功率分配结果'); 
xlabel('子载波序号'); 

Chow.m

function [bits_alloc, power_alloc,iter_count] = Chow(SNR,N_subc,gap,Rt) 
% ---------------------initialization------------------------- 
margin=0;         %初始化门限        
Max_count=10;     %最大迭代次数        
iter_count=0;     %初始化迭代                  
total_bits=0;     %总比特殊            
power_alloc=zeros(1,N_subc);  %初始化功率分配结果
bits_alloc=zeros(1,N_subc);   %初始化比特分配结果
temp_bits=zeros(1,N_subc);    %初始化中间变量
bit_round=zeros(1,N_subc);    %初始化取整比特
diff=zeros(1,N_subc);         %初始化差分比特
%-----------------------------bits allocation------------------------- 
while (total_bits~=Rt)&&(iter_count<Max_count) 
    iter_count=iter_count+1; 
    N_use=N_subc;                            %初始化啊可用子载波数目
    temp_bits=log2(1+SNR./(1+margin/gap));   %分配比特的初始结果
    bit_round=round(temp_bits);              %取整后的比特数
    diff=temp_bits-bit_round;                %比特数差值
    total_bits=sum(bit_round);               %取整后的比特总数
    if(total_bits==0) 
        disp('the channel is not be used');  %信道太差，无法使用
    end 
    N_notuse=length(find(bit_round==0));       %未使用的子载波数目
    N_use=N_subc-N_notuse;                      %使用的子载波数目 
    margin=margin+10*log10(2^((total_bits-Rt)/N_use));    %更新门限
end 
%------------------------------bits alteration-------------------------- 
while(total_bits>Rt)                %如果 bit_total>Rt，找到最小的bit_diff(i)
    use_ind=find(bit_round>0);      %相应地bit_round减1，bit_total=bit_total-1
    bit_diff=diff(use_ind);         %bit_diff(i)加1，直至bit_total=Rt
    id=find(bit_diff==min(bit_diff),1); 
    ind_alter=use_ind(id);  
    bit_round(ind_alter)=bit_round(ind_alter)-1; 
    diff(ind_alter)=diff(ind_alter)+1; 
    total_bits=sum(bit_round); 
end 
while(total_bits<Rt)               %如果 bit_total<Rt，找到最大的bit_diff(i)
    use_ind=find(bit_round~=0);     %相应地bit_round加1，bit_total=bit_total+1
    bit_diff=diff(use_ind);          %bit_diff(i)减1，直至bit_total=Rt
    id=find(bit_diff==max(bit_diff),1); 
    ind_alter=use_ind(id); 
    bit_round(ind_alter)=bit_round(ind_alter)+1; 
    diff(ind_alter)=diff(ind_alter)-1; 
    total_bits=sum(bit_round); 
end 
bits_alloc=bit_round; 
%--------------------------power allocation----------------------------- 
power_alloc=(2.^bits_alloc-1)./SNR; 
end