数学实验第二章课件.pptVIP

表达式和运算符 命令文件(程序文件) 程序设计中的流程控制 函数文件的编写格式;表达式——由变量、运算符、函数、数字组成 球坐标变换;例2.2 格林威治天文台建于1675年，其经度为零,纬度为北纬510。计算格林威治天文台在地心直角坐标系中的坐标并在球面上标记 ;例2.3 四边形柱面绘图原理和方法。 ;cylinder(R，N)创建单位高度绕Z轴旋转曲面，R为母线,N表示圆圈上点数,N缺省默认值为 20。 Cylinder %单位高度二十边形柱面 [x,y,z]=cylinder([1,1],4); mesh(x,y,z) %单位高度四边形柱面 [x,y,z]=cylinder([1,0],3); mesh(x,y,z) %单位高度的三角形锥面;+ 加 - 减 * 乘 .* 点乘 / 右除 ./ 点右除 \ 左除 .\ 点左除 ^ 矩阵幂 .^ 点幂;运算优先级;————命令文件————;命令文件——MATLAB命令的有序集合。 文件执行——对文件中命令进行批处理，即从第一 条命令开始按顺序执行，直到最后一条命令。如果 中间某条命令出错，则中断并输出错误信息;条件控制; 循环控制; 错误控制;终止运行控制;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的必威体育精装版位置 当 (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经差绘三维地球。;球函数的图形 ;function sphefun() t=linspace(0,pi,40);s=linspace(0,2*pi,60); [theta,fai]=meshgrid(t,s); x=cos(theta);w=sqrt(1-x.^2); P31=3*(5*x.^2-1)/2; YR=w.*P31.*cos(fai); X=YR.*sin(theta).*cos(fai); Y=YR.*sin(theta).*sin(fai); Z=YR.*cos(theta); surf(X,Y,Z),axis off;;问题分析：n张盘片，A、B、C三根柱子;if na