function [bits_alloc,energy_alloc] = BitLoad(subchan_gains,total_bits,num_subc,gap,noise,M)
%比特分配算法
%----------
%
%輸入： 
%   subchan_gains ----子載波增益
%   total_bits ----比特總數
%   num_subc ----子載波數
%   gap ----系統間隔
%   noise ----噪聲功率
%   M ----最大星座大小
%輸出：
%   bits_alloc ----分配給每個子信道德比特數
%   power_alloc ----分配的總功率
%--------------
%計算每個信道的SNR
SNR = ComputeSNR(subchan_gains,noise,gap);
%這個函數初始化系統，即根據Chow算法分配比特數和功率，這個算法要比Campello算法有效
[bits_alloc,energy_alloc] = chow_algo(SNR,num_subc,M);
%根據目前信道的增益形成能量增加表，以供Campello算法使用 
energytable = EnergyTableInit(SNR,M);
%使用Campello算法
[bits_alloc,energy_alloc] = campello_algo(bits_alloc,energy_alloc,energytable,total_bits,num_subc,M);