数学实验第二章.pptx

命令文件(程序文件) 表达式和运算符 程序设计中的流程控制 函数文件的编写格式;————命令文件————;命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息;例2.2 生肖问题——08年是鼠年还是狗年？;表达式——由变量、运算符、函数、数字组成 球坐标变换;例2.3 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算天文台经线和纬线数据，在球面上绘出通过天文台的经纬线。 ;+ 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂;运算优先级;条件控制; 循环控制; 错误控制;终止运行控制;1. for 循环 语法： for 变量 =初值：步长：终值 可执行语句 end;例2.11 程序功能如下：把边长为2以原点为中心的正方形旋转pi/24，将其压缩(r=0.89),重复24次并绘图 ;例 利用帧动画方法演示马鞍面旋转过程。 ;2. while循环 将循环体中语句循环执行直到条件不成立为止。 语法： while 条件表达式 可执行语句 end; 例2.14 3n + 1 问题. ;3. continue 命令;红、绿两队从相距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的最新位置 A?A+10*tk; B?B – 8*tk. 当 (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;function 返回变量列表＝函数名(输入变量列表) %注解说明 输入变量检测,输出变量检测 函数体;函数内所有变量是局部变量，既不影响其他M文件 中同名变量，也不被其他M文件中同名变量所影响 函数文件中的输出变量要等于某个确定的表达式;例2.19 设地球半径6400km,以150经差绘三维地球。;例21编写递归函数,计算并输出裴波拉奇数列的第n项。 ;;问题分析：n张盘片，A、B、C三根柱子;if nargin==1, A= A;B= B;C= C;end if n==1 disp(strcat(No,int2str(n),:,A, - ,C)) else hanoi(n-1,A,C,B); disp( strcat(No,int2str(n),:,A, - ,C)) hanoi(n-1,B,A,C); end;No1: A - C No2: A - B No1: C - B No3: A - C No1: B - A No2: B - C No1: A - C;数据文件的输入方法;使用WinXP的记事本建立：data0.txt数据文件 ;n=length(x); x(n+1)=x(1);y(n+1)=y(1); dx=diff(x);dy=diff(y); Ln=sum(sqrt(dx.^2+dy.^2)); %计算多边形周长 Ln=20*Ln/1000 %周长换算为公里 An=400*An/1000000