关于L_moment method

灰色节能君

不知道有没有人了解 L_moment method  estimation 就是用来估计雨滴谱伽马函数的参数的。大致就是基于初值，使得这个函数最小，得到参数。我用了matlab里的fminsearch函数，但是因为这里需要求不完全gamma函数，参数不能为负，运行总是报错。  不知道我的方法做的对不对，还是说 L_moment method有特定的迭代公式，求大佬指导一下。

伽蓝鸟

气象家园刚刚恢复。。额，听说水文学里经常用L moment来估计伽马分布的参数，我们大气里主要用最大似然来估计（矩法估计在形状参数alpha比较小的时候不太靠谱）。
Wilks的书上给出了伽马分布的最大似然近似公式，比较简单，不用迭代；在第四章4.4.3节。那本书我曾经发过。

灰色节能君

伽蓝鸟 发表于 2020-1-31 17:43
气象家园刚刚恢复。。额，听说水文学里经常用L moment来估计伽马分布的参数，我们大气里主要用最大似然来估 ...

谢谢！