|
z=c*sqrt（d-y.*y/b^2-x.*x/a^2）；u=1；%u=1，表示z要取正值
z1=real（z）；%取z的实部z1
fork=2：N-1；%以下7行程序的作用是取消z中含虚数的点
forj=2：N-1；%定义变量
ifimag（z（k，j））～=0z1（k，j）=0；
end
ifall（imag（z（[k-1：k+1]，[j-1：j+1]）））～=0z1（k，j）=NaN；
end
end
end
surf（x，y，z1），holdon%画空间曲面
ifu==1z2=-z1；surf（x，y，z2）；%u=1时加画负半面
axis（[-abs（a），abs（a），-abs（b），abs（b），-abs（c），abs（c）]）；%选择坐标轴的范围
end
xlabel（'x'），ylabel（'y'），zlabel（'z'）；%X，Y，Z轴的说明
holdoff%停止
运行结果（略）。
案例8、椭球面的绘制。
用MATLAB作图的程序为
ellipsoid（1，2，3，5，4，3，20）%椭球面绘制函数
运行结果（略）。
案例9、绘制双叶双曲面。
用MATLAB作图的程序为
a=5；b=4；c=10；
[x，y]=meshgrid（-20：1：20，-20：1：20）；
z=sqrt（c*（1+x.^2/a^2+y.^2/b^2））；
mesh（z）
holdon
mesh（-z）
运行结果（略）。
2.5二元函数的图象，二元函数的定义域案例10、作出函数z=xy的图形
用MATLAB作图的程序同案例2。
运行结果（略）。
2.6偏导数的几何意义案例11、函数z=2-2x2-y2在点（0.5，0.5）处的偏导数的几何意义。
用MATLAB作图的程序为
clear；x=-1：0.1：1；
y=-1：0.1：1；[X，Y]=meshgrid（x，y）；
Z=2-2*X.^2-Y.^2；mesh（X，Y，Z）；
运行结果（略）。
图1、2分别是函数z=2-2x2-y2在点（0.5，0.5）处对x、y的偏导数。
参考文献：
[1]于万波编著.基于MATLAB的计算机图形与动画技术[M].清华大学出版社，2007-6.
[2]赵静，但琦主编.数学建成模与数学试验[M].高教出版社，2009-1.
[3]薛定宇，陈阳泉著.高等应用数学问题的MATLAB求解[M].清华大学出版社，2008-3.
|
|
核心期刊网（www.hexinqk.com）秉承“诚以为基，信以为本”的宗旨，为广大学者老师提供投稿辅导、写作指导、核心期刊推荐等服务。 核心期刊网专业期刊发表机构，为学术研究工作者解决北大核心、CSSCI核心、统计源核心、EI核心等投稿辅导咨询与写作指导的问题。 投稿辅导咨询电话：18915033935 投稿辅导客服QQ： 投稿辅导投稿邮箱：1003158336@qq.com |
| 你好，欢迎来到！ 设为首页 ｜收藏本站 | ||
|
