第2章MATLAB程与作图.pptVIP

第二章 MATLAB编程与作图 MATLAB数学实验 第二章 MATLAB编程与作图 第二章 MATLAB编程与作图 2.1 程序设计 2.2 作图 2.3 在线帮助和文件管理 2.1 程序设计 循环语句 for 循环变量=初值:增量:终值, 语句；end while（条件式）, 语句；end 分支语句 if（条件式）, 语句; end if（条件式1), 语句1; elseif （条件式2), 语句 2; ……; else, 语句;end switch(分支变量)case(值1), 语句1;case(值2), 语句2;……;otherwise 语句;end 其它：pause, break, return, error 2.1 程序设计 循环语句 for 循环变量=初值:增量:终值, 语句；end while（条件式）, 语句；end 分支语句 2.1 程序设计 ? s=0; for n=1:100 , s=s+1/n/n;end；s ? clear;s=0;n=1; while n=100, s=s+1/n/n; n=n+1; end;s 强行中断：Ctrl+C 2.1 程序设计 M脚本文件 eg2_1 在Editor窗口 文件名一律以字母开头，以字母、数字或下划线组成，不要含有空格、减号等. M文件名一般都用小写字母 保存在当前目录(Current Directory) 2.1 程序设计 M函数文件 function 输出变量=函数名(输入变量) 语句; eg2_1f M函数必须给予输入参数值 M函数在edit窗口编写，在Command窗口调用 M函数是以该函数的磁盘文件主名调用，而不是文件中的函数名称 M函数中变量为局部变量 2.1 程序设计 函数句柄(handle) fname=@eg2_1f; feval Inline函数 fun=inline(‘expr’, arg1, arg2, ...) 匿名函数 fun=@(arg1,arg2, …) 2.1 程序设计 注释：%开头，对本行后面字符起作用，不参与运算。 ?对话：input，disp ?全程变量与局部变量 ? nargin、nargout和varargin ?子函数和嵌套函数 提高速度 2.1 程序设计 普通编程 function s=f(m) s=0; for n=1:m s=s+1/n/n; end 向量化编程 function s=f(m) n=1:m; s=sum(1./n.^2); 2.1 程序设计 ?例2.4 编一M函数,对任意输入的向量x, 可计算分段函数值构成的向量。 分量方式 eg2_4a，慢 向量方式 eg2_4b, eg2_4c, 快 数组预分配y=zeros(size(x)) 2.2 作图 曲线图 plot(x,y) 以数据(x(i), y(i))为节点的折线图, 其中x, y为同长度的向量 plot(x1,y1,x2,y2,...) 多组数据折线图 fplot(fun,[a,b])函数fun在区间[a,b]上的函数图 plot3(x,y,z) 空间曲线图, 其中x, y, z为同长度的向量 图形导出到word（菜单：File-Save as, Edit-Copy figure) 线型与标记 P31 表 eg2_5 曲线图y=x3-x-1和y=|x|0.2sin(5x) 2.2 作图 曲面图 [x,y]=meshgrid(xa,ya)当xa, ya分别为m维和n维行向量，得到x和y均为n行m列矩阵。meshgrid常用于生成X-Y平面上的网格数据。 mesh(x,y,z) 绘制网面图，是最基本的曲面图形命令, 其中x, y, z是同阶矩阵，表示曲面三维数据。 surf(x,y,z) 绘制曲面图，与mesh用法类似。 eg2_6 二元函数图 z = x·exp(-x2-y2 ) Meshgrid解释 xa=6:8;ya=1:4; [x,y]=meshgrid(xa,ya) %生成X-Y面上网格 z=x.^2+y.^2 %计算X-Y面上各网格点的z轴高度 x y z 6 7 8 1 1 1 37 50 65 6 7 8 2 2 2 40 53 68 6 7 8 3 3 3 45 58 73 6 7 8 4 4 4 52 65 80