求助高手，该如何画这个三维图形？

alexqxp

iceberg_jrh提出了一个好问题，weiyan图比较有意思。感谢两位的精彩讨论。
lyrock 和bainhome不打不相识。
iceberg_jrh抱歉，没有找到那篇文章。我只记得文章中点应力在空间的分布是椭圆球形的，有点象纺锤。

iceberg_jrh

谢谢alexqxp版主，感谢你为找那篇文章花费你不少的宝贵时间。真是非常感谢。至于我这个问题我会继续去深究的，希望能尽快找到答案。也感谢weiyan，lyrock和bainhome，你们的热心让我备受鼓舞。感谢各位回答我的问题，谢谢

weiyan

在（a1,a2,a3)画a1-a2=±2k or a1-a3=±2k or a2-a3=±2k这六个平面，它们所围的一个闭曲面就是一个六面柱体。

[X,Y] = meshgrid(-2:.125:2);  %在X,Z平面上投影是直线z=x+1的平面
Z=X+1;
surf(X,Y,Z);
hold on;
Z=X-1;
surf(X,Y,Z);

hold on;
[X,Y] = meshgrid(-2:.125:2);
Z=Y+1;
surf(X,Y,Z);
hold on;
Z=Y-1;
surf(X,Y,Z);

hold on;
[X,Z] = meshgrid(-2:.125:2);
Y=X+1;
surf(X,Y,Z);
hold on;
Y=X-1;
surf(X,Y,Z);

iceberg_jrh

谢谢weiyan，已经向成功迈出了一步了。但是我想问的是我只想要这个封闭六面体，如何去掉这多出来的部分呢？要是用PhotoShop处理似乎不好控制。还有就是我想再请教各位一个问题，画一个平面的圆是很好画的，对吧？用极坐标就可以了。但是我想画个空间的圆该如何画呢？比如说画这个函数的图形，x+y+z＝0，我想将它的分为定义成一个圆形的区域，matlab能否实现。

bainhome

这个问题的原理大家都知道,对图形裁减可以用NaN赋值来实现,但是刚才鼓捣了一阵,没有找到合适的约束条件,bzzz,nostalgica,你们如果有空给做做这个题目:如何将封闭六面体时的多余部分去掉,感觉有点意思.^_^
至于第二个问题,你参考这个程序,如果这个都能画,我想你的那个应该有相似的地方
function cone_plot(m,x0,y0,z0)
% (x0,y0,z0)---顶点坐标
% m---底圆在y轴的坐标位置.
x1=linspace(-5,-eps,100);
x2=linspace(eps,5,100);
x=[x1,x2];
z1=sqrt(25-x.^2);
z2=-sqrt(25-x.^2);
x0=x0*ones(1,length(z1));
y0=y0*ones(1,length(z1));
z0=z0*ones(1,length(z1));
y=m*ones(1,length(z1));
for i=1:length(z1)
plot3([x0(i),x(i)],[y0(i),y(i)],[z0(i),z1(i)])
hold on
plot3([x0(i),x(i)],[y0(i),y(i)],[z0(i),z2(i)])
end
grid on

iceberg_jrh

谢谢bainhome的热心回答，你让我见识到了matlab的强大功能。谢谢，也非常佩服你能把matlab用的这么好。我自己再琢磨琢磨。

mn5207

作者还在吗
我也在思考这个问题

refeihc

作者可能不在了，还有人感兴趣吗？这个问题其实不难。

mn5207

楼上的，你的联系方式是什么，我的qq67665265
请加我为好友