表达式和运算符命令文件(程序文件)程序设计中的流程控制函数文件的编写格式表达式——由变量、运算符、函数、数字组成球坐标变换例2.1格林威治天文台建于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.4成都地理位置:北纬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.7程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图例2.8空间四面体底面是一三角形,将三角形三个顶点与四面体顶部的顶点相连接就形成三角形锥面。利用MATLAB的帧动画方法演示一个四面体旋转过程。2.while循环将循环体中语句循环执行直到条件不成立为止。语法：while条件表达式可执行语句end例2.93n+1问题.3.continue命令例2.10海伦公式立方倍积问题实验:作一立方体,使其体积等于单位正方体的两倍.对边长取不同位数,显示计算结果误差。红、绿两队从相距100公里的地点同时出发相向行军红队速度为10(公里/小时)绿队速度为8(公里/小时)开始时,通讯员骑摩托从红队出发为行进中的两队传递消息。摩托车的速度为60(公里/小时)往返于两队之间.每遇一队,立即折回驶向另一队.当两队距离小于0.2公里时,摩托车停止。计算通讯员驾驶摩托车跑了多少趟(从一队驶向另一队为一趟)。将红队、绿队和摩托车假设为A、B、C三个点.A点初始位置A=0,速度va=10(运动向右);B点初始位置B=100,速度vb=8(运动向左);C点初始位置C=0,速度vc=60(f=±1表示运动方向)当C向右运动时,C、B相遇时间:tk=(B–A)/(vb+vc)当C向左运动时,A、C相遇时间:tk=(B–A)/(va+vc)利用相遇时间tk及时计算A,B的最新位置当(B–A)<0.2时程序结束。while(B-A)>0.2iff==1t(k)=(B-A)/(vb+vc);elset(k)=(B-A)/(vc+va);endA=A+va*t(k);B=B-vb*t(k);plot(A,0,‘r.',B,0,'g.'),pause(1)f=-f;k=k+1;endkfunction返回变量列表＝函数名(输入变量列表)%注解说明输入变量检测,输出变量检测函数体functionY=young(n)ifnargin==0,n=3;endY=eye(n);Y(:,1)=ones(n,1);fork=3:nY(k,2:k-1)=Y(k-1,1:k-2)+Y(k-1,2:k-1);end函数内所有变量是局部变量，既不影响其他M文件中同名变量，也不被其他M文件中同名变量所影响函数文件中的输出变量要等于某个确定的表达式例2.13设地球半径6400km,以150经纬差绘三维地球。采用常规算法和递归算法分别编写计算数n的阶乘的函数文件问题分析：n张盘片，A、B、C三根柱子ifnargin==1,A='A';B='B';C='C';endifn==1disp(strcat(‘No’,in