在特定点求解隐函数

Question

我正在尝试求解一个隐式函数，并将值写入一个.txt文件

当x= 0，1且u= - 0.5，1.5时，u=0.5+ sin(x - 0.75*u)

我必须在特定的x点处获得相应的u值(0 0.1 0.2 0.3...1.0)。有什么办法吗？

我正在绘制隐式函数，如下

bla = ezplot('y - sin(2*pi*x - 0.4*y)');
h = gca
get(h)
tmp = get(bla,'contourMatrix');
h = gca
get(h)
set(gca,'ylim',[-1.0 1.0])
set(gca,'xlim',[0 1.0])

Answer 1

可以使用interp1进行插值，以使用tmp矩阵查找给定X值的相应Y值。

XValues = 0:0.1:1;
YValues = interp1(tmp(1,:), tmp(2,:), XValues);

默认的插值方法是linear。您可以查看其他可能的方法here。