第三周 2014-3-10 曲面 编程.ppt

* 例1.4 作螺旋线x=sint，y=cost，z=t. 键入： t=0:pi/50:10*pi; plot3(sin(t),cos(t), sin(t).*cos(t)); 对于 t 的501个取值，绘制并连接501个空间点，形成空间曲线 1.5.2 三维图形 1.空间曲线 2. 空间曲面绘图 MATLAB要求我们对平面网格的每一个点{(x,y)}，定义立坐标z。 然后在mesh命令下联接各个空间点(x,y,z)，形成曲面。 例：x=[1 2 3]; y=[4 5]; [X Y]=meshgrid(x,y); % 平面点集的横坐标矩阵为X，纵坐标矩阵为Y X= 1 2 3 1 2 3 Y= 4 4 4 5 5 5 实现上述任务的命令为： 定义向量x,y; [X Y]=meshgrid(x,y); 由X,Y定义空间点的立坐标 Z, 即Z=f(X,Y) mesh(X,Y,Z); (或surf(X,Y,Z)) x y Z=X.^2+Y.^2; Z= 17 20 25 26 29 34 mesh(X,Y,Z) 例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文件编写 前面看到了一些用MATLAB语言编写的简单程序，要想实现更强的功能，需要用到循环控制、选择控制。几乎所有实用的程序都包含循环控制、选择控制，熟练使用MATLAB的循环结构、选择结构是编程的基本要求。MATLAB提供四种循环和选择控制结构，它们是：for循环，while循环，if-else-end结构和switch-end结构。 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