Minimize distance between curves

Question

Deepa Maheshvare 2022-7-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1753325-minimize-distance-between-curves

评论： Deepa Maheshvare 2022-7-7

I have a dataset with values of multiple curves. An example plot is shown below. I want to scale the curves so that all curves overlap or get closer to the target curve. The following is the plot of a sample dataset

I want to scale the function values in the curves so that all curves overlap. Say, green curve is the target (g) and I want to move the other curves (f) towards target curve by scaling (scale facotor is a).

The function that I want to minimize is:

Expected result:

Data:

data = {{{157, 60}, {140, 57}, {121, 52}, {103, 47}, {86, 39}, {67, 28}, {50,21}, {32, 15}}, {{159, 70}, {140, 67}, {122, 63}, {105, 58}, {85, 45}, {68, 42}, {50, 34}, {33, 29}}, {{157, 74}, {140, 71}, {121, 70}, {105, 65}, {86, 55}, {69, 52}, {52, 43}, {32, 37}}, {{158, 100}, {141, 100}, {123, 95}, {105, 90}, {88, 78}, {69, 72}, {51, 57}, {33, 46}}, {{160, 143}, {140, 146}, {123, 146}, {105, 145}, {86, 136}, {70, 121}, {51, 97}, {34, 60}}}

Could someone please help me with the function that has to be used for solving this problem in MATLAB ?

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

请先登录，再进行评论。

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

Answer 1

Torsten 2022-7-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1753325-minimize-distance-between-curves#answer_999960

编辑：Torsten 2022-7-4

在 MATLAB Online 中打开

Differentiating with respect to a, setting the derivative to 0 and solving for "a" gives

a = sum_i (f_i*g_i) / sum_i (f_i^2)

I leave it to you to insert the data points.

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

Torsten 2022-7-6

编辑：Torsten 2022-7-6

在 MATLAB Online 中打开

What is the "absolute error" ?

sum_i abs(a*f_i - g_i)

?

If this is the case, use "linprog" to solve

min: sum_i u_i
s.c.
f_i * a - u_i <= g_i
-f_i * a - u_i <= -g_i

with unknowns a and u_i (i=1,...,number of data points).

Deepa Maheshvare 2022-7-7

在 MATLAB Online 中打开

@Torsten

Thanks so much. Yes, I mean sum_i abs(a*f_i - g_i)

% Loss Function (absolute error loss)
%% input curves
scale =  1.5;
x1 = [0,4,6,10,15,20]*scale;
y1 =  [18,17.5,13,12,8,10];
x2 = [0,10.5,28]*scale;
y2= [18.2,10.6,10.3];
%% y2 is the target function and y1 is the function to be shifted
f = y1;
g = interp1(x2,y2,x1);
%% linprog inputs
% x = linprog(f,A,b)
% i = 1:length(x1)
A = [fi - 1; -fi -1];
x = [a ui];
b = [gi -gi];
% obj = sum_i u_i (function to minimize)
%% solve system