matlab计算spi代码的解释

帅帅f

其中这一步的意思是：A1=[];
for is=1:scale
A1=[A1,Data(is:length(Data)-scale+is)];  
end
XS=sum(A1,2)
Data是降雨量数据。如果计算尺度（scale=3），那么A1就变成了n行3列的数据，第一列是从第一个月降水开始，第二列是从第二个月降水开始，第三列是从第3个月降水开始，再使用求和工具sum，使矩阵A1变成n行1列的数据，也就是上面的XS。
for is=1:nseas
    tind=is:nseas:length(XS);%这里的nseas取12，当is从1开始取时，tind的矩阵包括了研究时期所有1月份的降雨量。当is取2时，以此类推。
    Xn=XS(tind);
    [zeroa]=find(Xn==0);%找出降雨量为0的位置
    Xn_nozero=Xn;Xn_nozero(zeroa)=[];
    q=length(zeroa)/length(Xn);%降雨量为0的概率
    parm=gamfit(Xn_nozero);%拟合pdf
    Gam_xs=q+(1-q)*gamcdf(Xn,parm(1),parm(2));%看过spi计算原理的都知道当考虑降雨的等于0的情况，伽马函数会作调整。
    Z(tind)=norminv(Gam_xs);%转换为标准正态分布
end


代码来自mathworks.cn 由 Taesam Lee编写

dongyi

谢谢您解释

帅帅f

dongyi 发表于 2019-11-18 08:10
谢谢您解释

一起学习

leijiajin

在matlab上打开你的程序后，怎么输入数据啊？能不能演示下

帅帅f

leijiajin 发表于 2020-3-11 20:44
在matlab上打开你的程序后，怎么输入数据啊？能不能演示下

这个很简单的  你把降水量搞成一列数据  在运用那个函数（matlab加载函数会把）就可以得出结果  我和个只是解释函数中的某些东西