人行走的最少能量消耗的迈步角度分布

 问题称述

一身高为H的人以一定的速度V行走,假设其重心在H/2的位置,并且腿长也为H/2,试分析该人迈多大的步伐,能使额外消耗的能量最少?

 

 

问题分析和模型建立

由于行走速度V是定值,所以不同的步伐会对应不同的频率.同时,不同的步伐只需要用两腿迈步时的交角θ就可以确定.所以本题可以用穷举θ∈［0，π］计算消耗的多余能量就可以解决。消耗的多余能量可以用人行走的重心上下移动而产生的重力作的功来计算:即E₀=mgΔh   其中, Δh=H-H^{`</sup>  H<sup>`}=f(θ, H) 再加上每一步所损失的生物能量E₁，总消耗能量为E=t *（E₁＋E₀ ）――其中t为人行走的频率.

 

 

求解方法

高为H，腿的张角为θ，人行走的时间为Δt，人的质量为m。

走的路长为S=V×Δt。

每步长为T＝sqrt（2×（H/2）²－cosθ×2×（H/2）²）。

上下移动距离为△h＝H/2×（1－cos(θ/2)）。

上下移动损失的能量为E₀=S/T×△h×m。

每一步所损失的生物能量E₁＝S/T×c(c为一常数)

总消耗能量为E=E₁＋E₀=S×（(△h×m)＋c）/T=[((1－cos(θ/2))×S×m)/(sqrt2×sqrt（1－cos(θ)）)]+S×c/（H/2）/(sqrt2×sqrt（1－cos(θ)）)。

依次计算出上述步骤的结果，使张角θ从0取到180，便可得到一个关于消耗能量E的函数图，得到一个最小值，即是本次实验所求结果。

四、计算机结果及分析

写出以上计算的源程序：

 

H=10;% 高为H

V=10;%V 速度

t1=10;%t1 为人行走的时间Δt

m=10;% 人的质量为m

c=10;%c 为任意值,表每一步生物能量损失

for a=1:180;%a 表示θ

E=(1-cos(a/2*(pi/180)))*V*t1*m/sqrt(2-2*cos(a*(pi/180)))+V*t1*c*2/H/sqrt(2-2*cos(a*(pi/180)));%E 为走路损失的总能量

    b(a)=E;%b 为存值的变量

end;

plot(b)

for o=1:180;% 求消耗能量最少的角度aa

    if b(o)==min(b)

        aa=o

    end;

end;

得到如下图：

当把（H=112;%高为H，V=132;%V速度，t1=12;%t1为人行走的时间Δt，m=145;%人的质量为m，c=145;%c为任意值,表每一步生物能量损失）这些量加以变化成此值后，发现函数图像又变成了另一种样子：

于是推断，此最值应还与上述几个量有关，经公式演算判断，此最值与速度V,时间t1无关。与H,m,c有关,而与H,m关系不是很大，而且H,m应该是一个人的定值，故此次建模取男人高为H=1.8; 质量为m=70公斤、女人高为H=1.5; 质量为m=45公斤为例。

 

继续编程――取男人高为H=1.8; 质量为m=70公斤计算得到θ-c图：

c 步长为1

对局部进行放大：（下图）

得到其jd()图（显示出对应的c取值时人行走消耗能量最少的角度θ的值）：图如下：

故此次建模取男人高为H=1.8; 质量为m=70公斤行走的结论如下：

高为H(m)	质量m（kg）	每步损失能量C(焦耳)	θ(度)
1.8	70	1	6
1.8	70	50	44
1.8	70	100	61
1.8	70	150	72