Matlab: точка пересечения кривых

Question

Скажем, например, у меня есть данные, которые формируют параболическую кривую y = x^2, и я хочу прочитать значение x для заданного значения y. Как мне это сделать в MATLAB?

Если бы это была прямая линия, я мог бы просто использовать уравнение линии наилучшего соответствия, чтобы легко вычислить, однако я не могу сделать это с изогнутой линией. Если я не смогу найти решение, я решит для корней

Заранее спасибо.

Answer 1

Если все данные являются массивами (не аналитические выражения), я обычно делаю, что нахождение минимальной абсолютной погрешности

x=some_array; 
[~,ind]=min(abs(x.^2-y0)) 

Здесь y0 является данное значение у

Если данные представлены в виде функции, вы можете использовать fsolve:

function y = myfun(x) 
    y=x^2-y0 

[x,fval] = fsolve(@myfun,x0,options) 

Для символьных вычислений, можно использовать solve

syms x 
solve(x^2 - y0) 

Answer 2

Предполагая, что ваши две кривые - это всего лишь два вектора данных, я бы предложил вам использовать Fast and Robust Curve Intersections из Файлового обмена. См. Также эти два похожих вопроса: how to find intersection points when lines are created from an array и Finding where plots may cross with octave/matlab.