Matlab 3 variable function plot

The equation is:

f(X1,X2,X3)=bt1.x1+bt2.X2+bt3.X3+b1

bt1 to bt3 and b1 are all constants.

and I want to plot it in 3D. I tried a couple of functions like:

figure
syms x1 x2 x3
fimplicit3(Beta(1)*x1+Beta(2)*x2+Beta(3)*x3+b(1))

but it retunrs nothing.

[Tommy]

fimplicit3 wants a function handle. Try this:

fimplicit3(@(x1,x2,x3) Beta(1)*x1+Beta(2)*x2+Beta(3)*x3+b(1))

[Zeynab Mousavikhamene]

Again returned nothing.

[Tommy]

Hmm, are you able to provide all of the code which you are running?

[Zeynab Mousavikhamene]

Beta is:

Beta(1)= -294449.131783462

Beta(2)=14.7170998874722

Beta(3)=-0.127560549560172

b(1) is

87293272725.0805

f=@(x1,x2,x3) Beta(1).*x1+Beta(2).*x2+Beta(3).*x3+b(1);
fimplicit3(f)

[Tommy]

I believe nothing shows because there are no solutions within the default interval [-5 5]. But yes, maybe I incorrectly assumed you were trying to plot solutions to f=0.

[Zeynab Mousavikhamene]

You are right I am plotting solution for f=0. I changed the interval and did not show anything.

Answer 1

[Ameer Hamza]

fimplicit3 is used to plot implicit equations with three variables. What you have is 4D data (3 input variables, 1 output variable). You need to use a 4D visualization function, like slice(), to visualize your function. See this example

Beta(1)= -294449.131783462;
Beta(2)=14.7170998874722;
Beta(3)=-0.127560549560172;
b(1) = 87293272725.0805;

f=@(x1,x2,x3) Beta(1).*x1+Beta(2).*x2+Beta(3).*x3+b(1);

[X1,X2,X3] = meshgrid(linspace(-1,1));
V = f(X1,X2,X3);
slice(X1, X2, X3, V, [-0.5 0.5], 0.3, 0)
colorbar
shading interp

[Zeynab Mousavikhamene]

@Ameer Hamza

Thanks for the explanation. It is 3D not 4D.I want to plot f=0.

[Ameer Hamza]

If you want to solve it for f=0, then your equation is linear, and you can write an expression of x3 in terms of x1 and x2 using a symbolic toolbox and then use it to plot a surface. For example

Beta(1)= -294449.131783462;
Beta(2)=14.7170998874722;
Beta(3)=-0.127560549560172;
b(1) = 87293272725.0805;

syms x1 x2 x3
f = Beta(1).*x1+Beta(2).*x2+Beta(3).*x3+b(1);
x3_sol = solve(f, x3);
x3_sol_f = matlabFunction(x3_sol, 'Vars', {x1 x2});

[X1, X2] = meshgrid(-5:0.2:5);
X3 = x3_sol_f(X1, X2);

surf(X1, X2, X3);
shading interp

[Zeynab Mousavikhamene]

when I use zlime even in the surf case, the plot show nothing.

[Ameer Hamza]

I haven't used zlim in my code in the comment. Have you tried running that code?

Answer 2

[Tommy]

You can pick any interval. Your plot will only show something if solutions to f=0 lie within the interval. [0, 0, -b1/Beta(3)] is a clear solution. -b1/Beta(3) is on the order of 1e11, so how about this:

Beta(1)= -294449.131783462;
Beta(2)=14.7170998874722;
Beta(3)=-0.127560549560172;
b(1)=87293272725.0805;
f=@(x1,x2,x3) Beta(1).*x1+Beta(2).*x2+Beta(3).*x3+b(1);
fimplicit3(f, [-5 5 -5 5 1e10 1e12])