【推荐】数学实验2.doc

数学实验.实验2

4.某日凌晨一住所发生一件凶杀案，警方于6时到达现场后测得尸温26℃，室温17℃，2小时后尸温下降了3℃

解：设t0为正常体温，y为死亡后人的体温，x为时间。差分方程由冷却定理得。

冷却定律：

T=C+(T0-C)*e^-kt

C是室温，T0是人体常温37℃

程序：

函数

functiony=swsj1(x)

%带参数求得k=0.2027，令a=k

t0=37;

a=0.2027

y=17+(t0-17)*exp(-a*x)

end

n=6;

x=0:n;

y=swsj1(x)

plot(x,y,r*)

gtext(死亡时间),

运行结果：

a=

0.2027

y=

37.000033.330530.334227.887725.890024.258922.9271

y=

37.000033.330530.334227.887725.890024.258922.9271

由结果知此人的死亡时间为2?:00.

5.据报道，某种山猫在较好、中等及较差的自然环境下，年平均增长率分别是1.68%，0.55%和-4.50%，假定开始时有100只山猫，按以下情况讨论山猫数量逐年变化过程及趋势：

（1）3种自然环境下25年的变化过程（作图）；

（2）如果每年捕获三只，会发生什么情况？山猫会灭绝吗？如果每年捕获一只呢？

（3）在较差的自然环境下，如果想使山猫数量稳定在60只左右，每年需要人工繁殖多少只？

（1）3种自然环境下25年的变化过程（作图）

程序:

functionx=cat(x0,n,r)

x=x0;

fork=1:n

x(k+1)=(1+r)*x(k);

end

k=(0:25);

y1=cat(100,25,0.0168);

y2=cat(100,25,0.0055);

y3=cat(100,25,-0.0450);

round([k,y1,y2,y3])

plot(k,y1,k,y2,r:,k,y3,g--),grid

gtext(较好环境),gtext(中等环境),gtext(较差环境)

结果：

ans=

0100100100

110210196

210310191

310510287

410710283

510910379

611110376

711210472

811410469

911610566

1011810663

1112010660

1212210758

1312410755

1412610852

1512810950

1613110948

1713311046

1813511044

1913711142

2014011240

2114211238

2214411336

2314711335

2414911433

2515211532

（2）如果每年捕获三只，会发生什么情况？山猫会灭绝吗？如果每年捕获一只呢？

（i）每年捕获三只

程序:

functionx=cat(x0,n,r)

x=x0;

fork=1:n

x(k+1)=(1+r)*x(k)-3;

end

command窗口程序同第一问

结果：

ans=

0100100100

1999893

2979585

3969378

4959072

5938866

6928560

7908354

8898049

9877743

108675