數(shù)學(xué)實(shí)驗(yàn)死亡時(shí)間推測(cè)ppt完整

1、死亡時(shí)間推測(cè)死亡時(shí)間推測(cè) 小組成員：小組成員： 死亡時(shí)間（TIME OF DEATH)死亡時(shí)間在法醫(yī)學(xué)上是指死后經(jīng)歷時(shí)間或死后間隔時(shí)間，即發(fā)現(xiàn)、檢查尸體距死亡發(fā)生時(shí)的時(shí)間間隔。注重尸表檢查、判定，具有實(shí)際價(jià)值。死亡時(shí)間的推測(cè)是法醫(yī)學(xué)鑒定中首要解決的問題，其意義在于確定作案時(shí)間，認(rèn)定和排除嫌疑人有無作案時(shí)間，從而劃定偵查范圍乃至案件的最終的偵破。特別是涉及到多個(gè)嫌疑人時(shí)，死亡時(shí)間的準(zhǔn)確推測(cè)顯得尤為重要。法醫(yī)往往根據(jù)死后尸體變化發(fā)生的規(guī)律推測(cè)死亡時(shí)間。依據(jù)尸體變化發(fā)生的先后及法醫(yī)學(xué)實(shí)踐的需求，將死亡時(shí)間分為三個(gè)時(shí)期：死后早期、死后晚期和白骨期。某公寓發(fā)生一起謀殺案，死者是下午某公寓發(fā)生一起謀殺案，

2、死者是下午19:3019:30被發(fā)現(xiàn)的，法被發(fā)現(xiàn)的，法醫(yī)醫(yī)20:2020:20趕到現(xiàn)場(chǎng)，法醫(yī)在趕到現(xiàn)場(chǎng)，法醫(yī)在20:2020:20時(shí)，測(cè)得死者體溫為時(shí)，測(cè)得死者體溫為32.632.6C C ，一小時(shí)后，死者被移走時(shí)，又測(cè)量了一下體溫為一小時(shí)后，死者被移走時(shí)，又測(cè)量了一下體溫為31.431.4C C，當(dāng)，當(dāng)時(shí)室內(nèi)溫度與時(shí)間的關(guān)系如下表所示。經(jīng)過調(diào)查，種種跡象時(shí)室內(nèi)溫度與時(shí)間的關(guān)系如下表所示。經(jīng)過調(diào)查，種種跡象表明，此案最大的嫌疑犯是其單位的張某，但有人證明，張表明，此案最大的嫌疑犯是其單位的張某，但有人證明，張某下午某下午17:0017:00之前一直在辦公室，之前一直在辦公室， 17:0017:

3、00時(shí)張某才匆匆離開，時(shí)張某才匆匆離開，從其辦公室到公寓需要從其辦公室到公寓需要1010分鐘，此能否證明張某絕對(duì)不在現(xiàn)分鐘，此能否證明張某絕對(duì)不在現(xiàn)場(chǎng)？場(chǎng)？當(dāng)?shù)貧庀蟛块T其對(duì)一天室內(nèi)溫度的詳細(xì)記錄表當(dāng)?shù)貧庀蟛块T其對(duì)一天室內(nèi)溫度的詳細(xì)記錄表注：上表是時(shí)間段注：上表是時(shí)間段1717：0000 2121：2020每隔十分鐘一次的溫度每隔十分鐘一次的溫度記錄。記錄。 22.5322.4722.4122.3522.2922.2322.1722.1122.0521.9921.9421.8821.8321.7721.7221.6621.6121.5621.5121.4621.4021.3521.3021.25

4、21.2121.1621.11實(shí)驗(yàn)分析：首先確定案發(fā)時(shí)間，若死者死亡時(shí)間在5點(diǎn)10分之前，則張某就不是嫌疑人，否則就不能排除他的嫌疑。人體體溫受大腦神經(jīng)中樞調(diào)節(jié)。人死亡后體溫調(diào)節(jié)的功能喪失，尸體的溫度受外界溫度的影響。由物理知識(shí)可得，溫度高的物理下降的速度與其自身溫度同外界溫度的差值成正比。實(shí)驗(yàn)原理：設(shè) T(t) 表示t時(shí)刻尸體溫度，20:20為t=0， ，死者死亡時(shí)的體溫為 ，其中 t 就是死者的受害時(shí)間。K 為比例系數(shù)，X 為室溫。根據(jù)牛頓冷卻定理，設(shè)尸體的溫度為T，t為時(shí)間，-k為比例系數(shù)，則 CTCT4 .31) 1 (,6 .32)0(CtT5 .36)()(XTkdtdTCT6 .

5、32)0(實(shí)驗(yàn)程序：（1）擬合室溫X隨時(shí)間變化的函數(shù)t=0:1:26;T=22.53,22.47,22.41,22.35,22.29,22.23,22.17,22.11,22.05,21.99,21.94,21.88,21.83,21.77,21.72,21.66,21.61,21.56,21.51,21.46,21.40,21.35,21.30,21.25,21.21,21.16,21.11;p1=polyfit(t,T,1) p2=polyfit(t,T,2)z1=polyval(p1,t);z2=polyval(p2,t);r1=T-z1;r2=T-z2;s1=r1*r1,s2=r2*r

6、2 figure(1),plot(t,T,*,t,z1)figure(2),plot(t,T,*,t,z2) p1-0.0547 22.5006p2 0.0003-0.062322.5325s1=0.0071s2=2.2631e-004根據(jù)上式有：s1s2 ，所以函數(shù)二次擬合比一次擬合好。 即：T的擬合函數(shù)：T=0.0003t2-0.0623t+22.5325擬合結(jié)果：（2）求常微分方程的解析解T=dsolve(DT=-k*(T-(0.0003t2-0.0623*t+22.5325), T(0)=32.6) 運(yùn)行結(jié)果： T=(225325*k2+623*k+6)/(10000*k2)+(3*t

7、2)/10000-(exp(-k*t)*(-100675*k2+623*k+6)/(10000*k2)-(t*(623*k+6)/(10000*k) （3）令t=1，求比例系數(shù)kk=solve(225325*k2+623*k+6)/(10000*k2)+(3*t2)/10000-(exp(-k*t)*(-100675*k2+623*k+6) /(10000*k2)-(t*(623*k+6)/(10000*k)=31.4) 運(yùn) 行 結(jié) 果：k=0.12649478275160304876265778327726 這一頁答案不正確，但是步驟是正確的，只需將t=1換成t=6。然后得出正確的k值即可。因?yàn)閠的間隔是10分鐘，一個(gè)小時(shí)需要t=6才對(duì)，后面算出來的t，也需要這樣轉(zhuǎn)換。（4）求解t（當(dāng)T=36.5時(shí)，有t個(gè)小時(shí)）t=solve(-11/200*t+11/200/0.12+45/2+1/200*exp(-0.12*t)*(-11+2020*0.12)/0.12=36.5) 運(yùn) 行 結(jié) 果：t =-2.68:20-2:36=5:445:10所以