%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear all %% raw_data=xlsread('1.xlsx','sheet1','A131:F296');%data depth=raw_data(:,1);%depth data=raw_data(:,4); data_guiyi=mapminmax(data',0,1)'; data_new=sum(data_guiyi,2); %% for i=1:length(data_new) for j=1:i A(j)=sum(data_new(1:j))-(j/i)*sum(data_new(1:i)); R(i)=max(A)-min(A); end S(i)=((1/(i))*(sum(data_new(1:i)))^2-((1/(i))*sum(data_new(1:i)))^2)^(1/2); C=R./S; end %% C(1:3)=[]; C=C'; %% y=1:length(C); y=y'; % fitType=fittype('a*(C)^b','independent','C','coefficients',{'a','b'}); P=fit(y,C,fitType); b=P.b; %% H=2-b %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% [imf,residual,info]=emd(data_new); %% [hs,f,t,imfinsf,imfinse] = hht(imf); %% % I=(length(depth)*ones(length(info),1)./info.NumExtrema)'; A=cell(1,length(info.NumIMF)); for i=1:length(info.NumIMF) A{1,i}=findpeaks(imf(:,i)); end for i=1:length(info.NumIMF) Q(i)=length(A{i}); end I=(ones(length(info.NumIMF),1)*length(depth)./Q')' %% x=info.NumIMF; fitType=fittype('a*(x)^b','independent','x','coefficients',{'a','b'}) P=fit(x,I',fitType) b=P.b %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% F=H/b