函数的表示跟人的自然的逻辑思维之间还是有差距
来源:互联网 发布:锐步椭圆机 知乎 编辑:程序博客网 时间:2024/05/06 07:44
ArcTan是个奇葩函数,Reduce, Simplify不处理。
RegionPlot3D[Sqrt[x^2+y^2+z^2]<=2 &&((0<ArcTan[z,Sqrt[x^2+y^2]]<\[Pi]&&(0<ArcTan[x,y]<\[Pi]||0<ArcTan[-x,-y]<\[Pi]/2))||(0<ArcTan[-z,Sqrt[x^2+y^2]]<\[Pi]/2&&0<ArcTan[x,-y]<\[Pi]/2))||(Sqrt[x^2+y^2+z^2]<=1),{x,-2,2},{y,-2,2},{z,-2,2},PlotPoints->200,Mesh->None,Boxed->False,Axes->False,PlotTheme->"Classic"]
RegionPlot3D[1<=Sqrt[x^2+y^2+z^2]<=2&&((0<ArcTan[z,Sqrt[x^2+y^2]]<\[Pi]&&(0<ArcTan[x,y]<\[Pi]||0<ArcTan[-x,-y]<\[Pi]/2))||(0<ArcTan[-z,Sqrt[x^2+y^2]]<\[Pi]/2&&0<ArcTan[x,-y]<\[Pi]/2)),{x,-2,2},{y,-2,2},{z,-2,2},PlotPoints->160,Mesh->None,Boxed->False,Axes->False,PlotTheme->"Classic"]
Mathematica中绘制3D模型的函数,典型的,比较符合人们思维方式的是 RegionPlot3D。 3D模型可以转stl之类格式的起码有两类,一类是solid model, 一类surface model。
估计 RegionPlot3D 是为了 solid model 而生的。这使得它在绘图的时候虽然方便理解,但是存在诸多缺陷。尤其是棱角的地方,如果PlotPoints选项用默认、或者设定太小,就明显会出现视觉上的缺陷、直至根本不能用;如果设定太大,一般的电脑内存和CPU的情况下,累死也画不出来。
这里上面的例子简单,RegionPlot3D 都能搞定。最近发现,碰上稍微繁琐复杂的情况时,还是要考虑用Surface model,或者ContourPlot3D。这就涉及进一步的,如何把对体solid model的描述转化成对面surface model的描述。
这里有一个独立于built-in-function的很好的例子,让人看了非常开心。改用Surface model之后,大大提高了效率,改进了模型出图的视觉效果,简直是网友们智慧结晶的精品自定义函数。函数名字 contourRegionPlot3D
,还是挺有意义的:
contourRegionPlot3D[region_,{x_,x0_,x1_},{y_,y0_,y1_},{z_,z0_,z1_},opts:OptionsPattern[]]:=Module[{reg,preds},reg=LogicalExpand[region&&x0<=x<=x1&&y0<=y<=y1&&z0<=z<=z1];preds=Union@Cases[reg,_Greater|_GreaterEqual|_Less|_LessEqual,-1];Show@Table[ContourPlot3D[Evaluate[Equal@@p],{x,x0,x1},{y,y0,y1},{z,z0,z1},RegionFunction->Function@@{{x,y,z},Refine[reg,p]&&Refine[!reg,!p]},opts],{p,preds}]]
放个例子:
contourRegionPlot3D[(x < 0 || y > 0 || z < 0) && 0.5 <= x^2 + y^2 + z^2 <= 0.99, {x, -1, 1}, {y, -1, 1}, {z, -1, 1}, Mesh -> None, BoundaryStyle -> None, PlotPoints -> 80, Boxed -> False, Axes -> False, PlotTheme -> "Classic"]
以及
contourRegionPlot3D[-2x+2y-2z+Sqrt[2]>=0&&-2x-2y-2z+Sqrt[2]>=0&&2x-2y-2z+Sqrt[2]>=0&&2x+2y-2z+Sqrt[2]>=0&&2x+2y+2z+Sqrt[2]>=0&&-2x+2y+2z+Sqrt[2]>=0&&2x-2y+2z+Sqrt[2]>=0&&-2x-2y+2z+Sqrt[2]>=0&&(-(1/Sqrt[2])-x)^2+y^2+z^2>1/4&&x^2+(1/Sqrt[2]-y)^2+z^2>1/4&&x^2+y^2+(-(1/Sqrt[2])-z)^2>1/4&&x^2+y^2+(1/Sqrt[2]-z)^2>1/4&&x^2+(-(1/Sqrt[2])-y)^2+z^2>1/4&&(1/Sqrt[2]-x)^2+y^2+z^2>1/4,{x,-1,1},{y,-1,1},{z,-1,1},PlotPoints->80,Mesh->None,PlotTheme->"Classic",BoundaryStyle->None,Axes->False,Boxed->False]
阅读全文
0 0
- 函数的表示跟人的自然的逻辑思维之间还是有差距
- 理想和现实还是有差距的
- 差距实在是太大了,我跟人的差别实在是太大了
- 记住企业与股票之间总是有差距的
- 国内外设计师之间的差距
- 豬跟人的區別
- 计算两个时间之间的差距
- Android_SDK各版本之间的差距
- 找一找自己和大神之间的差距
- 跟人聊天的总结:为什么选择SAP
- 自然二进制和格雷码之间的转换
- 现实和梦想的差距有多大
- 程序设计的逻辑思维
- 中国人的逻辑思维能力
- 逻辑思维的训练
- JS的逻辑思维练习
- 我的逻辑思维
- 传入时间与现在时间的差距及其表示
- 使用plupload压缩图片
- MySQL触发器使用详解
- python多线程下载器
- Java面向对象_4
- js一个页面嵌套另外一个页面
- 函数的表示跟人的自然的逻辑思维之间还是有差距
- PHP的内存限制 Allowed memory size of 134217728 bytes exhausted (tried to allocate 1099 bytes) in
- 利用Ajax发送请求(一)
- 发送Https请求(Java)
- vim使用
- Retrofit2-如何在请求时使用动态URL
- linux下eclipse运行c程序提示Launch failed.Binary not found
- String常见功能
- windows远程上传文件到ubuntu