Maple中曲面交线的绘制
给定两个曲面的方程:x^2 + y^3 = 1
z =x^2 + y^2
试绘制这两个曲面的交线.
要求如下:
一:
不能绘制出任何一个曲面,
二:
x=-2..2,y=-2..2,z=-2..2 2# liweicai990
看不见。 restart;
with(plots,intersectplot):
intersectplot(x^2+y^2=z, x^2+y^3-1=0, x=-2..2, y=-2..2, z=0..2, axes=box, thickness=3);
3# TBE_Legend
漂亮! 2# liweicai990
图片看不见. 该曲线用OpenFC（Forcal程序）绘制，代码不如maple简洁，但可以通过鼠标拖动变换任意角度来观察物体。
gleDrawScene,stop(); //设置场景绘制函数后退出
//让图形随鼠标移动（这段代码是固定的）。factor为比例因子，factor越大则移动幅度越大。
fcglFollowMouse(factor:x,y,xx,yy,xxx,yyy,static,fcglFollowMouse,xMouse,yMouse)=
{
if,
gleGetWindowSize[&xxx,&yyy],
if,gleAngle],
if,
glRotated, //使图象沿X轴旋转
glRotated, //使图象沿Y轴旋转
glRotated //使图象沿Z轴旋转
};
Scene(:x,y,z,t,k)=
{
glClear[], //清除屏幕以及深度缓存
glLoadIdentity[], //重置视图
glTranslated, //移动坐标，向屏幕里移动40个单元
fcglFollowMouse, //通过鼠标拖动变换任意角度来观察物体
glColor3d, //设置颜色为红色
fgMesh3d[-k,k,-k,k,0,2*k,-2,2,-2,2,-2,2], //绘制三维网格
glColor3d, //设置颜色为绿色
k=10, t=0.001, x=-2,
glBegin(GL_POINTS),
(x<=2).while{
y=(1-x*x)^(1/3), z=x*x+y*y,
if, //绘制点
x=x+t
},
glEnd()
};
OpenFC下载：http://xiazai.zol.com.cn/detail/27/262791.shtml 观察角度以及其他可观数目的可控项都可以在Maple的右键中完成,非常简单.所以初次绘图可以只是初具雏形即可. 6# wanglu
关于http://forum.simwe.com/thread-908387-1-1.html
中拟合问题,Maple可以给出一种基于视觉而不必基于最优化的调整最佳参数的方法.你可以通过自己调节k和b的值来自行判断拟合是否已经"合适"了.我随便弄一下,不一定是最好的,因为是基于个人观感的,但看上去还不错.
本帖最后由 liweicai990 于 2009-12-1 22:20 编辑
贴代码怎么很乱 6# wanglu
关于http://forum.simwe.com/thread-908387-1-1.html
中拟合问题,Maple可以给出一种基于视觉而不必基于最优化的调整最佳参数的方法.你可以通过自己调节k和b的值来自行判断拟合是否已经"合适"了.我随 ...
maplelab 发表于 2009-12-1 22:04 http://forum.simwe.com/images/common/back.gif
谢谢关注我的帖子。
这是一个好方法，特别在数值求解时很有用。
在另一个帖子里，也是用绘图的方法直观的得到了方程组的解的大致范围和解的个数。一并致谢！ 6# wanglu
关于http://forum.simwe.com/thread-908387-1-1.html
中拟合问题,Maple可以给出一种基于视觉而不必基于最优化的调整最佳参数的方法.你可以通过自己调节k和b的值来自行判断拟合是否已经"合适"了.我随 ...
maplelab 发表于 2009-12-1 22:04 http://forum.simwe.com/images/common/back.gif
能否帮忙用Maple绘出我的第二个问题的图：用我的优化参数带入微分方程组求解作为理论值，再绘制实验值进行比较。
祝贺Maple开坛，祝贺TBE_Legend, maplelab两位版主。 不错，学习了。。
页:
[1]