Code covered by the BSD License  

Highlights from
Multifractal detrended fluctuation analyses

image thumbnail

Multifractal detrended fluctuation analyses

by

 

The multifractal

plot1.m
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','1000','2000','3000','4000','5000','6000','7000','8000'},...
    '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

Contact us