Question

用matlab，已知x的值，和f(x,y)=0的隐函数，如何求y

已知x=linspace(-2.23,2.23,100);
和 隐函数sin(2^0.5*y)-0.5*sin(0.5*2^0.5*(x-y)).*sin(0.5*2^0.5*(x+y))=0，如何用matlab求与x相对应的y的值。

Answer 1

h=ezplot('sin(2^0.5*y)-0.5*sin(0.5*2^0.5*(x-y)).*sin(0.5*2^0.5*(x+y))=0',[-2.23,2.23 0 1]);
x1=get(h,'XData');
y1=get(h,'YData');
x=linspace(-2.23,2.23,100);
y=interp1(x1,y1,x);
hold on
plot(x,y,'r')

使用ezplot函数来实现吧。先画出图像，然后图像上取点进行插值。当然这样可能不够精确。

要想更精确就给定一个x，然后只用数值方法求y。

追问

请问数值方法怎么做？我用subs(solve())，matlab一直busy，不知道是我的电脑太慢了，还是方程太复杂了。

追答

是你的方法不好
f=@(x,y)sin(2^0.5*y)-0.5*sin(0.5*2^0.5*(x-y)).*sin(0.5*2^0.5*(x+y));
x1=linspace(-2.23,2.23,100);
for i=1:length(x1)
x=x1(i);
y(i)=fsolve(@(y) f(x,y),0);
end
plot(x1,y)

追问

我输入
f=@(x,y)sin(2^0.5*y)-0.5*sin(0.5*2^0.5*(x-y)).*sin(0.5*2^0.5*(x+y));
x1=linspace(-2.23,2.23,100);
for i=1:length(x1)
x=x1(i);
y(i)=fsolve(@(y) f(x,y),0);
end
后出现错误：
Optimization terminated: first-order optimality is less than options.TolFun.
请问该怎么办？

追答

你的matlab是什么版本啊？

追问

我看到你求的结果，y的个数是x的2倍。。。

追答

不可能的，你到底是怎么运行的啊，我都试过了的，肯定没问题