matlab隐函数求解问题

matlab小菜鸡寻求大家的帮助~
我的问题如下：

  • 首先，x、y是未知常数，0

img

  • 接下来，我需要把x、y、n_B_0代入到第二个等式，希望最终能画出n_B和x、y之间的三维函数图

img

由于式子比较复杂，我无法直接用x、y写n_B_0的表达式，想着只能通过隐函数来求解n_B_0。我自己尝试了一下，画出了n_B_0和x、y之间的函数关系图，但是不知道怎么将这个n_B_0代入第二个式子求解n_B和x、y之间的关系图

下面附上我画n_B_0和x、y三维函数图的代码，不知道有没有帮助

clear;
clc;
close all;

%% 设置变量和参数
delta = 0.00001;
step = 0.01;

%% 设置函数
f = @(x,y,n_B) (log(1-x)-log(1-y))*(5/27*((1-x).^(3-n_B))*y*((1-y).^(n_B-1))+...
    32/27*((1-x).^(3-n_B))*((1-y).^(n_B)))+32/27*(((1-y).^(n_B))*...
    log(1-y)-((1-x).^(3-n_B))*log(1-x))+5/27*y*((1-y).^(n_B-1))*log(1-y);

%% 定义取值范围
interval = [0+delta 1-delta 0+delta 1-delta 0 3];

%% 绘制函数图像
fimplicit3(f,interval,'EdgeColor','none','FaceAlpha',.5);

%% 绘制多个分界面
% 采用 y = x 进行分割
hold on;
% 绘制分界面
f1 = @(x,y,z) (x-y);
% 定义取值范围
interval = [0 1 0 1 1 3];
% 绘制函数图像
fimplicit3(f1,interval,'EdgeColor','none','FaceAlpha',.5,'FaceColor','k');

%% 设置图形参数
title('函数三维图形')
xlabel('x');
ylabel('y');
zlabel('n_B');

由于第二个式子比较复杂，也附上这个式子的代码~
%% 第二个隐含式求解

f = @(x,y,n_B) ((1-x).^(3-n_B))*((1-y).^n_B)*(log(1-x)-log(1-y))*(1-(1-x).^(3-n_B_0))*(y*((1-y).^(n_B_0-1))+32/27*(1-((1-y).^n_B_0)-y*((1-y).^(n_B_0-1))))+(1-(1-x).^3)*((1-x).^(3-n_B))*(((1-y).^(n_B-1))*(log(1-y)-log(1-x))*(32/27-y)+32/27*log(1-x))+((1-x).^(3-n_B))*y*((1-y).^(n_B-1))*(log(1-y)-log(1-x))+32/27*((1-x).^(3-n_B))*((1-y).^n_B)*(log(1-x)-log(1-y))+32/27*y*((1-y).^(n_B-1))+32/27*((1-x).^(3-n_B))*y*((1-y).^(n_B-1))*(log(1-x)-1)-32/27*((1-x).^(3-n_B))*log(1-x)-y*((1-y).^(n_B-1))*log(1-y)+32/27*((1-y).^n_B)*log(1-y)+(32/27-32/27*((1-y).^3)-5/27*y*((1-y).^2))*((1-y).^n_B)*(((1-x).^(3-n_B))*(log(1-y)-log(1-x))-log(1-y));

总结一下，大致意思就是0