Empirical Study on Linear, Non-linear, and System Dynamics: Oriented Human Motor Behavior

Chulwook Park^a^,b^,c; Dong-Wook Han^d^*
^aResearch Fellow, Institute of Sport Science, Seoul National University, Seoul,^bResearch Fellow, International Institute for Applied Systems Analysis (IIASA), Laxenburg,^cResearch Fellow, Okinawa Institute of Science and Technology (OIST), Okinawa,^dProfessor, Department of Sport Science, Jeonbuk National University, Jeonju-si,Correspondence: handw@jbnu.eduInternational Journal of Applied Sports Sciences, Vol.34, No.2, pp.160-182, 31 December 2022
https://doi.org/10.24985/ijass.2022.34.2.160

초록

Human movement and motor behavior science involve various theoretical frameworks and methodologies. This study describes the force control of motor phenomenon, processes of linear or nonlinear dynamics, and proposes a system approach as an additional potential aspect. Thus far, issues regarding these motor control and learning paradigms were critically examined, and their respective mechanisms were compared using descriptive analysis. Elaborate simulations based on the transitions and development flow of each component at issue contributed to the linear approach, laid concrete emphasis on the advantages of the nonlinear approach, and empirically derived the rationale for the indispensable application of system dynamics. Sports science can benefit from system dynamics associated with human motor behavior, as demonstrated in this study.

Abstract

keyword: force controlelementary coordinationsystem dynamicsnetworkagent-based model

Introduction

Humans generate, control, and learn motor behaviors. During various movements or motor states, such processes and principles can be explained from three perspectives or based on three different approaches: (1) describing mechanisms arising from a biological system as having direct linear causal relationships (neuroanatomy and robotics), (2) explaining the phenomena as nonlinear self-organizing structures (Sternad, 2008), and (3) observing quantitative and qualitative states (uncertainty and complexity) as processes that inevitably accompany a system (Smith & Thelen, 2003).

The first perspective views our movements as essential neural mechanisms, focusing on the role of the motor execution system. Multiple motor programs are stored in the upper center (i.e., brain), and motor control occurs as the nervous system is activated by retrieving and adjusting motor programs as needed for execution. Considering that control occurs in the lower corticospinal tracts in a network of the main or unique central nervous system (CNS), such control is considered a classic problem in cognitive neuroscience (Kandel et al., 2000; Purves et al., 2001). From this perspective, similar to a marionette moving through the manipulation of its strings (Rene Descartes, 1596–1650: Ghost in the machine or Homunculus or Tertium quid), our body is interconnected by numerous unit elements and moves by direct motor nerve commands, signals, or transmissions comprising the substrates of the pertinent system (i.e., Turing computation or Newtonian mechanics) according to kinematic and kinetic principles (Turvey, 1990). A framework for this process has become clearly formulated through studies describing phenomena as the product of linear cause-and-effect interactions that are determined based on strict laws. This attempts to break down all phenomena into simple components that can be modelled using linear equations (Lu et al., 2019).

The second perspective regards motion control because of interactions among three factors: the nervous system, body, and environment (Chiel & Beer, 1997). This approach is attentive to events in terms of co-participants in the relevant context rather than to individual interventions based on neural commands or actions that trigger motions (Gibson, 1979). Behavior is much more interesting than linearity. The notion of the linearity (static nature) of mechanical reason has been challenged because it does not solve many questions about how our movements are highly interconnected and nonlinear. According to Bernstein (1967), organisms have an excessive number of degrees of freedom. Moreover, if we consider that all parts are strongly defined by their connections within a certain context, a great complexity can be observed (Rosen, 1991). From a functional perspective, appropriate motions are controlled by elements, actors, and the surrounding events (Fowler & Turvey, 1978). It has been argued that motion is a function rather than a mechanism, and movement is not a simple control of force but an affordance-based coupling, which is the perceived action potential of specific functionally enabled elements (Reed, 1985).

The third perspective involves observing system dynamics that integrate various constituents (Abrams & Strogatz, 2003). The relationship between the individuals or parts refers to how the elements can function together to create a structure or adaptation. For example, particles that are connected in a process are distributed in space and autonomously change over time according to their surroundings (Bell, 1964). When certain substances are mixed in a certain container (i.e., shallow flat dish), spatial-temporal patterns spontaneously emerge (Winfree, 1987). This evidence provides images of the system dynamics as continua that spontaneously order themselves. It emphasizes an understanding of the phenomenon of activation and control at the interconnected macro-micro level (Haken, 1983). To examine the processes rather than the transient states of movement or behavior, a more professional in-depth understanding of the related mechanisms and algorithms has been acquired using essential tools, such as linear algebra, agent-based modeling, and scientific programming, (Davis & Yen, 2019). This idea was first introduced in academia in the mid-20th century as an interdisciplinary approach. Further, it has provided valuable insight into interlinked phenomena (agent-based models and networks) that cannot be solved from a single point of view or subject matter (Ashby, 1947; Schöner, 2002).

This study compares the aforementioned perspectives, evolving along different contextual directions, to evaluate them empirically through related simulations. Accordingly, the concept of a change in force control (linear dynamics) has been used to demonstrate that causal control alone cannot adequately explain various properties of movement. Thus, detailed descriptions of the interactional relationships between multidimensional forces and (nonlinear) elementary coordination dynamics with inherent physical dynamics are provided. Furthermore, the necessity of the system dynamics approach is highlighted by simulating the state and adaptation of these systematically interdependent variables of a macro-(environmental) and micro (biological) level system. This fresh perspective on refocusing the related model will help bridge traditional bias, and will be an inspiration such that a dynamic system, which may not always be complicated, may gain the ability to distill simple principles (Park, 2022a).

Searching for Paradigms

Part 1: Linear Dynamics

Control theories related to computational devices and engineering concepts have established the operating principles of objects such as robots. They have attained a quasi-perfect description and demonstration of their mechanical movements (Greene, 1972). However, the human body differs considerably from these artifacts. Whereas mechanical devices deal with linear input-output relationships, interactions between human muscles and the CNS are nonlinear (Haken, 1987; Kugler & Turvey, 1987). The signal processing used to produce a force is subject to an interdependent influence of numerous internal and external conditions (Newell & Simon, 1976). This limitation in the movement of biological systems and the significant differences in the flow characteristics are explained by the basic concept of force (Kim et al., 1999). To control the elementary variables of a structure with an incomplete and complex design, the CNS must be equipped with a mechanistically strong and elaborately designed structure to generate accurate movements and predict the results.

Feedback Loop

Force control for motion is induced by the CNS and results in movement through various processes. CNS nerves use command signals called control parameters, provided solely by the person who performs a motor behavior and cannot be changed under the influence of external factors or peripheral information. Here, a signal provided independently of its result is called a feedforward or open-loop control. For instance, when a player hits a ball in a ball sport, the brain sends commands to the muscles. Further, the movement thus generated can be executed regardless of whether the player hits the ball or allows a prediction of its result before confirmation.

By contrast, a signal transformed by the controller depending on the self-generated result is called feedback or closed-loop control. An essential element of this type of control is the comparator (or reference) (Adams, 1971), which allows a comparison between the actual and desired effects. For example, when cycling, the cyclists rely on sensory information to obtain information about their speed from the passing wind or external moving environments and adjust their command signals to decide whether to pedal faster or slower or use the brake. This type of control also allows the cyclists to maintain their preferred speed. This achieves the functional purpose of reducing or increasing errors using the ratio of change in an environmental variable (Δ/Δy) through a positive or negative feedback loop.

A typical function of a negative loop is to reduce the initial deviation, whereas that of a positive loop is to add a certain proportional number of control parameters to the deviation of the environmental parameter, which leads to the maintenance of a specific value of the output parameter. However, a significant delay can occur in unexpected situations. For example, in the case of perturbations by an external functional signal or feedback loop under a stable state, a complete compensation for the external functional perturbation needs to be achieved through sufficient error correction (Figure 1, left). Furthermore, its failure can cause delays, leading to unexpected results (Figure 1, right). This response delay is a significant drawback of feedback control; it has strong explanatory power for a task requiring a long time or high accuracy but weak explanatory power for an unexpected intervention or fast movement (Schmidt & Lee, 1999).

In addition, open- and closed-loop control mechanisms combine contrasting actional characteristics. When applying them to a basketball or soccer game, open-loop control allows a player to snatch the ball by predicting the next move of the opponent. Moreover, closed-loop control executes accurate movements based on the present ball movement, relying on sensory information. These control patterns are attributable to a servomechanism (see Appendix 1’s image on the left side for more detail) consisting of muscle spindles (alphaMNs) and a γ-system (gammaMNs), which is a neurophysiological muscle reflex mechanism that produces spontaneous movements from a control circuit combining different schemas (Matthews, 1972; Merton, 1953). Using an open-loop scheme, the controller sends signals set to static values to the servo loop (tonic stretch reflex or TSR loop), measures the difference between the actual and set values (Δ = X1–X2) with the aid of a feedback mechanism. Through this process, errors are indispensable for the servo function. While an experienced servomotor instantly corrects any deviation, a novice servomotor frequently lags and makes mistakes. This is attributable to the auto-control mechanism of the set value within a specified range. For example, in an isometric contraction, the muscle activity (X) is set to a level enabling continuous maintenance by the comparator, causing a spontaneous increase (X↑) or decrease (X↓) in activity when the difference between the actual and set loads (Δ = X1–X2) exceeds a certain level.

However, muscle activity control in spontaneous contractile movements may not be a hierarchical or independent activity in which a feedback mechanism corrects the deviations in an encoded signal transmitted from the center through an open-loop scheme. While maintaining the position of the elbow through bicep muscle activation to resist an external load, the electromyogram signals of the relevant muscle disappeared for a certain period when the load was abruptly removed (unloading reflex) (Vallbo, 1971). However, the muscle resisting a load through a CNS command should maintain its activation level according to the mechanism based on the above hypothesis. This experimental proof, i.e., extensive stretch reflex (Matthews, 1972), is a limitation of the force control hypothesis (i.e., feedback loop and servo hypothesis) in explaining the motor control performed by the human body. Thus, it has shifted the attention of researchers to a more attractive alternative hypothesis.

Equilibrium Point

A gymnast on a high bar can form a suspended and swaying motor structure using the elastic forces of the muscles interacting with the gravitational force.

(1)

U = \sum_{i = 1}^{n} {m g y}_{i} + \frac{1}{2} \sum_{i, j = 1}^{n} {k a}_{i j} {(| p_{i} - p_{j} | - l_{i j})}^{2}

Here, m = mass, g = gravity, k = muscle stiffness, p_i = position of mass defined as i in the x_i – y_i relationship, a_ij = value of the mass (i ) added by the intrinsic properties of the muscle (j ), and l_ij = length of the relaxed muscle. When the body hangs on a high bar, the state of the muscles (U ) is formed by the equilibrium configuration between the potential gravitational energy and the potential elastic energy stored in the muscle. This phenomenon (i.e., isometric = total potential energies of the system), which involves a coupling affecting the muscle length and force until reaching the equilibrium point (equilibrium configuration), is appropriate for CNS commands (elastic potential of the spring), and external loads (gravitational potential of the mass) are reached (Asatryan & Feldman, 1965; Feldman, 1966; Latash, 1993) (see Appendix 1’s image on the right side).

As illustrated through a feedback loop, muscle reflexes are utilized when muscles are activated to characterize the parameters of reflexes that change along with the CNS commands. Because this process is backed by the unrestricted control capacity of the reflex circuit, it can overcome the limitations of the feedback loop or servo control hypothesis mentioned above. Furthermore, it can be appropriately applied to all reflex effects of a spontaneous muscle activation. Moreover, the coupling of internal and external factors reflects the nonlinear form (invariant characteristic) of the equilibrium point hypothesis for muscle control mechanisms (Feldman, 1965, 1986; Feldman & Levin, 1995; Latash, 1993).

According to a detailed report related to a hypothesis quantifying the changing relationships between the elementary (force) and performance variables in a similar coupling context (Latash, 2000), motor control occurs synergistically among various elementary variables in achieving a goal.

(2)

Δ V = (V_{U C M} - V_{O R T}) / (V_{T O T} / n)

Elementary variables involved in achieving the goal of a given movement can be good or bad. In Equation (2), to index their relative amounts, V_UCM, i.e., the variance of the uncontrolled manifold (VUCM), is a variance that has no effect on the performance variables and corresponds to a coupling between elementary variables. In addition, V_ORT, the variance in the space orthogonal to the UCM space (VORT) induces its effect on the performance variables for problem posing and corresponds to performance errors between elementary variables. Moreover, V_TOT denotes the total variance (V_UCM+V_ORT), and n is the degree of freedom of the entire space (n=df_UCM+df_ORT). This hypothesis includes the possibility that the index reflecting the goal to be achieved (ΔV) is determined by V_UCM (good variance), which does not affect the performance variables (Scholz & Schöner, 1999). Therefore, control depends on the extent of synergy among the elementary variables to achieve the intended motor behavior.

The above hypothesis, which indexes goal-directed synergy considering neurophysiological mechanisms, demonstrates that the pattern of force required for movement cannot be quantified as a single structural output of the CNS. In addition, it suggests that the main issue should be expanded from an analysis of only elementary units to determine what controls motor behavior. This may be understood as suggesting that given that the force pattern (control) is combined with other fluid variables, the threshold of the invariant characteristics is nonlinear (Schöner, 2002). Moreover, there is a functional and anatomical collaboration between elementary variables, including balanced signal combinations of all spinal nerve forms (Turvey & Carello, 1996).

Part 2: Nonlinear Dynamics

Classically, the primary concern in movement and motor behavior theories is the concept of force or the mechanism of force generation. Force is mainly controlled at the level of the higher intelligence (i.e., frontal lobe), and the CNS is indispensable in explaining the correct judgment and mechanism of a force used. However, limitations regarding the explanation of the motor execution system are difficult to reveal or generalize when applied to various cases, highlighting the need for an alternative approach (Greene, 1972). Driven by an inquiry into whether neural mapping features remain constant across different situations and whether ambiguous neural schemas can be discerned by motor execution subsystems, the perspective shifted to viewing movement as a nonlinear, self-organizing phenomenon (Kugler et al., 1980).

Point Attractor

Research on nonlinear dynamics began in the 1960s using nonlinear equations (Haken, 1987). Returning to a gymnast swinging on a high bar, suppose the gymnast’s body sways significantly once and then less the next moment. The primary conditions, such as the stiffness of the arm muscles with the gymnasts’ hands hanging on the bar (with stiffness k), gymnast’s mass (m), and friction between the hands and bar (b), are generally equal. The gymnast’s cyclic sway varies depending on the extent to which the gymnast pulls and releases their body when hanging on the bar (position x, velocity x′, and acceleration x″). The harder the pull is, the longer the oscillation, whereas the weaker the pull is, the shorter the oscillation. In both cases, the initial pull eventually ceases. Despite the different durations and shapes, the phenomenon that applies to both cases is a reduction of the energy generated by the pull when swaying toward a nonlinear damping of the oscillation. In both cases, the initial pull eventually ceases. Despite the different durations and shapes, the phenomenon that applies to both cases is a reduction of the energy generated by the pull when swaying toward a nonlinear damping of the oscillation with a gradually diminishing amplitude.

(3)

f (t) = m x^{″} + b x^{'} + k x

This system phenomenon of converging to a static equilibrium point is termed a “point attractor.” (Kugler et al., 1980).

(4)

m x^{″} + b x^{'} + k x = 0

A movement has multiple spatiotemporal paths and trajectories, depending on its initial state or disturbance (Figure 2). However, the system soon converges to a specific preferred state (damping or pumping) according to the general tendency to maintain a state with low variability and high stability for energy efficiency (Saltzman, 1985). This general characteristic is termed nonlinear attractor dynamics.

For example, when we lift our fingers or limbs and move them in-phase (Figure 3, left) or anti-phase (Figure 3, center), we can comfortably maintain the phase. If we then move the oscillators (fingers or limbs) at a different pace (Figure 3, right, where one limb has an oscillation frequency of 1 Hz, and the other has a frequency of 0.25 Hz), the phase changes owing to the eigenvalue of each oscillator, and this out-of-phase mode is less comfortable than the in-phase or anti-phase mode. This is due to the increased coordination variability, with the energy consumption increasing to maintain the eigenvalue of each oscillator. Varying the phase coordination requirements signifies that our movement system converges to the preferred state, that is, in in-phase or anti-phase mode, as in the oscillation motion of our example. This is termed a point attractor (Pikovsky et al., 2001), which is a purely stable pattern that appears during the locomotion of the limbs, head, and body.Table 1

Elementary Coordination

This property of invariant dynamics is the limb oscillation coordination characteristic of human motor behavior and constitutes the basis for explaining the regular locomotion of the limbs (Kay et al., 1987). The phase of an oscillator can be quantified by adding the time-dependent phase angle [ω₁(t)] at the point where the oscillator motion begins [θ₁(0); see Appendix 2 for more details of the schematic depiction), where r is a constant. When similarly estimating the phase of the other oscillator (θ₂), [ω₂(t) + θ₂(0)], the relative phase of the limbs is stable when in-phase (θ₂ - θ₁ ≈ 0) or anti-phase (θ₂ - θ₁ ≈ π).

(5)

φ = ω_{2} - ω_{1} = Δ ω (t) + φ (0) = φ (t)

Our limb coordination system prefers in-phase (φ ≈ 0) mode (Figure 4, reconfirmed by the Haken-Kelso-Bunz (HKB) model) (Haken et al., 1985). Moreover, when the energy demand [movement speed = control parameter (V)] required for coordination stability, such as an increased speed in basic symmetrical motion, increases, the anti-phases of the upper and lower limbs show a relatively unstable pattern. This is because the energy demand increases up to a certain point where a sudden shift to the in-phase mode occurs. By contrast, this phase transition, according to the energy demand, does not always occur under in-phase mode. There is also no phase transition under the same phase because of the energy demand (Kelso, 1984).

(6)

\dot{φ} = - φ δ ω + α cos (φ) - b cos (2 φ)

Here, φ is the phase angle of each activated limb oscillator, - φδω is the difference in the eigenvalue of each oscillator, and α and b are the coefficients of the two limbs, indicating the coupling strength. Equation (6) describes how the relative phase (φ), which represents the difference in the phase angle between two limbs, is determined based on the relative coupling strength according to the ratio between the two coefficients (α/b) (Kugler & Turvey, 1987). It experimentally reflects the nonlinear preferred stability that appears in our coordination process, that is, the attractor of the potential function (see Appendix 3 for more details) (Haken et al., 1985).

In addition, the symmetry of movements derived from the empirically tested nonlinear “point attractor” is spatiotemporally broken depending on the structure and function of the body (Amazeen et al., 1998). This is because routine coordination (e.g., grasping, twisting, and manipulating) involves our limbs and occurs in various ways due to the voluntary (or imposed implementation of) numerous interactions (Collins & Stewart, 1993). Therefore, to explain the coordination pattern and transition characteristics, we extended the basic theoretical model, which can be implemented by breaking the symmetry through frequency competition (Δω=ω_L-ω_R) (Amazeen et al., 1998) in the interlimb coordination model (Haken et al., 1985).

It was also implemented without symmetry, breaking through the asymmetric coupling of the interlimb eigenvalues (=Δω) or coupling intensity (coefficient of b/a) by expanding the coordination asymmetry (preference (dominant arm), attention, and oscillator axis of rotation) (Treffner & Turvey, 1996). For instance, based on the interlimb coordination defined by Equation (6), condition-dependent differences in the coupling intensity of the additional variables under examination were quantified through symmetry breaking (|c and d|>0, |c and d| ≈ 0) to examine their relative differences in phase.

(7)

\dot{φ} = Δ ω - [α sin (φ) + 2 b s i n (2 φ) - [c s i n (φ) + 2 d s i n (2 φ)] + \sqrt{ϱ ξ_{t}}]

In the expansion (Equation (7)), ϕ̇ = change in the relative phase, Δω= the extent of disequilibrium through the eigenvalue of each oscillator, [a cos(ϕ)+ b cos (2ϕ)] is the extent of equilibrium defined at an oscillation value of 0 or π (= inter-oscillator coupling intensity, and ρ = Gaussian white noise, i.e., asymmetric coupling attractor (= [α sin(φ) + 2b sin(2φ)]) intentionally added to the basic interlimb coordination structure (Amazeen et al., 1997; Park et al., 2001). Such applications have enabled us to examine from various angles how essential coordination dynamics are microscopically manifested through asymmetric potentials (see Appendix 4 for more details) (Park, 2020). We can recognize the dynamic structure of a neurological or behavioral variable as a representative state, as expressed by the mean [φ_ave-φ₀ (rad)] and standard deviation [SDφ(rad)] of the relative phase and the coefficient of variation (CV) [SDφ(rad)/φ_ave-φ₀(rad)].

However, this approach cannot elucidate whether the outputs of a movement or motor behavior of individual elements continuously change within a time series, not to mention the situation itself (Horn & Newton, 2019). This has contributed significantly to our understanding of the nonlinear dynamics of coordination patterns and the structural stability of various movements. At the same time, limitations are inevitable under the methodological aspect of experimental research, such as individual differences, sampling, and time series, in objectively explaining the principles and changing modes of motor behavior (Narizuka & Yamazaki, 2018).

Part 3: System Dynamics

System dynamics examines a set of individual elements or the time-series process of a structure and then discusses a single element or causal relationship regarding a specific state (Rosen, 1985). While accepting both a linear approach for force control and a nonlinear approach for coordination dynamics, as discussed above, system dynamics presupposes the likelihood that individual elements may not show unidirectional (i.e., linear or nonlinear) movements or responses concerning internally or externally significant interactions, that is, a gain (i.e., payoff) or loss (i.e., failure) (Chen & Billings, 1992).

(8)

x \to f \vec{\Rightarrow} f (x_{i})

Here, x denotes an element with several entities, and f denotes a function pointing in a specific direction (→) that coincides with the target value. Because the value f of f (x_i ) includes specific entities (x_n ), in addition to other intrinsic entities, the output can be either linear or nonlinear. Moreover, the system assumes that each entity internally interacts with all external entities, which we call the environment (Schöner, 2002).

(9)

x_{t} \to f (x_{i}) \to y_{t}

Here, x_i can perform two types of interactions with the external world (environment). The environment exerts its effect (x) on (x_i) and is given an effect (y) from x_i at a given time (t). Accordingly, a system can be defined by several interactive individual elements inherent in coherent behavior (Juarrero, 1999). This is reflected in the mechanism by which each element can actively change its behavior. It also constitutes the basis for quantitative or qualitative observations of the specific change processes of actual variables (Ashby, 1947; Marchal, 1975).

Network-agent-based Evolution

As previously described, system dynamics can show how parts or components jointly shape a structure and adapt to environmental changes. A greater emphasis is placed on the emergence of stationarity than on the accuracy or variability of force control, stability, and coordination uncertainty. Its crucial variables are the actions of other individuals and the environment (i.e., neighbors) interacting with individuals (e.g., molecules, cells, or athletes). With the changes in its state at each given point in time (ȳ_i) being monitored within the time series (g_i(t)) and its components (x) reflecting autonomous agents (b_i(x)) acting in a bottom-up manner, we assign practical meaning to behavior and state (f (χ̄)) (Poggio & Bizzi, 2004; Reynolds, 1987).

(10)

\bar{f} (x, t) = \sum_{i}^{n} b_{i} (x) g_{i} (t) {\bar{y}}_{i}

Because it can also include mechanisms exploring and imitating the strategies of other constituent elements (i.e., training habits and a tendency to seek healthcare), it enables more realistic estimates of the change and adaptation process based on interactions (Oliva, 2016). Specifically, many individuals and various components, their organizational form (structure: see Table 2 and Figure 5), and basic properties with constituent elements such as a propagation of positive (payoff) and negative (failure) variables allow more detailed (i.e., simulation) observations of processes such as changes in motor behavior (Pastor-Satorras et al., 2015).

(11)

A = G [n, p, β], p a n d β \in (0, 1) \to A m \times n

The significance of this system dynamics approach (i.e., network and agent-based simulation) is emphasized in sports-related research (Clemente et al., 2016). Furthermore, its effectiveness has been demonstrated through myoneural activity (Ryan et al., 2018), sports performance transfer (Bock et al., 2012), and outcome prediction (Bunker & Thabtah, 2019; Oldham & Crooks, 2019). This implies that studies have examined complicated intra- and inter-system interactions between individuals rather than observing the relationship between a force and interlimb coordination (Glazier & Davids, 2010). We previously verified empirically that our movements or motor behaviors are formed by the relative hierarchy of the CNS or its components and arise from the relationships between organizational members, whereby the role of influential individuals can serve as the strategy of that organization (system) (Passos et al., 2017).

According to Newell (1986), who investigated the principles of human motor behavior in the same context, the action of dynamic systems, such as human movement and motor behavior, emanates from interdependent system constraints between individuals, tasks, and the environment. Performer limitations such as motivation, skill level, or physical and psychological characteristics; task limitations such as rules, methods, and performance goals; and environmental constraints such as stadium conditions, jet lag, and sociocultural differences are inevitable constraints constantly encountered by performers. Furthermore, information constraints, such as the awareness of a performer regarding an individually available situation under a specific spatiotemporal setting, also affect decision-making at any instant. When encountering these constraints, we must appropriately interact with the constraining factors under other states for a successful task execution (Araújo & Davids, 2016). Therefore, the systematic handling of these relationships provides a realistic and accurate understanding of the phenomenon of motor behavior. Specifically, it captures the complicated interactions between components and the centrality of patterns and movements, and then analyzes intra-system interactions (Korte & Lames, 2018; Sampaio & Leite, 2013).

The following findings are noteworthy: (i) Functional variables of numerous nerves and muscles interconnected to produce movements (Grillner, 1975) are characterized by systematically different network properties (i.e., visual and motor systems) (Georgopoulos et al., 1982).

(12.1)

V i s u a l s y s t e m : \bar{f} (x) = \sum_{i}^{n} A_{i} k_{i} (\bar{x})

(12.2)

M o t o r s y s t e m : \bar{f} (\bar{x}) = \sum_{i}^{n} {\bar{A}}_{i} k_{i} (\bar{x})

(12.3)

P a r a l l e l s y s t e m : \bar{f} ({\bar{x}}_{j}) = ψ \sum_{i = 1}^{n} {\bar{A}}_{j i} {\bar{x}}_{i} + θ_{j}

Table 3The visual system (Equation (12.1)) is a set of regulated unity activities [A = matrix(eigenvalue), k_i (χ̄ ) = vector(eigenvector)] and mostly derives a single output signal (f̄ (x )) (Figure 6, plot on the left side). By contrast, the motor system (Equation (12.2)) regulates the unit activity of the individual components, and various types of output signals appear as the output vector field (f̄ (χ̄ )) or activations in parallel (f̄ (χ̄ _j )) (Figure 6, plot of the middle side). However, regardless of the origin (i.e., top-down or bottom-up) and destination (i.e., ā→b̄→c̄, b̄→c̄→ā, ..., c̄→ā→b̄) of a signal stimulus and its output mode (single, as the visual system, or multiple and parallel, as the motor system), all possible interactions are formed as networks of the matrix (A ) and vector k_i (χ̄ )) (Wickelgran, 1969) (Figure 6, plot on the left and middle). Even if the action of each component of a different structure is perceived with different intensities (i.e., centrality) depending on the constraints, they are positively or negatively interconnected and affect each other. Thus, a problem in any component results in a structural stability loss in the functionality (Sussillo et al., 2015).

(ii) These network properties estimate the central tactical position during the game, based on the frequency and success rate of passes made in team sports such as soccer, basketball, and handball by expanding the scope of the application (Korte & Lames, 2018) (Figure 6, plot on the right).

(13)

x_{i j} = {\begin{array}{l} 1 i f t i e n_{i}, n_{j} \\ 0 o t h e r w i s e \end{array}

(iii) Athletes mimic, unconsciously and automatically, the momentary emotional expressions of their teammates (Barsade, 2002). Furthermore, the positive or negative mood of a person (athlete or coach) rapidly spreads to other team members, exerting significant effects on their collaboration, conflicts, and task performance (Barsade et al., 2018).

(iv) In addition to changes in the form and injuries, small and large incidents, such as significant changes in individual performance (sports streaks (Bar-Eli et al., 2006)) were observed to considerably impact the functionality of the system (Thacker et al., 2003). We must estimate such system dynamics in a myoneural network with a high interconnection of the components (Pacheco et al., 2006).

Because an individual can be identified by a specific arrangement (n × m) of interdependent relationships within the matrix of the network components to which they belong, more specific and direct content at issue is investigated, such as the connective structure, degree, hierarchy, and centrality. As empirically illustrated in Figure 6, a structure similar to a sports team that forms a low variance at an extremely high degree of interconnection has an exponentially faster and higher risk transmission rate (see Appendix 5 for more details) than a general social structure that does have such a formation (Cui et al., 2021). This system dynamics model has also been applied in situation to basic motor behavior (Park, 2020). Furthermore, the phenomenon of an emergent pattern of key debates of individual movements (i.e., velocity and displacement) was also empirically observed (Figure 7, right) by adding mechanical variables (i.e., imitation and exploration) that can be estimated from simple individual movement characteristics (Figure 7, left: movement).

Complex Behavior

Reverting to the action of forces for more direct relevance to the principles explaining movements and motor behaviors, system dynamics (Haken, 1987) is of particular importance. Inspired by nonlinear dynamics studied during the 1960s, system dynamics characterized the concept of force as an elemental low-energy force field that directly forms a motion and can be applied to all systems as the motion and direct time-dependency change. This is because a slight change in force under the initial state leads to a small change on a short time scale but results in a qualitatively dramatic change on a longer time scale, drawing attention to the change in energy flow in tandem with a change over time rather than the force itself (Kugler et al., 1980). Specifically, the force can be viewed as a dynamic phenomenon inherent in a system governed by a self-recursion rule, rather than simply as a static state, such as the position (x) or velocity (x′ (Rosen, 1987, 1991)).

(14)

(x_{0^{'}} x_{0}) \to (x_{0} + x_{0^{'}} d t, x_{0^{'}} + F d t)

Movement is not simply a change from one position or velocity to another (x_o, x_o _′), but rather a phase change (x_o +x_o _′ dt, x_o _′+Fdt) under a constraint in which individuals, tasks, and environments interact to a limited extent (Bohm, 1969). This represents the view that an individual cannot reach the goal intended by the motor behavior by individually producing all forces required for implementing a task or motion but rather contributes to the flow of the total energy by being involved in the event or direction of concern.

This conceptual or methodological application of system dynamics (i) can estimate the extent to which the phase or relative pattern of a motor behavior within a time series is based on the degree of linear or nonlinear stability, which gradually changes. For example, (ii) based on a linear mean or standard deviation, system uncertainty can be assessed using the probabilistic (p_i ) expected value [log (1/p_i )] of each variable (Shannon, 1948).

(15)

\begin{array}{l} H (X) & = - \sum_{i = 1}^{n} p (x_{i}) \times log p (x_{i}) \\ = \sum_{i = 1}^{n} p_{i} \times log (1 / p_{i}) \end{array}

This is also (iii) a query into structured nonlinear quantitative (Stergiou & Decker, 2011) information on a time-series analysis, including the complexity, which varies according to the sensitivity of the initial conditions of the analyzed object (variable r).

Figure 7 shows the dynamic character of the system (vector operation), which the distance of the basic system phase defined as a given initial property differs in a time series. It denotes the time-dependent attractors’ changes to the parameter of interest and the topographic branching that appears as it passes through a specific time point in the phase space, allowing us to discern the qualitatively different behavioral contributions included in various constraints or motor parameters. The emergent phenomenon that shows the distance in the time-series-defined phase space and the stability (i.e., coordination or cooperation) of the object of observation begins to be driven from the initial state (Figure 8).

Notably, this fundamental process indicates that the second quantity of v⃗ contains the trait of the individual (agent-self interactions, where agents can interact with themselves), based on whether the following condition holds:

(16)

\begin{array}{l} \vec{v} \to \vec{f} \to f (\vec{v}), f (\vec{v}) \\ = {\begin{array}{l} \vec{v} (+) i f f (‖ v_{i} ‖) < ‖ v_{a v g} ‖ \\ \vec{v} (-) i f f (‖ v_{i} ‖) > ‖ v_{a v g} ‖ \end{array} \end{array}

where the function assumes that the attribute of the component is conditional upon the value yield in the other direction (−) if the length of the magnitude ||v_i|| is greater than the other length (||v_avg||). This trait was implemented according to the following quantities:

(17)

\begin{array}{l} \vec{v_{i}} = ‖ v_{i} ‖ * \vec{d_{i}}, ‖ v_{i} ‖ = \sqrt{{v_{i}}_{x}^{2} + {v_{i}}_{y}^{2}}, \vec{d_{i}} \\ = \frac{(v_{i x}, v_{i y})}{\sqrt{{v_{i}}_{x}^{2} + {v_{i}}_{y}^{2}}} \end{array}

(18)

{\vec{v}}_{a v g} = \frac{1}{N} \sum_{i = 1}^{N} \vec{v_{i}} = ‖ v_{a v g} ‖ * {\vec{d}}_{a v g}

where v⃗_i is the velocity of the individual, represented by the length of the individual’s magnitude (||v_i|| ) with the direction of the individual (d⃗_i ). In addition, v⃗_avg is the average velocity, which includes the navigation of the entire population of N individuals. The results of v⃗_i and v⃗_avg produce a new quantity v⃗_inew, which is written simply as (v⃗→v⃗_inew = v⃗_i+v⃗_avg ). This application expands a relatively simple nonlinear equation to gauge the system complexity (Wolf et al., 1985). In addition, we can more elaborately determine the collective dynamic relationships of a group of periodically operated systems.

Conclusions

Human movement and motor behavioral mechanisms were explained using primary data and theories from multidisciplinary research. The concepts of neurophysiological dynamics and cognitive-psychological theories can explain these highly complex phenomena, adding to the breadth and depth of research. In this context, irrespective of theoretical frameworks and methodologies (i.e., equilibrium-point method, uncontrolled manifold method, and HKB model), system dynamics (Strogatz, 1994) may be considered a topic of interest for many researchers. Force control, CNS motor programs, and representation are important in explaining movement principles (Schmidt & Lee, 1999). However, the interactional relationships between multidimensional internal and external forces (Bernstein, 1967) must also be integrated into established frameworks and theories given that the various properties constituting a movement are adequately systematized in the context of an executed movement (Kelso, 1994). Linear or nonlinear perspectives (Turvey & Carello, 1996) imitate all phenomena involved.

Motion control mechanisms have long been explained in many ways (Newell & Simon, 1976). Applications of basic mechanical concepts such as wedges, levers, wheels, pulleys, and screws have made it possible for models of scientific phenomena during the 17th and 18th centuries to develop new types of explanations for the movement principles. A dynamics-based approach, which began in the 1960s with an emphasis on nonlinear methodology, has made it possible to explain in detail how functional behavioral changes occur in geometric methods—in particular, using terms such as orbit, attractor, and branching, establishing theoretical frameworks for various connections (synapses), schemas, and types of logic (programs) for functional structures and data processing processes (Ashby & Hall, 1952).

This approach attracted the attention of researchers during the 1980s and 90s as a suitable structure for understanding human movements (Beer, 1995). Algorithms such as evolution, individual-based, and network properties are beneficial for understanding complex changes for an emergent behavior and self-organization. By substantializing the observations of relatively stable functional parameters and interactive processes and providing methods to conceptualize systems subject to continuous change, we have another tool to explain a movement.

Practical Applications

To move, nerves, muscles, joints, and segments must work efficiently and jointly according to the limited context of the individual, task, and environment involved in the executed movement. Moreover, the conditions of various contexts must be met (Abernethy & Sparrow, 1992). Simultaneously, duly admitting the achievements of classical linear or nonlinear approaches (i.e., motor program or representation) in explaining the dynamics of motion, note that the perspective of multidimensional system dynamics is also essential (Strogatz, 1994). As concretely rationalized in the main body of this paper, the properties of a movement or motor behavior are dynamically interrelated (Kelso, 1994), which cannot be fully understood through causal links or linear/nonlinear descriptions (Iberall & Soodak, 1987).

This system dynamics approach enables movement researchers to observe qualitative motor behavioral changes; thus, it is worth adopting when examining human motor parameters and stable characteristics (Bernstein, 1967). Given that human movements are dynamic and accompanied by time-series structures such as relativity, environmental stability, a coupling between perception and motor parameters, and context dependency effects (Kelso, 1994; Kugler & Turvey, 1987), attention must be paid to all types of adaptive relationships that emerge in diverse and complex human movements (Amazeen, 2002; Latash et al., 1996; Park, 2022b; Zanone & Kelso, 1997).

Against this background, this study empirically examined the dominant perspectives in studies on human movement. Detailed descriptions of the linear force control model were provided, followed by a nonlinear model, including a typical system dynamics model (Enoka, 2008). We realistically proposed the possible role of a system dynamics approach, in addition to explaining complex movement (motor) control. Furthermore, the considerations and proposals made in this study serve as interesting, up-to-date topics conducive to improving the potential of motor behavior research.

Acknowledgments

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education [Grant Number: 2020R1l1A1A01056967] (PI: Chulwook Park).

Author Contributions

Design-Chulwook Park,; Writing and original draft preparation-Chulwook Park,; Review and editing-Chulwook Park,; Validation-Dongwook, Han. All authors have read and agreed to the published version of the manuscript.

References

Abernethy, B., & Sparrow, W. A. (1992) Approaches to the study of motor control and learning (Summers, JJ, Ed.) Amsterdam, Netherlands: North-Holland The rise and fall of dominant paradigms in motor behavior research, pp. 3-45

Abrams, D. M., & Strogatz, S. H. (2003) Modeling the dynamics of language death Nature, 424(6951), 900. .

Adams, J. A. (1971) Journal of Motor Behavior (Vol. 3) A closed-loop theory of motor learning, 2, pp. 111-149, https://doi.org/10.1080/00222895.1971.10734898

Amazeen, E. L., Amazeen, P. G., Treffner, P. J., & Turvey, M. T. (1997) Journal of Experimental Psychology: Human Perception and Performance (Vol. 23) Attention and handedness in bimanual coordination dynamics, 5, pp. 1552-1560, https://doi.org/10.1037/0096-1523.23.5.1552

Amazeen, P. G. (2002) Journal of Motor Behavior (Vol. 34) Is dynamics the content of a generalized motor program for rhythmic interlimb coordination?, 3, pp. 233-251, https://doi.org/10.1080/00222890209601943

Amazeen, P. G., Amazeen, E. L., & Turvey, M. T. (1998) Journal of Motor Behavior (Vol. 30) Breaking the reflectional symmetry of interlimb coordination dynamics, 3, pp. 199-216, https://doi.org/10.1080/00222899809601337

Araújo, D., & Davids, K. (2016) Team synergies in sport: Theory and measures Frontiers in Psychology, 7, 1449 .

Asatryan, D. G., & Feldman, A. G. (1965) Functional tuning of the nervous system with control of movement or maintenance of a steady posture: I. Mechanographic analysis of the work of the joint or execution of a postural task Biophysics, 10(5), 925-935 .

Ashby, R., & Hall, A. (1952) Design for a brain London, England: Chapman.

Ashby, W. R. (1947) Journal of General Psychology (Vol. 37) Principles of the self-organizing dynamic system, 2, pp. 125-128, https://doi.org/10.1080/00221309.1947.9918144

Bar-Eli, M., Avugos, S., & Raab, M. (2006) Psychology of Sport and Exercise (Vol. 7) Twenty years of “hot hand” research: Review and critique, 6, pp. 525-553, https://doi.org/10.1016/j.psychsport.2006.03.001

Barsade, S. G. (2002) Administrative Science Quarterly (Vol. 47) The ripple effect: Emotional contagion and its influence on group behavior, 4, pp. 644-675, https://doi.org/10.2307/3094912

Barsade, S. G., Coutifaris, C. G. V., & Pillemer, J. (2018) Emotional contagion in organizational life Research in Organizational Behavior, 38, 137-151 .

Beer, R. D. (1995) Artificial Intelligence (Vol. 72) A dynamical systems perspective on agent-environment interaction, 1–2, pp. 173-215, https://doi.org/10.1016/0004-3702(94)00005-L

Bell, J. S. (1964) Physics Physique Fizika (Vol. 1) On the einstein podolsky rosen paradox, 3, pp. 195-200, https://doi.org/10.1103/PhysicsPhysiqueFizika.1.195

Bernstein, N. A. (1967) The coordination and regulation of movements Oxford: Pergamon Press

Bock, J. R., Maewal, A., & Gough, D. A. (2012) PLOS ONE (Vol. 7) Hitting is contagious in baseball: Evidence from long hitting streaks, 12, p. e51367, https://doi.org/10.1371/journal.pone.0051367

Bohm, D. (1969) Toward theoretical biology (Waddington, CH, Ed., Vol. 2) Chicago, IL: Aldine Publishing Company Some remarks on the notion of order, pp. 18-40

Bunker, R. P., & Thabtah, F. (2019) Applied Computing and Informatics (Vol. 15) A machine learning framework for sport results prediction, 1, pp. 27-33, https://doi.org/10.1016/j.aci.2017.09.005

Chen, S., & Billings, S. A. (1992) International Journal of Control (Vol. 56) Neural networks for nonlinear dynamic system modelling and identification, 2, pp. 319-346, https://doi.org/10.1080/00207179208934317

Chiel, H. J., & Beer, R. D. (1997) Trends in Neurosciences (Vol. 20) The brain has a body: Adaptive behavior emerges from interactions of the nervous system, body and environment, 12, pp. 553-557, https://doi.org/10.1016/s0166-2236(97)01149-1

Clemente, F. M., Martins, F. M. L., & Mendes, R. S. (2016) Social network analysis applied to team sports analysis Chaam, NetherlandsInternational Publishing: Springer

Collins, J. J., & Stewart, I. N. (1993) Journal of Nonlinear Science (Vol. 3) Coupled nonlinear oscillators and the symmetries of animal gaits, 1, pp. 349-392, https://doi.org/10.1007/BF02429870

Cui, S., Gao, M., Xun, Y., Fung, S. F., Tan, Y., Zhang, Y., & Xiong, Y. (2021) Research on the structure and characteristics of the overall social network of professional athletes Complexity, 2021, 1-11 .

Davis, W. S., & Yen, D. C. (Eds.) (2019) The information system consultant’s handbook: Systems analysis and design Boca Raton, FL: CRC Press

Enoka, R. M. (2008) Neuromechanics of human movement Champaign, IL: Journal of Human Kinetics

Feldman, A. G. (1965) Functional tuning of the nervous system with control of movement or maintenance of a steady posture-II. Controllable parameters of the muscle Biophysical Journal, 11, 565-578 .

Feldman, A. G. (1966) Functional tuning of the nervous system with control of movement or maintenance of a steady posture--III. Mechanographic analysis of the execution by man of the simplest motor tasks Biophysics, 766775 .

Feldman, A. G. (1986) Journal of Motor Behavior (Vol. 18) Once more on the equilibrium-point hypothesis (λ model) for motor control, 1, pp. 17-54, https://doi.org/10.1080/00222895.1986.10735369

Feldman, A. G., & Levin, M. F. (1995) Behavioral and Brain Sciences (Vol. 18) The origin and use of positional frames of reference in motor control, 4, pp. 723-744, https://doi.org/10.1017/S0140525X0004070X

Fowler, C. A., & Turvey, M. T. (1978) Skill acquisition: An event approach with special reference to searching for the optimum of a function of several variables Information processing in motor control and learning, 1-40 .

Georgopoulos, A. P., Kalaska, J. F., Caminiti, R., & Massey, J. T. (1982) Journal of Neuroscience (Vol. 2) On the relations between the direction of two-dimensional arm movements and cell discharge in primate motor cortex, 11, pp. 1527-1537, https://doi.org/10.1523/JNEUROSCI.02-11-01527.1982

Gibson, J. J. (1979) The ecological approach to visual perception Hillsdale, NJ: Lawrence Erlbaum Associates

Glazier, P. S., & Davids, K. (2009) Sports Medicine (Vol. 39) Constraints on the complete optimization of human motion, 1, pp. 15-28, https://doi.org/10.2165/00007256-200939010-00002

Greene, P. H. (1972) Problems of organization of motor systems Progress in Theoretical Biology, 2, 123-145 .

Grillner, S. (1975) Physiological Reviews (Vol. 55) Locomotion in vertebrates: Central mechanisms and reflex interaction, 2, pp. 247-304, https://doi.org/10.1152/physrev.1975.55.2.247

Haken, H. (1983) Advanced synergetics Berlin, Germany: Springer

Haken, H. (1987) Self-organizing systems (Yates, FE, Ed.) Boston, MA: Springer Synergetics, pp. 417-434

Haken, H., Kelso, J. A., & Bunz, H. (1985) A theoretical model of phase transitions in human hand movements Biological Cybernetics, 51(5), 347-356 .

Horn, T. S., & Newton, J. L. (2019) Advances in sport and exercise psychology (Horn, TS, & Smith, AL, Eds., 4th ed) Champaign, IL: Human Kinetics Publishers Developmentally based perspectives on motivated behavior in sport and physical activity contexts, pp. 313-331

Iberall, A. S., & Soodak, H. (1987) Self-organizing systems (Yates, FE, Ed.) Boston, MA: Springer A physics for complex systems, pp. 499-520

Juarrero, A. (1999) Dynamics in action : Intentional behavior as a complex system Cambridge, MA: MIT Press

Kandel, E. R., Schwartz, J. H., Jessell, T. M., Siegelbaum, S., Hudspeth, A. J., & Mack, S. (Eds.) (2000) Principles of neural science (Vol. 4) New York, NY: McGraw-Hill pp. 1227-1246

Kay, B. A., Kelso, J. A., Saltzman, E. L., & Schöner, G. (1987) Journal of Experimental Psychology. Human Perception and Performance (Vol. 13) Space-time behavior of single and bimanual rhythmical movements: Data and limit cycle model, 2, pp. 178-192, https://doi.org/10.1037//0096-1523.13.2.178

Kelso, J. A. (1984) Phase transitions and critical behavior in human bimanual coordination American Journal of Physiology, 246(6), R1000-R1004 .

Kelso, J. A. S. (1994) Interlimb coordination: Neural, dynamical, and cognitive constraints (Swinnen, SP, Heuer, H, Massion, J, & Casaer, P, Eds.) San Diego, CA: Academic Press Elementary coordination dynamics

Kim, S., Carlton, L. G., Liu, Y. T., & Newell, K. M. (1999) Journal of Motor Behavior (Vol. 31) Impulse and movement space-time variability, 4, pp. 341-357, https://doi.org/10.1080/00222899909600999

Korte, F., & Lames, M. (2018) Characterizing different team sports using network analysis Current Issues in Sport Science, 3, 005 .

Kugler, P. N., & Turvey, M. T. (1987) Information, natural law, and the self-assembly of rhythmic movement Hillsdale, NJ: Erlbaum

Kugler, P. N., Scott Kelso, J. A., & Turvey, M. T. (1980) 1 on the concept of coordinative structures as dissipative structures: I. Theoretical lines of convergence Advances in Psychology, 1, 3-47 .

Latash, M. (2000) Motor Control (Vol. 4) There is no motor redundancy in human movements. There is motor abundance, 3, pp. 259-260, https://doi.org/10.1123/mcj.4.3.259

Latash, M. L. (1993) Control of human movement Champaign, IL: Human Kinetics

Latash, M. L., Turvey, M. T., & Bernshteĭn, N. A. (1996) Dexterity and its development Mahwah, NJ: Lawrence Erlbaum

Lu, S., He, Q., & Wang, J. (2019) A review of stochastic resonance in rotating machine fault detection Mechanical Systems and Signal Processing, 116, 230-260 .

Marchal, J. H. (1975) Philosophy of Science (Vol. 42) On the concept of a system, 4, pp. 448-468, https://doi.org/10.1086/288663

Matthews, P. B. (1972) Mammalian muscle receptors and their central actions American Journal of Physical Medicine & Rehabilitation, 53(3), 143-144 .

Merton, P. A. (1953) Ciba foundation symposium-The spinal cord (Malcolm, JL, & Gray, JAB, Eds.) Chichester, England: John Wiley & Sons Speculations on the servo-control of movement, pp. 247-260

Narizuka, T., & Yamazaki, Y. (2018) Characterization of the formation structure in team sports Proceedings of the Institute of Statistical Mathematics, 65(2), 299-307 .

Newell, A., & Simon, H. A. (1976) Communications of the ACM (Vol. 19) Computer science as empirical inquiry: Symbols and search, 3, pp. 113-126, https://doi.org/10.1145/360018.360022

Newell, K. (1986) Motor development in children: Aspects of coordination and control (Wade, MG, & Whiting, HTA, Eds.) Dordrecht, Netherlands: Springer Constraints on the development of coordination

Oldham, M., & Crooks, A. T. (2019) Drafting agent-based modeling into basketball analytics Proceedings of 2019 Spring Simulation Conference, 1-12 .

Oliva, R. (2016) System Dynamics Review (Vol. 32) Structural dominance analysis of large and stochastic models, 1, pp. 26-51, https://doi.org/10.1002/sdr.1549

Pacheco, J. M., Traulsen, A., & Nowak, M. A. (2006) Physical Review Letters (Vol. 97) Coevolution of strategy and structure in complex networks with dynamical linking, 25, p. 258103, https://doi.org/10.1103/PhysRevLett.97.258103

Park, C. (2020) Complex Adaptive Systems Modeling (Vol. 8) Evolutionary understanding of the conditions leading to estimation of behavioral properties through system dynamics, 1, p. 2, https://doi.org/10.1186/s40294-019-0066-x

Park, C. (2022a) Eigenproperties of perception (dynamic touch) and action (phase dynamic) out of diversities Human Movement Science, 85, 102999 .

Park, C. (2022b) Practical foundation of applied systems analysis Seoul: National University Press In press, 978-89-521-3098-3 (93560)

Park, C., & Kim, S. (2020) Network-agent systems dynamic modelling: A guide for sport science Korean Journal of Sport Science, 31(3), 514-533 .

Park, H., Collins, D. R., & Turvey, M. T. (2001) Journal of Experimental Psychology Human Perception and Performance (Vol. 27) Dissociation of muscular and spatial constraints on patterns of interlimb coordination, 1, pp. 32-47, https://doi.org/10.1037/0096-1523.27.1.32

Passos, P., Araújo, D., & Volossovitch, A. (2017) Performance analysis in team sports London, England: Routledge

Pastor-Satorras, R., Castellano, C., Van Mieghem, P., & Vespignani, A. (2015) Reviews of Modern Physics (Vol. 87) Epidemic processes in complex networks, 3, pp. 925-979, https://doi.org/10.1103/RevModPhys.87.925

Pikovsky, A., Rosenblum, M., & Kurths, J. (2001) A universal concept in non-linear sciences (Vol. 12) Cambridge, England: Cambridge University Press

Poggio, T., & Bizzi, E. (2004) Nature (Vol. 431) Generalization in vision and motor control, 7010, pp. 768-774, https://doi.org/10.1038/nature03014

Purves, D., Augustine, G. J., Fitzpatrick, D., Katz, L. C., Lamantia, A. S., McNamara, J. O., & Williams, S. M. (2001) Neuroscience (2nd ed) Sunderland, England: Sinauer

Reed, M. I. (1985) Redirections in organizational analysis (Vol. 293) London, England: Tavistock

Reynolds, C. W. (1987) ACM SIGGRAPH Computer Graphics (Vol. 21) Flocks, herds and schools: A distributed behavioral model, 4, pp. 25-34, https://doi.org/10.1145/37402.37406

Rosen, R. (1985) Anticipatory systems: Philosophical, mathematical, and methodological foundations (Rosen, R, Rosen, J, Kineman, JJ, & Nadin, M, Eds.) Oxford, England: Pergamon Press International federation for systems research. Anticipatory systems, pp. 339-399

Rosen, R. (1987) Quantum implications: Essays in honor of David Bohm (Hilley, BJ, & Platt, FD, Eds.) New York, NY: Routlege & Kegan Some epistemological issues in physics and biology, pp. 315-327

Rosen, R. (1991) Life itself New York, NY: Columbia University Press

Ryan, K., Lu, Z., & Meinertzhagen, I. A. (2018) Journal of Comparative Neurology (Vol. 526) The peripheral nervous system of the ascidian tadpole larva: Types of neurons and their synaptic networks, 4, pp. 583-608, https://doi.org/10.1002/cne.24353

Saltzman, E. (1985) Status report on speech research SR-84 Haskins Laboratories Task dynamic coordination of the speech articulators: A preliminary model, pp. 1-18

Sampaio, J., & Leite, N. (2013) Routledge handbook of sports performance analysis (McGarry, T, O’Donoghue, P, & Sampaio, AJ, Eds.) New York, NY: Routledge Performance indicators in game sports, pp. 133-144

Schmidt, R. A., & Lee, T. D. (1999) Motor control and learning Champaign, IL: Human Kinetics Publishers

Scholz, J. P., & Schöner, G. (1999) Experimental Brain Research (Vol. 126) The uncontrolled manifold concept: Identifying control variables for a functional task, 3, pp. 289-306, https://doi.org/10.1007/s002210050738

Schöner, G. (2002) International encyclopedia of the social and Behavioral Science (Wright, JD, Ed.) Amsterdam, Netherlands: Elsevier Dynamical systems approach to neural systems and behavior, pp. 10571-10575

Shannon, C. E. (1948) Bell System Technical Journal (Vol. 27) A mathematical theory of communication, 3, pp. 379-423, https://doi.org/10.1002/j.1538-7305.1948.tb01338.x

Smith, L. B., & Thelen, E. (2003) Trends in Cognitive Sciences (Vol. 7) Development as a dynamic system, 8, pp. 343-348, https://doi.org/10.1016/s1364-6613(03)00156-6

Stergiou, N., & Decker, L. M. (2011) Human Movement Science (Vol. 30) Human movement variability, nonlinear dynamics, and pathology: Is there a connection?, 5, pp. 869-888, https://doi.org/10.1016/j.humov.2011.06.002

Sternad, D. (2008) Coordination: Neural, behavioral and social dynamics (Fuchs, A, & Jirsa, VK, Eds.) Berlin, Germany: Springer Towards a unified theory of rhythmic and discrete movements—Behavioral, modelling, and imaging results, pp. 105-133, https://doi.org/10.1007/978-3-540-74479-5_6

Strogatz, S. H. (1994) Non-linear dynamics and chaos Reading, MA: Addison-Wiley

Sussillo, D., Churchland, M. M., Kaufman, M. T., & Shenoy, K. V. (2015) Nature Neuroscience (Vol. 18) A neural network that finds a naturalistic solution for the production of muscle activity, 7, pp. 1025-1033, https://doi.org/10.1038/nn.4042

Thacker, S. B., Stroup, D. F., Branche, C. M., & Gilchrist, J. (2003) Prevention of knee injuries in sports: A systemic review of the literature Journal of Sports Medicine and Physical Fitness, 43(2), 165-179 .

Treffner, P. J., & Turvey, M. T. (1996) Experimental Brain Research (Vol. 107) Symmetry, broken symmetry, and handedness in bimanual coordination dynamics, 3, pp. 463-478, https://doi.org/10.1007/BF00230426

Turvey, M. T. (1990) American Psychologist (Vol. 45) Coordination, 8, pp. 938-953, https://doi.org/10.1037//0003-066x.45.8.938

Turvey, M. T., & Carello, C. (1996) Dexterity and its development (Latash, M, & Turvey, MT, Eds.) Mahwah, NJ: Erlbaum Dynamics of Bernstein’s level of synergies, pp. 339-376

Vallbo, ÅB. (1971) Journal of Physiology (Vol. 218) Muscle spindle response at the onset of isometric voluntary contractions in man. Time difference between fusimotor and skeletomotor effects, 2, pp. 405-431, https://doi.org/10.1113/jphysiol.1971.sp009625

Wickelgran, W. A. (1969) Psychological Review (Vol. 76) Context-sensitive coding, associative memory, and serial order in (speech) behavior, 1, pp. 1-15, https://doi.org/10.1037/h0026823

Winfree, A. T. (1987) When time breaks down Princeton, NJ: Princeton University Press

Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. (1985) Physica D: Nonlinear Phenomena (Vol. 16) Determining lyapunov exponents from a time series, 3, pp. 285-317, https://doi.org/10.1016/0167-2789(85)90011-9

Zanone, P. G., & Kelso, J. A. (1997) Journal of Experimental Psychology. Human Perception and Performance (Vol. 23) Coordination dynamics of learning and transfer: Collective and component levels, 5, pp. 1454-1480, https://doi.org/10.1037//0096-1523.23.5.1454

Appendices

Appendix 1

Illustrations of the Servo Control and Equilibrium Point Hypothesis

The image of the left side, each circuit denotes a population of neurons (Ns), interneurons (INs), and motoneurons (alphaMNs, gammaMNs). The system in its stable configuration (the image of the right side), this realization appears to be equilibrium points (EP). The point of intersection between the load characteristic is the equilibrium point of the system corresponding to a combination of muscle length (φ) and muscle force (p). solid line: characteristic of the tonic stretch reflex (static muscle torque versus joint angle); : its threshold: dashed line: some load characteristic; p. φ: equilibrium torque and angle, respectively. Note: each level has a specific function in the illustration (Feldman, 1986).

Appendix 2

The mechanism of oscillation on two different pendulums

Corresponding to the model which introduced in background, the mechanism of oscillation on two different (left pendulum and right pendulum) but nearly identical process phases was defined by the following dynamic (φ = θ_L – θ_R ). Here, φ is the phase of the strength between the left hand (θ_L ) and the right hand (θ_R ). The degree of relative phases (0°–180°) depends on the difference in the two oscillators. If each θ was defined as a sine function, as follows $ω_{1} - \frac{α}{2} sin (θ_{1} - θ_{2}), θ_{2} = ω_{2} - \frac{α}{2} sin (θ_{2} - θ_{1})$

, the logic can simply be rewritten to φ as the same dynamic function (θ₁ – θ₂ = φ). Then, each θ can be specified as the following equation:

Data from both hands were dividing into left and right components. These are (real) numbers and computed from the above phase coordination mathematical model. It can be infinitely large data set corresponding time series, probability distribution (discrete state variable) or probability density (continuous state variable) describing certain aspects of a distribution (Amazeen et al., 1997; Treffner & Turvey, 1996; Park, 2022).

Appendix 3

Preferred movement frequencies of the individual oscillators

Simulation on the left denotes the energy of the function at each averaged relative phase (Blue line = the vertical-axis). 0 is the most stable value for the system. The horizontal-axis indicates averaged relative phase between two limbs from in-phase at 0 to the anti-phase at π= 180° (−π = −180°). At the local point of 0 and π (−π), the function is close to those minima (attractors; black discs), and at the local point around 90 (−90), the state is close to maxima (repellers: red discs). The green and red lines denote the variation of the potential (V) functions depending on the two components’ preferred frequencies. Simulation on the right: blue line indicates identical symmetry; green line indicates large different symmetry (function: b/a = 0.5, detuning = −0.5); and red line indicates even larger different symmetry (function: b/a = 0.5, detuning = −1.5) (Haken et al., 1985; Park, 2022).

Appendix 4

Experimental characteristics for the extended symmetry breaking forces (entropy)

Each simulation denotes the estimated entropy forces according to the time series for each experimental design. The simulation on the left side denotes the entropy forces of the normal condition (red line = 5:00, green line = 12:00, black line = 17:00). The simulation in the middle side denotes the entropy forces of the heat-based normal (red line = 5:00, black line = 17:00) vs. abnormal (green line = 5:00, yellow line = 17:00) conditions. The simulation on the right side denotes the entropy forces of the ice-based normal (red line = 5:00, black line = 17:00) vs. abnormal (green line = 5:00, yellow line = 17:00) conditions. Notice considering the experimental design, the delta T value was 0.06, and the number of realizations of the ensemble was N = 50 (averaged pendulum frequency per minute) (Park, 2020).

Appendix 5

Prototype risk propagation distributions

The left (low λχ) and right (high λχ) simulations show the number of failures for the time step (t = 10) with the corresponding area proportional to the values. As the failure propagates, the area between the failure (red) and non-failure (blue) increases in each time step. Notice that the network property was gathered from two different sports groups’ data (upper = social group, bottom = athlete group). The eigenvector centrality is a feature of an adjacency matrix in a network (Ax = λχ, λ denotes the eigenvalue). Further, a higher connection for a specific node in an organization increases the eigenvector centrality and less variation (low), and vice versa (see Table below) (Pacheco et al., 2006; Park & Kim, 2020).

χ (n_{i}) = \frac{1}{λ} \sum_{i \neq j} x_{i j} c_{E V} (n_{j}), C_{E V} (n_{j}) = e i g e n v e c t o r v a l u e o f n o d e n_{j}

View original

Table

Properties of a created network (N = node)

G(n, p)		N_0	N_1	N_2	N_3	N_4	N_5	N_6	N_7	N_8	N_9
n = 10, p = 0.2	degree	1	2	3	1	1	3	1	3	4	1
n = 10, p = 0.2	Ax = λx	0.148	0.312	0.427	0.174	0.164	0.402	0.207	0.362	0.509	0.207

n = 10, p = 0.9	degree	8	9	8	7	8	9	9	9	9	8
n = 10, p = 0.9	Ax = λx	0.302	0.334	0.305	0.269	0.302	0.334	0.334	0.334	0.334	0.305

Note: For a synthetic network produced using G(n, p). To quantify the probability that a node has a degree for all [0 < d < (n-1)], a node is considered to have a degree of zero when nothing is connected, and the expected averaged node degree is p(n-1) = 0.9 × 9 = 8.1 in case of the n = 10, p = 0.9.

Figures and Tables

View original

Figure 1

Representative experimental simulation of the estimated cumulative feedback function. The x-axis denotes the time steps, and the y-axis refers to the control function. The set of left plots denotes a uniform distribution with the function (+1) and its result. By repeating the same property with the perturb function [(+1)*p(−1)], the construct becomes a delayed and non-static result (shown in the plot on the right).

View original

Figure 2

Simulation of different attractors. Damped (or pumped) decay lines indicate an embedded invariant property in the relationship between system attractor and damping (or pumping) potential. Dotted line = reference point.

View original

Figure 3

Experimental representation of synchronous diagrams from two different oscillators (red = oscillator 1, blue = oscillator 2): Left = synchronized almost in-phase, with a phase difference; middle = anti-phase condition; right = unmatched initial input between objects (a.u. = arbitrary unit). Note that all the structures were created through a circular function to clearly show their differences.

View original

Figure 4

Schematic of the experimental conditions and results of different phases: Left = general experimental setting with apparatus; middle = plot; right = the observed relative phase or phase relation between two oscillators at φ≈0 deg (in-phase) or φ≈ ± 3.14 deg (anti-phase = π). The movement speed denotes the control parameter, and the arrow indicates the attractor. In the plot on the right, V is the intrinsic dynamics of the potential function, and the black balls symbolize stable states (attractors) and red balls correspond to unstable states (repellers), φ2–φ1 ≈ 0 denotes a condition of a nearly synchronized in-phase, and φ2–φ1 ≈ π indicates a condition of an anti-phase.

View original

Figure 5

Athlete network topologies from professional athletes in South Korea. The left and right plots depict random and regular structures, respectively. The middle plot depicts the network structure of male baseball players (n = 40, the actual level of support among team members = 5.5, and the actual level of support outside of team members = 0.18). These topologies were created using the parameters n = individual, p = regular connections within the team, and β = random connections outside the team. The questionnaire had a simple question regarding who would support the player (or contact) when having important matters to discuss (e.g., injury, advice, mistake, necessary decision-making, or poor behavior).

View original

Figure 6

Generalization architectures of motor systems. In the case of vision (left), a single output signal is a combination of tuned unit activities (with weighted edges of a node). In the case of motor control (middle), each component of the output field is a combination of tuned unit activities (with weighted edges of a node). The right side shows its adjacent matrix corresponding to the network connection scheme. The plot of the right set shows a created team sport passing the network property of football data (average value on successful passes). The tracking of the passing distribution at the FIFA World Cup 2014 tournament was provided in the official FIFA match reports on their website (WorldCup/archive). The plot of the football shows that a defensive midfielder is significantly more central (P < .05) than other positions (see Table 3 for more details). Note a set of nodes connected by ties that link them (xij defines a binary tie between nodes).

View original

Figure 7

Behavioral dynamics underlying the network structure (see Table 2) and simple movement rule. Based on the simulation, the plot on the left shows that a displacement separates the individuals at the initial time steps (T = 1), and the relative position structure is controlled by the initial setting. Blue dots represent their position in an x, y coordinate plane, and the red lines denote links. The middle plots show the dynamic influence based on the individual characteristics at T = 10 and T = 100 (middle right). The right plot represents its behavioral bias calculated by the density function (T = time steps). It refers to the fact that although the pattern of the individual behavior depends on a localized view of the initial conditions, a slight change in the individual movement characteristic (velocity) has a remarkably diverging (or converging) displacement.

View original

Figure 8

Schematic of the individual (vector = v⃗i, averaged = v⃗avg ) operation in R2.

View original

Table 1

Data collection for experiment

Participants (N)	Joint	Trial	Task/rest (min)
8	Wrist (in-phase)	6	1m / 5m
8	Wrist (anti-phase)	6	1m / 5m

Note. Data collection for the phase test: Group one = 8 participants, 6 trials, at wrist with in-phase (total number of datasets = 48) and wrist with anti-phase (total number of datasets = 48). Each trial lasted 1 min, with 5-min rest intervals between trials.

View original

Table 2

Athletes Network structure (Small-world property)

Define	Parameterize	Range
Number of individuals	N	∈(1, ∞)
Connection prob (regular)	P	∈(0, 1)
Connection prob (random)	B	∈(0, 1)
Linear dimension (matrix)	m × n	row, col

Note. The data structure includes three parameters, where m × n represents the adjacent matrix.

View original

Table 3

Descriptive network property and its adjacent matrix for football data

	DM	GK	LCD	LD	LF	LM	OM	RCD	RD	RF	RM
Centrality	0.11	0.06	0.09	0.09	0.08	0.09	0.10	0.10	0.10	0.08	0.10
Eigen_Cent	0.49	0.13	0.26	0.26	0.21	0.32	0.37	0.27	0.26	0.21	0.32