数学实验第一篇.pptVIP

一个M文件便保存在磁盘上了，便于修改、调用、运行和今后访问。 一个M文件可被别的M文件调用。这如同在电脑美术设计中我们绘制了一个冉冉升起的太阳，以后在别的电脑美术作品中，可随意调用、插入这个太阳。 例：x=[1 2 3]; y=[4 5]; [X Y]=meshgrid(x,y); % 平面点集的横坐标矩阵为X，纵坐标矩阵为Y 对每一个平面点(x,y)，即X，Y的对应分量定义函数值 z=f(x,y) 这些函数值其实是空间点的立坐标矩阵Z。 (4) 对空间点集{ ( X , Y , Z) }完成曲面绘图。 实现上述任务的命令为： 定义向量x,y; [X Y]=meshgrid(x,y); 由X,Y定义空间点的立坐标 Z, 即Z=f(X,Y) mesh(X,Y,Z); (或surf(X,Y,Z)) X= 1 2 3 1 2 3 Y= 4 4 4 5 5 5 Z=X.^2+Y.^2; mesh(X,Y,Z) meshgrid 平面点的横坐标构成X矩阵；平面点的纵坐标构成Y矩阵；空间点的是定义在平面点上 x向量为平面点的横坐标 y向量为平面点的纵坐标 图形示例： mesh 例1.5 作曲面z=f(x,y)的图形 键入： x=-7.5:0.5:7.5;y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2); Z=sin(R)./(R+eps) ; mesh(X,Y,Z); 课堂例题 x=0:pi/100:2*pi; y=x;[X Y]=meshgrid(x,y); Z=sin(X).*cos(Y); mesh(X,Y,Z) 等高线 MATLAB还提供了画二维和三维等高线图的函数contour和contour3. 例1.6 作出库函数peaks产生的二元函数的曲面及其等值线图。 键入： [X,Y,Z]=peaks(30); surf(X,Y,Z)； figure(2)； %打开另一个图形窗口 contour(X,Y,Z,16); figure(3)； contour3(X,Y,Z,16); 可将上述的画网格图的mesh命令改为surf(X，Y，Z)；则输出的曲面图，效果有所不同。 %产生网格细度为30的描述山峰特征的空间点的横坐标、纵坐标、立坐标矩阵 %16个等距的水平平面与山峰相交，交线投影在平面上 % 16个等距的水平平面与山峰相交，交线放置在原有空间位置 1.6 M文件编写 前面介绍M文件的概念和组成，看到了一些用MATLAB语言编写的简单程 序，要想实现更强的功能，需要用到循环控制、选择控制。几乎所有实用 的程序都包含循环控制、选择控制，熟练使用MATLAB的循环结构、选 择结构是编程的基本要求。MATLAB提供四种循环和选择控制结构，它们 是：for循环，while循环，if-else-end结构和switch-end结构。它们经常出现在M文件中。 1.6.1运算符 MATLAB的运算符可分为三类：算术运算符、关系运算符和逻辑运算 符。其中算术运算符的优先级最高，其次是关系运算符，再其次是逻辑运 算符。算术运算符 +、-、*、/、^、.*、./、.^ (后三个是矩阵的运算符) 前面已经介绍，这里再介绍关系运算符和逻辑运算符 关系运算符对于程序的流程控制非常有用.MATLAB共有六个关系运算符,它们分别是 小于；=小于等于；大于；=大于等于；==等于；~= 不等于。关系运算符用于比较同类型的算术表达式，运算结果为逻辑值：1(表示真即成立)，0(表示假即不成立)。 例如：3*52*3/4+1 x=9;y=4;z=5; 3*x+y-z/8=sin(z)/(1/x) 2*x==18  3*x~=5*z+2 3*y~=10关系运算符还可以比较同型矩阵，此时将生成一个0-1矩阵，当相应元素经关系运算为真时，对应位置上生成1，否则为0。 1.关系运算符 例如a=[1 2 3 ; 4 5 6 ; 6 4 2];b=[1 2.2 2.3 ; 4.1