%計算各站點SPI1，SPI3SPI6SPI9SPI12及季SPI、年SPI值
%改變line7—讀取excel數(shù)據(jù)（只包括站點和日期包圍的降雨數(shù)據(jù)）
%若計算其他尺度如SPI24，改變line9?（數(shù)據(jù)尺度）?；?改變line66-68
%改變line74—78（輸入到excel中）
%code?by?cc?2018-9-2
close?all;clc;
rainData?=?xlsread（‘E:\weihe‘1‘B5:P640‘）;?%讀取降雨數(shù)據(jù)
[MN]?=?size（rainData）;??%m總月數(shù)n站點數(shù)
years=M/12;??scale=[136912];
[spi1spi3spi6spi9spi12]=deal（zeros（MN））;?%136912各站點的spi值????
[spiSprspiSumspiAutspiYear]=deal（zeros（yearsN））;?spiWin=deal（zeros（years-1N））;?%四季各站點的spi值
time=0;
%%%%%%%%%%%%%%fList=zeros（M1）;?mList=zeros（M1）;?nList=zeros（M1）;%%%%%%%%%
for?station=1:N
????rainDataSta=rainData（:station）;???????%某一站降雨數(shù)據(jù)
????spiArrayT=zeros（Mlength（scale））;??????%某站5個尺度spi矩陣
????sheet=1;
????for?timeScale=scale????????????????????%分不同時間尺度
????????rainAdd=zeros（M-timeScale+11）;????%時間尺度下的累計雨量序列
????????num=length（rainAdd）;
????????spiColumn=zeros（M1）;?
????????sheet=sheet+1;
????????for?i=timeScale:M
????????????for?j=0:timeScale-1
????????????????rainAdd（i-timeScale+1）=rainAdd（i-timeScale+1）+rainDataSta（i-j）;
????????????end
????????end?
????????%#################################################
????????%計算12個序列下SPI值所需的參數(shù)和降雨為0的概率
????????qList=zeros（121）;para1=zeros（121）;para2=zeros（121）;
????????for?month=1:12?????%列年里12個同期序列
????????????rainSameT=rainAdd（month:12:num）;???%某時間尺度下同期降雨序列
????????????number=length（rainSameT）;
????????????%計算降雨量為0的概率
????????????m=0?;?n=0;
????????????for?i=1:number
????????????????if?rainSameT（i）==0
????????????????????m=m+1;
????????????????else
????????????????????n=n+1;
????????????????end
????????????end
????????????q=m/number;
????????????rainNot0=rainSameT（rainSameT~=0）;?