directory_name = uigetdir('D:\','打開(kāi)要處理的數(shù)據(jù)文件夾');
files = dir(directory_name);
files_length = length(files);
outpath = uigetdir('G:\','創(chuàng)建保存數(shù)據(jù)的文件夾');
for k = 3:files_length
    %%%%%%%%%%%%%%%%%%%%%%%%%  次雨量 單次大于等于12  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    numdata=xlsread([directory_name,'\',files(k).name],'次雨量 單次大于等于12');
    len=length(numdata)+1;
    y= xlsread([directory_name,'\',files(k).name],['l2:l',num2str(len)]);
    x= xlsread([directory_name,'\',files(k).name],['f2:f',num2str(len)]);
    myfunc=inline('beta(1).*(x.^beta(2))','beta','x');
    beta=nlinfit(x,y,myfunc,[0 0]);
    filename=str2num(files(k).name(1:5));
    xlswrite([outpath,'\','output'],filename,['A',num2str(k),':A',num2str(k)]);
    xlswrite([outpath,'\','output'],beta,['B',num2str(k),':C',num2str(k)]);
    
   %%%%%%%%%%%%%%%%%%%%%%%%%  日雨量 單次大于等于12  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    numdata=xlsread([directory_name,'\',files(k).name],'日雨量 單次大于等于12');
    len=length(numdata)+1;
    y= xlsread([directory_name,'\',files(k).name],['l2:l',num2str(len)]);
    x= xlsread([directory_name,'\',files(k).name],['f2:f',num2str(len)]);
    myfunc=inline('beta(1).*(x.^beta(2))','beta','x');
    beta=nlinfit(x,y,myfunc,[0 0]);
    filename=str2num(files(k).name(1:5));
    xlswrite([outpath,'\','output'],beta,['E',num2str(k),':F',num2str(k)]);

    
       %%%%%%%%%%%%%%%%%%%%%%%%%  日雨量單日大于等于12  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    numdata=xlsread([directory_name,'\',files(k).name],'日雨量單日大于等于12');
    len=length(numdata)+1;
    y= xlsread([directory_name,'\',files(k).name],['l2:l',num2str(len)]);
    x= xlsread([directory_name,'\',files(k).name],['f2:f',num2str(len)]);
    myfunc=inline('beta(1).*(x.^beta(2))','beta','x');
    beta=nlinfit(x,y,myfunc,[0 0]);
    filename=str2num(files(k).name(1:5));
    xlswrite([outpath,'\','output'],beta,['H',num2str(k),':I',num2str(k)]);
end;
%以下為測(cè)試
% xx=min(x):max(x);
% yy=beta(1).*xx.^beta(2);
% plot(x,y,'o',xx,yy,'r');