FInding the breakpoint of a frequency spectra and its associated error

Question

Mahith Madhana Kumar 2021-2-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/732218-finding-the-breakpoint-of-a-frequency-spectra-and-its-associated-error

评论： Mahith Madhana Kumar 2021-2-1

Hi all,

I was wondering whether somebody could help me in finding the first breakpoint of a frequency spectra and its associated error. The spectra is as shown in the attached figure (spectrum.jpg). The spectrum should be fitted as shown in spectrum_fitted.jpg inorder to find the first breakpoint (labelled as A) so as to get the corresponding frequency which is about 1.2Hz (in figure spectrum_fitted.jpg). The fit might not be accurate and hence there would be an associated error in the value obtained for the frequency at A. So I would need both the value for frequency at A( first breakpoint) as well as its associated error. Wondering if that is possible.

The data used to get the frequency plot has been attached as data.mat

Thank you.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Mathieu NOE 2021-2-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/732218-finding-the-breakpoint-of-a-frequency-spectra-and-its-associated-error#answer_611413

在 MATLAB Online 中打开

hello

I combinde this

Piecewise linear least square fit - File Exchange - MATLAB Central (mathworks.com)

and fminsearch to create that code - seems to work finding the optimal values for break points

% load x and y data
load('data.mat');
xdata = xdata(:);
ydata = ydata(:);
% Init of fminsearch 
a = max(xdata)/4;
b = max(xdata)*3/4;
global xdata ydata
x0 = [a,b];
x = fminsearch(@obj_fun, x0);
a_sol = x(1);
b_sol = x(2);
XI = [min(xdata),a_sol,b_sol,max(xdata)]; % vector of 1-D look-up table "x" points
YI = lsq_lut_piecewise( xdata, ydata, XI ); % obtain vector of 1-D look-up table "y" points
% plot fit
plot(xdata,ydata,'.',XI,YI,'+-')
legend('experimental data (x,y(x))','LUT points (XI,YI)')
title('Piecewise 1-D look-up table least square estimation')
function err = obj_fun(x)
global xdata ydata
    XI = [min(xdata),x(1),x(2),max(xdata)];  % vector of 1-D look-up table "x" points
    YI = lsq_lut_piecewise( xdata, ydata, XI ); % obtain vector of 1-D look-up table "y" points
    Yint = interp1(XI,YI,xdata);
    err = sum((Yint-ydata).^2);
end