最新数学建模第三次作业.docx

  • 格式:docx
  • 大小:184.40 KB
  • 文档页数:11

精品文档院系:数学学院专业:信息与计算科学年级:2014 级学生姓名:王继禹学号:201401050335教师姓名:徐霞6.6 习题3.一个慢跑者在平面上沿着他喜欢的路径跑步,突然一只狗攻击他,这只狗以恒定速率跑向慢跑者,狗的运动方向始终指向慢跑者,计算并画出狗的轨迹。

解:(1)模型分析建立:狗的轨迹:在任意时刻,狗的速度向量都指向它的目标慢跑者。

假设 1:慢跑者在某路径上跑步,他的运动由两个函数 X(t)和 Y(t)描述。

假设 2:当 t=0 时,狗是在点 (x0,y0)处,在时刻 t 时,它的位置是 (x(t),y(t)) 那么下列方程成立:222(1)狗以恒定速率跑:X’+y’=w(2)狗的速度向量平行于慢跑者与狗的位置的差向量:将上述方程带入等式:,可得:再将λ代入第二个方程,可得狗的轨迹的微分方程:(2)程序及结果dog 函数[dog.m]function[zs,isterminal,direction] = dog(t,z,flag)global w; % w=speed of the dogX=jogger(t);h = X-z;nh=norm(h);if nargin<3 || isempty(flag)zs=(w/nh)*h;elseswitch (flag)case 'events'zs = nh-1e-3;isterminal = 1;direction = 0;otherwiseerror(['Unknow flag:'flag]);endend慢跑者的运动轨迹方程,水平向右[jogger.m]function s = jogger(t);s = [8*t;0];标记的函数[cross.m]function cross(Cx,Cy,v)Kx = [Cx Cx Cx Cx-v Cx+v];Ky = [Cy Cy+2.5*v Cy+1.5*v Cy+1.5*v Cy+1.5*v]plot(Kx,Ky);plot(Cx,Cy,'o' );主程序:静态显示[main1.m]global wy0 = [60;70];w=10;options = odeset('RelTol',1e-5,'Events', 'on' ); [t,Y] = ode23('dog' ,[0,20],y0,options);clf;hold on ;axis([-10,100,-10,70]);plot(Y(:,1),Y(:,2));J=[];for h=1:length(t),w = jogger(t(h));J=[J;w'];endplot(J(:,1),J(:,2),':');p = max(size(Y));cross(Y(p,1),Y(p,2),2)hold off ;动态显示[main2.m]global w;y0=[60;70];w=10;options = odeset('RelTol',1e-5,'Events', 'on' ); [t,Y]=ode23('dog' ,[0,20],y0,options); J=[];for h=1:length(t);w= jogger(t(h));J=[J;w'];endxmin = min(min(Y(:,1)),min(J(:,1)));xmax = max(max(Y(:,1)),max(J(:,1)));ymin = min(min(Y(:,2)),min(J(:,2)));ymax = max(max(Y(:,2)),max(J(:,2)));clf;hold on ;axis([xmin-10 xmax ymin-10 ymax]);title('The jogger and the Dog');for h = 1:length(t)-1,plot([Y(h,1),Y(h+1,1)],[Y(h,2),Y(h+1,2)],'-', 'Color', 'red', 'EraseMode ' , 'none');plot([J(h,1),J(h+1,1)],[J(h,2),J(h+1,2)],'-', 'Color', 'green', 'EraseMo de', 'none');drawnow;pause(0.1);endplot(J(:,1),J(:,2),':' );p = max(size(Y));cross(Y(p,1),Y(p,2),2)hold off;结果t=12.2761812635281,在 12.27 秒后狗追上慢跑者。

慢跑者轨迹是椭圆轨迹[jogger2.m]function s=jogger2(t)s=[10+20*cos(t) 20+15*sin(t)];狗的微分方程[dog.m]function[zs,isterminal,direction] = dog(t,z,flag)global w; % w=speed of the dogX=jogger2(t);h = X-z;nh=norm(h);if nargin<3 || isempty(flag)zs=(w/nh)*h;elseswitch (flag)case 'events'zs = nh-1e-3;isterminal = 1;direction = 0;otherwiseerror(['Unknow flag:'flag]);endend主程序[main3.m]global w;y0=[60;70];w=10;options = odeset('RelTol',1e-5,'Events', 'on' ); [t,Y]=ode23('dog' ,[0,20],y0,options); J=[];for h=1:length(t);w= jogger2(t(h));J=[J;w'];endxmin = min(min(Y(:,1)),min(J(:,1)));xmax = max(max(Y(:,1)),max(J(:,1)));ymin = min(min(Y(:,2)),min(J(:,2)));ymax = max(max(Y(:,2)),max(J(:,2)));clf;hold on ;axis([xmin-10 xmax ymin-10 ymax]);title('The jogger and the Dog');for h = 1:length(t)-1,plot([Y(h,1),Y(h+1,1)],[Y(h,2),Y(h+1,2)],'-', 'Color', 'red', 'EraseMode ' , 'none');plot([J(h,1),J(h+1,1)],[J(h,2),J(h+1,2)],'-', 'Color', 'green', 'EraseMo de', 'none');drawnow;pause(0.1);endplot(J(:,1),J(:,2),':' );p = max(size(Y));cross(Y(p,1),Y(p,2),2)hold off;结果取 w=25 有t=4.017776368842910,经过 4 秒左右狗追上慢跑者。

8.平面上有n(n>=2)个,任何两个都相交但无 3 个共点。

n 个把平面划分成多少个不通的区域?解:∵一个将平面分 2 份两个相交将平面分4=2+2 份,三个相交将平面分8=2+2+4 份,四个相交将平面分14=2+2+4+6 份,⋯平面内 n 个,其中每两个都相交于两点,且任意三个不相交于同一点,n 个分平面区域数 f ( n) =2+( n-1) n=n2-n+2明:( 1)当 n=1 ,一个把平面分成两个区域,而12-1+2=2 ,命成立.(2)假 n=k( k≥ 1),命成立,即 k 个把平面分成k2-k+2 个区域.当 n=k+1 ,第 k+1 个与原有的 k 个有 2k 个交点,些交点把第 k+1 个分成了2k 段弧,而其中的每一段弧都把它所在的区域分成了两部分,因此增加了2k 个区域,共有 k2-k+2+2k=( k+1)2-( k+1) +2 个区域.∴n=k+1 ,命也成立.由( 1)、( 2)知,任意的 n∈ N* ,命都成立.9.某人有n元,他每天一次物品,每次物品的品种很,或者一元的甲物品,或者二元的乙物品,他花完n 元有多少不同的方式?解: an 表示花完n 元的方案种数,若n=1,只能甲,有一种方法,故a1=1,若n=2,可以 2 个甲,或者 1 个乙或 1 个丙,即 a2 =3,当 n≥3 ,花的方式由甲和乙丙的种数之和构成,即 a n=a n-1+a n-2+a n-2=a n-1+2a n-2当 n≥ 3 , a n +a n-1=2( a n-1+a n-2),即{a n+1 +a n}是公比 q=2 的等比数列,首a2 +a1=1+3=4,a n+1+a n=4?2n-1=2n+1,n∴a n+a n-1=2 ,两式相减得a n+1-a n-1=2n+1-2-=2- ,(n≥ 2),若n 是奇数, a n=2n-1+2n-3+⋯ +22+a1=(2n+1-1)/3n-1 n-33n+1若 n 是偶数, a n=2 +2+⋯ +2 +a2=(2+1)/3.7.6 习题1.在化工生产中常常需要知道丙烷在各种温度T 和压力 P 下的导热系数K 。

下面是实验得到的一组数据:T /C686887871061061401409.798113.3249.007813.3559.791814.2779.656312.463P /103KPaK0.08480.08970.07620.08070.06960.07530.06110.0651试求 T =99 C 和P=10.3x103KPa下的 K。

解:找出温度T 相等时,导热系数 K 与压力 P 的关系。