Finding Max value on peak point of graph (logistics growth model)

Question

Nirmal Kishore Janakiraman 2021-3-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/784986-finding-max-value-on-peak-point-of-graph-logistics-growth-model

评论： Nirmal Kishore Janakiraman 2021-4-1

Is there a method in matlab where i can find the value of the point in the graph (red arrow) before it plateues. I was thinking of using the gradient function. But unable to get any. Some help please !

2 个评论
显示无隐藏无

dpb 2021-3-27

Have more success if attach the data file, probably, so folks have something to work from.

save the variable and attach the .mat file

Nirmal Kishore Janakiraman 2021-3-30

编辑：Nirmal Kishore Janakiraman 2021-3-30

Thank you sir. I have attached below and what is it that im trying to find.

请先登录，再进行评论。

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

Answer 1

Ananya Tewari 2021-3-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/784986-finding-max-value-on-peak-point-of-graph-logistics-growth-model#answer_662319

在 MATLAB Online 中打开

I understand you want the coordinates of the point after which the graph plateaus. Assuming you are having the data points of the graph, using movmax we can get the desired result.

maxi = movmax(data,[0 2]);  % Taking 3 values at a time and finding maximum of them
idx = 0;  % this will capture the index after which the graph plateaus
for k = 1:n  % where n is number of data points
    % when the middle value of the 3 elements is the largest we break the loop
    if(maxi(k) == data(k+1))  
        idx = k;
        break;
    end
end

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

Ananya Tewari 2021-3-31

在 MATLAB Online 中打开

Hi Nirmal,

The ischange() finds abrupt changes in a given data and this function can be leveraged to find the approximate point after which the graph plateaus.A possible workaround is as follows:

data = biomass_conc(1,:) %data-points
x=1:5000
% Find change points
changeIndices = ischange(biomass_conc,'linear','MaxNumChanges',20,...
    'SamplePoints',x);
% Display results
clf
plot(x,biomass_conc,'Color',[109 185 226]/255,'DisplayName','Input data')
hold on
% Plot change points
x2 = repelem(x(changeIndices),3);
y = repmat([ylim(gca) NaN]',nnz(changeIndices),1);
plot(x2,y,'Color',[51 160 44]/255,'LineWidth',1,'DisplayName','Change points')
title(['Number of change points: ' num2str(nnz(changeIndices))])
hold off
% The last value of x2 contains the approximate x-coordinate after which
% the graph plateaus
x_dash = x2(end)   % x-coordinate
y_dash = biomass_conc(x_dash)  % corresponding y-coordinate