ESAIM: P ROCEEDINGS , VOL .

9, 2000, 65-92
AUTOMATIQUE,BIOLOGIE ET SANT É c EDPSciences, SMAI 2000
[Link]

Le système cardio-vasculaire et sa régulation par le

système nerveux autonome : modélisation et mesures
J. Bestel ∗ J. Clairambault † C. Médigue∗ A. Monti∗ M. Sorine∗
Mai 1999

Résumé
Le système cardio-vasculaire (SCV) est modélisé dans son ensemble, avec ses différents comparti-
ments : pompe cardiaque, petite et grande circulation, circuit artériel et circuit veineux. Les problèmes de
modélisation sont abordés sous l’angle de la commande, c’est-à-dire en fonction d’une mission à remplir
pour le SCV et son contrôleur à court terme, le système nerveux autonome (SNA). Deux aspects du modèle
sont développés ici : d’une part, la pompe cardiaque, intégrant un modèle non linéaire de muscle, et son
insertion dans la circulation sanguine, d’autre part l’arc baroréflexe, principale boucle de régulation à court
terme de la pression artérielle. Des applications cliniques en cours analysant le gain de cet arc baroréflexe
sont également mentionnées.

1 Le système cardio-vasculaire vu comme un système contrôlé

(Résumé de [1], consultable en annexe ; pour plus de détails, on pourra se référer à [2, 3, 4])

1.1 Le cœur et les vaisseaux

La principale mission du système cardio-vasculaire (SCV), la fourniture d’oxygène et de nutri-
ments aux différents tissus de l’organisme, est assurée par une double circulation sanguine, systé-
mique (du cœur vers les tissus par les artères et retour par les veines) et pulmonaire (du cœur vers
le poumon et retour), chacun de ces deux compartiments étant mis en mouvement par la pompe car-
diaque correspondante (ventricule gauche ou droit). Cette double pompe est contrôlée, en fréquence
et en volume éjecté à chaque battement, par le système nerveux autonome (SNA).
La fréquence cardiaque, principal facteur d’adaptation du débit, est donnée par un pacemaker
naturel soumis au contrôle du SNA : le nœud sino-atrial, situé dans l’oreillette droite (il existe des
modèles mathématiques au niveau de la cellule, dérivés du modèle nerveux de Hodgkin et Huxley, de
ce pacemaker ; ils ne seront pas rappelés ici, une revue sommaire pouvant en être trouvée dans [5]).
Le volume de sang éjecté dans la circulation systémique à chaque battement, ou volume d’éjection,
dépend de 3 facteurs : a/ la précharge, ou volume de sang présent dans le ventricule gauche juste avant
l’éjection (volume dépendant de la pression du retour veineux au cœur, dite pression veineuse centrale,
∗
INRIA Rocquencourt
†
Université Paris VIII et INRIA Rocquencourt

65
Article published by EDP Sciences and available at [Link] or [Link]
PVC) ; b/ la contractilité du myocarde, sous la dépendance du SNA ; c/ la postcharge, ensemble des
facteurs qui s’opposent au travail du cœur : résistance pariétale du ventricule, impédance de l’aorte,
et résistances périphériques, ces dernières étant également sous la dépendance du SNA.
Le circuit vasculaire n’est pas seulement pourvu de résistances, il est aussi pourvu de capacités.
Les vaisseaux sanguins ne sont en effet pas rigides, mais compliants, et les veines systémiques, no-
tamment, constituent un important réservoir de sang. Les quatre compartiments du circuit vasculaire :
artériel et veineux systémique, artériel et veineux pulmonaire sont ainsi modélisables en première in-
tention chacun par une cellule du premier ordre (ou éventuellement du second ordre, si l’on prend en
compte les effets d’inertie du sang). La principale variable d’état de chacun de ces quatre comparti-
ments est la pression sanguine, et celle qui joue le rôle le plus important est la plus élevé : la pression
artérielle systémique.

1.2 Contrôle de la pression systémique par le SNA

Les barorécepteurs artériels sont sensibles à l’étirement de la paroi ; Ils transmettent des messages
nerveux, sous la forme de trains d’impulsions électriques, codés en fréquence, par des fibres afférentes
(= de la périphérie vers le centre) qui ont leur terminaison dans le tronc cérébral.
Les barorécepteurs du circuit veineux sont situés dans l’oreillette droite ; ils envoient aux centres
du tronc cérébral des messages qui ont pour effet d’adapter le contrôleur (le SNA) à la pression
veineuse de retour et aux variations de la pression intrathoracique dues à la respiration.
Du tronc cérébral partent les deux branches efférentes du SNA, sympathique et parasympathique,
qui ont des effets antagonistes sur la pression artérielle systémique : une stimulation sympathique
augmente la fréquence cardiaque et la contractilité des fibres myocardiques, augmente aussi les résis-
tances périphériques, et au total la pression, tandis qu’une stimulation parasympathique, principale-
ment à destination du pacemaker sinusal, ralentit le cœur et fait chuter la pression.
Ces voies ont des temps d’action très différents : la voie sympathique, cardioaccélatrice, est effi-
cace en une dizaine de battements cardiaques, tandis que le délai d’action du vague transmis par la
voie parasympathique est très court, de l’ordre du battement cardiaque (on parle de “coup de frein
vagal”). Leur action est de trois types : action sur la fréquence cardiaque, soit par modification du
rythme du pacemaker, soit par modification de la vitesse de transmission de l’influx électrique dans
la paroi du cœur ; action sur la puissance de la pompe, soit directement sur la vitesse de contraction
des fibres musculaires, soit sur la tension et la relaxation de ces fibres ; action enfin sur les résistances
périphériques : vasoconstriction ou vasodilatation.

2 Principes de modélisation de la circulation sanguine

On présente dans ce qui suit un modèle mathématique de la pompe cardiaque inséré dans la circu-
lation sanguine. L’ensemble est vu comme un système dynamique mécanique dont on étudie l’activité
périodique, initiée par une commande chimique elle-même reliée au pacemaker électrique du nœud
sinusal.
Un modèle simple de la géométrie du ventricule permet de relier les variables contrainte et défor-
mation de la fibre aux variables pression et volume du ventricule. Le reste de la circulation est un
ensemble de compartiments caractérisés par des constantes de temps de remplissage et de vidage.

66
Ce modèle, parce qu’il repose sur des considérations physiques, rend compte de phénomènes qualita-
tifs clefs.

2.1 Modèle de fibre musculaire

On suppose un montage de type Hill-Maxwell ([3], [4]) : les éléments contractiles (EC), actifs,
sont en série avec des éléments élastiques linéaires passifs qui “absorbent” les déformations quand la
fibre ne peut pas se déformer dans son ensemble (phases isométriques) ; le tout est en parallèle avec
des éléments élastiques non linéaires (EP) passifs, qui assurent plutôt la cohésion de l’ensemble.
La contrainte totale σ dans la direction des fibres est la somme des contraintes de EC et de EP :

σ = σC + σP

On note εC la déformation de EC et ε celle de EP (aussi égale à la déformation totale).

2.1.1 Lois de comportement

a/ Partie contractile :
Les mesures de l’activité du muscle en général et du myocarde en particulier montrent que la force
développée n’est pas proportionnelle à la déformation de la fibre. Nous proposons (voir [6]) le mo-
dèle :

σ̇C = −k1 σC |ε̇C | + k2 ε̇C + kC u (1)

où k1 ,k2 , et kC sont des raideurs.

Si on fait u = 0 dans (1) on obtient une relation contrainte-déformation élasto-plastique, qui fournit
une loi de comportement satisfaisante pour le muscle (voir [7]).
De plus, les parties contractiles (sarcomères) sont actives : elles développent une force en réponse à
une commande chimique (modification de la concentration en ions calcium intracellulaire). Le terme
en u en rend compte dans (1), il permet en particulier de “démarrer” la contraction à partir du repos.
Ce modèle vaut pour la fibre entière (mise en série de sarcomères) qui, dans le cas du myocarde,
est sollicitée périodiquement ; nous proposons en première intention pour u un “créneau” (fonction
constante par morceaux, ne prenant que deux valeurs), périodique à la période cardiaque, pouvant
correspondre par un intermédiaire chimique aux deux temps électriques (dépolarisation / repolarisa-
tion) de la membrane de la cellule excitable cardiaque.

b/ Éléments passifs :
L’élément passif parallèle contribue à augmenter la contrainte au-delà d’une déformation seuil (pour
protéger le cœur des trop grands remplissages : il concourt à la fermeture de la valve quand la défor-
mation devient trop grande). On propose la loi de comportement : σP (ε) = kP |ε − ε̄|+ , où kP est une
raideur, et ε̄ représente une déformation seuil.

2.1.2 Équation du mouvement

Des mesures réelles montrent que lors de la contraction et de la relaxation isovolumiques la paroi
endocavitaire est soumise à une accélération importante qui a une signification physiologique (voir

67
[8]). On en rend compte en écrivant la relation de la dynamique pour les éléments contractiles ; on
propose l’équation (2) dans laquelle on a ajouté un terme d’amortissement (frottement visqueux) et
un autre de rappel élastique :
ε̈C + 2ζω ε̇C + ω 2 εC = −α σC + β ε (2)
Les paramètres ζ,ω,α,β sont reliés à ceux, physiques, caractérisant la fibre (raideur, masse, longueur
au repos).

2.2 Modèle de ventricule

On suppose un seul degré de liberté pour la déformation du cœur. Par exemple, on assimile le
cœur à un cylindre de hauteur a, de rayon b et d’épaisseur h. On suppose que la déformation a lieu
uniquement suivant a : a = a0 (1 + ε) et b = b0 .
V est le volume du cylindre ; on a V = V0 (1 + ε) et donc V̇ = V0 ε̇, V0 étant le volume au repos.
Ce modèle fournit donc une relation simple entre la déformation de la fibre et celle de l’objet en trois
dimensions.
On suppose dans un premier temps que l’épaisseur du ventricule est petite et constante : h = h0 .
La contrainte σ est alors reliée à la pression dans le ventricule P par l’intermédiaire de la relation de
Laplace :
2h0
P = σ (3)
b0
La pression interne est donc donnée, en notant V̄ = V0 (1 + ε̄) par (4) :
 
2h0 kP
P (V,σC ) = σC + V − V̄ + (4)
b0 V0

2.3 Le ventricule inséré dans la circulation

On couple les équations (1) et (2) à l’équation de conservation du volume ventriculaire ; par
exemple, on obtient pour le ventricule gauche l’équation (5) :
1 1
V̇ = |Ppv − P (V,σC )|+ − |P (V,σC ) − Psa |+ (5)
Rpv Rsa
où P (V,σC ) est donnée par (4).
Rpv et Rsa sont les résistances à l’écoulement du sang à travers les valves cardiaques respectivement
veineuse (mitrale) et artérielle (sigmoïdes aortiques).
D’autre part les pressions Ppv (t) et Psa (t) désignent respectivement les pressions pulmonaire veineuse
et systémique artérielle. Ce sont des variables d’état du système complet ; par exemple, l’évolution de
la pression artérielle systémique est régie par l’équation différentielle (6) :
1 1
Csa Ṗsa = |P (V,σC ) − Psa |+ − (Psa − Psv ) (6)
Rsa RS
où Psv (t) est la pression systémique veineuse ; Csa est la compliance artérielle systémique et RS est
la résistance totale systémique.
On pourra se référer à [9] pour la vision schématisée (par compartiments) de la boucle complète de
circulation du sang.

68
2.4 Comportements qualitatifs
On présente une liste non exhaustive de comportements qui servent à tester qualitativement le
modèle.

2.4.1 Contraction isotonique

Par définition de la contraction isotonique, σ̇C = 0, avec ε̇C < 0, alors (1) donne

kC
ε̇C = − u (7)
k1 σC + k2
Cette équation rend compte de résultats expérimentaux selon lesquels le muscle se raccourcit d’autant
moins vite qu’il doit équilibrer une charge plus importante (voir p. ex. [3]). Si cette charge est nulle,
kC
le muscle se raccourcit à la vitesse maximale vmax qui, d’après (7), vaut − u.
k2
À l’échelle du sarcomère la physiologie relie vmax aux taux d’hydrolyse de l’ATP par le complexe
Actine-Myosine ([4]) ; ceci donne donc une information sur l’amplitude de u durant sa phase positive.

2.4.2 Fonctionnement de la pompe et loi de Starling

Quand le cœur se remplit, son volume V augmente, ce qui se traduit par une augmentation de
ε. Donc, d’après (4), à partir d’un certain volume de remplissage, P augmente, jusqu’à dépasser la
pression dans le compartiment veineux, ce qui provoque la fermeture de la valve (voir (5)). À cet
instant P est égale à la précharge. Puis P continue d’augmenter, maintenant grâce à σC , jusqu’à
atteindre la pression dans le compartiment artériel : postcharge (voir [2] pour une introduction à ces
notions, et [4] pour une description détaillée).
Si on élargit le cadre au contexte non pacemaker, on peut donc constater que : plus le cœur aura
admis de sang - sous l’effet d’une pression veineuse importante et/ou d’une pression intraventriculaire
de début de remplissage faible - plus la composante passive de P contribuera à l’éjection, i.e. à
atteindre la pression dans le compartiment artériel. Comme le débit de sortie est fonction de cette
différence de pression de part et d’autre de la valve (voir (5)), le volume éjecté en sera d’autant plus
important. À cet effet peut s’ajouter celui d’une contractilité augmentée (inotropisme positif) à travers
la commande chimique sur σC .

2.4.3 Contraction et relaxation isovolumiques

En situation isométrique les équations (1) et (2) forment un système isolé de la circulation exté-
rieure, que l’on peut raccrocher à des grandeurs discrètes caractéristiques d’un cycle cardiaque :
pendant la contraction isovolumique le volume intraventriculaire V est constant et vaut le volume
télédiastolique ; il reflète donc la précharge. Pendant la relaxation isovolumique V est égal au volume
télésystolique, qui est l’image de la postcharge.

2.5 Modèle de la circulation complète

On donne les équations du système dynamique complet, dont chacun des paramètres a été introduit
précédemment pour le ventricule gauche.

69
2.5.1 Ventricule gauche


 σ̇CG = −k1G σCG |ε̇CG | − k2G ε̇CG + kCG u






 VG − V0G
ε̈CG = − 2ξG ωG ε̇CG − ωG 2 εCG − γG σCG + δG
 V0G (8)





 1 1
 V̇G = |Ppv − PG (VG ,σCG )|+ − |PG (VG ,σCG ) − Psa |+
Rpv Rsa
 
2h0G kPG
où PG (VG ,σCG ) = σCG + VG − V̄G +
b0G V0G

2.5.2 Ventricule droit


 σ̇CD = −k1D σCD |ε̇CD | − k2D ε̇CD + kCD u






 VD − V0D
ε̈CD = − 2ξD ωD ε̇CD − ωD 2 εCD − γD σCD + δD
 V0D (9)





 1 1
 V̇D = |Psv − PD (VD ,σCD )|+ − |PD (VD ,σCD ) − Ppa |+
Rsv Rpa
 
2h0D kPD
où PD (VD ,σCD ) = σCD + VD − V̄D +
b0D V0D

2.5.3 Circuit vasculaire


 1 1

 C sa Ṗ sa = |P G − P sa | − (Psa − Psv )

 R +
R
 sa
1
S
1
Csv Ṗsv = − |Psv − PD |+ + (Psa − Psv ) (10)

 Rsv RS

 1 1

 Cpa Ṗpa = |PD − Ppa |+ − (Ppa − Ppv )
Rpa RP

avec Cpv Ppv (t) = Csa Psa (0) + Csv Psv (0) + Cpa Ppa (0) + Cpv Ppv (0)
− [Csa Psa (t) + Csv Psv (t) + Cpa Ppa (t)] + VT
− [VG (t) + VD (t) + Vsa (0) + Vsv (0) + Vpa (0) + Vpv (0)]

On a exprimé ici que sur les échelles de temps considérées le volume total VT est constant.

70
3 Modélisation du baroréflexe et mesures chez l’homme
Le système nerveux autonome joue un rôle pivot dans la régulation à court terme de la pression
artérielle via l’arc baroréflexe. Celui-ci, de façon très simplifiée, peut être modélisé comme sur la
fig. 1, qui met en évidence les principaux mécanismes de contrôle de la pression artérielle par le
SNA.
SNA
voies efférentes

EFFET CHRONOTROPIQUE

SYMPATHIQUE régulation de la fréquence cardiaque

EFFET CHRONOTROPIQUE SYSTEME

PARASYMPATHIQUE

CARDIO

EFFET INOTROPIQUE régulation de la contractilité

SYMPATHIQUE VASCULAIRE
du myocarde

EFFET PERIPHERIQUE régulation au niveau du

SYMPATHIQUE système vasculaire périphérique

voies afférentes pression

BARORECEPTEURS
artérielle

F IG . 1 – Arc baroréflexe: schéma de la régulation de la pression artérielle. Le SNA peut être vu com-
ment le contrôleur. Quatre différents mécanismes de contrôle, divisés en parasympathique et sympa-
thiques, ont été inclus dans le modèle.(voir le texte pour une brève introduction ; pour une description
complète du modèle voir [10] ou [11]).

Le modèle qui a été développé est un modèle en boucle fermée caractérisé par les éléments clas-
siquement présents dans une boucle de contrôle :
– les capteurs et les voies afférentes
dans le modèle, ils ont été réduits simplement aux barorécepteurs à haute pression. On ne consi-
dère ni les récepteurs à basse pression ni l’intégration des informations au niveau central.
– mécanismes de contrôle de la pression artérielle
• réflexe chronotropique : contrôle par les branches efférentes antagonistes du SNA du rythme
cardiaque ;
• réflexe inotropique : contrôle par la voie sympathique de la contractilité du myocarde ;
• réflexe vasculaire périphérique : contrôle par la voie sympathique du calibre des vaisseaux
périphériques ;

71
 – organes effecteurs
• le pacemaker cardiaque : le SNA, agissant sur le nœud sinusal, est capable de régler la
fréquence de contraction du muscle cardiaque;
• la contractilité du myocarde : le SNA peut agir sur la puissance de la pompe cardiaque ;
• les résistances périphériques : l’augmentation de l’activité du nerf sympathique sur le sys-
tème vasculaire périphérique induit une vasoconstriction qui se traduit par une augmenta-
tion de ces résistances.
Le modèle a été développé dans [10] et repris dans [11] ; il reproduit les principales caractéristiques du
système cardiovasculaire et ses réponses aux stimuli extérieurs : en particulier il reproduit l’Arythmie
Sinusale Respiratoire (ASR) et l’apparition d’une composante en basse fréquence grâce à la présence
d’un comportement résonnant autour de 0.1 Hz (ondes de Mayer). Des développements récents du
modèle visent à rendre mieux compte d’autres phénomènes : entraı̂nement des basses fréquences par
la respiration, réponse du système cardiovasculaire à un changement de posture (tilt-up test), évalua-
tion de la fonction cardiaque (en particulier dans l’insuffisance cardiaque).

3.1 Arc baroréflexe : pression artérielle - rythme cardiaque

Le but est de décrire, par de simples indicateurs numériques, l’interaction entre la pression arté-
rielle (PA) et le rythme cardiaque (signal RR, suite des délais entre battements cardiaques successifs,
détectés sur l’électrocardiogramme). De tels indicateurs doivent permettre d’évaluer le fonctionne-
ment du SNA, à l’etat physiologique comme dans les états pathologiques. De ce point de vue, il
est licite de simplifier le modèle, de façon à mettre en évidence dans l’arc baroréflexe seulement le
contrôle du RR par le SNA (fig. 2). Sur cette figure est représentée l’interaction entre le centre respi-

CENTRE
RESPIRATOIRE
effets
mécaniques
effets
centraux
fréquence système
SNA
cardiaque cardiovasculaire

pression
barorécepteurs
artérielle

F IG . 2 – Schéma simplifié de l’arc de contrôle du rythme cardiaque (RR), et des interactions aux
deux niveaux, central et mécanique, de la respiration et du système cardio-vasculaire.

ratoire et le système cardiovasculaire : d’un côté une interaction au niveau central (interaction directe
entre les centres autonomes bulbaires qui contrôlent le cœur et le centre respiratoire bulbaire), de
l’autre, une interaction située plutôt au niveau mécanique (intrathoracique), due à la modification du
retour veineux à l’inspiration.

72
3.1.1 Point de vue entrée-sortie

Pour évaluer les influences de la PA sur le RR, le système étant en boucle fermée, on doit supposer
faibles les interactions centrales entre le centre respiratoire et le centre autonome qui règle le rythme
cardiaque. Sous cette hypothèse, on se place dans un contexte de relation entrée - sortie, où l’entrée
est la PA et la sortie le RR. Il devient naturel d’étudier conjointement les variations des deux signaux
pour bien interpréter le fonctionnement du contrôleur.

3.1.2 Le rôle de la respiration

De nombreuses recherches ont démontré que la fréquence et le volume respiratoires ont une large
influence sur la variabilité du RR et de la PA. Pour que les fluctuations de l’intervalle RR puissent
être considérées comme des indicateurs fiables du contrôle nerveux, en particulier quand l’ASR est
utilisée pour quantifier le tonus vagal, il faut donc faire en sorte que les interactions entre la respiration
et le système nerveux autonome soient constantes durant chaque phase de mesures cliniques. Ainsi
dans l’expérience qui sera décrite ci-dessous, nous avons fixé la fréquence respiratoire (imposant aux
sujets de rythmer leur respiration à l’aide d’un métronome), de façon à obtenir une bande bien définie
en haute fréquence (HF), la même pour tous les sujets.

3.1.3 Analyse spectrale

L’analyse spectrale de la variabilité du RR et de la PA, et leur étude conjointe par l’ananalyse

interspectrale constituent des outils complémentaires d’analyse du SNA. Si xi et yi sont les i-ème
valeurs du RR et de la PA (systolique, par exemple) on peut calculer les variables suivantes :

PN −1
Transformée de Fourier Discrète X(k∆f ) = i=0 xi wi e−j2πk∆f iTo

P M2
Densité spectrale d’énergie D̄xx (k∆f ) = 1 1
M +1 T k=− M
|X(k∆f )|2
2

1 1
P M2
Densité spectrale d’énergie croisée D̄xy (k∆f ) = M +1 T k=− M
X̃(k∆f ) Y (k∆f )
2

|D̄xy (k∆f )|
Module de la fonction de transfert |Hxy (k∆f )| =
D̄xx (k∆f )
 
−1 imag[D̄xy (k∆f )]
Phase de la fonction de transfert Φxy (k∆f ) = tan
real[D̄xy (k∆f )]

2 |D̄xy (k∆f )|2

Fonction de Cohérence γxy (k∆f ) =
D̄xx (k∆f )D̄yy (k∆f )
où X̃(k∆f ) est le conjugué de X(k∆f ), To est le pas d’échantillonage en temps, ∆f est le pas
d’échantillonnage en fréquence et wi est la i-ème composante de la fenêtre de lissage. Pour évaluer
l’état du système nerveux autonome, on calcule les paramètres suivants dans chaque fenêtre fréquen-

73
tielle, en basse (description des ondes de Mayer) et en haute fréquence (liée à l’activité respiratoire) :

P
Énergie x∗F = k∆f ∈∗F D̄xx (k∆f )

Cohérence moyenne dans la fenêtre coh∗F = 1 P γ 2 (k∆f )

L∗F k∆f ∈∗F xy
P
Gain moyen dans la fenêtre g∗F = L1 k∆f ∈∗F |Hxy (k∆f )|
∗F

où la bande de fréquence ∗F peut être HF (haute fréquence) ou LF (low frequency, i.e. basse fré-
quence), L∗F est la largeur de la fenêtre considérée. Chaque signal cardiovasculaire a donc été carac-
térisé par son contenu énergétique en haute et en basse fréquence.

3.1.4 Proposition d’indicateurs pour l’évaluation de l’arc baroréflexe

On introduit ici deux autres indicateurs, déduits d’une analyse du modèle de l’arc baroréflexe, et
qui vont décrire l’état du fonctionnement du SNA :

Gain normalisé : ḡ∗F = g∗F xmean

ymean
ḡ
Index de sensibilité : ḡ LF
HF

Ces indicateurs peuvent être calculés pour chaque bande de fréquence (ce qui à notre connaissance
n’avait encore jamais été fait), l’index de sensibilité permettant alors la comparaison directe entre le
comportement du contrôleur en haute et en basse fréquence séparément. Le gain a été normalisé par
les valeurs moyennes du RR et de la PA, ce qui permet de prendre en compte le point de consigne (en
pression) du système cardiovasculaire.

3.1.5 Résultats expérimentaux

Les changements du SNA provoqués par un stimulus extérieur (tilt-test) ont été évalués en exami-
nant les signaux cardio-vasculaires et respiratoire. On s’est limité à l’étude des périodes stationnaires,
et on a donc exclu de l’analyse les périodes qui suivent des changements de respiration (spontanée
versus contrôlée) ou de posture (allongée versus debout).
Le protocole est représenté sur la figure 3, il se décompose en 3 étapes successives : 5 minutes
(mn) de respiration spontanée en position allongée ; 5 mn de respiration contrôlée à 0.25 Hz (on
contrôle seulement la fréquence respiratoire, pas le volume) et finalement, après un changement de
posture (tilt-test passif), 5 mn en gardant la fréquence respiratoire contrôlée à 0.25 Hz .
On constate (fig. 4) en basse fréquence (LF) une augmentation importante de l’énergie pour la
PAS, et une augmentation moins évidente pour le RR. L’énergie en haute fréquence (HF) disparaı̂t
pour le RR (inhibition parasympathique) tandis qu’elle augmente pour la PAS (effet mécanique de la
respiration plus important).

74
 Tilt 60o
Position couchée Position debout

Respiration Respiration Respiration

spontanée controlée à 0.25Hz controlée à 0.25Hz

0 5 10 15

périodes sélectionées pour l’analyse des signaux

F IG . 3 – protocole expérimental utilisé pour étudier le SNA et la réponse du système cardiovasculaire

à des changements de position.

RR PAS Respiration

1
0 11
00
00
11
0
1
0
1 00
11
0
1 00
11
00
11
0
1
0
1 00
11
0000000000
1111111111
0000000000
1111111111
0.1 0.25
0000000000
1111111111
0000000000
1111111111
0.1 0.25
1111111111
0000000000
0000000000
1111111111
0.25

1
0
0
1
0
1
0
1
0
1
1
0 0
1
0
1
0
1 0
1
0
1 0
1
0
1 0
1
0
1
0
1 0
1
0
1
0
1 0
1
0
1 0
1
0
1 0
1
0
1
0
1
111111111
000000000 111111111
000000000 111111111
000000000
000000000
111111111
0.1 0.25 0.1 0.25 000000000
111111111
000000000
111111111
0.25

F IG . 4 – Variations des composantes spectrales des signaux cardio-vasculaires induites par tilt-test.
En haut en position couché, en bas en position debout. En abscisses sont représentées les fréquences
et en ordonnées la puissance spectrale [ms2 /Hz] pour le rythme cardiaque (RR), [mmHg 2 /Hz]
pour la pression artérielle systolique (PAS), et pour la respiration [en unités arbitraires]. L’énergie
en basse fréquence est représentée par une barre verticale noire, celle en haute fréquence par une
barre grise.

75
 Comme on peut le remarquer à l’aide d’une analyse temps-fréquence, ou en comparant les ré-
sultats obtenus par DFT (Transformation de Fourier Discrète) sur plusieurs fenêtres d’analyse, les
composantes en basse fréquence du RR et de la PAS présentent un comportement en bouffées, bien
visible sur la figure 5 (cette figure inclut toutes les périodes, sans limitation aux périodes station-
naires).
Densité spectrale
de puissance 10 LF
[ms2 /Hz]

RR

chute de la composante HF après tilt

1 bouffées de la composante LF avant et après tilt
HF

TILT Temps (mn)

Densité spectrale
1.6
de puissance
[mmHg2 /Hz]

0.1

PAS
augmentation de la composante HF après tilt LF

augmentation de la composante LF après tilt

HF

TILT Temps (mn)

F IG . 5 – Évolution temporelle de l’énergie en basse (LF) et haute fréquence (HF) des signaux rythme
cardiaque (RR) et pression systolique (PAS) avant et après tilt. La composante spectrale en basse
fréquence présente des bouffées, tandis que celle en haute fréquence est plus stable, mais subit en
revanche un grand changement lors du tilt-test. Les spectres sont obtenus par DFT, sur des périodes
de 2 minutes, avec recouvrement de 50 %.

Même si la DFT n’est pas nécessairement le meilleur outil pour étudier le système cardiovascu-
laire, puisque son emploi suppose la stationnarité des signaux (et il est possible que les 5 mn suivant le
changement de posture ne soient pas suffisantes pour atteindre une situation stationnaire), on peut ten-
ter d’interpréter l’instabilité de la composante en basse fréquence qui se manifeste ici par l’apparition
de bouffées comme le reflet de la mise en jeu du régulateur : le SNA.
La figure 6 représente l’évolution temporelle des gains et de l’index de sensibilité proposé. Le gain
en basse fréquence gLF évolue par bouffées et sans différences significatives avant et après tilt. Le gain
en haute fréquence gHF décroı̂t significativement, entraı̂nant une chute de l’index de sensibilité.
Cette étude confirme la nécessité d’analyser conjointement les composantes en haute et basse
fréquence du rythme cardiaque et de la pression artérielle, et surtout de décrire le contrôleur par des
indicateurs liés aux gains et aux points de consigne du système cardiovasculaire.

76
 Gains
normalisés 3

[nu]

1.5
gLF

gHF
0
TILT Temps (mn)
Index de
sensibilité 2.5

[nu] 1
0
0
1
1
0
0
1
1
0 00
11
00
11
0
1
1
0
0
1

1
0
0
1
gHF
gLF
1
0 00
11
0 0
1 1
0 11
1 00
11
00
0
Temps (mn)

F IG . 6 – Évolution temporelle des gains normalisés en basse (gLF ) et en haute (gHF ) fréquence et de
l’index de sensibilité après un stimulus postural. Les spectres sont obtenus par DFT, par périodes de
2 minutes avec recouvrement de 50%.

77
4 Perspectives de la modélisation et mesures en clinique
On a présenté ici deux aspects du travail de modélisation en cours : le muscle cardiaque vu comme
un système dynamique commandé de manière périodique, d’une part, et l’arc baroréflexe, vu comme
une fonction de transfert non linéaire, d’autre part. Beaucoup reste à faire pour ajouter à cet édifice
d’autres sous-systèmes, et enfin pour les intégrer en un modèle simple, dont les paramètres soient
identifiables et utilisables en routine par le clinicien.

4.1 Nouvelles tâches de modélisation

Dans les tâches entreprises, ou à entreprendre, pour obtenir un modèle du système cardio-vasculaire
et respiratoire intégrant les divers sous-systèmes impliqués, on peut relever principalement :
– la modélisation simplifiée de la commande chimique musculaire à partir de la dépolarisation
électrique membranaire des cellules ventriculaires ;
– l’introduction d’un couplage électrique entre cellules pacemaker (du nœud sinusal) et non pa-
cemaker (ventriculaires) pour prendre en compte le délai (variable) de transmission de l’influx
électrique à l’intérieur de la paroi cardiaque (de tels modèles existent, cf. par exemple [12]) ;
– la prise en compte des effets d’inertie du sang dans les compartiments vasculaires, en particulier
dans la circulation veineuse systémique : actuellement, ces compartiments sont conçus comme
des cellules du premier ordre, un ordre 2 peut s’avérer nécessaire, en particulier lors de la mise
en jeu du SCV à l’exercice physique ;
– la modélisation mécanique des barorécepteurs (transducteurs mécano-électriques) sur lesquels
on ne possède à notre connaissance actuellement que des résultats fonctionnels, conduisant à
des modèles du type p + αp0 = k(ν + βν 0 ), où p est la pression artérielle et ν la fréquence de
décharge du nerf efférent [13], mais dans lesquels la base physiologique est questionnable ;
– l’inclusion d’un modèle d’oscillateur respiratoire, et, sur le même principe, d’un modèle de
centre nerveux autonome, et de leurs interactions centrales : de ce point de vue, les travaux de
J. Champagnat et J.-F. Vibert exposés dans le présent volume sont une avancée importante dans
cette voie ;
– à chacun des étages où la régulation par le SNA est mise jeu, l’introduction du contrôle nerveux
comme une commande venant remplacer un paramètre considéré comme fixe dans le sous-
système cible ; en particulier :
– modélisation des effets chronotropes : contrôle par le SNA de la période de la cellule
pacemaker du nœud sinusal (voir p. ex. [14]) ;
– modélisation des effets dromotropes : contrôle par le SNA du temps de conduction élec-
trique intracardiaque (retard variable) ;
– modélisation des effets inotropes : contrôle par le SNA des paramètres de la commande
chimique du muscle cardiaque, en particulier de la concentration en ions calcium intracel-
lulaires ;
– modélisation des muscles circulaires lisses des résistances périphériques et de leur com-
mande électrique par les rameaux sympathiques (modèle vraisemblablement du même
type que celui du muscle cardiaque, avec bien sûr des paramètres différents).

78
4.2 Identification de paramètres par de nouvelles mesures cliniques
Parallèlement à ce travail de modélisation, un travail de traitement et d’exploitation de mesures
est en cours, dans le but immédiat d’identifier des paramètres des sous-modèles du SCV (et du SNA)
considérés. Ces mesures, invasives ou non, sont réalisées par les équipes médicales avec lesquelles
nous collaborons.

4.3 Proposition de nouveaux indicateurs de fonctionnement du SCV et du SNA

L’objectif de ce travail de modélisation, qui s’appuie pour l’identification des paramètres du mo-
dèle sur des campagnes de mesures cliniques, est en effet à terme de fournir aux praticiens des indi-
cateurs de fonctionnement du SCV et de son contrôleur, le SNA :
– des indicateurs de la fonction cardiaque peuvent être obtenus à partir de mesures (a priori inva-
sives) réalisées pendant les phases de contraction et de relaxation isovolumiques, en particulier
par identification de caractéristiques de la commande u et des paramètres du modèle de la fibre
musculaire (cf. ci-dessus l’équation (1) dans “Modèle de la fibre musculaire”) ; dans ces phases,
le sous-système mécanique considéré est en effet isolé et d’identification plus facile ;
– des indicateurs globaux du SCV considéré comme un système dynamique sont attendus des
mesures portant sur les sorties telles que pression artérielle, rythme cardiaque et respiration ;
– enfin des indicateurs portant sur le contrôleur à court terme (le SNA) devraient fournir aux cli-
niciens des renseignements permettant d’évaluer par exemple la stabilité de ce contrôleur, et
éventuellement de diagnostiquer et de corriger certains troubles tels que hypotension orthosta-
tique et syncope vasovagale, et certaines formes d’hypertension ; un gain équivalent a déjà été
proposé (cf. ci-dessus “Modélisation du baroréflexe et mesures chez l’homme”) : il s’agirait ici
notamment d’analyser l’adaptation du point de consigne du contrôleur.

Références
[1] Clairambault, J., Médigue, C., Bestel, J. ( 1997): Le système cardiovasculaire et sa régulation à
court terme par le système nerveux autonome. In : Journées d’Étude “Automatique et Santé” du
club EEA (organisateurs : E. Dombre et A. Fournier). ISIM, Montpellier, juin 1997.
(inclus en annexe à la fin de cet article)
[2] S. Silbernagl and A. Despopoulos. Atlas de poche de physiologie. Flammarion, Médecine-
sciences, 1995.
[3] G. Chauvet. Traité de Physiologie théorique : De la cellule à l’homme. Vol. 2, Masson, 1987.
[4] A. M. Katz. Physiology of the Heart. Raven Press, 1992.
[5] Clairambault, J., Claude, D. Systèmes dynamiques et biologie : quelques exemples. Séminaire
T. I. P. E. sur les systèmes dynamiques. Journée de formation (ENSTA-UPS) pour les professeurs
de Mathématiques Spéciales (organisateur : Max Bezard). ENSTA, Paris, mai 1996.
[6] Bestel, J., Sorine, M. Modélisation mathématique du système cardio-vasculaire. Rapport de Re-
cherche INRIA, à paraı̂tre.
[7] I. Mirsky and W. W. Parmley. CARDIAC MECHANICS: Physiological, Clinical, and Mathe-
matical Considerations. Chap. 4. 1974.

79
 [8] SORIN-BIOMEDICA. Documentation technique sur la sonde BEST.
[9] F. C. Hoppensteadt and C. S. Peskin. Mathematics in Medicine and the Life Sciences. Chap. 5.
Springer Texts in Applied Mathematics, 1992.
[10] Vermeiren, C. Analyse et modélisation du système cardio-vasculaire et sa régulation à court
terme par le système nerveux autonome. Thèse de doctorat en génie biologique et médical,
Université Paris Val de Marne, décembre 1996.
[11] Monti, A. Parameter assessment in a model for the control of cardiovascular functions. Thèse de
fin d’études, Politecnico di Milano, octobre 1998.
[12] Landau, M., Lorente, P., Henry, J., Canu, S. Hysteresis phenomena between periodic and statio-
nary solutions in a model of pacemaker and nonpacemaker coupled cardiac cells. J. Math. Biol.,
25:491-509, 1987.
[13] Taher, M., Cecchini, A.B.P., Allen, M.A., Gobran, S.R., Gorman, R.C., Guthrie, B.L., Lingen-
felter, K.A., Rabbany, S.Y., Rolchigo, P.M., Melbin, J., Noordergraaf, A. Baroreceptor responses
derived from a fundamental concept. Ann. Biomed. Eng., 16: 429-443, 1988.
[14] Clairambault, J. A model of the autonomic control of heart rate at the pacemaker cell level
through G-proteins. 17e IEEE-EMBS, 1379-1380, Montréal, 1995.

80
 Annexe : Le système cardiovasculaire et sa régulation
à court terme par le système nerveux autonome
J. Clairambault, C. Médigue, J. Bestel

Journées Automatique et Santé, Montpellier, mai 1997, Réf. [1] du présent article (1999)

Résumé
On présente le système cardiovasculaire (SCV) dans son ensemble, avec ses compartiments : pompe
cardiaque, petite et grande circulation, circuit artériel et circuit veineux, les variables observables en cli-
nique, les mécanismes de contrôle, dont le système nerveux autonome (SNA), et les interactions du SCV
avec la respiration. Les problèmes de modélisation sont abordés sous l’angle de la commande, c’est-à-
dire en fonction d’une mission à remplir pour le SCV et son contrôleur. Les possibilités de commande
pharmacologique du SNA sont également examinées, et des applications à la clinique envisagées.

1 Présentation du système cardiovasculaire (SCV)

La principale mission du SCV est la fourniture d’oxygène et de nutriments aux différents tissus
de l’organisme, en respectant une hiérarchie dans les distributions régionales : cerveau, puis rein, ter-
ritoire splanchnique (celui de la digestion...), et membres. La sélection d’un territoire aux dépens des
autres se fait par vasoconstriction (réduction du calibre des vaisseaux), dans les territoires négligés.

1.1 Moyens mis en œuvre

Les moyens mis en œuvre pour assurer cette mission sont une pompe : le cœur, et un circuit
hydraulique qu’on peut considérer, du moins à court et même à moyen terme (c’est-à-dire indépen-
damment de la réabsorption rénale), à volume constant : le système vasculaire.

1.1.1 Le cœur en tant que pompe

La pompe cardiaque est un générateur de pression, ou plus exactement de débit (la pression étant
le produit du débit par la résistance opposée au travail de la pompe : si on assimile la pression à un
voltage et le débit à une intensité de courant, il ne s’agit ni plus ni moins que de la loi d’Ohm).
Cette pompe est en fait elle-même séparée de manière étanche (séparation achevée à la naissance)
en deux pompes, l’une à haute pression vers l’ensemble de l’organisme (le ventricule gauche), et
l’autre à basse pression vers le poumon (le ventricule droit). L’égalisation des débits droit et gauche
est réalisée par un mécanisme d’adaptation purement mécanique : l’effet Starling (voir ci-dessous :
cahier des charges physiologiques). Cette double pompe a une fréquence et un volume d’éjection
(les deux facteurs du débit cardiaque) variables ; ces deux facteurs sont en permanence contrôlés, en
particulier par le système nerveux autonome (SNA).

1.1.2 Pacemaker et conduction électrique

La fréquence cardiaque est le principal facteur d’adaptation du débit. Elle est donnée par un pa-
cemaker naturel, situé dans l’oreillette droite : le nœud sinusal, qui est soumis au contrôle du SNA.

81
À partir du nœud sinusal l’influx électrique se propage (à vitesse variable, et cette vitesse est aussi
soumise au contrôle du SNA) vers le myocarde ventriculaire, lieu du couplage excitation-contraction.

1.1.3 Volume d’éjection

Le volume d’éjection dépend de 3 facteurs : a/ la précharge, ou volume de sang présent dans le
ventricule gauche juste avant l’éjection (volume proportionnel à la pression du retour veineux au cœur,
dite pression veineuse centrale, PVC) ; b/ la contractilité du myocarde, sous la dépendance du SNA ;
c/ la postcharge, ensemble des facteurs qui s’opposent au travail du cœur : résistance pariétale du
ventricule, impédance de l’aorte, et résistances périphériques.

1.1.4 Les deux circulations

Le circuit vasculaire est divisé en deux parties : la circulation dite systémique, ou grande circula-
tion, qui apporte oxygène et nutriments aux tissus, et la circulation pulmonaire, ou petite circulation,
qui assure la navette entre le cœur, générateur de débit, et le poumon, générateur d’oxygène. L’oxy-
gène est la matière première de la chaı̂ne de phosphorylation oxydative de la mitochondrie, centrale
energétique de la cellule, et on peut définir la puissance consommée par un organe comme le produit
du débit sanguin le traversant par la différence artérioveineuse, différence des concentrations en oxy-
gène dans le sang artériel (à l’entrée de l’organe, à haute pression) et dans le sang veineux (à sa sortie,
à basse pression). Le maintien d’une pression artérielle efficace, permettant à la circulation d’assurer
une bonne oxygénation des tissus, est assuré par deux mécanismes complémentaires : la régulation du
débit cardiaque et celle des résistances périphériques.

1.1.5 Résistances périphériques

L’ensemble des résistances périphériques, encore appelées shunts artérioveineux, constitue la ré-
sistance systémique, qu’on peut définir par la loi d’Ohm comme le quotient d’une différence de pres-
sion artérioveineuse globale (mesurée entre la sortie du cœur au niveau des valvules sigmoı̈des aor-
tiques et le retour veineux au niveau de l’abouchement de la veine cave) par le débit cardiaque. Ces
résistances périphériques sont soumises principalement à deux types de régulations. L’une, dite auto-
régulation, est purement locale et a pour but de maintenir une fourniture en oxygène, c’est-à-dire un
débit artériel local, efficace même lorsque la pression chute ; un modèle simple consiste à rendre cette
résistance tissulaire locale proportionnelle (grâce à des chémorécepteurs) à la concentration veineuse
en oxygène. L’autre est sous la dépendance du même contrôleur que le débit cardiaque (le SNA).

1.1.6 Compliance vasculaire

Le circuit vasculaire n’est pas seulement pourvu de résistances, il est aussi pourvu de capacités.
Les vaisseaux sanguins ne sont en effet pas rigides, mais compliants, et une augmentation de pression
se traduira par une augmentation de leur calibre. Les veines systémiques, notamment, constituent un
important réservoir de sang (environ 3 litres, contre moins d’un litre pour les artères systémiques),
et le circuit veineux, circuit à basse pression, est aussi appelé circuit capacitif, par opposition au
circuit artériel, à haute pression, dit circuit résistif. Mais les artères présentent elles aussi une certaine
compliance ; en particulier la compliance de l’aorte se manifeste par un aspect bifide (“dicrote” : à

82
deux sommets) de l’onde de pression aortique, due à la réflexion de l’onde de propagation sur les
résistances périphériques, lorsque les valvules sigmoı̈des sont fermées.

1.2 Critères à optimiser par la commande du SCV

Il est légitime d’imposer au contrôleur de minimiser le travail de la pompe cardiaque, c’est-à-dire
essentiellement celui du ventricule gauche, la plus puissante des deux pompes, celle qui débite à haute
pression. Mais les deux pompes sont en fait synchrones, et le débit de la pompe droite est égal à celui
de la pompe gauche ; on peut donc définir ce travail comme une valeur moyenne au cours du temps du
produit du débit cardiaque par la différence de pression artérioveineuse globale (entre la veine cave à
l’entrée et l’aorte à la sortie).
La principale variable régulée est la pression artérielle. Normalement comprise entre 70 (valeur
minimale = diastolique) et 140 mm Hg (maximale = systolique), elle maintient un débit cérébral (le
territoire cérébral est prioritaire, son approvisionnement se faisant éventuellement aux dépens des
autres) satisfaisant jusqu’à environ 60 mm Hg de pression systolique, après quoi le débit s’effondre.

1.3 Mécanismes de contrôle de la pression artérielle

1.3.1 La régulation à long terme est réalisée :
– soit par modification du volume plasmatique au niveau rénal ; il s’agit d’une régulation pour
l’essentiel hormonale : système rénine-angiotensine-aldostérone, ADH (hormone antidiurétique,
alias vasopressine), facteur atrial natriurétique, dont la cible est rénale, mais aussi également
d’une régulation passive, conséquence directe des variations de pression artérielle dans le rein ;
– soit par modification des résistances périphériques : régulation hormonale par action vasocons-
trictrice directe de l’angiotensine II.

1.3.2 La régulation à moyen terme

C’est une régulation hormonale par les catécholamines circulantes (noradrénaline surtout), qui
sont sécrétées par la médullosurrénale sous l’influence de centres cérébraux (zones pressogènes),
activés par la chute du débit sanguin local (ischémie cérébrale), mais aussi en cas de stress, et ont une
action vasoconstrictrice qui peut atteindre jusqu’à 6 fois celle obtenue par le contrôle nerveux.

1.3.3 La régulation à court terme

C’est une régulation uniquement nerveuse, assurée par une partie du système nerveux qui fonc-
tionne à l’état basal de manière indépendante de la vie de relation : le système nerveux végétatif, ou
système nerveux autonome (SNA).
Les centres du SNA reçoivent des influences externes en provenance de l’hypothalamus (tempé-
rature), et, au-delà, du système limbique ou du cortex (émotions, à effet parfois quasi-immédiat sur le
rythme cardiaque). De plus, pendant le sommeil paradoxal, la substance réticulée ascendante activa-
trice du mésencéphale exerce sur ces centres d’importantes perturbations. On se contentera de décrire
ci-dessous le fonctionnement du SNA à l’état de base, ce qui suppose un sujet vigile, ou en état de
sommeil profond, et aussi déconnecté que possible de perturbations extérieures éventuelles.

83
2 Le contrôleur à court terme : le système nerveux autonome
2.1 Capteurs et voies afférentes
Des barorécepteurs, sensibles à la distension, sont situés dans la paroi des artères au niveau de
la crosse de l’aorte et du sinus carotidien (élargissement de la lumière de l’artère carotide, situé à la
hauteur de la bifurcation entre carotides externe et interne), et aussi sur le circuit veineux de retour,
au niveau de l’oreillette droite. Les barorécepteurs artériels sont sensibles à des variations de pression
minimes, de l’ordre de quelques mm Hg. Ces barorécepteurs transmettent des messages nerveux,
sous la forme de trains d’impulsions électriques, codés en fréquence, par des fibres afférentes (= de
la périphérie vers le centre) qui rejoignent les IXe et Xe paires crâniennes pour se terminer dans le
noyau du tractus solitaire (NTS, sous le plancher du 4e ventricule, dans le tronc cérébral).
Les barorécepteurs de l’oreillette droite, ainsi que des mécanorécepteurs pulmonaires stimulés
comme eux pendant l’inspiration, envoient des messages au NTS qui ont pour effet d’adapter le SNA
à la pression veineuse de retour et aux variations de la pression intrathoracique dues à la respiration.

2.2 Centres de contrôle

À l’intérieur du NTS, des neurones secondaires déchargent suivant une fréquence nettement plus
basse que ceux provenant directement des barorécepteurs, et répondent par une augmentation de leur
fréquence de décharge à une augmentation de pression au niveau des barorécepteurs. Cette dépen-
dance peut être modélisée comme intégrant à la fois une pression moyenne et la dérivée de la pression.
C’est vraisemblablement dans le NTS qu’on peut trouver la première non-linéarité de l’arc baro-
réflexe artériel (= la boucle de contrôle de la pression artérielle par les barorécepteurs et le système
nerveux autonome). Cette non-linéarité, placée entre une entrée constituée des trains d’impulsions
électriques en provenance des barorécepteurs, et une sortie à destination du noyau ambigu et de l’aire
pressive bulbaire, est assimilable à un relais tout-ou-rien, sans doute avec hystérésis.
Le NTS est relié à deux centres bulbaires, d’une part au noyau ambigu, d’où partent des effé-
rences parasympathiques (= vagales) à destination du cœur, d’autre part à l’aire pressive bulbaire,
d’où partent des efférences sympathiques à destination du cœur et des vaisseaux. Il existe des inhibi-
tions intercentrales réciproques, qui permettent un couplage dès l’origine de ces deux voies.

2.3 Voies efférentes : la voie (ortho)sympathique et la voie parasympathique

Les deux branches efférentes du SNA ont des effets antagonistes : une stimulation sympathique
augmente la fréquence cardiaque et la contractilité des fibres myocardiques, augmente aussi les ré-
sistances périphériques, tandis qu’une stimulation parasympathique, principalement à destination du
pacemaker sinusal, ralentit le cœur.
Ces voies ont des temps d’action très différents : la voie sympathique, cardioaccélatrice, est effi-
cace en en une dizaine de battements cardiaques, tandis que le délai d’action du coup de frein vagal
transmis par la voie parasympathique est quasi instantané, de l’ordre du battement. En dehors d’une
raison téléologique immédiatement compréhensible par tout conducteur de véhicule à moteur, il y a à
ce fait plusieurs explications d’ordre anatomique : la voie parasympathique part directement du centre
encéphalique pour ne faire relais que dans la paroi même du cœur (le neurone postganglionnaire est
ainsi très court), tandis que la voie sympathique fait relais très tôt, dans les ganglions sympathiques

84
paravertébraux, donc loin du cœur. Or les fibres préganglionnaires sont myélinisées, et transmettent
les impulsions électriques beaucoup plus rapidement que les fibres postganglionnaires, qui ne le sont
pas ou peu. D’autre part, la synapse ganglionnaire sympathique est à simple entrée pour des sorties
multiples (environ 1 pour 20, d’où un codage en fréquence, mais aussi en amplitude si on considère
la sommation des potentiels d’action sur plusieurs voies parallèles), contre un codage (1 pour 1), en
fréquence pure pour la voie parasympathique, d’où une possibilité de perte en ligne sur la voie sym-
pathique. Enfin, la synapse neuro-effectrice cardiaque a un neuromédiateur différent dans les deux
voies : acétylcholine pour la parasympathique, noradrénaline pour la sympathique, et la synapse no-
radrénergique serait plus lente que la synapse acétylcholinergique.
Il résulte de ces considérations anatomiques qu’on peut considérer la voie sympathique comme
une voie à retard par rapport à la voie parasympathique.

2.4 Organes effecteurs du SNA

2.4.1 Le pacemaker cardiaque, principal effecteur du SNA
Le nœud sinusal est constitué de plusieurs milliers de cellules qui sont synchronisées, mais sont
toutes douées d’automatisme. Beaucoup de modèles à base d’équations différentielles (non linéaires)
en ont été proposés, qui tous reposent sur des variantes du modèle de Hodgkin et Huxley (1952) pour
le potentiel d’action de la cellule nerveuse (on en trouvera une brève revue dans [6]). Ces modèles
ont l’inconvénient d’être gourmands en nombre de variables, rendant compte de courants ioniques
différents, et de ne pas être observables facilement à l’échelle d’un individu.
Le mode d’action des deux voies antagonistes du SNA au niveau de la cellule pacemaker car-
diaque est assez bien connu ; la modulation par la voie sympathique (neuromédiateur : noradrénaline)
se traduit essentiellement par un redressement de la pente (sensiblement constante) de la phase de
dépolarisation diastolique lente du potentiel d’action de cette cellule, réduisant la durée du cycle car-
diaque ; la modulation parasympathique (neuromédiateur : acétylcholine) a l’effet opposé. Des inhibi-
tions réciproques mettent en place un couplage local sympathique / parasympathique au niveau de la
jonction neuroeffectrice, qui vient compléter le couplage central. Cet effet du SNA est dit chronotrope.
Il existe aussi un effet dit bathmotrope, d’action sur l’excitabilité des cellules du tissu de conduction
par variation de la période réfractaire, mais cet effet semble peu important au niveau macroscopique.
Le nœud sinusal bat, hors de toute influence du SNA, à un rythme propre d’environ 105 battements
par minute (bpm) ; mais il est en permanence soumis à des tonus, sympathique et parasympathique,
qui ramènent ce rythme autour de 70 bpm (il y a normalement un tonus vagal dominant).

2.4.2 La conduction auriculoventriculaire

Les cellules du nœud sinusal ne sont pas les seules cellules cardiaques qui soient douées d’automa-
tisme, mais ce sont les plus rapides, et elles imposent leur rythme au cœur. L’influx électrique né dans
le nœud sinusal se propage dans l’oreillette droite, puis dans la paroi du septum interventriculaire,
pour gagner par le faisceau de His les cellules contractiles du myocarde ventriculaire ; la dépolarisa-
tion de l’oreillette se produit bien avant celle du ventricule et n’a que très peu d’effets mécaniques. La
vitesse de conduction auriculoventriculaire (appréciable sur l’électrocardiogramme de surface par la
durée de l’intervalle PR), déterminée par des résistances (électriques) variables entre cellules du tissu
de conduction, est aussi sous la dépendance du SNA : une stimulation sympathique augmente cette
vitesse, tandis qu’une stimulation parasympathique la diminue. C’est l’effet dit dromotrope du SNA.

85
2.4.3 La contractilité du myocarde
Le couplage excitation-contraction dans la cellule ventriculaire est encore mal connu, mais le
principal facteur de la puissance de la pompe cardiaque est la vitesse de raccourcissement des fibres
contractiles. Là encore le SNA joue un rôle : le sympathique augmente la contractilité (c’est l’effet
inotrope positif du sympathique) ; en revanche il y a peu, ou pas, de terminaisons nerveuses parasym-
pathiques ventriculaires, et donc peu d’effet inotrope direct du parasympathique.

2.4.4 La tension et la relaxation du myocarde

De même, le sympathique augmente la tension au repos et la relaxation (= la distensibilité, i.e. la
longueur de la fibre au repos) de la fibre myocardique, et augmente par là la puissance de la pompe.
Ces effets sont dits tonotrope (tension), et lusitrope (relaxation). Ils augmentent l’effet Starling (dit
aussi “loi fondamentale du cœur”), indépendant par ailleurs du SNA, qui permet au ventricule de se
contracter d’autant plus efficacement qu’il est plus rempli, l’élongation de la fibre myocardique au
repos étant plus important.

2.4.5 Les résistances périphériques

Elles sont sous l’influence du sympathique seul, principalement par l’intermédiaire de récepteurs
noradrénergiques de type α (à la différence des récepteurs cardiaques, qui sont surtout de type β [2]).
L’effet du sympathique sur les résistances vasculaires est une vasoconstriction qui se traduit par une
augmentation des résistances. Cet effet est sans doute celui qui est mis en place le plus rapidement
dans le passage à l’orthostatisme (= la position debout), permettant de maintenir une pression artérielle
efficace avant même que l’effet de l’augmentation du débit cardiaque ait pu se faire sentir.

2.5 Mise en jeu physiologique et physiopathologique du SNA

L’analyse spectrale du rythme cardiaque (identifié à la suite des intervalles de temps entre batte-
ments) et de la pression artérielle systolique, chez l’homme comme chez l’animal, met en évidence
deux pics principaux : un pic, autour de 0.25 Hz, correspond à la fréquence respiratoire, et est parti-
culièrement net chez le sujet couché, en particulier en sommeil profond, et un autre, autour de 0.1 Hz,
représentatif des “ondes de Mayer”, est un rythme endogène (on ne lui a pas trouvé d’origine phy-
siologique extérieure au SCV) de basse fréquence bien visible en particulier sur un enregistrement
effectué après passage à l’orthostatisme. Ces pics disparaissent en cas de blocage pharmacologique
(électif : injection de β-bloquants pour la voie sympathique, d’atropine pour la voie parasympathique)
des récepteurs du SNA. Ces modulations visibles sur un enregistrement de rythme cardiaque sont ac-
tuellement en clinique le principal moyen non invasif d’accès au SNA.

2.5.1 L’arythmie sinusale respiratoire

Elle n’a rien de pathologique, et est au contraire un élément du bon fonctionnement du SNA.
Elle est définie par une augmentation de la fréquence cardiaque à l’inspiration, et une diminution à
l’expiration. Il existe à ce mécanisme essentiellement deux origines, centrale et mécanique.
L’explication centrale est une interaction directe entre les centres autonomes qui contrôlent le
cœur et le centre respiratoire bulbaire : le noyau dorsal du vague (alias cardiopneumoentérique).

86
 L’explication mécanique est fortement liée à la modification de la pression veineuse de retour due
à l’inspiration : l’abaissement de la coupole diaphragmatique induit un pompage de sang abdominal
dans la cage thoracique, avec distension de mécanorécepteurs (barorécepteurs du circuit veineux)
situés dans l’oreillette droite, qui transmettent par la voie du nerf vague un message aux centres
autonomes, provoquant une inhibition parasympathique, ce qui revient à une stimulation sympathique
immédiate (réflexe de Bainbridge). Il faut y ajouter une stimulation parallèle de mécanorécepteurs
pulmonaires à l’inspiration, qui produit les mêmes effets.
Enfin une explication mécanique locale directe, par étirement des fibres du nœud sinusal à l’inspi-
ration (donc indépendamment du SNA) existe aussi, mais son importance est sans doute négligeable.

2.5.2 Le passage à l’orthostatisme

Au cours du passage à l’orthostatisme, il y a chute brutale de sang veineux vers les membres
inférieurs, diminution de la pression veineuse de retour, diminution du débit cardiaque et donc de
la pression artérielle au niveau des barorécepteurs artériels. Cette stimulation du baroréflexe artériel
entraı̂ne une élévation de la fréquence cardiaque, sans effets immédiats néanmoins, car la chute du
retour veineux ne permet pas au débit cardiaque d’augmenter assez rapidement. C’est pourquoi la
compensation de la pression artérielle lors du passage à l’orthostatisme se fait essentiellement par
l’intermédiaire de la vasoconstriction périphérique, par une inhibition parasympathique consécutive à
la chute de pression au niveau des barorécepteurs.

2.5.3 L’exercice physique

Les mêmes phénomènes que lors du passage à l’orthostatisme existent ; il s’y ajoute aussi l’inter-
vention d’un chémoréflexe musculaire (sensible aux pressions partielles en oxygène et en gaz carbo-
nique, et un déplacement de la courbe stimulus-réponse du baroréflexe (réajustement d’une grandeur
de consigne sur la pression artérielle) [15].
Mais chez les sujets sportifs entraı̂nés, le tonus vagal est tel qu’il n’y a que très peu de réponse du
sympathique à une chute de pression ; au contraire les cas d’hypotension orthostatique à l’arrêt trop
brusque de l’effort sont fréquents.

2.5.4 L’insuffisance cardiaque

Dans l’insuffisance cardiaque, en revanche, le débit cardiaque étant insuffisant, les résistances
périphériques sont stimulées par la voie sympathique au maximum de leurs capacités, avec (régulation
à moyen terme, hormonale) des taux très élevés de catécholamines circulantes.

2.6 Commande pharmacologique du SNA

Les médicaments du SNA ont une action soit centrale, soit périphérique. Ils stimulent ou bloquent
les récepteurs, α ou β, à la noradrénaline, ainsi que les récepteurs, nicotiniques ou muscariniques
(les récepteurs cardiaques sont muscariniques) à l’acétylcholine. Le plus connu des parasympatholy-
tiques est l’atropine. On les utilise, expérimentalement ou en clinique, pour bloquer électivement une
branche ou l’autre du SNA. Il existe aussi des α- et β-mimétiques, et des parasympathomimétiques.

87
3 Principes de la modélisation du SCV sous le contrôle du SNA
3.1 Problèmes généraux
La construction d’un modèle du SCV sous la dépendance du SNA se heurte à des contraintes très
fortes. Les variables mesurables (autant que possible de manière non invasive, c’est-à-dire par des
instruments de mesure externes) sont bien peu nombreuses : pression artérielle en continu par bras-
sard avec contre-pression (FinapresTM ), électrocardiogramme, respiration, électroencéphalogramme
dans l’étude du sommeil. La pression veineuse centrale est par définition invasive, mais il peut être
nécessaire d’en disposer sur des enregistrements animaux, par exemple, pour estimer le retour vei-
neux. Quelles variables d’état considérer ? Une variable d’“état sympathique” et une autre, d’“état
parasympathique” sont de bons candidats pour décrire l’état interne du contrôleur.
Les conflits entre description microscopique (la cellule) et macroscopique (l’organisme) dans une
modélisation à base physiologique sont représentatifs de l’opposition entre ce qu’on sait de la phy-
siologie, souvent dans des détails très fins au niveau cellulaire, et ce qu’on peut en observer au niveau
intégré du corps entier.
En dehors de l’observabilité, des problèmes d’identifiabilité sont à prévoir, et la construction d’un
modèle doit en tenir compte, si on veut pouvoir s’en servir en clinique ; mais suivant le problème
physiopathologique considéré, tel ou tel aspect du mécanisme de commande peut être privilégié, aux
dépens des autres, afin de rendre un modèle identifiable.
Enfin, signalons que l’intérêt des médecins et des physiologistes, dans notre pays, pour la modéli-
sation mathématique est très variable ; si certains y sont réfractaires, d’autres sont prêts à monter des
expérimentations animales ou humaines pour tester une hypothèse.

3.2 Cahier des charges physiologiques et applications à la modélisation

3.2.1 Pompe cardiaque

La “loi fondamentale du cœur” (l’effet Starling) est que toute augmentation de la précharge (vo-
lume, ou pression au repos, avant la contraction) augmente la force de contraction du ventricule, non
par augmentation de la vitesse de contraction (contractilité), mais par accroissement de la longueur
de la fibre myocardique au repos. Ce mécanisme a pour effet d’égaliser les débits droit et gauche
et d’éviter toute stase, et tout pompage inutile, dans la circulation pulmonaire. En revanche, toute
augmentation de la postcharge diminue la contractilité du myocarde.
Un effet positif direct d’une augmentation de la fréquence cardiaque sur la contractilité, dit effet
(de l’escalier de) Bowditch, est aussi à prendre en compte. Il est totalement indépendant et de l’effet
Starling et du SNA, et est lié au largage intracellulaire de calcium : plus le cœur bat vite, et plus la
fraction de temps pendant laquelle les canaux calciques membranaires sont ouverts est élevée ; or le
calcium intracellulaire est un déterminant essentiel de la contractilité musculaire.
L’intervention du SNA doit être modélisée à tous les niveaux du cœur. Sur le pacemaker (effet
chronotrope), ce peut être par un oscillateur de type Hodgkin-Huxley à paramètres dépendant du SNA
[18, 5], ou plus simplement par un oscillateur harmonique de fréquence dépendant directement du
SNA [21]. Les effets bathmotrope (excitabilité générale du myocarde) et dromotrope (sorte d’“avance
à l’allumage”) peuvent être négligés dans un premier temps. L’effet inotrope est aussi essentiel à
prendre en compte, la puissance de la pompe dépendant fortement d’une stimulation sympathique.

88
3.2.2 Circulation artérielle systémique
L’impédance et la compliance des grosses artères doivent aussi être prises en compte dans la
modélisation de la pression artérielle systémique, si on veut la mettre en rapport de manière réaliste
avec le débit cardiaque et la résistance systémique.

3.2.3 Baroréflexe
La modélisation des non-linéarités est un aspect crucial. La voie sympathique est vraisemblable-
ment une voie à retard par rapport à la voie parasympathique, et ce seul retard autorise l’apparition
du rythme endogène à basse fréquence ; mais il faut noter que ce n’est pas la seule manière d’obtenir
ce résultat. En particulier, modéliser le baroréflexe central comme un relais avec hystérésis (le dé-
clenchement, l’“onset”, du réflexe se produisant pour un point de consigne en pression plus élevé que
son arrêt, l’“offset”) permet, ainsi que le montre une analyse par la fonction de description, d’obtenir
une stabilisation de la pression au prix de l’apparition d’un cycle limite (le “0.1 Hz”) sans qu’il soit
nécessaire d’introduire un retard. Ce rythme endogène de basse fréquence est-il transitoire ou non ?
On l’a décrit après passage à l’orthostatisme, et aussi dans certaines situations qui peuvent être liées à
des perturbations passagères du baroréflexe. Si un tel phénomène hystérétique est avéré (il ne semble
en particulier pas exister chez le chien au repos [11], et de même le 0.1 Hz est souvent totalement
absent chez l’homme en position couchée au repos), sa quantification (par exemple par la largeur de
la fenêtre d’hystérésis) peut être un élément précieux pour évaluer la sensibilité du baroréflexe.

3.2.4 Résistances périphériques

Elles doivent être décrites en incluant et l’autorégulation (qui, isolée, est dangereuse, puisqu’elle
produit une ouverture des shunts périphériques en cas de diminution de la pression artérielle), et leur
contrôle par le SNA (effets α-adrénergiques).

3.2.5 Circulation veineuse systémique

La modélisation doit être celle d’une circulation à très basse pression, aidée par a/ une succion par
la pompe cardiaque (“vis a fronte”) ; b/ une pression résiduelle ayant passé la barrière artérioveineuse
(“vis a tergo”) ; c/ une pression régionale (muscles) relayée par des valvules veineuses (“vis a latere”).
Les interactions mécaniques régionales avec les mouvements respiratoires (veine cave supérieure et
oreillette droite) doivent prendre en compte la très grande compliance des veines, et leur fonction
fondamentale de réservoir de sang.

3.2.6 Petite circulation et mouvements respiratoires

La petite circulation représente moins de 10 % du volume sanguin circulant ; la réserve de sang
pulmonaire est d’environ 6 systoles (6 fois le volume d’éjection moyen), et cette marge permet au
SCV de s’adapter dans l’urgence à un effondrement de la pression veineuse de retour.
Si on prend en compte les concentrations artérielle et veineuse en oxygène, leur dépendance de la
fréquence et de l’amplitude des mouvements respiratoires est à intégrer au modèle ; mais il faut surtout
rendre compte de l’arythmie sinusale respiratoire, du moins dans sa partie mécanique en rendant la
pression veineuse centrale dépendante de ces mouvements.

89
3.2.7 États physiopathologiques complexes
Il s’agit de l’exercice physique extrême (chémoréflexes : [O2 ], [CO2 ], pH), de pathologies impli-
quant des perturbations importantes du SNA ([Link]. diabète avec neuropathie), et du sommeil para-
doxal. Le sommeil paradoxal est un état physiologique du système nerveux central au cours duquel
des interactions très complexes peuvent survenir entre notamment la réticulée et les centres cardiores-
piratoires du bulbe, avec déconnection concomitante du contrôle par l’hypothalamus et le cortex, mais
surtout une interaction moindre que dans le sommeil profond entre centre respiratoire et noyau am-
bigu (donc une diminution de l’arythmie sinusale respiratoire) [1]. C’est pourquoi une identification
des paramètres du SNA doit prendre en compte les états de vigilance et de comportement.

3.3 État de l’art

Une description de base du SCV peut être trouvée dans [16]. Un point de vue plus général sur la
physiologie intégrative peut être trouvé dans [4], et le petit atlas de Silbernagl et Despopoulos [19]
est une référence précieuse pour la description des grandes régulations physiologiques. Toujours du
côté de la physiologie, la revue de Melchior [13] traite de la réponse à l’orthostatisme. Eckberg [7]
a mené de nombreuses expériences décrivant le fonctionnement du baroréflexe ; Taher [20] traite en
particulier du problème de l’hystérésis dans le baroréflexe. Le livre de Levy et Schwartz [12] fait
un point récent complet sur le contrôle vagal, et aussi accessoirement sur le contrôle sympathique.
Schwaber [17] étudie le fonctionnement des mécanismes centraux du baroréflexe (au niveau du NTS).
Les interactions entre mathématiciens, automaticiens, et physiologistes ne datent pas d’hier. En
particulier, Guyton [8] a développé dès les années 60 un modèle analogique complet du SCV, incluant
aussi le rein (ce qui est nécessaire pour étudier l’hypertension artérielle) ; de même, Milhorn [14]
propose des principes de modélisation par schémas-blocs du SCV, et pas seulement par des outils
linéaires. Kitney [10] a dès 1979 proposé l’analyse du SCV par la fonction de description. Plus ré-
cemment, à l’aide de modèles à base d’équations différentielles, Vermeiren [21], par la méthode de la
fonction de description, Seidel [18], et Cavalcanti [3], par une analyse de continuation-bifurcation sui-
vant un paramètre “retard sur la voie du baroréflexe”, ont rendu compte du rythme endogène de basse
fréquence. Le débit cardiaque et la pression artérielle, avec la résistance systémique, sont étudiés dans
[22]. Un modèle du SCV récent et complet, mais non pulsatile, est celui de Kappel [9].

3.4 Diverses pathologies et disciplines médicales concernées

En cardiologie, il a été montré que sur l’ensemble d’une population, un mauvais état du SNA
augmentait les risques de mort subite chez les survivants après infarctus du myocarde.
En neurologie, on étudie des pathologies à l’origine d’une atteinte du SNA (p. ex. neuropathie
diabétique) ou liées à un dysfonctionnement du SNA (syndrome des apnées du sommeil d’origine
centrale, syndrome des myoclonies nocturnes, notamment).
En pédiatrie, un défaut de maturation du SNA a été proposé comme facteur de risque de mort
subite inexpliquée du nourrisson.
En obstétrique, il est connu de longue date qu’une chute de la fréquence cardiaque fœtale suivie
d’un plateau (perte de variabilité) est un signe de souffrance fœtale aiguë.
Ces disciplines, ainsi que d’autres (anesthésie-réanimation pour surveillance-diagnostic au bloc
opératoire ou en unité de soins intensifs) pourraient bénéficier d’un modèle du SCV et de son contrôle

90
par le SNA qui leur donnerait accès aux paramètres internes du SCV.

Références
[1] Benoit, O., Foret, J. Le sommeil humain, bases expérimentales et physiopathologiques. Masson,
Paris, 1992.
[2] Bioulac, B. La neurotransmission dans le système nerveux végétatif. La revue du praticien, 31
(37): 2605-2631, 1981
[3] Cavalcanti, S., Belardinelli, E. Modeling of Cardiovascular Variability Using a Differential
Equation. IEEE-BME, 34 (10): 982-988, 1996.
[4] Chauvet, G. Traité de physiologie théorique (3 tomes). Masson, Paris, 1987.
[5] Clairambault, J. A model of the autonomic control of heart rate at the pacemaker cell level
through G-proteins. Actes de la 17e conférence internationale IEEE-EMBS, 1379-1380, Mont-
réal, septembre 1995.
[6] Clairambault, J., Claude, D. Systèmes dynamiques et biologie : quelques exemples. Séminaire
T. I. P. E. sur les systèmes dynamiques. Journée de formation (ENSTA-UPS) pour les professeurs
de Mathématiques Spéciales (organisateur : Max Bezard). ENSTA, Paris, mai 1996.
[7] Eckberg, D.L., Rea, R.F., Andersson, O.K., Hedner, T., Prnow, J., Lundberg, J.M., Wallin, B.G.
Baroreflex modulation of sympathetic activity and sympathetic neurotransmitters in humans.
Acta Physiol. Scand., 133: 221-231, 1988.
[8] Guyton. A.C. Arterial Pressure and Hypertension. Saunders, Philadelphie, 1980.
[9] Kappel, F., Peer, R.O. A mathematical model for fundamental regulation processes in the car-
diovascular system. J. Math. Biol., 31: 611-631, 1993.
[10] Kitney, R.I. A nonlinear model for studying oscillations in the blood pressure control system. J.
Biomed. Eng., 1 (2): 89-98, 1979.
[11] Levison, W.H., Barnett, G.O., Jackson, W.D. Nonlinear analysis of the baroreceptor reflex sys-
tem. Circ. Res., 18: 673-682, 1966.
[12] Levy, M.N., Scwartz, P.J., Eds. Vagal control of the heart. Futura, New-York, 1994
[13] Melchior, F.M., Sirnivasan, R.S., Charles, J.B. Mathematical modeling of human cardiovascular
system for simulation of orthostatic response. Am. J. Physiol., 262: H1920-H1933, 1992.
[14] Milhorn, H.T. The application of control theory to physiological systems. Saunders, Philadel-
phie, 1966.
[15] Papelier, Y. Adaptation cardiovasculaire à l’exercice physique chez le sportif : rôle de l’interac-
tion entre baroréflexes et chémoréflexe musculaire. Thèse de doctorat en sciences et techniques
des activités physiques et sportives, Université Paris Sud, juin 1994.
[16] Hoppensteadt, F.C., Peskin, C.S. Mathematics in Medicine and the Life Sciences. Springer Texts
in Applied Mathematics, 1992.
[17] Schwaber, J.S., Graves, E.B., Paton, J.F.R. Computational modeling of neuronal dynamics for
systems analysis: application to neurons of the cardiorespiratory NTS in the rat. Brain Res., 604:
126-141, 1993.
[18] Seidel, H., Herzel, H. Modelling Heart Rate Variability Due to Respiration and Baroreflex. In :
E. Mosekilde, O.G. Mouritsen, Eds. Modelling the Dynamics of Biological Systems. Springer
Series in Synergetics, 1995.

91
[19] SilBernagl, S., Despopoulos, A. Atlas de poche de physiologie. Médecine-Sciences, Flamma-
rion, Paris, (4e de la 2e édition française) 1995.
[20] Taher, M., Cecchini, A.B.P., Allen, M.A., Gobran, S.R., Gorman, R.C., Guthrie, B.L., Lingen-
felter, K.A., Rabbany, S.Y., Rolchigo, P.M., Melbin, J., Noordergraaf, A. Baroreceptor responses
derived from a fundamental concept. Ann. Biomed. Eng., 16: 429-443, 1988.
[21] Vermeiren, C. Analyse et modélisation du système cardio-vasculaire et sa régulation à court
terme par le système nerveux autonome. Thèse de doctorat en génie biologique et médical,
Université Paris Val de Marne, décembre 1996.
[22] Wesseling, K.H., Jansen, J.R.C., Settels, J.J., Schreuder, J.J. The computation of aortic flow from
pressure in humans using a nonlinear, three-element model. J. Appl. Physiol., 74 (5): 2566-2573,
1993.

92