function?[allocatedpower?waterlv]?=?waterfill2（transpowernoisearr）

????%?WATERFILL2?returns?waterfill?algorithm?results?of?the?given?SNR?paras.
????%?注水算法，返回各個發(fā)射天線上分配的信號功率
????%?伴隨有完整算法注釋。為避免和MATLAB內(nèi)嵌函數(shù)重名，故加上2后綴
????%?Author:?Jumper_Le?（Southeast?University）?20171125?Version.1
????
????[arisenoisedt]=sort（noisearr）;
????%?升序排列noise方差序列
????for?p=length（arisenoise）:-1:1
????????%?開始降序號循環(huán)
????????waterlevel=（transpower+sum（arisenoise（1:p）））/p;
????????%?計算當(dāng)前前p個噪聲水平算出的水平面值
????????pouringpower=waterlevel-arisenoise;
????????%?計算以當(dāng)前水平面，需要注入的水量
????????pp=pouringpower（1:p）;
????????%?復(fù)制注入水量數(shù)組的前p個值到新變量pp
????????if（pp（:）>=0）
????????????break;
????????????%?如果當(dāng)前情況下，注入水量均為正值（有意義），那么退出循環(huán)，確認(rèn)注水結(jié)束
????????end
????????%?如果注入水量有負(fù)值，需要舍棄當(dāng)前序列噪聲最大項，重新進(jìn)行剩余項的注水分配
????end
????allocatedpower=zeros（1length（arisenoise））;
????allocatedpower（dt（1:p））=pp;
????%?返回p長度的信號功率分配數(shù)組，剩余的項均為0，代表超出水平面（不考慮）
????waterlv=waterlevel;
????%?返回waterlevel
????
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1190??2017-11-26?02:12??waterfill2.m

?????文件???????1054??2017-11-26?12:08??waterfill3.m

-----------?---------??----------?-----??----

?????????????????2244????????????????????2