Mode identification on a Graph.

Question

mpho bosupeng 2021-4-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/812190-mode-identification-on-a-graph

回答： Clayton Gotberg 2021-4-24

download.jpg

Hello everyone,

I am a beginner with MATLAB and i need some assistance. I need to produce a graph similar to this one attached. I want to include the vertical dotted lines and the labelling on the graph (HF, IF, AND Trend). My question is how to I do this, I only know that this will need some threshold.

Thanks.

3 个评论
显示 1更早的评论隐藏 1更早的评论

mpho bosupeng 2021-4-24

编辑：mpho bosupeng 2021-4-24

Thanks @Clayton Gotberg. The code is

plot(time,signal)

xlabel('year')

ylabel('PIL')

IMFH=IMF_R(1:end,2*L+1:end-2*L);

NumIMFH=size(IMFH,1);

sumIMFH=0;

for j=1:NumIMFH

sumIMFH=IMFH(j,:)+sumIMFH;

STH(j)=nanmean(sumIMFH)/std(sumIMFH);

indiH(j)=j;

end

figure

plot(indiH(1:end-1),STH(1:end-1),'ok')

ylabel('\mu/\sigma')

xlabel('c_j')

I need to analyze the points to find the regions

Adam Danz 2021-4-24

You can also add the vertical lines using xline()

> I need to analyze the points to find the regions

Ok, let us know exactly what you need help with.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Clayton Gotberg 2021-4-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/812190-mode-identification-on-a-graph#answer_683880

在 MATLAB Online 中打开

It looks like you already have the code to find out where the divisions should be. In that case, all that remains to be solved is:

The plot is missing the vertical lines to separate the regions
The plot doesn't have the region labels

HF_IF_divider = ; % this value should be the point in x at which the HF and IF regions touch
IF_Trend_divider = ; % this value should be the point in x where IF and Trend meet
plot(data) % replace with all of the plotting you've done so far
hold on
% Add the vertical lines
xline(HF_IF_divider,'k--') % Thanks @Adam Danz for mentioning this great function
xline(IF_Trend_divider,'k--') % 'k' means black lines, '--' means dashed lines
label_x = [1/2*HF_IF_divider, 1/2*(IF_Trend_divider+HF_IF_divider),...
    1/2*(length(signal)+IF_Trend_divider)]; % x positions of labels
label_y = [0.5 0.5 0.5]; % Change this to modify the height of the label
zone_names = ["HF" "IF" "Trend"];
text(label_x,label_y,zone_names,'HorizontalAlignment','center','VerticalAlignment','middle','FontSize',12);