METHOD FOR DETERMINING THE METABOLIC PARAMETERS OF A SUBJECT

Information

  • Patent Application
  • 20230157581
  • Publication Number
    20230157581
  • Date Filed
    November 23, 2021
    2 years ago
  • Date Published
    May 25, 2023
    a year ago
Abstract
The present invention relates to a method of determining metabolic parameters of one or several muscular motor units by coupling muscular response to ventilatory response of a subject, comprising, from a step of performing efforts (i) using a test device by the subject moving on the test device; (ii)-a step of instantaneously measuring the movement speed of one or several muscular motor units of the subject v(t), and the oxygen flow rate φO2t); (iii)-a step of plotting the COT curve of the oxygen flow rate of the subject against the movement speed of one or several muscular motor units of the subject v(t), (iv)-a step of determining the following metabolic parameters from the COT curve: Basal power B which is the vertical asymptote at the origin, Flow resistance RM which is the oblique asymptote at high intensities.
Description
FIELD OF THE INVENTION

The present invention relates to a method of, and a device for characterizing the metabolic response of a subject or an individual to exercise.


State of the Art


The study of the response to exercise is an old subject covering many fields of researches, from physiology, and more specifically sports physiology, to biomechanics.


Energy expenditure has been known and studied since the end of the 19th century, including from a thermodynamic point of view via a calorimetric approach, and then a metabolic and nutritional approach.


Studies on ventilation and ventilatory adaptation are also widespread and are integrated into exercise testing.


The book by the authors Monod, Flandrois, and Vandewalle (Physiologie du sport 6th edition Masson (2007)), makes an exhaustive review of these issues for the humans. Equivalent review can be found, for example, for the horse, in the book “The Athletic Horse, Principles and Practice of Equine Sports Medicine”, edited by David R. Hodgson, Kenneth Harrington McKeever and Catherine M. McGowan.


To date, there is no single protocol for obtaining this characterization directly. The existing measurements and protocols are carried out separately and do not allow a synthetic mapping of the response to exercise of an individual or a subject. It should be noted that the main protocols for identifying an individual's energy potential are based on the following measurements requiring specific equipment:

    • Ventilation, measured by oxygen Cost of Transport (COT) curves;
    • Heart rate and rhythm (ECG), and blood gases;
    • Force-movement speed response of muscles (associated with the Hill 1938 and 1964 model);
    • Electromyograms.


Moreover, producing a synthetic mapping is only accessible at the cost of multiple protocols on different equipment,


In addition, the publication (Goupil et al. 2019. New J. Phys. 21, 023021 “Thermodynamics of metabolic energy conversion under muscle load”) builds the thermodynamic foundation for the response of a muscle, or a set of muscles, under stress.


It demonstrates that two metabolic parameters can be estimated on the basis of the classical measurement of the force-movement speed response of the muscle: the isometric force Fiso, on the one hand, and a feedback resistance, Rfb, on the other hand.


DISCLOSURE OF THE INVENTION

One purpose of the invention is to allow the quantification of the health condition and performance in terms of energetic potential of an individual from an experimental protocol comprised of two tests. The key innovative character of this invention emerges from the thermodynamic description of the response to exercise, that provide the possibility to bring together the results of these tests in a single descriptive framework.


The method and device of the invention allow for the mapping of a subject's performance and energetic potential in a single test protocol by directly identifying muscular response coupled with ventilatory response of the subject.


To date, there is no single protocol for obtaining this characterization directly. The existing measurements and protocols are carried out separately and do not allow a synthetic mapping of the response to exercise of an individual or a subject.


At the end of the proposed protocol carried out with the device, the synthetic mapping involves 5 traits that are easily calculated:

  • Basal power, noted B
  • Threshold metabolic intensity, noted IT
  • Metabolic resistance, noted RM
  • Feedback resistance, noted Rfb
  • Isometric force, noted Fiso


All these parameters together give value of metabolic merit factor, Fm=(Fiso Rfb IT)/(B RM), which is the main indicator of the individual's potentiality.


As such, the proposed device can be used both for monitoring sports performance and for functional rehabilitation.


More specifically, the present invention presents a method of determining metabolic parameters of one or several muscle motor units of any subject by coupling muscle response with ventilatory response of the subject.


Comprising the following steps:

  • (i)+(ii)-a step of acquiring movement speed of one or several muscle motor units of the subject as a function of time v(t), and an oxygen flow rate (φO2(t) of the subject;
  • (iii)-a step of plotting the COT curve: the oxygen flow rate against the movement speed of one or several muscle motor units v(t),


the COT curve being coupled to the muscle response by the following equation:






COT
=



φ
02

v

=


a

0



+


R
m


v

+

B
v







with B the basal power of one or several muscle motor units


RM the flow resistance of one or several muscle motor units

  • (iv)-a step of determining the following metabolic parameters from the COT curve:
    • Basal power B which is the vertical asymptote at the origin;
    • Flow resistance RM which is the oblique asymptote;
    • a0 which is the y-intercept of the oblique asymptote.


Indeed, the COT curve has a vertical asymptote at the origin and an oblique asymptote also called the oblique asymptote at high intensities. Alternatively, the step of acquiring a movement speed value of one or several muscular motor units of a subject v(t) and an oxygen flow rate value φO2(t) can be replaced by a step of performing exercise (i) using a test device and a step of instantaneously measuring (ii) the movement speed of one or several muscular motor units of the subject v(t) and the oxygen flow rate φO2(t). Alternatively, the step of acquiring may be preceded by a step of instantaneously measuring the movement speed of one or several muscle motor units of the subject v(t), and of the oxygen flow rate φO2(t) of the subject. In addition, these values can be measured instantaneously during the performance of exercise by the subject performing an exercise test or a stress test on a test device.


According to Other Optional Features of the Method:



  • the step (i)+(ii) acquiring also a force F(t),

  • the step (iii) plotting also a second curve: the force F(t) against movement speed v(t) of one or several muscle motor units,



governed by the following equation from Hill's muscle model:





(F+a) (v+b)=c

  • in step (iv), the following metabolic parameters are determined from the second curve:
    • Isometric force Fiso corresponding to the force at zero movement speed,
    • Feedback resistance Rfb corresponding to Fiso/vx with vx equal to the movement speed at zero force,
    • Threshold metabolic intensity IT,


      with the following equations:










I
T

=
b



1.












R
fb

=



a
0


I
T


=


a
0

b





2.












F
iso

=


c

I
T


-

a
0





3.






Alternatively, the step of instantaneously measuring (ii) the movement speed of one or more muscular motor units of the subject v(t), and of the oxygen flow rate φO2(t) may be completed by a step of acquiring a value of developed power P(t) and/or force F(t).

  • in step (iv), the parameters are determined at different times
  • a figure of merit fm is calculated as fm=(Rfb* IT* Fiso)/(RM*B)
  • the movement speed of one or several muscle motor units was measured from a test device corresponding to a treadmill or a bicycle
  • the movement speed of one or several muscle motor units v(t), and the oxygen flow rate φO2(t) were obtained when the subject was subjected to a variable load Zload(t) which is an instantaneous time-modulated load;
  • the step (i)+(ii) acquiring also as a function of time a force F(t);
  • in the step (iv) of determining metabolic parameters of one or several muscle motor units, by:
    • calculating Laplace transforms F(p) and V(p) from the force F(t) and movement speed v(t);
    • calculating impedances Z(p)=F(p)/V(p);
    • calculating, from the impedances Z(p), the following metabolic parameters:
      • Isometric force Fiso,
      • A resistance impedance of organism of the subject ZIM(p)


The organism being represented by the isometric force and the resistance impedance of the organism such that:






F(p)=Fiso*Z(p)/(ZIM(p)+Z(p))






V(p)=Fiso/(ZIM(p)+Z(p))

    • Calculating from ZIM(p) the real part which is the metabolic parameter R(IM) and then the following metabolic parameters:
      • Feedback resistance Rfb
      • Threshold metabolic intensity IT
      • Basal power B
      • Flow resistance RM


Alternatively, the step of acquisition as a function of time may correspond to or be preceded by a step of performing exercise, with the subject moving and being subjected to a variable load which is an instantaneous time-modulated load, as well as to a step of instantaneous measurement, measuring the movement speed of one or several assemblies of muscular motor units, of the developed power or force and of the oxygen flow.

  • a metabolic power resulting from the difference between input and output powers is modeled,
  • total mechanical resistance of the output is calculated as






R
in
=R
M
+R
fb
=R
M+(Fiso+a0)/(vT+v),


the measured values being the mechanical force FM and the movement speed v, and


with the force Fiso=FB constant, Rin is determined by modulating a load RL (imposed effort)

  • the test device has real-time controlled braking to produce time-modulated braking over a frequency range of about 10−2 Hz to 10 Hz, to achieve the load RL
  • the braking is performed by a motor, mechanically coupled to a driven wheel of the test device, and which acts under active load.
  • the shape of the time modulation signals of the load RLis:
  • a time-harmonic, variable-frequency load;
  • an impulse load;
  • an indexed load.


According to another embodiment, the invention relates to a device for determining metabolic parameters of one or several muscle motor units coupling muscle response to ventilatory response of a subject, comprising:

  • means for acquiring efforts coupled to a test device, to measure as a function of time: at least the movement speed of one or several muscular motor units of the subject v(t), and oxygen flow rate φO2(t);
  • processor for plotting a COT curve: the oxygen flow rate of the subject against the movement speed of one or several muscle motor units v(t),


the COT curve being coupled to the muscle response by the following equation:






COT
=



φ
02

v

=


a

0



+


R
m


v

+

B
v







with B the basal power of one or several muscle motor units;


RM the flow resistance of one or several muscle motor units;


the processor being configured to determine the following metabolic parameters from the COT curve:

    • Basal power B which is the vertical asymptote at the origin;
    • Flow resistance RM which is the oblique asymptote;
    • a0 which is the y-intercept of the oblique asymptote.


According to other optional features of the device:

  • the means for acquiring are configured to measure a force F(t),
  • the processor is configured to plot a second curve: the force F(t) against movement speed v(t) of the one or several muscle motor units, governed by the following equation from Hill's muscle model:





(F+a) (v+b)=c

  • the processor is configured to determine the following metabolic parameters from the second curve:
    • Isometric force Fiso corresponding to the force at zero movement speed,
    • Feedback resistance Rfb corresponding to Fiso/vx with vx equal to the movement speed at zero force,
    • Threshold metabolic intensity IT,


      With the following equations:










I
T

=
b



1.












R
fb

=



a
0


I
T


=


a
0

b





2.












F
iso

=


c

I
T


-

a
0





3.






  • the processor is configured to determine a figure of merit fm equal to fm=(Rfb*IT*Fiso)/(RM*B)

  • the test device is of the treadmill or bicycle type,

  • the means for acquiring are configured to measure a force F(t),

  • the test device is configured to expose the subject to a variable load Zload(t) which is an instantaneous time-modulated load;

  • the processor is configured to determine metabolic parameters of one or several muscle motor units, by:
    • calculating the Laplace transforms F(p) and V(p) from the force f(t) and movement speed v(t);
    • calculating the impedances Z(p)=F(p)/V(p);
    • calculating, from the impedances Z(p), the following metabolic parameters:
      • Isometric force Fiso,
      • A resistance impedance of organism of the subject ZIM(p)


        the organism being represented by the isometric force and the resistance impedance of the organism such that:







F(p)=Fiso*Z(p)/(ZIM(p)+Z(p))






V(p)=Fiso/(ZIM(p)+Z(p))

    • calculating from ZIM(p) the real part which is the metabolic parameter R(IM) and then the following metabolic parameters:
      • Feedback resistance Rfb
      • Threshold metabolic intensity IT
      • Basal power B
      • Flow resistance RM
    • the test device has real-time controlled braking to produce time-modulated braking over a frequency range of about 10−2 Hz to 10 Hz, for example to achieve a load RL
    • the braking is performed by a motor, mechanically coupled to a driven wheel of the test device, and which acts under active load.





DESCRIPTION OF THE FIGURES

Other objectives, features and advantages become apparent from the following detailed description with reference to the drawings given for illustrative and non-exhaustive purposes:



FIG. 1 is an example of a force-movement speed response (top), and a power-movement speed response (bottom) for three different muscles: slow muscles (dashed line), medium muscles (dotted line), fast muscles (solid line);



FIG. 2 shows a plot of ϕ− and COT as a function of metabolic intensity for different bodies;



FIGS. 3.a and 3.b show the principle of the invention which allows the determination of metabolic parameters from COT measurement (FIG. 3.a) and from Hill's parametrization (FIG. 3.b);



FIG. 4 shows a diagram of a device according to the invention;



FIG. 5 shows the schematization of the balance between power a) and the associated nodal formalism b);



FIG. 6 shows a synthetic mapping of the metabolic parameters of two individuals or subjects.





DESCRIPTION OF THE INVENTION

The present invention relates to the determination of intrinsic and extrinsic metabolic parameters of one or several muscular motor units consisting of one or several muscles defining an identity card of a subject's physical condition. The determination is based on:

    • A thermodynamic model defining the parameters unambiguously.
    • An exercise measurement device according to an appropriate protocol and parameters extraction.


Each muscular motor unit thus consists of a muscle, and the metabolic parameters are generally determined on several motor units considered.


Preferably, the method according to the invention is implemented by at least one computing means 3, such as a processor 3 associated with a memory. With the one or more processors being adapted to, preferably configured to plot a COT curve :the subject's oxygen flow rate against the subject's movement speed v(t), and to determine metabolic parameters from the curve. In addition, the one or more processors are capable of, preferably configured to execute the various steps of a method according to the invention.


Furthermore, by “means for acquiring” 2a & 2b, is meant, within the meaning of the invention, a test device 1 and/or measurement means. The terms “subject” and “individual” in the present invention, refer to humans but also to any animal species, domestic or wild. In the present invention, the terms “subject” and “individual” are used without distinction.


THEORETICAL SUMMARY

The core of the modeling is based on metabolic power production by a body (or an organism) and is derived from the publication (Christophe Goupil et al. 2019 New J. Phys. 21, 023021 “Thermodynamics of metabolic energy conversion under muscle load”).


A metabolic power can be broken down into a metabolic force and a metabolic intensity, the product of which gives the metabolic power.


It is considered that the metabolic power is transmitted, via the limbs, in the form of mechanical power and movement speed which are the quantities actually measured experimentally.


The thermodynamic model developed is in perfect agreement with Hill's phenomenological modeling, for which it provides the thermodynamic foundations.


Metabolic power production results from the difference in input, Φ+, and output, Φ−, powers of the organism.


Φ+ corresponds to nutrient supply (mainly glucose), while Φ− is a waste stream that includes excess heat and especially the production of lactic acid-like reaction metabolites.


The point is key, because it follows that the flow Φ− can be considered, to a large extent, proportional to the flow of oxygen consumed by the organism.


The theoretical model leads to the following expressions:







Φ
+

=





F
iso



I
M

/

η
c


+
B



I
T

+

I
M





I
T









Φ
-

=



R
M



I
M
2


+





R
fb



I
M
2


+



F
iso

(


1

η
c


-
1

)



I
M


+
B



I
T

+

I
M





I
T







With:

    • IM: metabolic intensity (considered proportional to the movement speed of movement of the limb, which is measured);
    • Fiso: isometric force (extracted from the exercise curves);
    • nc: maximum thermodynamic efficiency;
    • B: basal power (extracted from the exercise curves);
    • IT: threshold metabolic intensity (extracted from the exercise curves).


These concepts were introduced in the publication (Christophe Goupil et al. 2019 New J. Phys.21 023021 “Thermodynamics of metabolic energy conversion under muscle load”).


The metabolic power results from the difference between the input and the output powers, and is written:






P
=
 



Φ
+

-

Φ
-


=



F
M



I
M


=


[


F
iso

-


(


R
M

+



F
iso

+


R
fb



I
T





I
T

+

I
M




)



I
M



]



I
M








The metabolic power involves two additional terms compared to the previous ones, they are RM and Rfb, with:

    • RM: flow resistance (extracted from the exercise curves) also called the viscous resistance;
    • Rfb: feedback resistance (extracted from the exercise curves).


The five-parameter parameterization thus contains both metabolic (IT, B) and mechanical (Fiso,RM) parameters. As regards the term Rfb, it has a metabolic origin but it translates mechanically.


All the parameters can be combined into a metabolic merit factor which is, in a way, the summary scalar of the subject's health condition,







f
M

=



R
fb



F
iso



I
T




R
M


B






All the parameters are now established, there remains to specify some relationships and experimental methods.


Hill's Parameterization


As mentioned above, the metabolic model provided the thermodynamic basis for Hill's muscle model.


Hill proposes a phenomenological expression of the force-movement speed response of one or several muscle motor units, according to:





(F+a) (v+b)=c


With:

F mechanical force (measured);


v: movement speed of contraction (measured).


The theoretical model developed resulted in the following correspondences:






a=R
fbIT+RMIM=a0+RMIM≈a0=cte





b=IT=cte






c=(Fiso+RfbIT) IT=cte


The typical force-movement speed response curve can be traced back to the values of the three parameters, since the overall shape is defined by two extreme points, Fiso and the movement speed at zero force.







v
X





F
iso


R
fb


.





The characteristic dip of the curve defines IT.


The following parameters are therefore extracted:

    • Fiso
    • Rfb
    • IT


Balance of the Parameters


We have five parameters:

    • Fiso: isometric force (extracted from the exercise curves);
    • B: basal power (extracted from the exercise curves);
    • IT: threshold metabolic intensity (extracted from the exercise curves);
    • RM: flow resistance (extracted from the exercise curves);
    • Rfb: feedback resistance (extracted from the exercise curves).


And we have a force-movement speed experiment (FIG. 1 and FIGS. 3.a and 3.b) that determines three of them:










I
T

=
b



1.












R
fb

=



a
0


I
T


=


a
0

b





2.












F
iso

=


c

I
T


-

a
0





3.






There are still two to be determined. This is the role of the COT measurement.


COI and COT


Using the proportion between the oxygen flow rate and ϕ−, the missing terms can be extracted.


The Cost of Oxygen Index (COI) is defined as:






COI
=



Φ
-


I
M


=


a
0

+


R
M



I
M


+

B

I
M








The shape of the plot is shown in FIG. 2.


It should be noted that the curve has a vertical asymptote at the origin, and an oblique asymptote at high intensities. These two asymptotes allow the extraction of, respectively:

    • B
    • RM


It should also be noted that that the y-intercept also makes it possible to find Hill's parameter a0. All five parameters characterizing the organism are therefore accessible.


In the case of movement, and using the proportion between IM and movement speed, the oxygen Cost of Transport (COT) can be defined, which is a quantity widely used in the literature.


In fact, from a dimensional point of view, the COT corresponds to the amount of energy waste produced to move body of the subject over a unit length.


Its expression is:






COT
=



φ
02

v

=


a

0



+


R
m


v

+

B
v







The form of the COT is therefore isomorphic to that of the COI.


Conclusion: by means of an exercise test measuring powers, movement speeds and oxygen flows, it is possible to calculate the five metabolic parameters that define the mapping of a subject's condition and performance. In addition, the merit factor provides a summary of the subject's state.


Experimental Embodiment

From an experimental point of view, a treadmill or bicycle type exercise test device 1 should be considered, as shown in FIG. 4.


The quantities to be measured will be:

  • Walking or pedaling movement speed;
  • Power developed (the force or torque produced will be extracted, knowing the movement speed);
  • Oxygen consumption.


Protocol


Principle


From a practical point of view, the measurements are equivalent to extracting the mechanical output impedance of the organism.


Indeed, the equivalent schematic according to a nodal approach is illustrated in FIG. 5. The one or several limbs, or more precisely the one or several acting muscle units, are assimilated to a source of mechanical power delivering a force, here represented by a generator. In series with the generator are the resistances RM and Rfb which induce a power loss by viscous dissipation.


The actual mechanical power available is at the end of the chain. The power results from the product of the force available at the end of the chain, multiplied by the movement speed of movement of the limbs.


The mechanical output resistance is written






R
in
=R+R
H(v)=RM+(Fiso+a0)/(vT+v)


Knowing that the measured quantities are the mechanical force FM and the movement speed v, on the one hand, and that the force Fiso=FB is constant, on the other hand, the value Rin can be extracted by modulating the load RL. Depending on the type of protocol envisaged, the load may be a braking force that prevents movement, or a mass that opposes movement, as in the case of so-called “leg-press” devices.)


Embodiment

The one skilled in the art will understand that the experimental embodiment can be carried out by means of various devices (treadmill, bicycle, leg press), depending on whether the determination of the metabolic parameters relates to a limb, a portion or all of the organism.


Let us consider the case of a training bike type device, as shown in FIG. 4. The load RL is experimentally achieved by real-time controlled braking to produce time-modulated braking over a frequency range of about 10−2 Hz-10 Hz. Braking may be performed by a motor, mechanically coupled to the driven wheel, and which acts under active load.


The experimental procedure is equivalent to a harmonic load as encountered in impedance spectroscopy measurements. The electrical analog of the mechanical protocol is the impedance spectroscopy technique used in electrochemistry to determine the performance of electric batteries.


The shape of the modulation signals of the load RL can be of different types, following the tradition of impedance spectroscopy analyses:

  • a time-harmonic, variable-frequency load;
  • an impulse load, the load being equivalent to abruptly changing the braking and immediately canceling it (not recommended due to risk of injury);
  • an indexed load, the load being equivalent to abruptly changing the braking and maintaining it (a technique known as “Wingate” in English), but its exploitation is generally not made beyond the analysis of the shape of the curves).


Thus, at the end of the protocol carried out on the equipment, the synthetic mapping is composed of the following elements:

  • Basal power;
  • Chemical energy release rate;
  • Metabolic, static and dynamic viscosity;
  • Mechanical, static and dynamic viscosity;
  • Isometric force.


Measured under stationary stress, these quantities have real values.


On the other hand, when measured under transient stress, the viscosity parameters have an imaginary part which accounts for the elastic and inertial terms.


All of these parameters together are used to obtain the value of the metabolic merit factor, which is the main indicator of the individual's potentiality. As such, the proposed equipment can be used both for monitoring sports performance and for functional rehabilitation. It makes it possible to obtain a synthetic mapping of the response of an individual to exercise, which is not possible to obtain using the state of the art where the measurements and the existing protocols are carried out separately.


As shown in FIG. 6, two individuals have different values of metabolic parameters.


Individual I1 has a fairly low basal power, high threshold isometric strength and metabolic intensity. These last two quantities show that the individual I1 is capable of subjecting organism of the subject to a significant physical effort. Low basal power of the subject characterizes a dominant of slow type muscles. It is thus an individual adapted to an endurance type effort, rather than an explosive one. The fairly high values of the flow and feedback resistances show that there is room for improvement both physically, through training (lowering of the flow resistance), and possibly metabolically (lowering of the feedback resistance) through better metabolization of nutrients.


Individual I2 has high basal power, low threshold intensity and low isometric strength. The resistances are also of low value. In this case it is about an individual having presumably a very good potential for explosive effort (important Basal), but who is strongly reduced by the weakness of the threshold intensity. This contradiction can be explained by the presence of an injury that does not allow a normally performing body to express its full potential. It can therefore be observed that the collection of the five parameters makes it possible to distinguish between the measured performance and an individual's potentialities.

Claims
  • 1. A method of determining metabolic parameters of one or several muscle motor units of a subject by coupling muscle response to ventilatory response of said subject, comprising the following steps: (i)+(ii)-a step of acquiring movement speed of one or several muscle motor units of the subject as a function of time v(t), and an oxygen flow rate φO2(t) of the subject;(iii)-a step of plotting a COT curve of the oxygen flow rate against the movement speed of one or several muscle motor units v(t), the COT curve being coupled to the muscle response by the following equation:
  • 2. The method according to claim 1, wherein: the step (i)+(ii) comprises acquiring also a developed force F(t),the step (iii) comprises plotting also a second curve of the force F(t) against the movement speed v(t) of one or several muscle motor units, governed by the following equation from Hill's muscle model: (F+a) (v+b=c4in step (iv), the following metabolic parameters are determined from said second curve: Isometric force Fiso corresponding to the force at zero movement speed,Feedback resistance Rfb corresponding to Fiso/vx with vx equal to the movement speed at zero force,Threshold metabolic intensity IT,
  • 3. The method according to claim 1, wherein in step (iv) the parameters are determined at different times.
  • 4. The method according to claim 2, wherein a factor of merit fm is calculated as fm=(Rfb*IT*Fiso)/(RM*B).
  • 5. The method according to claim 1, wherein the movement speed of one or several muscle motor units was measured from a test device corresponding to a treadmill or a bicycle.
  • 6. The method according to claim 2, wherein: the movement speed of one or several muscle motor units v(t), and the oxygen flow rate φO2(t) were obtained when the subject was subjected to a variable load Zload(t) which is an instantaneous time-modulated load;the step (i)+(ii) comprises acquiring a developed force F(t);the step (iv) of determining metabolic parameters of one or several muscle motor units, comprises: calculating Laplace transforms F(p) and V(p) from the force F(t) and movement speed v(t);calculating impedances Z(p)=F(p)/V(p);calculating, from the impedances Z(p), the following metabolic parameters: Isometric force Fiso,A resistance impedance of organism of the subject ZIM(p);the organism being represented by the isometric force and the resistance impedance of the organism of the subject such that: F(p)=Fiso*Z(p)/(ZIM(p)+Z(p))V(p)=Fiso/(ZIm(p)+Z(p))calculating from ZIm(p) the real part which is the metabolic parameter R(IM) and then the following metabolic parameters: Feedback resistance Rfb, Threshold metabolic intensity IT,Basal power B,Flow resistance RM.
  • 7. The method according to claim 2, wherein: a metabolic power resulting from the difference between input and output powers is modeled,total mechanical resistance of the output is calculated as Rin=RM+Rfb=RM+(Fiso+a0)/(vT+v),the measured values being the mechanical force FM and the movement speed v, andwith the force Fiso=FB constant, Rin is determined by modulating a load RL.
  • 8. The method according to claim 7, wherein the movement speed of one or more several muscle motor units was measured from a test device corresponding to a treadmill or a bicycle, said the test device having a real-time controlled braking to produce time-modulated braking over a frequency range of about 10−2 Hz to 10 Hz, to achieve the load RL.
  • 9. The method according to claim 8, wherein the braking is performed by a motor, mechanically coupled to a driven wheel of the test device, and which acts under active load.
  • 10. The method according to claim 9, wherein the shape of time modulation signals of the load RL is: a time-harmonic, variable-frequency load;an impulse load; oran indexed load.
  • 11. A device for determining metabolic parameters of one or several muscle motor units of a subject by coupling muscle response to ventilatory response of said subject, comprising: means for acquiring efforts coupled to a test device, to measure as a function of time: at least movement speed of one or several muscular motor units of the subject v(t), and oxygen flow rate φO2(t);a processor for plotting a COT curve of the oxygen flow rate of the subject against the movement speed of one or several muscle motor units v(t), the COT curve being coupled to the muscle response by the following equation:
  • 12. The device according to claim 11, wherein: the means for acquiring are configured to measure a developed force F(t),the processor is configured to plot a second curve: of the developed force F(t) against movement speed v(t) curve of the one or several muscle motor units, governed by the following equation from Hill's muscle model: (F+a) (v+b)=c the processor is configured to determine the following metabolic parameters from the second curve: Isometric force Fiso corresponding to the force at zero movement speed,Feedback resistance Rfb corresponding to Fiso/vx with vx equal to the movement speed at zero force,Threshold metabolic intensity IT,With the following equations:
  • 13. The device according to claim 11, wherein the processor is configured to determine a factor of merit fm equal to fm=(Rfb*IT*Fiso)/(RM*B).
  • 14. The device according to claim 11, wherein the test device is of the treadmill or bicycle type.
  • 15. The device according to claim 11, wherein: the means for acquiring are configured to measure a force F(t),the test device is configured to expose the subject to a variable load Zload(t) which is an instantaneous time-modulated load;the processor is configured to determine metabolic parameters of one or several muscle motor units, by: calculating Laplace transforms F(p) and V(p) from the force F(t) and movement speed v(t);calculating impedances Z(p)=F(p)/V(p);calculating, from the impedances Z(p), the following metabolic parameters: Isometric force Fiso,A resistance impedance of organism of the subject ZIM(p)the organism of the subject being represented by the isometric force and the resistance impedance of the organism such that: F(p)=Fiso*Z(p)/(ZIM(p)+Z(p))V(p)=Fiso/(ZIM(p)+Z(p))calculating from ZIM(p) the real part which is the metabolic parameter R(IM) and then the following metabolic parameters: Feedback resistance Rfb;Threshold metabolic intensity IT;Basal power B;Flow resistance RM.
  • 16. The device according to claim 15, wherein the test device has real-time controlled braking to produce a time-modulated braking over a frequency range of about 10−2 Hz to 10 Hz.
  • 17. The device according to claim 16, wherein the braking is performed by a motor, mechanically coupled to a driven wheel of the test device, and which acts under active load.