数学建模-MATLAB用法.ppt

数学建模培训—— MATLAB软件 刘盼萍 中北大学数学系ping8191@ 用计算机解决数值分析问题 1.1 Matlab 的特点与功能 1.2 Matlab 的集成环境 Matlab 的工作界面 M 文件的建立、打开、运行 3.2 函数文件 函数文件如下： function [s,p]=fcircle(r) %CIRCLE calculate the area and perimeter of a circle of radii r %r 圆半径 %s 圆面积 %p 圆周长 ? s=pi*r*r; p=2*pi*r; switch 语句 switch 语句 例 已知y的表达式 ，当n=100时，求y的值。 图形的基本属性 ezplot 作图 法2、二维曲线的绘制： ezplot ezplot(f(x,y),[a,b,c,d]) 绘f(x,y)=0在区间 axb,cyd 上的图形 ezplot(f,g,[a,b]) 绘x=f(t),y=g(t)在区间 atb 上的图形 ezplot(f(x,y),[a,b]) 绘f(x,y)=0在区间 ax,yb 上的图形 缺省的绘图区间为 [-2?, 2?] 缺省的绘图区间为 [0, 2?] ezplot(f,[a,b]) 绘 y=f(x)在区间 axb 上的图形 缺省的绘图区间为 [-2?, 2?] 例： 4.2 三维曲线 ： plot3 绘制一条三维曲线：x=x(t),y=y(t),z=z(t),则： plot3的用法与 plot 类似 plot3(x,y,z,’选项’) plot3(x1,y1,z1,’选项1’，x2,y2,z2,’选项2’ …,xn,yn,zn,’选项n’) 绘制多条三维曲线： 三维曲线的简捷绘制： ezplot3 缺省的绘图区间为 [0, 2?] ezplot3(x,y,z,[a,b]) 绘空间曲线x=x(t),y=y(t),z=z(t)在 atb 上的图形 4.3 三维曲面 绘制由函数 z=z(x,y) 确定的曲面时，步骤如下： setp1: x=a:dx:b;y=c:dy:d; [X,Y]=meshgrid(x,y); setp2: 再计算函数z=z(x,y)在各网格点上的值 Z=…… ;注意此处采用点运算 setp3: 其中，X, Y 是网格划分后得到的网格矩阵 若 x = y, 则可简写为 [X,Y]= meshgrid(x) mesh(X,Y,Z,C) 或 surf(X,Y,Z,C) 或 其它绘图函数 绘制由矩阵 X,Y,Z 所确定的曲面网格图(或着色的表面图),矩阵 C 用于确定网格颜色，省略时 C=Z。 例：“墨西哥帽子” ( –a x a, -a y a ) a=8 时的曲面图形 x=-8:0.5:8; y=-8:0.5:8; [X,Y]=meshgrid(x,y); r=sqrt(X.^2+Y.^2)+eps; Z=sin(r)./r; subplot(1,2,1), mesh(X,Y,Z) subplot(1,2,2), surf(X,Y,Z) ＋eps？ mesh 与 surf 的比较 三维曲面的简捷绘制： ezmesh、ezsurf ezmesh(z,[a,b,c,d]) ezmesh(x,y,z,[a,b,c,d]) ezsurf 的用法 与 ezmesh 相同 标题标注 title(‘text’, ‘属性1’, 属性1的值, ‘ 属性2’, 属性2的值, ...) 属性: linewidth, markersize, fontsize,  fontweight, fontname, … 4.4 辅助操作（适合二维和三维图形） 坐标轴标注 xlabel 、ylabel、zlabel 添加图例 legend 在指定位置（x,y)处添加文本 text(x,y,string1,string2, ...) 显示网格 grid on/off 画/不画网格线 grid 在两种状态之间切换 axis([xmin, xmax, ymin, ymax, zmin, zmax]) 作用：控制坐标轴的显示范围 坐标轴控制命令： axis 划分绘图区域 subplot(m,n,p) 将一个绘图窗口分割成 m*n 个子区域，并 按行 从左至右 依次编号 。并指定第