MATLAB单简程序设计.ppt

MATLAB 程序设计方法 ;表达式——由变量、运算符、函数、数字组成 球坐标变换;例2.1 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记 ;+ 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂;运算优先级;命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息;例2.2 成都地理位置:北纬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.3 程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图 ;例2.4空间四面体底面是一三角形，将三角形三个顶点与四面体顶部的顶点相连接就形成三角形锥面。利用MATLAB的帧动画方法演示一个四面体旋转过程。 ;2. while循环 将循环体中语句循环执行直到条件不成立为止。 语法： while 条件表达式 可执行语句 end; 例2.5 3n + 1 问题. ;3. continue 命令;如果一个正整数 n 只能被1和它自身整除，则称这个数为素数(或质数)。例如：2、3、5、7、11、······;例2.7海伦公式;function 返回变量列表＝函数名(输入变量列表) %注解说明 输入变量检测,输出变量检测 函数体;红、绿两队从相距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)/(8+60) 当C向左运动时,A、C相遇时间: tk= (B – A)/(60+10) 利用相遇时间tk、及时计算A, B的必威体育精装版位置 当 (B– A)0.2 时程序结束。;while (B-A)0.2 if f==1 tk=(B-A)/(vb+vc); else tk=(B-A)/(vc+va); end A=A+va*tk;B=B-vb*tk; plot(A,0,R.,B,0,g.),pause(1) f=-f;k=k+1; end k,A,B;练习与思考题