实现方法来自论文:

Quantification of scaling exponents and crossover phenomena in nonstationary
heartbeat time series
C.‐K. Peng , Shlomo Havlin , H. Eugene Stanley , and Ary L. Goldberger
输入:一个数组B,长度为len
输出:斜率b(2)
第一步:计算y(k),其中Bave位B的均值。

第二步:计算yn(k),n是能被len整除的数,比如我这是3500个数,那么n就可以取1,2,5,10,35等。

根据n的大小划分,比如我取n=500,那么y就被划分成了7块,对应了1-500;501-1000;1001-1500;1501-2000;2001-2500;2501-3000;3001-3500。

然后计算每一块的regress(即利用回归分析得到预测值)。

第三步:根据如下公式计算F(n)

第四步:得到n和F(n),以log10为底,然后画图,再次用regress得到斜率b。

function out=dfa(a)avg=mean(a);sum=0.0;len=length(a);for i=1:3500sum=sum+a(i)-avg;y(i)=sum;end%x=1:3500;%X=[ ones(length(y),1), x' ];%b=regress(y',X)index=1;for n=10:3500if(rem(3500,n)==0)    count=0;k=3500/n;for i=1:kbegin=(i-1)*n+1;finish=i*n;q=y(begin:finish);p=begin:finish;X=[ ones(length(q),1), p' ];b=regress(q',X); %回归 y=b1+b2*x;for j=begin:finishcc=b(1)+b(2)*j;op=y(j)-cc;count=count+op*op;endendf=(count/3500)^0.5;K(index)=log10(n);F(index)=log10(f);index=index+1;endendX=[ ones(length(F),1), K' ];[b,bint,r,rint,stats]=regress(F',X);out=b(2);z=b(1)+b(2)*K;figure,plot(K,F,'rp',K,z,'b');title('散点图及回归线');
end

matlab 计算DFN(去趋势波动分析) regress回归分析相关推荐

  1. 【数据分析】多重分形去趋势波动分析附matlab代码

    ✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信.

  2. 预处理丨去趋势(Matlab和C++)

    预处理丨去趋势(Matlab和C++) 预处理丨去趋势(Matlab和C++) 介绍 原理 MATLAB 实现 C++ 实现 本科没学过信号处理,对采集的数据需要做预处理,如去趋势,本文介绍一下去趋势 ...

  3. 关于CDO中去趋势的一点理解

    碎碎念 1.去趋势基本概念理解 维基百科给出的解释: 在随机过程, 混沌理论和时间序列分析中, 去趋势波动分析(英文:Detrended Fluctuation Analysis, DFA)是一种判断 ...

  4. matlab去趋势,[转载]使用Matlab对数据进行去趋势(detrend)

    使用Matlab对数据进行去趋势(detrend) 介绍 去趋势(detrend)处理可以消除传感器在获取数据时产生的偏移对后期计算产生的影响.从数据中删除趋势可以将分析集中在数据趋势本身的波动上.但 ...

  5. matlab 趋势分析,科学网—SEN趋势度分析及其MATLAB实现 - 杨建华的博文

    SEN趋势度分析及其MATLAB实现 一.博文概述 1.目的 学习理解SEN趋势度分析方法及其应用,并编写简易的SEN趋势度分析函数 2.时间 2016年12月24日 3.关键词 SEN趋势度分析   ...

  6. 股票市场周期波动分析的MATLAB实现(1)

    [摘要]本文主要基于MATLAB结合霍华德·马克斯的<周期>一文实现对上证指数的分析.<周期>一书着重强调了股市周期的存在,而在周期中的买进和卖出位置就决定了投资收益.因而利用 ...

  7. Matlab 去趋势处理

    1.介绍 去趋势处理可以去除传感器获取数据时产生的偏移,可以将数据分析集中在波动上 2.原理 对数据减去一条最优(拟合)的直线,平面或者曲面,使得处理后的数据均值为零 3.函数 3.1连续线性趋势 y ...

  8. matlab趋势面分析,第三章趋势面分析.ppt

    <第三章趋势面分析.ppt>由会员分享,可在线阅读,更多相关<第三章趋势面分析.ppt(61页珍藏版)>请在人人文库网上搜索. 1.第三章 地理学中的经典统计分析方法,第6节 ...

  9. 风险管理matlab,风险管理KMV模型Matlab计算—-实例分析 - Ariszheng

    风险管理KMV模型Matlab计算--实例分析 编写自己编写的 KMVOptSearch 添加两个KMV模型文档 2009-6-5 2009-8-14 程序内容将作为 金融数量分析的一节进行仔细讲解 ...

最新文章

  1. 一口气搞懂《虚函数和纯虚函数》
  2. SWIG,C#沟通C++的桥梁
  3. The table(CF226D)
  4. Zedboard安装桌面系统ubuntu及opencv(2)
  5. 【git】Git 提示fatal: remote origin already exists 错误解决办法
  6. python连接redis集群如何释放内存_python 连接redis集群
  7. matlab波形反白,基于MATLAB的海岸污染物浓度扩散实验分析
  8. IP65防水步进电机可以在什么环境使用?
  9. idea运行maven:No URLs will be polled as dynamic configuration sources
  10. MySQL memo优化_mysql memo
  11. Puppet自动化Nginx+Mongrel负载均衡配置
  12. 新版itunes添加铃声
  13. 如何免费使用阿里云服务器?【一篇文章教会你,真香】
  14. 硬件模块应用之超声波测距模块SRF05应用
  15. gRPC-拦截器简单使用
  16. [论文阅读笔记14]Nested named entity recognition revisited
  17. 枸杞子泡水喝有哪些功效
  18. jQuery基础文档(持续更新)
  19. Redis 存储List对象
  20. 计算机图形学学习记录 系统及图形硬件

热门文章

  1. restore和recover 区别
  2. 交互设计要学什么,先弄懂这52个概念再说
  3. Beam Search 及5种优化方法
  4. Apache Drill简介
  5. ISME Commun:RNA-seq助力发现菌群收敛新机制
  6. python操作word题注样式_题注在论文图表编号中的使用
  7. sorthern杂交_您必须在动物杂交上做的10件事:新视野
  8. 你们猜猜,哪个省的女生最爱丰胸?文末有彩蛋!
  9. 重生大玩家-王者归来
  10. 客厅窗帘要安装纱帘吗?怎么选择纱帘?-好佳居窗帘十大品牌