WSDI的计算MATLAB程序

herbking

最近在进行各种干旱指数的计算，今天先在这里分享一下water storage deficit index (WSDI)的计算程序，希望大家多多交流指正。

1. clear;
   
2. TWSK=load('test.txt');%2003/01-2015/12,a time series of water storage anomaly
   
3. for i=1:13
   
4.    yK=TWSK(1+(i-1)*12:12+(i-1)*12);%y every year data
   
5.    rawK(i,:)=yK;%
   
6. end
   
7. for j=1:12
   
8.     damK=rawK(:,j);%13 year every day data
   
9.     meanK(j)=mean(damK);
   
10. end
    
11. for i=1:13
    
12.     for j=1:12
    
13.         WSDK(i,j)=rawK(i,j)-meanK(j);
    
14.     end
    
15. end
    
16. for i=1:13
    
17.     WSDK1(1+(i-1)*12:12+(i-1)*12)=WSDK(i,:);
    
18. end
    
19. WSDK1=WSDK1';
    
20. WSDIK=zscore(WSDK1);
    
21. save('WSDI_text.txt','WSDIK','-ascii');

复制代码