Is there a way to plot everything together using the surf command in matlab?

Question

Gregory Drake 2019-4-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/456306-is-there-a-way-to-plot-everything-together-using-the-surf-command-in-matlab

评论： Star Strider 2019-4-14

Is there a way to plot everything together using the surf command in matlab? right now I can only one z value at a time (i.e.

surf(R,Rc,N_diesel(:,:,1)); , surf(R,Rc,N_diesel(:,:,2));, surf(R,Rc,N_diesel(:,:,3)); ... and so forth)

Can someone point me in the right direction?

My Current Code:

--------------------------------------------------------------------------------------------

T = linspace(250,1000,16); 
T_size = length(T);
Rc = linspace(1.001,4.001,7);
rc_size = length(Rc);
R = linspace(12,24,25);
r_size = length(R);
N_diesel = zeros(r_size,rc_size,T_size);
for ki=1:1:T_size
    for j=1:1:rc_size
        for i=1:1:r_size
        k = get_k_air(T(ki));
        Calc_Diesel = efficiency(R(i),Rc(j),k);
        N_diesel(i,j,ki) = Calc_Diesel;
        end
    end
end
[Rc,R] = meshgrid(Rc,R);
surf(R,Rc,N_diesel(:,:,1));
x = array2table(round(N_diesel(:,:,1),3,'significant'));
function k = get_k_air(T1)
    keys = linspace(250,1000,16);
    values = [1.401,1.4,1.398,1.395,1.391,1.387,1.381,1.376,1.37,1.364,1.359,1.354,1.349,1.344,1.34,1.336];
    lookup_table = containers.Map(keys,values);
    k = lookup_table(T1);
end
function Calc_Diesel = efficiency(R,Rc,k)
    N1 = 1/(R.^(k - 1));
    N2 = (Rc.^k - 1)/(k.*(Rc - 1));
    Calc_Diesel = 1 - (N1*N2);
end

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-4-14

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/456306-is-there-a-way-to-plot-everything-together-using-the-surf-command-in-matlab#answer_370584

在 MATLAB Online 中打开

Use the hold function, and plot the surfaces in a loop.

Try this:

figure
hold all
for k1 = 1:size(N_diesel,3)
    surf(R,Rc,N_diesel(:,:,k1));
    x{k1} = array2table(round(N_diesel(:,:,k1),3,'significant'));
end
hold off
view(30,30)
grid on

Full context:

T = linspace(250,1000,16); 
T_size = length(T);
Rc = linspace(1.001,4.001,7);
rc_size = length(Rc);
R = linspace(12,24,25);
r_size = length(R);
N_diesel = zeros(r_size,rc_size,T_size);
for ki=1:1:T_size
    for j=1:1:rc_size
        for i=1:1:r_size
        k = get_k_air(T(ki));
        Calc_Diesel = efficiency(R(i),Rc(j),k);
        N_diesel(i,j,ki) = Calc_Diesel;
        end
    end
end
[Rc,R] = meshgrid(Rc,R);
figure
hold all
for k1 = 1:size(N_diesel,3)
    surf(R,Rc,N_diesel(:,:,k1));
    x{k1} = array2table(round(N_diesel(:,:,k1),3,'significant'));
end
hold off
view(30,30)
grid on
function k = get_k_air(T1)
    keys = linspace(250,1000,16);
    values = [1.401,1.4,1.398,1.395,1.391,1.387,1.381,1.376,1.37,1.364,1.359,1.354,1.349,1.344,1.34,1.336];
    lookup_table = containers.Map(keys,values);
    k = lookup_table(T1);
end
function Calc_Diesel = efficiency(R,Rc,k)
    N1 = 1/(R.^(k - 1));
    N2 = (Rc.^k - 1)/(k.*(Rc - 1));
    Calc_Diesel = 1 - (N1*N2);
end