MATLAB--第二章.ppt

*/37 表达式和运算符 命令文件(程序文件) 程序设计中的流程控制 函数文件的编写格式 ? ? ? ? 《数学实验》第二章 表达式——由变量、运算符、函数、数字组成 球坐标变换 ————表达式和运算符———— 用赋值语句创建变量 变量名=数据 变量名=表达式 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); th=linspace(-pi/2,pi/2,100); x=cos(th);y=zeros(size(x));z=sin(th); plot3(x,y,z,b,x0,y0,z0,ro,linewidth,2) view(124,30) + 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂 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 R=6400; S0=4*pi*R*R; d=10000; S=2*pi*R*R*d/(R+d); S/S0*100 命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息 M文件分为命令文件和函数文件两种 ①在编辑窗口中编写; ②保存并对文件命名； ③命令窗口键入文件名运行； ④观察运行结果； ······························· 例2.3直线平行于Z轴沿x-y平面上的四边形移动，形成四边形柱面。利用矩阵方法绘制四边形柱面. 设四边形顶点为(-1, -1), (1, -1), (1, 1), (-1, 1).由于四边形是封闭的图形，将第五个点设为第一个点。设柱面高为1,创建Z坐标矩阵 创建柱面的X坐标和Y矩阵 x=[-1,1,1,-1,-1]; %确定四边形X坐标 y=[-1,-1,1,1,-1]; %确定四边形Y坐标 X=[1;1]*x ;Y=[1;1]*y; %创建X和Y坐标矩阵 Z=[0;1]*ones(1,5); %创建柱面Z坐标矩阵 mesh(X,Y,Z) %绘图 colormap([0 0 0]) %设定黑色 axis off %去掉坐标框架 按顺序录入以下语句创建程序文件cylind1.m 比较cylinder([1,1],4)? 例2.7