清华数学实验第二章MATLAB程序设计方法 II36.pptVIP

*/19 MATLAB程序设计方法 II 摩托车问题与数学模型 函数文件的编写格式 数据文件的输入方法 牟合方盖的几何图形 ? ? ? ? 红、绿两队从相距100公里的地点同时出发相向行军 红队速度为10(公里/小时)绿队速度为8(公里/小时) 开始时,通讯员骑摩托从红队出发为行进中的两队 传递消息。摩托车的速度为60(公里/小时)往返于两队 之间. 每遇一队,立即折回驶向另一队.当两队距离小 于0.2公里时,摩托车停止。计算通讯员驾驶摩托车跑 了多少趟(从一队驶向另一队为一趟)。 例2.15 摩托车问题 将红队、绿队和摩托车假设为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 who Your variables are: A f tk vb B k va vc A=0; B=100; va=10;vb=8;vc=60; f=1;k=0; plot(A,0,ro,B,0,go),hold on function 返回变量列表＝函数名(输入变量列表) %注解说明 输入变量检测,输出变量检测 函数体 ————函数文件的编写———— function Y=young(n) if nargin==0,n=3;end Y=eye(n);Y(:,1)=ones(n,1); for k=3:n Y(k,2:k-1)=Y(k-1,1:k-2)+Y(k-1,2:k-1); end 例2.18 杨辉三角形 函数内所有变量是局部变量，既不影响其他M文件 中同名变量，也不被其他M文件中同名变量所影响 函数文件中的输出变量要等于某个确定的表达式 函数文件的第一行必须按特定格式书写 ————函数文件的编写———— 输入/输出 变量检测命令：nargin 、 nargout 当函数文件被用户调用时(程序执行时) nargin返回函数被调用时输入变量的个数。 nargout返回当函数被调用时输出变量的个数。 例2.19 地球半径6400km,输入经差 绘三维地球。 function earthface(Dtheta) if nargin==0,Dtheta=15;end R=6400; theta=(-180:Dtheta:180)*pi/180; fai=(-90:Dtheta:90)*pi/180; X=R*cos(theta)*cos(fai); Y=R*cos(theta)*sin(fai); Z=R*sin(theta)*ones(size(fai)); colormap([0 0 1]) mesh(X,Y,Z),axis off 特例之一：有输入无输出变量的函数文件 ————函数文件的编写———— 设计网格点矩阵绘三维地球。 R=6400; [theta,fai]=meshgrid(-90:15:90,-180:15:180); theta=theta*pi/180; fai=fai*pi/180; X=R*cos(theta).*cos(fai); Y=R*cos(theta).*sin(fai); Z=R*sin(theta); colormap([0 0 1]) mesh(X,Y,Z),axis off 例2.21编写判素数子函数ispri(n)，主函数myprime调用该子函数并产生100以内的所有素数 function p=myprime() p=2;q=3:2:100; for k=1:49 n=q(k);T=ispri(n); if T,p=[p,n];end end function T=ispri(n) for k=2:n-1 if mod(