scholarly journals Rendering neuronal state equations compatible with the principle of stationary action

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Erik D. Fagerholm ◽  
W. M. C. Foulkes ◽  
Karl J. Friston ◽  
Rosalyn J. Moran ◽  
Robert Leech

AbstractThe principle of stationary action is a cornerstone of modern physics, providing a powerful framework for investigating dynamical systems found in classical mechanics through to quantum field theory. However, computational neuroscience, despite its heavy reliance on concepts in physics, is anomalous in this regard as its main equations of motion are not compatible with a Lagrangian formulation and hence with the principle of stationary action. Taking the Dynamic Causal Modelling (DCM) neuronal state equation as an instructive archetype of the first-order linear differential equations commonly found in computational neuroscience, we show that it is possible to make certain modifications to this equation to render it compatible with the principle of stationary action. Specifically, we show that a Lagrangian formulation of the DCM neuronal state equation is facilitated using a complex dependent variable, an oscillatory solution, and a Hermitian intrinsic connectivity matrix. We first demonstrate proof of principle by using Bayesian model inversion to show that both the original and modified models can be correctly identified via in silico data generated directly from their respective equations of motion. We then provide motivation for adopting the modified models in neuroscience by using three different types of publicly available in vivo neuroimaging datasets, together with open source MATLAB code, to show that the modified (oscillatory) model provides a more parsimonious explanation for some of these empirical timeseries. It is our hope that this work will, in combination with existing techniques, allow people to explore the symmetries and associated conservation laws within neural systems – and to exploit the computational expediency facilitated by direct variational techniques.

2016 ◽  
Vol 31 (07) ◽  
pp. 1650027 ◽  
Author(s):  
Mahdi Hajihashemi ◽  
Ahmad Shirzad

Following the Poincaré algebra, in the Hamiltonian approach, for a free spinning particle and using the Casimirs of the algebra, we construct systematically a set of Lagrangians for the relativistic spinning particle which includes the Lagrangian given in the literature. We analyze the dynamics of this generalized system in the Lagrangian formulation and show that the equations of motion support an oscillatory solution corresponding to the spinning nature of the system. Then we analyze the canonical structure of the system and present the correct gauge suitable for the spinning motion of the system.


From the general principles of quantum mechanics it is deduced that the wave equation of a particle can always be written as a linear differential equation of the first order with matrix coefficients. The principle of relativity and the elementary nature of the particle then impose certain restrictions on these coefficient matrices. A general theory for an elementary particle is set up under certain assumptions regarding these matrices. Besides, two physical assumptions concerning the particle are made, namely, (i) that it satisfies the usual second-order wave equation with a fixed value of the rest mass, and (ii) either the total charge or the total energy for the particle-field is positive definite. It is shown that in consequence of (ii) the theory can be quantized in the interaction free case. On introducing electromagnetic interaction it is found that the particle exhibits a pure magnetic moment in the non-relativistic approximation. The well-known equations for the electron and the meson are included as special cases in the present scheme. As a further illustration of the theory the coefficient matrices corresponding to a new elementary particle are constructed. This particle is shown to have states of spin both 3/2 and 1/2. In a certain sense it exhibits an inner structure in addition to the spin. In the non-relativistic approximation the behaviour of this particle in an electromagnetic field is the same as that of the Dirac electron. Finally, the transition from the particle to the wave form of the equations of motion is effected and the field equations are given in terms of tensors and spinors.


2013 ◽  
Vol 28 (01) ◽  
pp. 1250234 ◽  
Author(s):  
A. A. DERIGLAZOV

We propose Lagrangian formulation for the particle with value of spin fixed within the classical theory. The Lagrangian is invariant under non-Abelian group of local symmetries. On this reason, all the initial spin variables turn out to be unobservable quantities. As the gauge-invariant variables for description of spin we can take either the Frenkel tensor or the Bargmann–Michel–Telegdi (BMT) vector. Fixation of spin within the classical theory implies O(ℏ)-corrections to the corresponding equations of motion.


Author(s):  
C Mellace ◽  
A P Lai ◽  
A Gugliotta ◽  
N Bosso ◽  
T Sinokrot ◽  
...  

One of the important issues associated with the use of trajectory coordinates in railroad vehicle dynamic algorithms is the ability of such coordinates to deal with braking and traction scenarios. In these algorithms, track coordinate systems that travel with constant speeds are introduced. As a result of using a prescribed motion for these track coordinate systems, the simulation of braking and/or traction scenarios becomes difficult or even impossible. The assumption of the prescribed motion of the track coordinate systems can be relaxed, thereby allowing the trajectory coordinates to be effectively used in modelling braking and traction dynamics. One of the objectives of this investigation is to demonstrate that by using track coordinate systems that can have an arbitrary motion, the trajectory coordinates can be used as the basis for developing computer algorithms for modelling braking and traction conditions. To this end, a set of six generalized trajectory coordinates is used to define the configuration of each rigid body in the railroad vehicle system. This set of coordinates consists of an arc length that represents the distance travelled by the body, and five relative coordinates that define the configuration of the body with respect to its track coordinate system. The independent non-linear state equations of motion associated with the trajectory coordinates are identified and integrated forward in time in order to determine the trajectory coordinates and velocities. The results obtained in this study show that when the track coordinate systems are allowed to have an arbitrary motion, the resulting set of trajectory coordinates can be used effectively in the study of braking and traction conditions. The results obtained using the trajectory coordinates are compared with the results obtained using the absolute Cartesian-coordinate-based formulations, which allow modelling braking and traction dynamics. In addition to this numerical validation of the trajectory coordinate formulation in braking scenarios, an experimental validation is also conducted using a roller test rig. The comparison presented in this study shows a good agreement between the obtained experimental and numerical results.


2021 ◽  
Author(s):  
Friedrich Pfeiffer

Abstract Constraints in multibody systems are usually treated by a Lagrange I - method resulting in equations of motion together with the constraint forces. Going from non-minimal coordinates to minimal ones opens the possibility to project the original equations directly to the minimal ones, thus eliminating the constraint forces. The necessary procedure is described, a general example of combined machine-process dynamics discussed and a specific example given. For a n-link robot tracking a path the equations of motion are projected onto this path resulting in quadratic form linear differential equations. They define the space of allowed motion, which is generated by a polygon-system.


Author(s):  
Jinhao Qiu ◽  
Junji Tani

Abstract Equations of motion for multi-layer piezoelectric cylindrical shells and the equations of the integrated piezoelectric sensors are derived. The state equation is obtained by solving the equations of motion with modal expansion method. The feedback control, feedforward control, and their combination are applied in the control of forced vibration of the piezoelectric cylindrical shell with integrated sensor and actuators. The simulation and experimental results show that good control effectiveness can be obtained by using the integrated piezoelectric sensor and actuators in conjunction with the combination of feedback and feedforward control methods.


Author(s):  
T. Hatch ◽  
A. P. Pisano

Abstract A two-degree-of-freedom (2-DOF), analytical model of a hydraulic valve lifter is derived. Special features of the model include the effects of bulk oil compressibility, multi-mode behavior due to plunger check valve modeling, and provision for the inclusion of third and fourth body displacements to aid In the use of the model in extended, multi-DOF systems. It is shown that motion of the lifter plunger and body must satisfy a coupled system of third-order, non-linear differential equations of motion. It is also shown that the special cases of zero oil compressibility and/or 1-DOF motion of lifter plunger can be obtained from the general third-order equations. For the case of zero oil compressibility, using Newtonian fluid assumptions, the equations of motion are shown to reduce to a system of second-order, linear differential equations. The differential equations are numerically integrated in five scenarios designed to test various aspects of the model. A modal analysis of the 2-DOF, compressible model with an external contact spring is performed and is shown to be in excellent agreement with simulation results.


Author(s):  
H K Kim ◽  
Y-S Park

An efficient state-space method is presented to determine time domain forced responses of a structure using the Lagrange multiplier based sub-structure technique. Compared with the conventional mode synthesis methods, the suggested method can be particularly effective for the forced response analysis of a structure subjected to parameter changes with time, such as a missile launch system, and/or having localized non-linearities, because this method does not need to construct the governing equations of the combined whole structure. Both the loaded interface free-free modes and free interface modes can be employed as the modal bases of each sub-structure. The sub-structure equations of motion are derived using Lagrange multipliers and recurrence discrete-time state equations based upon the concept of the state transition matrix are formulated for transient response analysis. The suggested method is tested with two example structures, a simple lumped mass model with a non-linear joint and an abruptly parameter changing structure. The test results show that the suggested method is very accurate and efficient in calculating forced responses and in comparing it with the direct numerical integration method.


1987 ◽  
Vol 177 ◽  
pp. 1-18 ◽  
Author(s):  
D. L. R. Oliver ◽  
J. N. Chung

The steady-state equations of motion are solved for a fluid sphere translating in a quiescent medium. A semi-analytical series truncation method is employed in conjunction with a cubic finite-element scheme. The range of Reynolds numbers investigated is from 0.5 to 50. The range of viscosity ratios is from 0 (gas bubble) to 107 (solid sphere). The flow structure and the drag coefficients agree closely with the limited available experimental measurements and also compare favourably with published finite-difference solutions. The strength of the internal circulation was found to increase with increasing Reynolds number. The flow patterns and the drag coefficient show little variation with the interior Reynolds number. Based on the numerical results, predictive equations for drag coefficients are recommended for both moderate- and low-Reynolds-number flows.


Author(s):  
Seyed Hasan Miri Roknabadi ◽  
Mohamad Fakhari Mehrjardi ◽  
Mehran Mirshams

This paper presents an optimal attitude maneuver by Reaction Wheels to achieve desired attitude for a Satellite. At first, Dynamic Equations of motion for a satellite with just three Reaction Wheels of its active actuators are educed, and then State Equations of this system are obtained. An optimal attitude control with the LQR method has exerted for a distinct satellite by its Reaction Wheels. As a result simulation has presented an optimal effort by calculated Gain matrix to achieve desired attitude for chosen Satellite. It shows that satellite becomes stable in desired attitude with a low energy and time consumption. Furthermore equations derivation, coupling of electrical Reaction Wheel equations with dynamic equations of satellite motion, linearizes them and Reaction wheel saturation avoidance approaches are innovative. Simulation results, accuracy of achieving desired attitude and satellite stability support this statement.


Sign in / Sign up

Export Citation Format

Share Document