explict (pointwise) function from numerical solution to implicit equation

Question

Laurence Kranich 2020-6-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543146-explict-pointwise-function-from-numerical-solution-to-implicit-equation

评论： Ameer Hamza 2020-6-11

I can plot the numerical solutions to a complicated expression f(x,y)=0 over intervals for x and y using fimplicit(f,interval). However, how do I then define/specify an explicit relation y=g(x) pointwise using the numerical values computed from f(x,y)=0 for use in a separate expression? (There is a functional relation.)

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

请先登录，再进行评论。

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

Answer 1

Ameer Hamza 2020-6-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/543146-explict-pointwise-function-from-numerical-solution-to-implicit-equation#answer_447176

编辑：Ameer Hamza 2020-6-6

在 MATLAB Online 中打开

Try interp1 to create an explicit relation for 'y' in term of 'x' using the datapoints from fimplicit

f = @(x, y) y - x.^2;
fp = fimplicit(f, [-2 2 0 10]);
x = fp.XData;
y = fp.YData;
g = @(xq) interp1(x, y, xq);

Then run

>> g(1)
ans =
    1.0002
>> g(0.5)
ans =
    0.2501

Alternate solution:

Use fsolve() directly. The advantage is that you can control the spacing of x-points

f = @(x, y) y - x.^2;
x = -2:0.01:2;
y = fsolve(@(y) f(x, y), rand(size(x)));
g = @(xq) interp1(x, y, xq);

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

Ameer Hamza 2020-6-11

在 MATLAB Online 中打开

Laurence's comment posted as answer moved here:

Ameer,

As I said, your first suggestion worked. (Thank you again!) However, I am now trying to generalize this method verbatim to the case of three variables and am having trouble. My attempt is included below. Even if the 3-D commands are not well-defined, it is clear what I am attempting to do. That is, I am trying to solve f(x,y,z)=0 pointwise for an explicit function z=g(x,y), for which I can then store the data and use in the expression for dtdp, below. Do you have any suggestions?

Laurence

f = @(x,y,z)(2.*sqrt(z).*(1+sqrt(z)).*((y+sqrt(z).^2).*((8.*(1+x))./((1-x).*(3+x).^2))+((8./9).*(1-y).^2)+(2.*(y+sqrt(z)).*(1-y).*(8./(3-x).^2))))-((1-z).*(1-y).*((y+sqrt(z)).*(8.*(1+x)./((1-x).*(3+x).^2))-((1-y).*(8./9))+((1-y).*(8./((3-x).^2)))-((y+sqrt(z)).*(8./((3-x).^2)))));
interval = [0.0001 .9999 0.0001 .9999 0.0001 .9999];
fp = fimplicit3(f, interval); 
x = fp.XData ; 
y = fp.YData ; 
z = fp.ZData ; 
g = @(xq, yq) interp2(x, y, z, Xq, Yq ); 
dtdp = (((1-z).*(1-y))-((1-z).*(y+sqrt(z)))-(4.*sqrt(z)).*(1+sqrt(z)).*(y+sqrt(z)).*(((8.*(1+x))./((1-x).*((3+x).^2)))-(8./((3-x).^2)))-(2.*(1-z).*(1-y))-((1-z)+(4.*sqrt(z)).*(1+sqrt(z)).*(1-y)).*((8./((3-x).^2))-(8./9)))./((((1+sqrt(z))./(sqrt(z))).*((y+sqrt(z)).^2).*((8.*(1+x))./((1-x).*((3+x).^2)))+((y+sqrt(z)).^2).*((8.*(1+x))./((1-x).*((3+x).^2)))+((2.*(1+sqrt(z)).*((y+sqrt(z))).*((8.*(1+x))./((1-x).*((3+x).^2)))+((1-y)).*(y+sqrt(z)).*((8.*(1+x))./((1-x).*((3+x).^2)))-(((1-z)./(2.*(sqrt(z))))).*(1-y).*((8.*(1+x))./((1-x).*((3+x).^2)))+(((1+sqrt(z))./(sqrt(z))).*(1-y).^2).*(8./9)+(((1+sqrt(z))./(sqrt(z)))).*2.*(y+sqrt(z)).*(1-y).*(8./((3-x).^2))+((y+sqrt(z)).*(1-y))*(8./((3-x).^2))+(2.*(1+sqrt(z)).*(1-y))*(8./((3-x).^2))+((y+sqrt(z)).*(1-y))*(8./((3-x).^2))+((1-y).^2).*(8./((3-x).^2))+((2.*(1+sqrt(z)).*((y+sqrt(z))).*(8./((3-x).^2))))))));
plot3 (x,y,dtdp);
xlabel('alpha')
ylabel('po')
zlabel('dtdp')

Laurence Kranich 2020-6-11

That worked. Thanks again!!

Ameer Hamza 2020-6-11

I am glad to be of help!

请先登录，再进行评论。

explict (pointwise) function from numerical solution to implicit equation

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

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

explict (pointwise) function from numerical solution to implicit equation

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

采纳的回答

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

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

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

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