add additional parameter to contour plot to generate 3d plot

Question

Richard 2013-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/63931-add-additional-parameter-to-contour-plot-to-generate-3d-plot

I have an equation which has 3 unknown variables V, s, and T. I want to produce a contour plot that shows what happens when each of these variables changes. Currently I can do this by creating a vector for both V and T but keep s as a constant:

V = 1:250; 
T = linspace(1,25,length(V));
s = 30;
for j = 1:length(T)
    D(j,:) = (3.*V.*(s^2))./(T(j).*((T(j).^2)-...
        (3.*T(j).*s)+(3.*(s.^(2)))));
end

Then by plotting D as a contour plot:

pcolor(T,V,D');shading interp

we can see that when V is large and T is small, D will be large. From here, I would like to generate a 3d plot where I can also show the influence of changing 's'. At the moment 's' is a constant, so I can show how s affects the outcome by producing say 4 different figures for 4 different values of s. However, it would be better if I could generate one plot with 3 axis showing how D varies with V, T, and s. Could anyone provide some information about the best way of doing this?

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

请先登录，再进行评论。

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

Answer 1

José-Luis 2013-2-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/63931-add-additional-parameter-to-contour-plot-to-generate-3d-plot#answer_75490

在 MATLAB Online 中打开

What you are asking is basically how to plot four-dimensional data in a two dimensional plane. It is not easy, and tends to be messy, IMO. They way you do it sounds reasonable to me, but if you really want to have everything in a single plot, you could always use the slice() function:

V = 1:250; 
T = linspace(1,25,length(V));
ii = 10:10:40;
all_data = nan(250,250,4);
counter = 1;
for s = ii
    for j = 1:length(T)
        D(j,:) = (3.*V.*(s^2))./(T(j).*((T(j).^2)-...
            (3.*T(j).*s)+(3.*(s.^(2))))); %This could be vectorized
    end
    all_data(:,:,counter) = D';
    counter = counter + 1;
end
[x y z] = meshgrid(V,T,ii);
sH = slice(x,y,z,all_data,[],[],ii);
set(sH,'EdgeColor','none');