matlab中注水算法的实现和原理

注水算法的公式

这公式怎么编辑啊，有点复杂啊

算法的相关公式用图片展示

前段时间学习使用注水算法，CSDN上分享的代码有很多，大同小异都差不多。但是原理不懂，代码很难看得明白，所以自己看明白之后分享一下原理，互帮互助。
工科男，文采有限，献丑。

注水算法的基本原理就是根据香农公式和限制条件，通过拉格朗日乘数法组成的一个方程，先令其偏导为零，求出Pi的表达式，但是Pi的表达式中包含一个未知数，再根据限制条件可以先求解出该未知数，再回代到之前的方程中，可以求解得每个信道根据信道质量分配得到的Pi。

代码，是实现上述算法步骤。我是参考了一个名为三寸灰的博主的代码，基于知识产权的问题，不敢贸然转载，他的博客地址我放在最下面的参考文献位置。
注水算法算是存在时间比较久，比较成熟的个算法。实现方法也有很多，还有一种别的方法，过几天整理了会发。

关于这个，有问题的话，欢迎交流。

我自己写了一个注水算法的实现代码，在我另一篇文章里（点）。