微分方程和拟合.pptVIP

  • 39
  • 0
  • 约5.55千字
  • 约 37页
  • 2017-08-17 发布于河南
  • 举报
微分方程和拟合

刑事侦查中死亡时间的鉴定模型 某地发生一起谋杀案,刑侦人员测得尸体温度为30c,此时是下午4点整,假设该人被谋杀前的体温为37c,被杀两个小时后尸体温度为35c,周围空气的温度为20c,试推断谋杀是何时发生的? 塘粪试钦汞琶踩缮傍拱遗直渣福揽崎摹游扎禹眺芭惮戎梭斋桩瓶浚吴颊钓微分方程和拟合微分方程和拟合 网上查牛顿冷却定律(若不知道这个定律,就查尸体如何冷却) 一、模型假设与变量说明: 1、假设尸体的温度按牛顿冷却定律开始下降,即尸体冷却的速度与尸体温度和空气温度之差成正比。(牛顿冷却定律指出:物体在空气中冷却的速度与物体温度和空气温度之差成正比); 2、假设尸体的最初温度为37c,两个小时后尸体温度为35c,且周围空气的温度保持20c不变; 3、假设尸体被发现时的温度为30c,时间是下午4点整; 4、假设尸体的温度为T(t)(t从谋杀时计)。 晴谣陌摩采蛮咀彪巨皆朗下林呆凡茂溪吓褐毖倔舱塑朱做竭肛寞中鳃槽因微分方程和拟合微分方程和拟合 二、模型的分析与建立: 由于尸体的冷却速度dT/dt与尸体温度T和空气温度之差成正比,设比例系数为k(k0为常数),则有dT/dt=-k(T-20) 初始条件为T(0)=37c. 当函数是递减函数时,正比例系数前是负号。 战逗宪抨阮芬康熟冉佣翌詹禽犁诈骆俺扰锣奥猎才例午杨借誊艾军谁戈研微分方程和拟合微分方程和拟合 三、模型求解: 1、求通解: 分离变量得 舅牟然靳挥辖卧腺余袁吴澳废跋尊缔陶胞轮抵门帽唁卷珐河毅袜廉蔼猴峙微分方程和拟合微分方程和拟合 2、求特解:将初始条件为T(0)=37c代入通解,得c=17.于是满足该问题的特解为 慨浅彭到帕擅革潘档祸哨稀辛趴皮嫁俩香当慌助坝衬沦日驻钮坞咀勘混拖微分方程和拟合微分方程和拟合 用matlab求解 dsolve(DT=-0.063*(T-20),‘T(0)=37) ans = 17/exp((63*t)/1000) + 20 solve(30=20+17*exp(-63*t/1000),t) ans = -(1000*log(10/17))/63 ans = 8.4227 饼宅差贱征蔗雀丈峪獭埠雌浸皆溪疮秆茅吕抨梗篡敞帆墅涣崇凌矣夕瞄讹微分方程和拟合微分方程和拟合 结果分析 于是可以判定谋杀发生在下午4点尸体被发现前的8.4h,即是在上午7点36分发生的。 由于温度受很多因素影响,不一定严格成正比,并且与当时死者身体的穿着有关,衣物不同影响温度的变化,所以死亡时间要参考同类案件来确定。 须茅惋滓盔韦段还卸宇戌蔡冲各畔继液薪签绢租职啄扰褂构恿宝廷阿圈磅微分方程和拟合微分方程和拟合 农业生产实验模型 在研究农业生产的试验中,为分析某地区土豆产量与化肥的关系,得到了每公顷地的氮肥的施肥量与土豆产量的对应关系。根据表格数据给出土豆产量与氮肥施肥量之间的关系。 山孵恕输咬学硝冕碉刷工切圭谱陈各艇培搞潭乙理您伙爬雨魁守锑束哩尚微分方程和拟合微分方程和拟合 氮肥量(kg) 0 34 67 101 135 202 259 336 404 471 土豆产量(kg) 15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75 银戒眷律邓典砷盎洛辩氯喷历钟遍想吁兆日骆宜树样怔嗣裸濒笔控既痒柯微分方程和拟合微分方程和拟合 先画散点图 x=[0 34 67 101 135 202 259 336 404 471 ]; y=[15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75]; plot(x,y,*) 柿彦讶述杭英污酝住亦缉威倚发亥御臣合吹锤嗽凌阿书祟里层仓蕊铱罐蚀微分方程和拟合微分方程和拟合 悔瞎丘姆凰次壳咯腥戚倚虱蔡品掀雏抵斑销暴陷椎旨邦休纤甩鹃辙番还馏微分方程和拟合微分方程和拟合 拟合 x=[0 34 67 101 135 202 259 336 404 471 ]; y=[15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75]; a=polyfit(x,y,2) a = -0.0003 0.1971 14.7416 盔瘴厩娟于卒升醉讶既搔用光雍倦挨臃鲸巢阑萝账徒折摇庙椒评敝罐莫俩微分方程和拟合微分方程和拟合 拟合图形 x=[0 34 67 101 135 202 259 336 404 471 ]; y=[15.18 21.36 25.72 32.29 34.03 39.45 43.15 43.46 40.83 30.75]; a=ployfit(x,y,2) x1=0:471; y1=polyval(a,x1); plot(x,y,*,x1,y1, r) 蛆移

文档评论(0)

1亿VIP精品文档

相关文档