一类时间序列的多重分形分析

%%
%作者：**
%邮箱：ydlrpch@126.com  QQ: 442700198
%完成时间：2015年6月
%参考论文：一类时间序列的多重分形分析，丁丽娟，华东师范大学，硕士学位论文
%功能描述：主要研究多重分形的应用（1）ｑ阶矩结构分割函数法（qth-order moment structure partiction function
%           methods）是进行时间序列多重分形分析较为直观的方法。它是用多重分形谱来
%           描述归一化后的分形时问序列在不同标度区间的特征，即从系统的局部出发来研
%           究其最终的特征．它不仅使计算复杂性大大降低，而且通过对多重分形谱的计
%           算，对分形测度在支集上的分布进行了刻画，从而扩大了分形的应用范围．
%
%           这种方法降低了计算的复杂性，且扩大了分形的应用范围，在物理、
%           材料、自然图象的处理等方面都有广泛应用．能检测出数据是否具有多重分形性．Ding—Shun Ito等运用一些统计工具对1987—2002年间台湾每日股价指数
%           （TSPI）进行了分析．发现其呈右偏态频率分布，通过自相关和功率谱分析分别粗
%           略检测到了长程相关性和标度变化性．分析证实了所研究时间序列中多重分形
%           性质的存在．王树珏、田华也用此方法对深圳证券市场进行了多重分形分析，
%           结果表明中国股市存在多重分形特征：同时通过多重分形谱的统计物理计算，得
%           出中国股市多重分形特征较弱的结论．在医学上，Plamen Ch则用此方法分
%           析了心跳动力系统的多重分形性．并对健康和非健康的对象进行了比较，结果表
%           明，非健康者系统的多重分形性弱些．
%% 清理环境变量，导入数据
clc ;
clear ;
%load 20150520脉搏谱分析实验第二组.mat;value = value ;
load mbl22015-06-04220443334;value = mbl2;

%% 归一化
vMax = sum（value）;
valuen = value/vMax;
%% 计算盒概率
L = length（valuen）;
p = 0 ;
sl = [] ;
for s = 6 :fix（L/4）;%把信号分成V个S长的小块
 sl = [sl  s]; %记录不同S的值
 Ns = L/s;
 Pv=[];
 for v = 1:N