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

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

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

更多Matlab仿真内容点击👇

智能优化算法  神经网络预测 雷达通信  无线传感器

信号处理 图像处理 路径规划 元胞自动机 无人机

⛄ 内容介绍

Fractal structures are found in biomedical time series from a wide range of physiological

phenomena. The multifractal spectrum identififies the deviations in fractal structure within

time periods with large and small flfluctuations. The present tutorial is an introduction to

multifractal detrended flfluctuation analysis (MFDFA) that estimates the multifractal spec

trum of biomedical time series. The tutorial presents MFDFA step-by-step in an interactive

Matlab session. All Matlab tools needed are available in Introduction to MFDFA folder at the

website www.ntnu.edu/inm/geri/software. MFDFA are introduced in Matlab code boxes

where the reader can employ pieces of, or the entire MFDFA to example time series. After

introducing MFDFA, the tutorial discusses the best practice of MFDFA in biomedical signal

processing. The main aim of the tutorial is to give the reader a simple self-sustained guide

to the implementation of MFDFA and interpretation of the resulting multifractal spectra.

⛄ 部分代码

YMatrix1=[multifractal.*30,RW3];

YMatrix2=[monofractal.*30,RW2];

YMatrix3=[whitenoise.*30,RW1];

X1=2600:3600;

Y1=RW3(2600:3600);

% Create figure

figure1 = figure('PaperType','a4letter','PaperSize',[20.98 29.68],...

    'Color',[1 1 1]);

% Create axes

axes1 = axes('Parent',figure1,'YTickLabel',{'0','200','400','600'},...

    'YTick',[0 200 400 600],...

    'XTickLabel',{},...

    'XTick',zeros(1,0),...

    'Position',[0.13 0.6545 0.7745 0.2705],...

    'LineWidth',2,...

    'FontSize',14);

% Uncomment the following line to preserve the Y-limits of the axes

ylim(axes1,[-220 700]);

hold(axes1,'all');

% Create multiple lines using matrix input to plot

pplot1 = plot(YMatrix1,'Parent',axes1);

set(pplot1(2),'LineWidth',2,'Color',[1 0 0]);

% Create axes

axes2 = axes('Parent',figure1,'YTickLabel',{},'YTick',zeros(1,0),...

    'XTick',zeros(1,0),...

    'Position',[0.6641 0.8264 0.1875 0.1359],...

    'LineWidth',2);

% Uncomment the following line to preserve the Y-limits of the axes

ylim(axes2,[370 570]);

box(axes2,'on');

hold(axes2,'all');

% Create plot

plot(X1,Y1,'Parent',axes2,'LineWidth',2,'Color',[1 0 0]);

% Create axes

axes3 = axes('Parent',figure1,'YTickLabel',{'0','200','400','600'},...

    'YTick',[0 200 400 600],...

    'XTickLabel',{},...

    'XTick',zeros(1,0),...

    'Position',[0.13 0.3833 0.7745 0.272],...

    'LineWidth',2,...

    'FontSize',14);

% Uncomment the following line to preserve the Y-limits of the axes

ylim(axes3,[-220 700]);

hold(axes3,'all');

% Create multiple lines using matrix input to plot

pplot2 = plot(YMatrix2,'Parent',axes3);

set(pplot2(1),'DisplayName','Noise like time series');

set(pplot2(2),'LineWidth',2,'Color',[1 0 0],...

    'DisplayName','Random walk like time series');

% Create ylabel

ylabel('Amplitude (measurement units)','FontSize',14);

% Create axes

axes4 = axes('Parent',figure1,'YTickLabel',{'-200','0','200','400','600'},...

    'YTick',[0 200 400 600],...

    'XTickLabel',{'0','1?00','2?00','3?00','4?00','5?00','6?00','7?00','8?00'},...

    'XTick',[0 1000 2000 3000 4000 5000 6000 7000 8000],...

    'Position',[0.13 0.1108 0.7745 0.2723],...

    'LineWidth',2,...

    'FontSize',14);

ylim(axes4,[-220 700]);

hold(axes4,'all');

% Create multiple lines using matrix input to plot

pplot3 = plot(YMatrix3,'Parent',axes4);

set(pplot3(2),'LineWidth',2,'Color',[1 0 0]);

% Create xlabel

xlabel('time (sample number)','FontSize',14);

% Create legend

legend1 = legend(axes3,'show');

set(legend1,'LineWidth',2);

% Create textbox

annotation(figure1,'textbox',[0.2362 0.6476 0.1934 0.04038],...

    'String',{'Monofractal time series'},...

    'FontSize',14,...

    'FitBoxToText','off',...

    'LineStyle','none');

% Create textbox

annotation(figure1,'textbox',[0.2362 0.9246 0.1822 0.04038],...

    'String',{'Multifractal time series'},...

    'FontSize',14,...

    'FitBoxToText','off',...

    'LineStyle','none');

% Create textbox

annotation(figure1,'textbox',[0.238 0.2453 0.1057 0.04038],...

    'String',{'White noise'},...

    'FontSize',14,...

    'FitBoxToText','off',...

    'LineStyle','none');

% Create line

annotation(figure1,'line',[0.3828 0.3828],[0.8978 0.8278]);

% Create line

annotation(figure1,'line',[0.3828 0.4714],[0.8281 0.8291]);

% Create line

annotation(figure1,'line',[0.4713 0.4713],[0.8294 0.9004]);

% Create line

annotation(figure1,'line',[0.3828 0.4705],[0.8995 0.9004]);

% Create line

annotation(figure1,'line',[0.4705 0.6632],[0.8291 0.8277]);

% Create line

annotation(figure1,'line',[0.4696 0.6667],[0.9017 0.9637]);

% Create line

annotation(figure1,'line',[0.1302 0.1293],[0.9327 0.8977],'LineWidth',4,...

    'Color',[1 1 1]);

clear pplot1 pplot2 pplot3 legend1 axes1 figure1 axes2 figure2 axes3 figure3 ans axes4 X1 Y1 YMatrix1 YMatrix2 YMatrix3

⛄ 运行结果

正在上传…重新上传取消

⛄ 参考文献

[1]Espen Alexander Fürst E.A.F.I. Ihlen. Introduction to Multifractal Detrended Fluctuation Analysis in Matlab[J]. Frontiers, 2012.

❤️ 关注我领取海量matlab电子书和数学建模资料

❤️部分理论引用网络文献，若有侵权联系博主删除