注水法功率分配问题.m

 
% Mohammad Ismail Hossain  
% Jacobs University Bremen 
% Waterfilling Alogorithm 
 
clc 
clear all; 
Trans_Power=10;  
Noise_Power=[2 3 4 1 3 4 3 2]; 
Number_Channel= length（Noise_Power） ;  
[S_Number dt]=sort（Noise_Power）; 
sum（Noise_Power） 
for p=length（S_Number）:-1:1 
    T_P=（Trans_Power+sum（S_Number（1:p）））/p; 
    Input_Power=T_P-S_Number; 
    Pt=Input_Power（1:p）; 
    if（Pt（:）>=0） 
        break 
    end 
end