清华杨顶辉数学实验四作业.pdf

第3题 【模型建立】建立a和v 的函数关系，求出一系列时间点的v和a，然后利用数值积分的方 法求出一系列时间点的x. 函数 functiona=f(t,v) ift=60 a=(176.4*t+18280-0.4*v^2)/(1400-18*t); else a=-9.8-v^2/800; end 这是一个微分方程模型，初始条件是v0=0 【模型求解】 求出一系列时间点的v ts=0:0.01:71.4; v0=0; [t,v]=ode45(@f,ts,v0); [t,v] 得到关键数据 59.6300 267.2137 71.0000 0.6908 求出一系列时间点的a fori=1:1:7141 a(i,1)=f(t(i),v(i)); end; [t,a] 得关键数据 59.6400 0.9210 71.0000 -9.8006 求出一系列时间点的x fori=1:1:7141 ifi==1 x(i,1)=0; else x(i,1)=x(i-1,1)+0.01*v(i); end; end; x(6001) ans= 1.2189e+004 x(7141) ans= 1.3052e+004 【结论】引擎关闭瞬间火箭的高度是12.2km，速度为267.2m/s，加速度为0.92m/s^2；火箭 到达最高点时的高度为13.1km，加速度为-9.8m/s^2. 【心得】使用代码[t,v]=ode(@f,ts,v0)的时候，在定义函数时，要写成dv=f(t,v)的形式，而不 能写成dv=f(v,t)的形式，否则能计算并且结果是错的；ode45还是有一定误差，可能和我的 步长有关，理论上最大速度在60s取到，但是matlab 的计算结果是说在59.6s左右取到；关 闭引擎瞬间的加速度可能会和答案差别很大，因为在那个小区域上，加速度变化很快，数量 级在很快变化，一定的小误差就会导致和答案有数量级的差别。 第5题 (1)首先，将各物理量化为标准单位 水深 300*0.3048=91.44m 重力 2344.6N 质量 239.245kg 阻力系数 1.1813 临界速度 12.192m/s 浮力 2090.74N 由牛顿第二定律，重力-浮力-水的阻力=质量*加速度，即 2344.6-2090.74-1.1813v=239.245a 即a=dv/dt=1.061-0.00494v,v(0)=0 (2) 【模型建立】假设海洋够深，找出圆桶第一次达到40lbf/s 的速度时，下沉的深度，如果 这个深度大于300ft，则工程师赢，否则原子能委员会赢。 用数值方法 建立函数 functiona=f(t,v) a=-0.00494*v+1.061; end 这是一个微分方程模型 【模型求解】 求出一系列v 的值 ts=0:0.01:11.83; v0=0; [t,v]=ode45(@f,ts,v0); [t,v]得到关键数据 11.8300 12.1919，这说明在t=11.8300s时，圆桶刚刚达到12.1919m/s的速度（假设海洋 够深，圆桶能达到这样的速度） 求这个时刻，圆桶下降的深度 fori=1:1:1184 ifi==1 x(i,1)=0; else x(i,1)=x(i-1,1)+0.01*v(i); end; end; [t,x] 得到关键数据 11.8300 72.8785 由于72.878591.44，所以工程师赢 用解析方法 设v(t)=C(e^λt-1)，代入方程解得C=-214.8，λ=-0.00494 令v=0解得t=11.83s 对v(t)积分得x(t)=C((e^λt-1)/λ-t)，令t=11.83s，得x=72.8391.44 故工程师赢 【结论】假设海洋够深，那么圆桶在速度达到12.1919这样的临界速度时，只下沉到海底 72.83m处，由于圆桶的速度不断增加直到最后平稳，所以当圆桶继续下降，直到91.44m 的 海底时，就可能发生破裂，从而工程师赢。 【感想】这道题的数值解法与解析解法答案很吻合，证实了解析解法的准确性；本题题目全 部是用一套陌生的单位写出，我发现，新的单位总是标准单位的倍数，例如ft是m 的倍数， lbf是N 的倍数等等，所以我本来以为所有的计算都用新单位，然后得到了错误的答案。然 后我意识到这个