只贴基本的适合小白的Matlab实现代码，深入的研究除了需要改进算法，我建议好好研究一下混沌与分形，不说让你抓住趋势，至少不会大亏，这个资金盈亏回调我以前研究过。

function [line_H,RS]=dfaX_0322(x)
%% 初始化数据处理
% x=importdata('IF.txt');[line_H,RS]=dfaX_0322(x);line_H
%x=1000*rand(4000,1);[line_H,RS]=dfaX_0322(x);line_H
x=x(:);x=x+2;%防止x中有0
x1=log([x(2:end);10])-log(x);x=x1(1:end-1);%消除序列的短期自相关性，以满足R/S分析对观测对象独立的要求
% x=cumsum(x-mean(x));
%x=log(x(:));
N=length(x);
xulie=[8:floor(N/2)];%表示将x分成floor(N/xulie(i))段

%% 处理input

len_xulie=length(xulie);
RS=zeros(1,len_xulie);%初始化序列的重标极差

%% 求RS过程
for flag_i=1:len_xulie
m=xulie(flag_i);
junzhi_ij=zeros(1,floor(N/xulie(flag_i)));%初始化均值
biaozhuencha_ij=zeros(1,floor(N/xulie(flag_i)));%初始化标准差
jicha_ij=zeros(1,floor(N/xulie(flag_i)));%初始化极差
for flag_j=1:floor(N/xulie(flag_i))
yuansu_ij=x([m*(flag_j-1)+1:m*flag_j]);
junzhi_ij(1,flag_j)=mean(yuansu_ij);
biaozhuencha_ij(1,flag_j)=sqrt(mean((yuansu_ij-junzhi_ij(1,flag_j)).^2));
leijilicha_ij=cumsum(yuansu_ij-junzhi_ij(1,flag_j));
jicha_ij(1,flag_j)=max(leijilicha_ij)-min(leijilicha_ij);
end
RS(1,flag_i)=mean(jicha_ij./biaozhuencha_ij);
end

%% 拟合分段xulie和RS序列，求出Hurst指数
temp1=find(isnan(RS));RS(temp1)=[];xulie(temp1)=[];
line_H=polyfit(log10(xulie),log10(RS),1);

%% 作图
% subplot(1,2,1);plot(x);
% subplot(1,2,2);
% plot(log10(N./xulie(1:len_xulie)),log10(RS(1:len_xulie)),'ro');hold on;
% plot(log10(N./xulie(1:len_xulie)),polyval(line_H,log10(N./xulie(1:len_xulie))),'-b*');
% title(['拟合直线斜率',num2str(line_H(1))]);
% Log_N_RS=[log10(N./xulie(1:len_xulie);log10(RS(1:len_xulie))];
%

%% Hurst指数和相关性的关系
% (1)当0.5<h<1时，说明时间序列具有长程相关性，呈现出趋势不断增强的状态，即在某一时间段是递增(递减)趋势，下一个时间段也会是递增(递减)趋势，且h越接近于1，相关性越强。
% (2)当h=0.5时，说明时间序列不相关，是一个独立的随机过程，即当前状态不会影响将来状态。
% (3)当0<h<0.5时，说明径流时间序列只存在负的相关性，呈现反持久性的状态，即时间序列在某一个时间段是递增(递减)的趋势，则在下一个时间段是递减(递增)的趋势。

% 赫斯特指数有三种形式：
% 　　1．如果H=0.5，表明时间序列可以用随机游走来描述；
% 　　2．如果0.5<H≤1，表明黑噪声(持续性)即暗示长期记忆的时间序列；
% 　　3．如果0≤H<0.5，表明粉红噪声(反持续性)即均值回复过程。
% 　　也就是说，只要H ≠0.5，就可以用有偏的布朗运动(分形布朗运动)来描述该时间序列数据。

%% polyfit(拟合)和polyval（插值）应用例子
% x = (0: 0.1: 2.5)';
% y = erf(x);
% p = polyfit(x,y,6);
% y1=polyval(p,x);
% subplot(1,3,1);plot(x,y,'-b*');
% subplot(1,3,2);plot(x,y1,'-ro');
% subplot(1,3,3);plot(x,y,'-b*');hold on;plot(x,y1,'-ro');

%% n的输入
% 2.^[2:15]'
% 4
% 8
% 16
% 32
% 64
% 128
% 256
% 512
% 1024
% 2048
% 4096
% 8192
% 16384
% 32768
end