第二章程序设计1.pptVIP

*/28 MATLAB 程序设计方法 I 表达式和运算符 命令文件(程序文件) 程序设计中的流程控制 练习与思考题 ? ? ? ? 变量命名规则 变量名的第一个字符必须是一个英文字母，最多 可包括31个字符； x1, x2, x3,alpha, · · · · · · · · · 变量名可由英文字母、数字和下划线混合组成； 变量名不得包含空格和标点,但可以有下连字符 变量名和函数名对字母的大小写敏感； X，x, Z, z 是不同变量 MATLAB的变量和表达式 局部变量：具有函数内部作用范围 全局变量：可在几个函数和基本工作空间 之间共享，具有全局作用范围 全局变量的声明： global 变量名1 变量名2 … 局部变量和全局变量 永久变量 ans; eps; pi; i; j; inf; NaN; 表达式——由变量、运算符、函数、数字组成 球坐标变换 ————表达式和运算符———— 用赋值语句创建变量 变量名=数据 变量名=表达式 x=R*cos(theta)*cos(fai); y=R*cos(theta)*sin(fai); z=R*sin(theta); ——纬度 ——经度 例2.2 格林威治天文台建于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.3 四边形柱面绘图原理和方法。 设四边形顶点为(-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.1 将地球模型取为半径为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命令的有序集合。