2. MATLAB程序设计方法.ppt

* */32 * 数学科学学院：汪小平 wxiaoping325@163.com 第二章 MATLAB程序设计方法 表达式和运算符 命令文件(程序文件) 程序设计中的流程控制 函数文件的编写格式 ? ? ? ? 表达式——由变量、运算符、函数、数字组成 球坐标变换 ————表达式和运算符———— 用赋值语句创建变量 变量名=数据 变量名=表达式 x=R*cos(theta)*cos(fai); y=R*cos(theta)*sin(fai); z=R*sin(theta); ——纬度 ——经度 例2.1 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记 [X,Y,Z]=sphere(24); mesh(X,Y,Z),hold on fai=0; theta=51*pi/180; x0=cos(theta)*cos(fai); y0=cos(theta)*sin(fai); z0=sin(theta); plot3(x0,y0,z0,ro,linewidth,2) view(124,30) 例2.2 四边形柱面绘图原理和方法。 设四边形顶点为(-1, -1), (1, -1), (1, 1), (-1, 1). 设柱面高为1,创建Z坐标矩阵 创建柱面的X坐标和Y矩阵 cylinder(R，N)创建单位高度绕Z轴旋转曲面，R为母线,N表示圆圈上点数,N缺省默认值为 20。 Cylinder %单位高度二十边形柱面 [x,y,z]=cylinder([1,1],4); mesh(x,y,z) %单位高度四边形柱面 [x,y,z]=cylinder([1,0],3); mesh(x,y,z) %单位高度的三角形锥面 使用cylinder绘柱面 + 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂 1 .算术运算符 2.关系运算符 ? 小于 ? = 小于等于 ? 大于 ? = 大于等于 ? == 等于 ? ~= 不等于 3.逻辑运算符 与 (and) | 或 (or) ~ 非 (not) ————表达式和运算符———— 运算优先级 ①函数运算 ②算术运算 ③关系运算 ④逻辑运算 exp()、log()、sin()、abs()、fix()、··· .^、^、.*、./、*、/、+、–、； 、=、、=、==、~= 、| 表达式 (1+fix(pi))*mod(2,4)+2*3^2结果为ans=26； 三角形任意两边之和大于第三边的逻辑表达式 (A) a+b=c | a+c=b | b+c=a (B) a+b=c | a+c=b | b+c=a (C) a+bc | a+cb | b+ca (D) a+bc a+cb b+ca ————命令文件———— 计算机程序主要功能——接收数据和处理数据,并将处理后的数据完整有效地提供给用户 例2.3 将地球模型取为半径为R=6400(km)的球体,计算高度为1万公里的地球同步卫星对地面的覆盖率 在编辑窗口编写文件 planet.m 在命令窗口键入 planet 结果显示: ans = 30.4878 d=10000; R=6400; S0=4*pi*R*R; S=2*pi*R*R*d/(R+d); S/S0*100 命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息 M文件分为命令文件和函数文件两种 ①在编辑窗口中编写; ②保存并对文件命名； ③命令窗口键入文件名运行； ④观察运行结果； ······························· 条件控制; 循环控制; 错误控制;终止运行控制 条件控制——有选择地运行程序块 1. if / elseif 语句 例2 .8 判润年程序 ————程序设计中的流程控制———— year=input(input year:=); n1=year/4;