MATLAB程序设计方法I表达式——由变量、运算符、函数、数字组成球坐标变换例2.2格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记+加-减*乘.*点乘/右除./点右除\左除.\点左除^矩阵幂.^点幂运算优先级————命令文件————命令文件——MATLAB命令的有序集合。文件执行——对文件中命令进行批处理，即从第一条命令开始按顺序执行，直到最后一条命令。如果中间某条命令出错，则中断并输出错误信息例2.3直线平行于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([000])%设定黑色axisoff%去掉坐标框架例2.7成都地理位置:北纬30o,东经104o;北京地理位置:北纬40o,东经116o。计算两城市距离近似值。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))%计算并显示距离条件控制;循环控制;错误控制;终止运行控制1.for循环语法：for变量=初值：步长：终值可执行语句end例2.11程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图例2.12空间四面体底面是一三角形，将三角形三个顶点与四面体顶部的顶点相连接就形成三角形锥面。利用MATLAB的帧动画方法演示一个四面体旋转过程。2.while循环将循环体中语句循环执行直到条件不成立为止。语法：while条件表达式可执行语句end例2.143n+1问题.3.continue命令如果一个正整数n只能被1和它自身整除，则称这个数为素数(或质数)。例如：2、3、5、7、11、······例2.17海伦公式立方倍积问题实验。作一立方体,使其体积等于单位正方体的两倍.对边长取不同位数,显示计算结果误差。练习与思考题