MATLAB简单程序设计.ppt

*/22 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); 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 命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息 M文件分为命令文件和函数文件两种 ①在编辑窗口中编写; ②保存并对文件命名； ③命令窗口键入文件名运行； ④观察运行结果； ······························· 例2.2 成都地理位置:北纬30o,东经104o;北京地理位置:北纬40o,东经116o。计算两城市距离近似值。 由两城市经纬度可计算出地心直角坐标系的坐标 P1(x1，y1，z1)，P2(x2，y2，z2) 向径之间的角度 球面短程线计算公式 ①输入经纬度数据和地球半径； ②转换两城市的经纬度为地心直角坐标数据； ③提取两个点的向径坐标； ④计算向径间的夹角和短程线长度并输出计算结果。 datas=[30,104;40,116]; %经纬度数据 R=6400; theta=datas(:,1)*pi/180; %弧度制 fai=datas(:,2)*pi/180; x=R*cos(theta).*cos(fai); %直角坐标转换 y=R*cos(theta).*sin(fai); z=R*sin(theta); Op1=[x(1),y(1),z(1)]; %第一城市坐标 Op2=[x(2),y(2),z(2)]; %第二城市坐标 d=R*acos(Op1*Op2/(R*R)) %计算并显示距离 按顺序录入程序文件(文件名:distance.m) 在命令窗口运行程序distance，屏幕将显示： d =1563.42 条件控制; 循环控制; 错误控制;终止运行控制 条件控制——有选择地运行程序块 1. if / elseif 语句 例2 .8 判润年程序 ————程序设计中的流程控制———— year=input(input year:=); n1=year/4; n2=year/100; n3=year/400;