当前位置：文档之家 › matlab画二元函数

matlab画二元函数

  • 格式：pdf
  • 大小：84.89 KB
  • 文档页数：4

下载文档原格式

  / 4

合集下载

Matlab用法简介

Matlab用法简介

量对应元素的运算
x.*y=[x1*y1 x2*y2 x3*y3] (乘积) x./y=[x1/y1 x2/y2 x3/y3] (右除，右 边的y做分母) x.\y=[y1/x1 y2/x2 y3/x3] (左除，左 边的x做分母) x.^5=[x1^5 x2^5 x3^5] (乘幂) 2.^x=[2^x1 2^x2 2^x3] x.^y=[x1^y1 x2^y2 x3^y3]
左除“ \ ”：
求矩阵方程AX=B的解；（ A 、B的行要保持一致） 解为 X=A\B； 当A为方阵且可逆时有X=A\B=inv(A)*B；
右除“ / ”：
求矩阵方程XA=B的解 （A 、B的列要保持一致） 解为 X=B/A ， 当A为方阵且可逆时有X=B/A=B*inv(A)
矩阵的运算(左除和右除)
1. meshgrid指令：生成网格点
观察meshgrid指令的效果。 程序：
a=-3;b=3;c=-3;d=3;n=10; x=linspace(a,b,n); y=linspace(c,d,n); [X,Y]=meshgrid(x,y); plot(X,Y,'+')
观察结果
2. 计算函数值
3. 绘图指令

两个坐标系，用于绘制不同尺度的函数。
绘制平面曲线（线型、点形和颜色的控制）
plot(x,y,‘颜色＋线型＋点形’) plot(x,y,‘颜色＋线型＋点形’,x,y,‘颜 色＋线型＋点形’,… ) 句柄图形和 set 命令改变属性值，可套 用：


h=plot(x,y), set(h,‘属性’,属性值,‘属性’,属性值,…) 或plot(x,y,'属性','属性值')设置图形对象 的属性。

MATLAB图形工具箱

MATLAB图形工具箱

一、二维曲线的绘制
对于图形坐标的控制，请看图片自 己在上机实习的时候学会！
一、二维曲线的绘制
绘制函数的图形: 语法: fplot(‘fun’,lim,’s’):绘制函数fun的图形. 说明:fun为要绘制图形的函数,lim为图形的 范围,如:[0,2*pi],s为图形所用的点型,线 型和颜色的设定,与plot命令相同.
二、三维图形的绘制
2‘ 三维曲面图 语法：surf(z) surf(x,y,z) 说明：意义与mesh命令相同，只是图形不 是网线而是曲面。
三、特殊图形的绘制
1、条形图 语法：bar(x,y,width) %画条形图 bar3(x,y,width,’参数’) %画三维条形图 说明：x是横坐标向量，省略时默认值是1:m，m为y的向 量长度；y是纵坐标，当y是向量时，每个元素对应一 个竖条，当y是矩阵时，将画出m组竖条，每组包含n 条；width是竖条的宽度。默认时为0.8，（把bar命令 改为barh命令时将绘制出水平的条形图）。 三维的条形图中的参数可以detached, grouped,stacked三个,分别表示不同的显示方式.
一、二维曲线的绘制
数据点形 实点标记：. 圆圈标记：。 等等，其他更多的可以查相关的资料。 例如:>>x=0:0.1:2*pi; >>plot(x,sin(x),’r-.’) :用红色点划线 画出曲线 >>plot(x,cos(x),’b:*’)
一、二维曲线的绘制
6、交互式图形命令 （1）ginput命令 ginput命令与其他图形命令的原理不同， 不是把数据表现在图形上，而是从图上 获取数据。因此，ginput命令在数值优化、 工程设计中十分有用，仅适用于二维图 形。

matlab基础测习题(带答案)

matlab基础测习题(带答案)

精心整理Matlab基础练习题常量、变量、表达式1、MATLAB中，下面哪些变量名是合法的？（）2、3、4、5、）6、小数7、8、在MATLAB中，a=1，b=i，则a占_8__个字节，b占_16_个字节，c占________字节。

9、在MATLAB中，inf的含义是__无穷大__，nan的含义是__非数（结果不定）___。

数组1、在MATLAB中，X是一个一维数值数组，现在要把数组X中的所有元素按原来次序的逆序排列输出，应该使用下面的（）指令。

（A）X[end:1] （B）X[end:-1:1] （C）X(end:-1:1)（D）X(end:1) 2、在MATLAB中，A是一个字二维数组，要获取A的行数和列数，应该使用的MATLAB的命令是（）。

（A）class(A) （B）sizeof(A) （C）size(A)（D）isa(A)3、4、这5、6、___linspace(0,5,50)___7、在MATLAB中，A=[0:1/2:2]*pi，那么sin(A)=___[010-10]_____。

8、在MATLAB中，A=[1,2,3;4,5,6;7,8,0]，B=[2,1,6;8,5,2;14,2,1]。

写出下面MATLAB语句执行的结果：(为节省篇幅，把矩阵写成mat2str的形式)（1）A==B [000;010;000]（2）A.*B [2218;322512;98160]（3）A(:)' [147258360]（4）A(1,:)*B(:,3) 139、在MATLAB中，写出下面MATLAB语句执行的结果：（1）c lear,A=ones(2,6)A=（5）B=A（:,end:-1:1)B=211713951231915117310、请编写一段matlab程序，完成以下功能：（1） 生成一个100行，200列的二维随机数组；A=rand(100,200);（2） 找出数组A 中所有大于0.49且小于0.51的元素的单下标；Idx=find(A(:)>0.49&A(:)<0.51)（3） 数组A 中满足（2）中的条件的元素有多少个？11、12、 22cos 0.7sin (1)x y x x ⎡⎤=+⎢⎥+⎣⎦，编写一段matlab 程序，要求如下（1） 在[0,2]π区间，每隔0.01取一x 数值，计算出相应的y 的函数值；x=0:0.01:2*pi;y=(0.7+2*cos(x)./(1+x.^2)).*sin(x)（2）根据MATLAB计算出的数据，找出在[0,2]π内该函数的极小值的坐标。

matlab梯度法求二元函数最小值

matlab梯度法求二元函数最小值

梯度法是一种常见的数值优化方法，用于求解无约束优化问题。

以下是使用Matlab 梯度法求二元函数最小值的步骤：定义目标函数：首先需要定义要求解的目标函数。

比如对于一个二元函数f(x,y)，可以在Matlab 中定义为一个function handle 的形式：f = @(x) x(1)^2 + x(2)^2;计算梯度：使用gradient 函数计算目标函数的梯度，得到每个变量的偏导数。

例如，在Matlab 中可以这样计算f(x,y) 的梯度：grad = @(x) [2*x(1), 2*x(2)];初始化参数：确定起始点，即需要进行优化的初始值，例如：x0 = [1; 2];设置迭代次数和步长：根据实际情况，设置迭代次数和步长，例如：max_iter = 1000;alpha = 0.01;进行迭代：使用for 循环进行多次迭代，每次沿着梯度方向更新参数值，直至达到最小值或者迭代次数满足停止条件。

具体来说，可以实现如下代码：for i = 1:max_iterx1 = x0 - alpha*grad(x0);if abs(f(x1)-f(x0)) < 1e-6break;endx0 = x1;end其中，1e-6 表示允许的最小误差值。

输出结果：经过多次迭代之后，x0 中存储的即为所求的最小值点。

可以使用fprintf 函数将其输出：fprintf('The minimum point is (%f, %f)\n', x0(1), x0(2));这样就完成了在Matlab 中使用梯度法求解二元函数最小值的过程。

需要注意的是，在实际应用中，需要根据具体问题和实验数据调整步长和迭代次数等参数，以获得更好的优化效果。

1. Matlab使用入门

1. Matlab使用入门

1 2 3 4 5 6 7 8
1 2 3 4 5 6 7 8
24/32
绘网面命令mesh()使用格式: mesh(x,y,z) 或 mesh(z) 例1.12 绘二元函数 z = x exp( –x2 – y2)的图形。 [x,y]=meshgrid(-2:0.2:2); z=x.*exp(-x.^2-y.^2); mesh(x,y,z) colormap([0 0 1]) 注记: x,y是维数相同的矩阵; 表达式中“.*”和“.^”运算 使得z是与x,y维数相同矩阵。 功能相同的绘图命令:surf, meshc, meshz,· · · · · · · ·
11/32
三角函数与双曲函数
sin asin cos acos tan atan cot acot sec asec csc acsc 正弦函数 反正弦函数 余弦函数 反余弦函数 正切函数 反正切函数 余切函数 反余切函数 正割函数 反正割函数 余割函数 反余割函数 sinh asinh cosh acosh tanh atanh sech asech csch acsch coth acoth 双曲正弦函数 反双曲正弦函数 双曲余弦函数 反双曲余弦函数 双曲正切函数 反双曲正切函数 双曲正割函数 反双曲正割函数 双曲余割函数 反双曲余割函数 双曲余切函数 反双曲余切函数
1/3 1/4 1/5
19/32
————矩阵创建常用方法————
1.直接输入法; 2.特殊矩阵函数法; 3. 数据文件输入 注意事项 矩阵元素必须在方括号[ ]之内; 同一行相邻元素间用逗号或空格分隔; 矩阵的行与行之间用分号分隔. A=hilb(3) B=invhilb(3) A*B
掌握MATLAB程序设计方法
以MATLAB为操作平台完成实验作业

matlab二元函数已知函数值求自变量

matlab二元函数已知函数值求自变量

matlab二元函数已知函数值求自变量MATLAB是一种常用的数学软件，可以用于解决各种数学问题，包括求解二元函数。

在本文中，我们将探讨如何使用MATLAB求解二元函数的自变量。

二元函数是指具有两个自变量的函数，可以表示为f(x, y)。

在MATLAB中，我们可以通过给定函数值来求解对应的自变量。

具体而言，我们可以通过给定f(x, y)的函数值来求解x和y的取值。

我们需要明确给定的函数值。

假设我们已知二元函数f(x, y)在某些特定点的函数值，我们的目标是求解x和y的取值。

为了使用MATLAB求解这个问题，我们可以使用插值方法。

插值是一种通过已知数据点来推断未知数据点的方法。

在MATLAB中，有多种插值方法可供选择，如线性插值、多项式插值和样条插值等。

下面我们将以一个具体的例子来说明如何使用MATLAB求解二元函数的自变量。

假设我们有一个二元函数f(x, y)，我们已知它在以下四个点的函数值：f(1, 1) = 2f(1, 2) = 4f(2, 1) = 3f(2, 2) = 6我们的目标是求解x和y的取值。

我们可以使用线性插值方法来推断其他点的函数值。

线性插值是一种简单但有效的插值方法，它假设函数在已知点之间是线性的。

在MATLAB中，我们可以使用interp2函数来进行线性插值。

interp2函数的用法如下：z = interp2(x, y, v, xi, yi)其中，x和y是已知点的横纵坐标，v是已知点的函数值，xi和yi 是待求函数值的横纵坐标，z是待求函数值。

根据我们的例子，我们可以将已知点的横纵坐标和函数值输入interp2函数来求解其他点的函数值。

具体代码如下：x = [1 1 2 2];y = [1 2 1 2];v = [2 4 3 6];xi = [1.5];yi = [1.5];z = interp2(x, y, v, xi, yi)运行上述代码，我们可以得到z的值为4.5。

MATLAB课后习题及练习7.6

7.6 课后作业1．分别绘制下列函数图形：(1) r=3(1-cosθ) (极坐标)>> t=0:.01:2*pi;>> figure>> polar(t,abs(3*(1-cos(t))));运行结果：(2) y (t)=1.25e-0.25t +cos(3t)>> x=0:pi/180:2*pi;>> y1=1.25*exp(-0.25*x);y2=cos(3*x);>> plot(x,y1+y2)运行结果：2 ．绘制函数y (t)=1-2e-tsint (0 ≤t≤8)的图形，且在x 轴上标注“Time”，y 轴上标注“Amplitude”，图形的标题为“Decaying Oscillating Exponential”。

>> x=0:0.01:8;>> y=1-2*exp(-1*x).*sin(x);>> figure>> plot(x,y)>> xlabel('Time');ylabel('Amplitude');>> title('Decaying Oscillating Exponential');运行结果：3 ．在同一图中绘制下列两条曲线(x ∈[0,25]内) ：(1) y 1(t)=2.6e(-0.5x) cos(0.6x)+0.8 ；(2) y 2(t)=1.6cos(3x)+sin(x)要求用不同的颜色和线型分别表示y 1(t) 和y 2(t) ，并给图形加注解。

>> x=0:0.1:25;>> y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.8;>> y2=1.6*cos(3*x)+sin(x);>> figure>> plot(x,y1,'r-X',x,y2,'b--.');>> legend('y1=2.6*exp(-0.5*x).*cos(0.6*x)+0.8','y2=1.6*cos(3*x)+sin(x)')4 ．在一个图形窗口下绘制两个子图，分别显示下列曲线：(1)y =sin2xcos3x ；(2) y =0.4 x要求给x 轴、y 轴加标注，每个子图加标题。

matlab中mesh命令

matlab中mesh命令mesh命令是MATLAB中用于绘制三维曲面的函数，它可以将离散的数据点连接成平滑的曲面。

本文将介绍mesh命令的使用方法和一些相关知识。

我们来了解一下mesh命令的基本语法。

在MATLAB命令窗口中输入"help mesh"即可查看详细的使用说明。

一般情况下，mesh命令的语法如下：mesh(X, Y, Z)其中，X、Y和Z是三个相同大小的矩阵，分别表示曲面上各个点的横坐标、纵坐标和高度。

通过这些点的连接，mesh命令可以绘制出一个三维曲面。

下面我们来看一个具体的例子。

假设我们希望绘制一个二元函数z = f(x, y) = sin(x) + cos(y)，我们可以先生成一个横坐标和纵坐标的网格矩阵，然后计算对应的高度矩阵。

具体的代码如下：```matlabx = linspace(-pi, pi, 100);y = linspace(-pi, pi, 100);[X, Y] = meshgrid(x, y);Z = sin(X) + cos(Y);mesh(X, Y, Z);运行上述代码后，MATLAB会自动打开一个新的窗口，显示出绘制的三维曲面。

通过鼠标的拖动和缩放，我们可以自由地查看曲面的各个部分。

除了基本的绘制功能，mesh命令还提供了一些可选参数，用于进一步定制曲面的外观。

例如，我们可以使用以下代码将曲面的颜色设置为红色，并添加网格线：```matlabmesh(X, Y, Z, 'FaceColor', 'red', 'EdgeColor', 'black');```mesh命令还可以与其他绘图命令结合使用，实现更加复杂的图形效果。

例如，我们可以使用以下代码将曲面绘制在一个已有的三维坐标系上：```matlabfigure;hold on;plot3(X(:), Y(:), Z(:), 'o', 'MarkerFaceColor', 'blue', 'MarkerEdgeColor', 'blue');mesh(X, Y, Z);hold off;运行上述代码后，我们会得到一个包含散点和曲面的图形，其中散点表示曲面上的各个数据点。

matlab 求解二元方程

matlab 求解二元方程
Matlab是一种强大的计算软件，可以用来解决各种数学问题，包括求解二元方程。

二元方程是指只有两个未知数的方程，通常表示为ax + by = c和dx + ey = f。

如果我们知道a、b、c、d、e和f 的值，我们可以使用Matlab求解这个方程组，找到x和y的值。

要使用Matlab求解二元方程，需要使用线性代数中的求解方法。

具体步骤如下：
1. 将方程表示为矩阵形式。

将系数和常数项分别表示为矩阵A 和矩阵B，得到方程组AX = B。

2. 使用Matlab中的“linsolve”函数求解方程组。

该函数接受两个参数：A和B。

它返回一个向量X，其中包含方程组的解。

3. 使用向量X计算x和y的值。

根据方程ax + by = c和dx + ey = f，我们可以得到x和y的计算公式。

使用Matlab求解二元方程非常简单，只需几行代码即可完成。

以下是一个示例：
% 定义方程的系数和常数项
A = [2 3; 4 5];
B = [8; 14];
% 求解方程组
X = linsolve(A, B);
% 计算x和y的值
x = X(1);
y = X(2);
disp(['x = ', num2str(x)]);
disp(['y = ', num2str(y)]);
这将输出x = 1，y = 2，这就是方程的解。

MATLAB绘图与图形处理

MATLAB绘图与图形处理人们很难从一大堆原始的数据中发现它们的含义，而数据图形恰能使视觉感官直接感受到数据的许多内在本质，发现数据的内在联系。

MATLAB可以表达出数据的二维，三维，甚至四维的图形。

通过图形的线型，立面，色彩，光线，视角等属性的控制，可把数据的内在特征表现得淋漓尽致。

下面我们分别介绍图形的命令。

7.1 二维图形7.1.1 基本平面图形命令命令1 plot功能线性二维图。

在线条多于一条时，若用户没有指定使用颜色，则plot循环使用由当前坐标轴颜色顺序属性（current axes ColorOrder property）定义的颜色，以区别不同的线条。

在用完上述属性值后，plot又循环使用由坐标轴线型顺序属性（axes LineStyleOrder property）定义的线型，以区别不同的线条。

用法plot(X,Y) 当X,Y均为实数向量，且为同维向量（可以不是同型向量），X=[x(i)]，Y=[y(i)]，则plot(X,Y)先描出点(x(i)，y(i))，然后用直线依次相连；若X，Y为复数向量，则不考虑虚数部分。

若X，Y均为同维同型实数矩阵，X = [X(i)]，Y = [Y(i)]，其中X(i),Y(i)为列向量，则plot(X,Y)依次画出plot(X(i),Y(i))，矩阵有几列就有几条线；若X，Y中一个为向量，另一个为矩阵，且向量的维数等于矩阵的行数或者列数，则矩阵按向量的方向分解成几个向量，再与向量配对分别画出，矩阵可分解成几个向量就有几条线；在上述的几种使用形式中，若有复数出现，则复数的虚数部分将不被考虑。

plot(Y) 若Y为实数向量，Y的维数为m，则plot(Y)等价于plot(X,Y)，其中x=1：m；若y 为实数矩阵，则把y按列的方向分解成几个列向量，而y 的行数为n，则plot(Y)等价于plot(X,Y)其中x=[1;2;…;n]；在上述的几种使用形式中，若有复数出现，则复数的虚数部分将不被考虑。

  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

y = x sin x y = 2 x − 4 和 在 x ∈ [0,6π ] 上 z = x cos x z = 4 x + 2
的图象． 解 打开Ｍ文件编辑窗口，在其中输入下面命令集： x=0:pi/50:6*pi; y1=x.*sin(x); z1=x.*cos(x); y2=2*x-1; z2=4*x-8; plot3(x,y1,z1,x,y2,z2); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') legend('圆锥螺线','空间直线') 取名为 exa1 保存，再在命令窗口中输入命令 exa1，程序运行 结果如图 4-4． 说明：若要作单根曲线，和前面 plot 函数一样，只是多了一个立坐标． 二、空间曲面 1．三维网格图 调用格式一： [x,y]=meshgrid(x1:h1:x2,y1:h2:y2) z=f(x,y) plot3(x,y,z) 其中 x、y 和 z 为长度相同的向量，分别用于存储 x 坐标、y 坐标和 z 坐标数据． 例 2 作出空间平面 4 x − 2 y + z = 3 的图象． 解 打开Ｍ文件编辑窗口，在其中输入下面命令集： [x,y]=meshgrid(-3:0.1:3,-2:0.1:3); %产生一个 x,y 平面上 51× 61 的网格 z=3-4*x+2*y; plot3(x,y,z); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') title('4x-2y+z=3') 取名为 exa2 保存，再在命令窗口中输入命令 exa2，程序运行 结果如图 4-5． 例 3 作出空间曲面 z = 图 4-5
图 4-4
sin x 2 + y 2 x2 + y2
在区域
{( x, y ) || x |≤ 8, | y |≤ 8} 上的图象．
解 打开Ｍ文件编辑窗口，在其中输入下面命令集： %产生对应的 x,y 坐标用点幂 %产生对应于 x,y 的 z 坐标用点除 [x,y]=meshgrid(-8:0.1:8); r=sqrt(x.^2+y.^2); z=sin(r)./r; plot3(x,y,z);
图 4-8
xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') shading interp %插值处理 title('z=x^2-y^2 的插值处理图') 取名为 exa5 保存，再在命令窗口中输入命令 exa5，程序运行结果如图 4-8． 3．柱面图 调用格式一： [x,y,z]=cylinder(R,N) mesh(x,y,z) 其中， R 是一个向量， 其坐标分量依次对应柱面各横截面的半径， Ｎ是多边形的边数． 默 认值是Ｒ＝[1,1]，Ｎ＝20． 例 6 作出正六棱柱在Ｒ＝[２,２]的图象． 解 打开Ｍ文件编辑窗口，在其中输入下面命令集：
2
− y2
的图象．
[x,y]=meshgrid(-4:0.1:3,-3:.1:2); %产生一个 x,y 平面上 51×71 的网格 z=exp(-x.^2-y.^2); mesh(x,y,z); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') 取名为 exa2 保存，再在命令窗口中输入命令 exa2，程序运行 结果如图 4-7． 2．三维曲面图 调用格式： [x,y]=meshgrid(x1:h1:x2,y1:h2:y2) z=f(x,y) surf(x,y,z)
[x,y,z]=cylinder([2,2],6); mesh(x,y,z); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') hidden off %显示隐含线，若不想显示隐含线，则取消 图 4-9 此命令或设置为 hidden on 取名为 exa6 保存，再在命令窗口中输入命令 exa6，程序运行结果如图 4-9． 调用格式二： [x,y,z]=cylinder(R,N) surf(x,y,z) 其中，R、Ｎ的意义和上面一样． 例 7 作出正十棱台在Ｒ＝[1,2]的图象． 解 打开Ｍ文件编辑窗口，在其中输入下面命令集：
z = x 2 − y 2 在区域 {( x, y ) || x |≤ 1, | y |≤ 2}
图和经插值处理的图． 解 打开Ｍ文件编辑窗口，在其中输入下
面命令集： [x,y]=meshgrid(-1:0.1:1,-2:.1:2); z=x.^2-y.^2; subplot(1,2,1), surf(x,y,z); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') title('z=x^2-y^2'); subplot(1,2,2), surf(x,y,z);
数学实验四 用 Matlab 软件作二元函数的图象
一、空间曲线 调用格式： plot3(x1,y1,z1,x2,y2,z2,…,xn,yn,zn) 其中 xi、yi 和 zi 为分别是长度相同的向量，用于存储 xi 坐标、yi 坐标和 zi 坐标数据， i=1,2,…,n． 例 1 在同一坐标系中分别作出三维曲线
图 4-10
其中，N（N>2）是球面的边数．默认值是Ｎ＝20．sphere(N)不返回坐标，直接绘出球 面． 例 8 以子图形式分别作出 N=5、20、40 和 40 并作平滑处理的球． 解 打开Ｍ文件
sphere(5); title('子图 1') subplot(2,2,2),sphere; title('子图 2') subplot(2,2,3) [x,y,z]=sphere(40); surf(x,y,z); title('子图 3') subplot(2,2,4) [x,y,z]=sphere(40); surf(x,y,z); shading flat %对球面作平滑处理 title('子图 4') 图 4-11 取名为 exa8 保存，再在命令窗口中输入命令 exa8，程序运行结果如图 4-11． 三、上机实验 1．用 help 命令查看函数 plot3，mesh 和 surf 等的用法． 2．上机验证上面各例． 3．作相关小节练习中空间曲线和曲面的图象．
图 4-7
说明： surf 函数可借助 shading 函数作平滑和插值处理，而 shading 函数有三个参数， 分别为 flat（作平滑处理） 、interp（去掉连接线条，在各片之间使用颜色插值）和 faceted （默认值，对前面两种参数之一的作用进行还 原） ． 例 5 用子图分别作出马鞍面
xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') 取名为 exa3 保存，再在命令窗口中输入命令 exa3，程序运行结果如图 4-6． 调用格式二： [x,y]=meshgrid(x1:h1:x2,y1:h2:y2) z=f(x,y) mesh(x,y,z) 例 4 在 x、y 平面内选取一个区域，作出空间曲面 z = e − x 解 打开Ｍ文件编辑窗口，在其中输入下面命令集：
[x,y,z]=cylinder([1,2],10); surf(x,y,z); xlabel('x 轴'),ylabel('y 轴'),zlabel('z 轴') 取名为 exa7 保存，再在命令窗口中输入命令 exa7，程序运行 结果如图 4-10． 说明：若要作正棱锥，只须取Ｒ的某个分量为 0 即可． 4．球面图 调用格式： [x,y,z]=sphere(N) surf(x,y,z)