matlab中标准化降水指数程序,标准化降水指数（SPI）计算程序

之前遇到的一个客户需求就是计算标准化降水指数(SPI)。

主要参考的论文为：洪兴骏，等 “标准化降水指数SPI分布函数的适用性研究”。原文可以很容易搜索到。

论文将计算方法其实写的比较清楚了，本文主要提供其C++实现。

头文件

//SPI_Calculator.h

#include

#include

#include

#include

#include

using namespace std;

#ifndef SPICALCULATOR_H

#define SPICALCULATOR_H

class SPI_Calculator{

public:

/**

path is the file that provides the data of rain fall

file format:

yearamount of precipitation

yearamount of precipitation

.

.

.

*/

SPI_Calculator(char* path);

/**

parameter is the year to get the SPI

return SPI

*/

double getSPI(int year);

double f_probabilityDensity(double x);

double f_gamma(double z, unsigned Ntermsx);

double integral(double a,double b);</