Как автоматически определить, где график меняет уклон, чтобы вы могли соответствовать -Matlab - PullRequest
Новая
       49

Как автоматически определить, где график меняет уклон, чтобы вы могли соответствовать -Matlab

1 голос
/ 13 марта 2020

Следующий код вычисляет и строит график соответствия графика. Проблема в том, что каждый раз, когда я должен найти, где наклон меняется самостоятельно. См. Следующий график:

enter image description here

Есть ли способ найти его автоматически? 

start=50;             rd1start=1;
ending=100;             rd1end=100;

relative=zeros(1,100);
W_esc_half=linspace(1,100,100);                         % x axis values
relative(1:50)=W_esc_half(1:50).^(1.2);       % y axis values  fof smaller times

relative(50:100)=W_esc_half(50:100).^(1.5);      % y axis values for greater times

figure(1)
fitResults1 = polyfit(log10(W_esc_half(start:ending)),log10(relative(start:ending)),1);
pol=polyval(fitResults1,log10(W_esc_half(start:ending)));                                 % Compute the fit coefficients
a=fitResults1(1);
b=fitResults1(2);
polyfit_str = ['<(?r)^2> ~ ? ^{'  (num2str(a)) ' } '] ;


fit1=W_esc_half(start:ending).^a*10^(b);


hold on
loglog((W_esc_half(rd1start:rd1end)),(relative(rd1start:rd1end)),'blue-','LineWidth',2)  % draw the original function
loglog((W_esc_half(start:(length(pol)+start-1))),fit1,'cyan--','LineWidth',2);            % draw the fit

1 Ответ

2 голосов
/ 13 марта 2020

Мне обычно нравится делать это быстро, отыскивая в diff сигнал 

subplot(2,1,1)
plot(diff(relative),'.')
subplot(2,1,2)
findpeaks(diff(relative))

Но у вас есть медленно увеличивающийся наклон в дополнение к резкому изменению. Поэтому, если ваша функция не так хороша, вам может понадобиться немного настроить функцию findpeaks. С другой стороны, если у вас постоянные уклоны и смена аббатта, вы также можете go с find(abs(diff(relative)) > 1).

...