函数.zip

函数.zip

%该函数用于不同省份站点不同月份的取值
function Ka = calculateKa（stnIDMon）
%ka为四川、山东、山西、内蒙、宁夏、青海、陕西、河南、甘肃9省（行）12个月（列）的Ka值
ka = [1.0 	1.0 	1.1 	1.2 	1.0 	1.2 	1.2 	1.2 	1.0 	1.0 	1.0 	1.0 ;
      0.4 	0.8 	1.0 	1.2 	1.2 	1.2 	1.2 	1.0 	1.0 	0.8 	0.6 	0.4 ;
      0.4 	0.8 	1.0 	1.2 	1.2 	1.2 	1.2 	1.0 	1.0 	0.8 	0.6 	0.4 ;
      0.0 	0.0 	0.0 	0.6 	1.0 	1.2 	1.2 	1.0 	0.9 	0.4 	0.0 	0.0 ;
      0.4 	0.8 	1.0 	1.0 	1.0 	1.2 	1.2 	1.0 	0.9 	0.8 	0.6 	0.4 ;
      0.0 	0.0 	0.0 	0.6 	1.0 	1.2 	1.2 	1.0 	0.9 	0.4 	0.0 	0.0 ;
      0.4 	0.8 	1.0 	1.2 	1.2 	1.2 	1.2 	1.0 	1.0 	0.8 	0.6 	0.4 ;
      0.4 	0.8 	1.0 	1.2 	1.2 	1.2 	1.2 	1.0 	1.0 	0.8 	0.6 	0.4 ;
      0.4 	0.8 	1.0 	1.2 	1.2 	1.2 	1.2 	1.0 	1.0 	0.8 	0.6 	0.4 ];
%stnid为四川、山东、山西、内蒙、宁夏、青海、陕西、河南、甘肃各站点的代号
stnid = {56172;						
         [54725	54823	54843	54909	54916];		
         [53487	53673	53772	53863	53959];		
         [53463	53336];					
         [53614	53723];					
         [52836	52866	52955	56029	56046];		
         [53646	57127];					
         [53898	57083];					
         [52787	52884	53915	56096	57014	52996	52983]};
%求站点所在的行，对应ka所在的行
[mn] = find（cellfun（@（x） ismember（stnIDx）stnid））;
%取对应站点对应月份的Ka值
Ka = ka（mMon）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-05-04 08:13  函数\
     文件        1318  2020-04-30 17:09  函数\calculateKa.m
     文件         188  2020-04-28 23:43  函数\calculateMI.m
     文件        2720  2020-04-29 09:21  函数\CalculatePET.m
     文件         344  2020-04-30 12:15  函数\calculatePn.m
     文件         690  2020-04-30 12:54  函数\CalculateSPI.m
     文件         572  2020-04-30 16:32  函数\calculateSPIW60.m
     文件        4434  2020-04-30 17:23  函数\main.m
     文件         610  2020-04-30 17:12  函数\timemake.m