Avatar robotics involves remotely controlling a robot to interact with the physical environment on behalf of a human operator, enabling them to virtually embody the robot and perform actions as if physically present. This transformative technology extends human presence to remote or hazardous locations, with applications spanning space exploration, disaster response, remote inspection, telemedicine, and diverse domains. Leveraging progress in robotics, teleoperation systems, sensory feedback interfaces, and communication networks, avatar robotics enhances human capabilities, ensures safer operations, and broadens human presence and expertise in various fields.

Furthermore, cybernetic control functions as a regulatory system utilizing feedback mechanisms to uphold stability and achieve desired outcomes. The incorporation of feedback loops is central to cybernetic control systems, continuously monitoring a system’s behavior, comparing it to a reference state, and generating corrective actions to address any deviations. This iterative feedback process facilitates self-regulation and goal attainment within the system. The application domains of cybernetic control span engineering, biology, and psychology, with the goal of enabling robots to interact with humans in more intuitive ways. This involves adapting their actions and responses based on human feedback and behavior, cultivating a more seamless and responsive human-robot interaction.

Cybernetic biorobotics, at its core, is an interdisciplinary frontier that harmonizes principles from cybernetics, biology, and robotics. Its primary mission is the exploration and development of robots or robotic systems intricately inspired by the marvels of biological organisms. This field is driven by the ambition to conceive robots capable of mimicking and integrating the sophisticated principles and behaviors observed in living entities. Researchers draw inspiration from the intricate control systems of biological organisms and this creative synthesis results in the creation of robots characterized by adaptive and intelligent behaviors, thus mirroring the intricacies found in the natural world. Bioinspired robotics, a central focus within this discipline, involves distilling the fundamental principles and behaviors intrinsic to biological entities and skillfully incorporating them into the design and control of robotic systems. It has the potential to advance the development of robots endowed with locomotion and manipulation capabilities akin to animals, as well as robots capable of adapting to dynamic environments or interacting with humans in more natural and intuitive ways. Moreover, research in cybernetic biorobotics can offer valuable insights into comprehending biological systems, fostering advancements in disciplines like neuroscience and biomechanics. Furthermore, remote cybernetic robots may rely on haptic systems as essential interfaces. A haptic system, characterized by its ability to provide users with a sense of touch or tactile feedback through force, vibration, or other mechanical means, comprises a haptic interface and a haptic rendering system. Collaboratively, these components simulate touch sensations, enabling users to engage with virtual or remote environments in a tactile manner.

This research introduces a control and sensing architecture that integrates a cybernetic scheme based on the recognition of electromyographic control signals, governing a range of locomotive behaviors in a robotic fish. Conceptually, the human operator receives feedback signals from the sensors of the bio-robotic avatar, conveying information about its remote environment. The proposed approach stands out due to its key features and contributions, which include:

Section 2 provides a comprehensive discussion on the comparative analysis of the current state of the art. Section 3 provides a detailed description of the proposed architecture cybernetic system model. In Section 4 presents an approach for filtering electromyography (EMG) data and delves into an in-depth discussion of a classifier based on deep ANN for the recognition of hand-motion EMG stimuli patterns. Section 5 presents the development of a fuzzy-based oscillation patterns generator. Section 6 details the robot’s mechanism parts and its dynamic model. Section 7 focuses on the development of a feedback control for the fish’s ballasting system. Finally, Section 8 provides the concluding remarks of the research study.

This section syntheses the relevant literature and provides insights into the current state of the art. Further, it aims to examine and evaluate existing research and advancements in the field. This brief analysis identifies and compares different aspects, providing a comprehensive overview including relevant research and advancements in the field about methodologies, and outcomes.

Multiple basic concepts of cybernetics [1] in the intersection of physics and control theory and molecular systems were presented in [2] speed-gradient approach to modeling the dynamics of physical systems is discussed. A novel research approach namely Ethorobotics proposes the use and development of advanced bioinspired robotic replicas as a method for investigating animal behaviour [3]. In the domain of telepresence and teleoperation, diverse systems and methodologies have been devised to facilitate remote control of robots [4]. One such system is the multi-robot teleoperation system based on a brain-computer interface, as documented by [6]. This system aims to enable individuals with severe neuromuscular deficiencies to operate multiple robots solely through their brain activity, thus offering telepresence via a thought-based interaction mode. A comprehensive review addressing existing teleoperation methods and techniques for enhancing the control of mobile robots has been presented by [12]. This review critically analyzes, categorizes, and summarizes existing teleoperation methods for mobile robots while highlighting various enhancement techniques that have been employed. It makes clear the relative advantages and disadvantages associated with these methods and techniques. The field of telepresence and teleoperation robotics has witnessed substantial attention and interest over the past decade [14], finding extensive applications in healthcare, education, surveillance, disaster recovery, and corporate/government sectors. In the specific context of underwater robots, gesture recognition-based teleoperation systems have been developed to enable users to control the swimming behavior of these robots. Such systems foster direct interaction between onlookers and the robotic fish, thereby enhancing the intuitive experience of human-robot interaction. Furthermore, efforts have been made to enhance the consistency and quality of robotic fish tails through improved fabrication processes, and target tracking algorithms have been developed to enhance the tracking capabilities of these robots [11]. The work [10] developed teleoperation for remote control of a robotic fish by hand gestures recognition. It allowed direct interaction between onlookers and the biorobot. Another notable system is the assistive telepresence system employing augmented reality in conjunction with a physical robot, as detailed in the work by [7]. This system leverages an optimal non-iterative alignment solver to determine the optimally aligned pose of the 3D human model with the robot, resulting in faster computations compared to baseline solvers and delivering comparable or superior pose alignments. The review presented in [20] analyses the progress of robot skin as multimodal sensing and machine perception for sensory feedback in feeling proximity, pressure, temperature for collaborative robot applications considering immersive teleoperation and affective interaction. The work [24] reported an advanced robotic avatar system designed for immersive teleoperation, having some key functions such as human-like manipulation and communication capabilities, immersive 3D visualization and transparent force-feedback telemanipulation. Suitable human-robot collaboration in medical application has been reported [21], where force perception is augmented to human operator during needle insertion on the soft tissue. Telepresence of mobile robotic systems may incorporate remote video transmission to steer the robot by seeing through its eyes remotely. The work [22] presented an overview including social application domains. Research has been conducted on the utilization of neural circuits to contribute to limb locomotion [15] in the presence of uncertainty. By optimizing data, it showed the combination of circuits necessary for efficient locomotion. A review has also been conducted on central pattern generators (CPGs) employed for locomotion control in robots [16]. This review encompasses neurobiological observations, numerical models, and robotic applications of CPGs. The work [27] describes an extended mathematical model of the CPG supported by two neurophysiological studies: identification of a two-layered CPG neural circuitry and a specific neural model for generating different patterns. The CPG model is used as the low-level controller of a robot to generate walking patterns. Inclusion of ANN as a layer of the CPG, to produce rhythmic and non-rhythmic motion patterns. The work in [18] presented a review of bionic robotic fish, tackling major concepts on kinematics and control, learning, hydrodynamic forces and critical concepts on locomotion coordination. The research presented in [19] reviews the human manual control of devices in cybernetics using mathematical models and advances of theory and applications, from linear time-invariant modeling of stationary conditions to methods and analysis of adaptive and time-varying of cybernetics-human in control tasks. New foundations for cybernetics will emerge and impact numerous domains involving humans in manual and neuromuscular systems modeling control.

Building upon the preceding analysis regarding relevant literature, the subsequent table (Table 1) encapsulates the primary distinctions articulated in this study in relation to the most pertinent literature identified.

As delineated in Table 1, the present study introduces distinctive elements that set it apart from the recognized relevant literature. However, it is noteworthy to acknowledge that various multidisciplinary domains may exhibit commonalities. Across these diverse topics, shared elements encompass robotic avatars, teleoperation, telepresence, immersive human-robot interfaces, as well as haptic or cybernetic systems in different application domains. In this research, the fundamental principle of a robotic avatar entails controlling its swimming response to biological stimuli from the human operator. The human controller is able to gain insight into the surrounding world of the robotic fish avatar through a haptic interface. This interface allows the human operator to yield biological electromyography stimuli as the result of their visual and skin impressions (e.g. pressure, temperature, heading vibrations). The biorobotic fish generates its swimming locomotive behavior that is governed by EMG stimuli yielded in real-time in the human. Through a neuro-fuzzy controller, the neuronal part (cybernetic observer) classifies the type of human EMG reaction, and the fuzzy part determines the swimming behavior.

This section encompasses a comprehensive framework that highlight the integration of various components to propose a cohesive cybernetic robotic model. In addition, this section outlines the key concepts and elucidates their interactions within the system.

Figure 1 presents an overview of the key components constituting the proposed system architecture. This manuscript thoroughly explores the modeling of four integral elements: i) the cybernetic human controller, employing ANN classification of EMG signals; $ii$) a fuzzy-based locomotion pattern generator; $iii$) an underactuated bioinspired robot fish; and $iv$) the robot’s sensory system, contributing feedback for the haptic system. While we will discuss the relevance and impact of the latter item within the architecture, it is important to note that the detailed exploration of topics related to haptic development and wearable technological devices goes beyond the scope of this paper and will be addressed in future work. Nevertheless, we deduce the observable variables that serve as crucial inputs for the haptic system.

Essentially, there are six haptic feedback sensory inputs of interest for the human, representing the observable state of the avatar robot: Eulerian variables, including angular and linear displacements and their higher-order derivatives; biomechanical caudal motion; hydraulic pressure; scenario temperature; and passive vision. Figure 2-left provides an illustration of the geometric distribution of the sensing devices.

The instrumented robot is an embodiment of the human submerged in water, featuring an undulatory swimming mechanical body imbued with muscles possessing underactuated characteristics. These features empowers the bio-robotic avatar to execute movements and swim its aquatic surroundings.

The observation models aim to provide insights into how these sensory perceptions are conveyed to the haptic helmet with haptic devices, including wheel reaction mechanism. A comprehensive schema emerges wherein the haptic nexus, bolstered by pivotal human’s biosensorial components including gravireceptors, Ruffini corpuscles, Paccinian receptors, and retinal photoreceptors, converges to interface with the sensory substrate of the human operator. Such a convergence engenders a cascading sequence wherein biological input stimuli coalesce to yield discernible encephalographic activity, the primary layer of subsequent electromyographic outputs. These consequential EMG outputs, undergo processing a swim oscillatory pattern generator, thereby embodiment control of biomechanical cybernetic governance.

In accordance with Figure 1, it is noteworthy that among the various haptic input variables such as temperature, pressure, and the visual camera images represent direct sensory measurements transmitted from the robotic fish to the components of the haptic interface. Conversely, the robotic avatar takes on the role of a thermosensory adept, in order for the human to assess the ambient thermal landscape. Thus, from this discernment, a surrogate thermal approach is projected onto the tactile realm of the human operator, through modulation of thermally responsive plates enmeshed within the haptic interface. Therefore, a crafted replica of the temperature patterns detected by the robotic aquatic entity is seamlessly integrated into the human sensory experience. This intertwining of thermal emulation is reached by the network of Ruffini corpuscles, intricately nestled within the human skin, thereby enhancing the experiential authenticity of this multisensory convergence. As for interaction through the haptic functions, the Paccinian corpuscles function as discerning receptors, proficiently registering subtle haptic pressures. It finds its origin in the dynamic tactile signals inherent to the aquatic habitat, intricately associated with the underwater depth traversed by the robotic avatar. Integral to the comprehensive sensory scheme, the optical sensors housed within the robotic entity acquires visual data. This visual data is subsequently channeled to the human’s cognition, through the haptic interface’s perceptual canvas. Within this, the human sensory apparatus assumes the role of an engaged receptor, duly transducing these visual envoys through the lattice of retinal photoreceptors.

Embedded within the robotic fish’s body, several Inertial Measurement Units (IMU) plays a pivotal role in quantifying Eulerian inclinations intrinsic to the aquatic environment. These intricate angular displacements are subsequently channeled to the human operator, thereby initiating an engagement of the reaction wheel mechanism. As a consequential outcome of this interplay, a synchronized emulation of tilting motion is induced, mirroring the nuanced cranial adjustments executed by the human operator. Any alignment of movements assumes perceptible form, relayed through the humans’s network of gravireceptors nestled within the internal auditory apparatus. The Euler angular and linear speeds are not directly measured; instead, they must be integrated using various sensor fusion approaches to enhance the avatar’s fault tolerance in reading its environment. For example, the angular observations of the robot fish are obtained through numerical integro-differential equations, which are solved online as measurements are acquired. Let’s introduce the following notation for inclinometers (${\alpha }_i$) and accelerometers ($a_{\mathrm{a}}$), with the singular direct sensor measurement ${\dot{\alpha }}_g$ derived from the gyroscopes. The observation for the fish’s roll velocity combining the three inertial sensors is while the pitch velocity is modelled by and the yaw velocity is obtained by

Within this context, the tangential accelerations experienced by the robot body are denoted as $a_{\alpha ,\beta ,\gamma }$ [m/s${ }^2$]. Additionally, the angular velocities measured by the gyroscopes are represented by $\dot{\alpha },\dot{\beta },\dot{\gamma }$ [rad/s${ }^2$]. Correspondingly, the inclinometers provide angle measurements denoted as $\alpha ,\beta ,\gamma$ [rad]. These measurements collectively contribute to the comprehensive observability and characterization of the robot’s dynamic behavior and spatial orientation.

Furthermore, the oscillations of the caudal tail are reflections of the dynamics of the underactuated spine. These dynamics are captured by quantifying encoder pulses, denoted as ${\eta }_t$, which provide precise angular positions for each vertebra. Given that real-time angular measurements of the vertebrae are desired, higher-order data is prioritized. Consequently, derivatives are computed by initiating from the Taylor series to approximate the angle of each vertebral element with respect to time, denoted as t. thus, rearranging Math notation and trunking up to the first-derivative, dropping ${\varphi }_i^{\left(1\right)}$ off as a state variable, the first-order derivative (${\varphi }^{\left(1\right)}\left(t\right)\equiv \dot{\varphi }\left(t\right)$) is given by and assuming that a vertebra’s angular measurement model and in terms of the encoder’s pulses $\eta$ with resolution R, then it is stated that by substituting the pulses encoder model into the angular speed function for the first vertebra, as for the second vertebra, and

The preliminary sensing models serve as a comprehensive representation, strategically integrated into the control models as crucial feedback terms. A detailed exploration of this integration is elucidated in Section 6 and Section 7.

This section details the experimental acquisition of EMG data, its spatial filtering, and pattern extraction achieved through the statistical combination of linear envelopes. Additionally, an adaptive method for class separation and data dispersion reduction is described. The section also covers the structure of a deep neural network, presenting its classification output results from mapping input EMG stimuli.

The related work reported a system for automatic pattern generation for neurosimulation in [25], where a neurointerface was used as a neuro-protocol for outputting fingers deflection and nerves stimulation. In the present research, numerous experiments were carried out to pinpoint the optimal electrode placement and achieve precise electromyographic readings for each predefined movement in the experiment. The positions of the electrodes were systematically adjusted, and the results from each trial were compared. Upon data analysis, it was discerned that the most effective electrode placement is on the Ulnar nerve, situated amidst the muscles Flexor Digitorium Superficialis, Flexor Digitorium Profundus, and Flexor Carpi Ulnaris. A series of more than ten experiments was executed for each planned stimulus or action involving hands, allowing a 2s interval between each action, including the opening and closing of hands, as well as the extension and flexion of the thumb, index, middle, ring and little fingers. The data was measured by a g.MOBIlab+ device with two-channel electrodes and quantified in microvolts per second [$\mu v/s$], as depicted in Figure 3.

The data acquired from the electromyogram often exhibit substantial noise attributed to both the inherent nature of the signal and external vibrational factors. To refine the data quality by mitigating this noise, a filtering process is essential.

In this context, a second-order Notch filter was utilized. This filter is tailored to target specific frequencies linked to noise, proving particularly effective in eliminating electrical interferences and other forms of stationary noise [26]. A Notch filter is a band-rejection filter to greatly reduce interference caused by a specific frequency component or a narrow band signal. Hence, in the Laplace space, the second-order filter is represented by the analog Laplace domain transfer function: where ${\omega }_0$ signifies the cut angular frequency targeted for elimination, and $2\xi$ signifies the damping factor or filter quality, determining the bandwidth. Consequently, solving to obtain its solution in the physical variable space. The bilinear transformation relates the variable s from the Laplace domain to the variable $z_k$ in the $\mathbb{Z}$ domain, considering T as the sampling period, and is defined as follows, upon substituting the previous expression into the transfer function of the analog Notch filter and algebraically simplifying, the following transfer function in the $\mathbb{Z}$ domain is obtained, redefining the notation as ${\upsilon }_t=z_k$ just to meet equivalence with the physical variable, The 2^nd-order Notch filter $h\left(\upsilon \right)$ was employed on raw EMG data to alleviate noise resulting from electrical impedance and vibrational electrode interference, with parameters set at ${\omega }_0=256$Hz and $\xi =0.1$ and results depicted in Figure 4.

Following this, leveraging the outcomes of the Notch filter, the data undergoes processing through three distinct filters or linear envelopes. This serves as a secondary spatial filter and functions as a pattern extraction mechanism. These include a filter for average variability, one for linear variability, and another for average dispersion. Each filter serves a specific purpose, enabling the analysis of different aspects of the signal. Consider n as the number of measurements constituting a single experimental stimulus, and let N represent the entire sampled data space obtained from multiple measurements related to the same stimulus. Furthermore, denote ${\widehat{v}}_i$ as the $i^{th}$ EMG measurement of an upper limb, measured in microvolts ($\mu V$). From such statements, the following Propositions 1, 2 and 3 are introduced as new data patterns.

Hence, let the vector $\overrightarrow{\mathbf{\delta }}\in {\mathbb{R}}^3$ such that ${\mathbf{\delta }}_k={({\Omega }_k,{\gamma }_k,{\lambda }_k)}^{\top }$ and represent filtered data points in the $\Omega \gamma \lambda$-space. This work considered to include a brief data preprocessing as a method to improve multi-class separability and data scattering reduction in pattern extraction. Three distinctive patterns—$\gamma ,\lambda ,\Omega$—captivatingly converge in Figure 5. This illustration exclusively features patterns associated with sequences of muscular stimuli from both the right and left hands. For supplementary stimulus plots, refer to the appendix in Appendix A at the end of this manuscript.

From numerous laboratory experiments, over $75\%$ of the sampled raw data fall within the range of one standard deviation. Consider the vector $\overrightarrow{\mathbf{\sigma }}\in {\mathbb{R}}^3$ such that the standard deviations vector $\overrightarrow{\mathbf{\sigma }}={({\sigma }_{\Omega },{\sigma }_{\gamma },{\sigma }_{\lambda })}^{\top }$ encompasses the three spatial components by its norm. Building upon the preceding statement, we can formulate an adaptive discrimination criterion as elucidated in Definition 1.

Hence, considering the recent Definition 1, and in the current scenario with ${\kappa }_d=1.0$, which serves as a tuning discrimination factor. Therefore, the norm $l_h$ represents the distance between the frame origin and any class in the $\Omega \gamma \lambda$-space. This distance is adaptively calculated based on the statistics of each EMG class. where, the coefficients ${\kappa }_{\Omega }$, ${\kappa }_{\gamma }$, ${\kappa }_{\lambda }$ are smooth adjustment parameters to set separability along axes. Hence, relocating each class center to a new position is stated by Proposition 4.

Thus, by following the step-by-step method outlined earlier, Figure 6 showcases the extracted features of the EMG data, representing diverse experimental muscular stimuli. These results hold notable significance in the research, as they successfully achieve the desired class separability and data scattering, serving as crucial inputs for the multilayer ANN.

Henceforth, the focus lies in identifying and interpreting the EMG patterns projected in the $\gamma \lambda \Omega$-space, as illustrated in Figure 6. The subsequent part of the section delve into the architecture and structure of the deep ANN employed as a classifier, providing a detailed account of the training process. Additionally, the section highlights the performance metrics and results achieved by the classifier, offering insights into its effectiveness. Despite the challenges posed by nonlinearity, multidimensionality, and extensive datasets, various neural network structures were configured and experimented with. These configurations involved exploring different combinations of hidden layers, neurons, and the number of outputs in the ANN. To achieve the highest success rate in accurate data classification through experimentation, the final ANN was designed with perceptron units, as depited in Figure 7. It featured three inputs corresponding to the three EMG patterns $\gamma ,\lambda ,\Omega$ and included 12 hidden layers, each with 20 neurons. The supervised training process, conducted on a standard-capability computer, took approximately 20–30 minutes, resulting in nearly 1% error in pattern classification.

However, in the initial stages of the classification computations, with some mistuned adaptive parameters, the classification error was notably higher, even with much deeper ANN structures, such as 100 hidden layers with 99 neurons per layer. To facilitate the implementation, this work utilized the C/C++ library Fast Artificial Neural Networks (FANN), generating extremely fast binary code once the ANN was trained. In the training process of this research, about 50 datasets from separate experiments for each type of muscular stimulus were collectively stored, each comprising several thousand muscular repetitions. A distinct classification label was assigned a priori for each class type within the patterns space. To demonstrate the reliability of the approach, 16 different stimuli per ANN were established for classification and recognition, resulting in the ANN having 4 combinatory outputs, each with two possible states. Figure 8 depicts mixed sequences encompassing all types of EMG stimuli, with the ANN achieving a 100% correct classification rate.

Moreover, Table 2 delineates the mapping relationship between the ANN’s input, represented by the EMG stimuli, and the ANN’s output linked to a swimming behavior for controlling the robotic avatar.

This section delineates the methodology utilized to produce electric oscillatory signals, essential for stimulating the inputs of electromagnetic devices (solenoids) embedded within the mechanized oscillator of the bio-robotic fish. The outlined approach for generating periodic electric signals encompasses three key components: a) the implementation of a fuzzy controller; b) the incorporation of a set of periodic functions dictating angular oscillations to achieve desired behaviors in the caudal undulation of the fish; and c) the integration of a transformation model capable of adapting caudal oscillation patterns into step signals, facilitating the operation of the dual-coil electromagnetic oscillator.

Another study [28] reported a different approach, a neuro-fuzzy-topological biodynamical controller for muscular-like joint actuators. In the present research, an innovative strategy suggested for the fuzzy controller involves the combination of three distinct input fuzzy sets: the artificial neural network outputs transformed into Crisp sets, and the linear and angular velocities of the robot derived from sensor measurements. Simultaneously, the fuzzy outputs correspond to magnitudes representing the periods of time utilized to regulate the frequency and periodicity of the caudal oscillations. This comprehensive integration enables the fuzzy controller to effectively process both neural network-derived information and real-time sensor data, dynamically adjusting the temporal parameters that govern the fish’s undulatory motions.

The depiction of the outputs from the EMG pattern recognition neural network is outlined in Table 2. Each binary output in the table is linked to its respective crisp-type input fuzzy sets when represented in the decimal numerical base $y_C$, as illustrated in Figure 9a. Moreover, Definition 2 details the parametric nature of the input sets.

Therefore, the reasoning rules articulated in the inference engine have been devised by applying inputs derived from Table 2, specifically tailored to generate desired outputs that align with the oscillation periods T [s] of the fish’s undulation frequency (see Figure 10).

Hence, the output fuzzy sets delineate the tail undulation speeds of the robot fish across the caudal oscillation period T [s], as illustrated in Figure 10. Notably, three identical output sets with distinct concurrent values correspond to the periods of three distinct periodic functions forward: (f), right (r) and left (l); as subsequently defined by equations (22a), (22b), and (22c).

The inference engine’s rules dictate the application of fuzzy operators to assess the terms involved in fuzzy decision-making. As for the thrusting velocity fuzzy sets, where $v=any$ was previously established, and by applying Definition 3 the fuzzy operator is described by Likewise, the angular velocity fuzzy sets, where previously $\omega =any$ was stated, by applying the second part of Definition 3, the following fuzzy operator is described by Therefore, according to premise ($v=any$ or $\omega =any$), the following fuzzy operator applies, Essentially, the fuzzification process applies strictly similar for the rest of the inference rules, according to the following Proposition,

Executing the previously mentioned proposition, Remark Section 5 provides a demonstration of its application.

Moving forward, during the defuzzification process, the primary objective is to attain an inverse solution. The three output categories for periods T include ‘forward’, ‘right’, and ‘left’, all sharing identical output fuzzy sets of T (Figure 10). Nevertheless, the output fuzzy sets consist of two categories of distributions: Gauss and sigmoid distribution sets, as outlined in Definition 4. Regarding the Gauss distributions, their functional form is specified by:

In accordance with Definition 4, any ${\mu }_k$ possesses a normalized membership outcome within the interval ${\mu }_k\in [0,\cdots ,1]$. The inverse sigmoid membership function, denoted as $T_k\in \mathbb{R}$∀${\mu }_k$, is determined by the general inverse expression:Similarly, the inverse Gaussian membership function, where $T_k\in \mathbb{R}$∀${\mu }_k$, is defined by the inverse function:

Hence, exclusively for the $j^{th}$ category among the output fuzzy sets affected by the fuzzy inference rule essential for estimating the value of T, the Centroid method is employed for deffuzification through the following expression: or more specifically,

For terms $T\left({\mu }_f\right)$ and $T\left({\mu }_{ts}\right)$, the inverse membership function (18) is applicable, whereas for the remaining sets in the $j^{th}$ category, the inverse membership (19) is applied.

The work [29] reported a CPG model to control a robot fish’s motion in swimming and crawling, and let perform different motions influenced by sensory input from light, water, and touch sensors. Oscillators and Central Pattern Generators (CPGs) are closely related concepts. Oscillators are mathematical or physical systems exhibiting periodic behavior and are characterized by the oscillation around a stable equilibrium point (limit cycle). In the context of CPGs, these are neural networks that utilize oscillators that create positive and negative feedback loops, allowing for self-sustaining oscillations and the generation of rhythmic patterns, particularly implemented in numerous robotic systems [30]. CPGs are neural networks found in the central nervous system of animals (e.g. fish swimming [31]), that generate rhythmic patterns of motor activity and are responsible for generating and coordinating optimized [32] repetitive movements.

The present research proposes a different approach from the basic CPG model, and as a difference from other wire-driven robot fish’s motion approaches [33], this study introduces three fundamental undulation functions: forward, right-turning, and left-turning. These functions are derived from empirical measurements of the robot’s caudal fin oscillation angles. However, a distinctive behavioral undulation swim is achieved by blending these three oscillation functions, each incorporating corresponding estimation magnitudes derived from the fuzzy controller outputs. The formulation of each function involves fitting empirical data through Fourier series. As a difference from other approaches on CPG parameters adjustment[34], the preceding fuzzy outputs obtained from (21) to estimate the time periods $T_f,T_r,T_l$ play a pivotal role in parameterizing the time periods for the periodic oscillation functions, as outlined in Proposition 6.

The approaches to forward, right-turn, and left-turn based on the findings of Proposition 6 are illustrated in Figure 11. Additionally, a novel combined oscillation pattern emerges by blending these three patterns (23), each assigned distinct numerical weights through the neuro-fuzzy controller. The proposed robotic mechanism features a multilink-based propulsive spine, driven by an electromagnetic oscillator composed of a pair of antagonistic solenoids that necessitate a synchronized sequence of electric pulses (see Figure 12b). The amplitudes generated by $\psi (\varphi ,T_f,T_r,T_l)$ in Equation 23 essentially represent the desired undulation pattern for the robotic fish’s caudal fin. However, these oscillations are not directly suitable for the inputs of the coils. To address this, our work introduces a decomposition of $\psi$ into two step signals centered around a stable equilibrium point (limit cycle), one for the right coil (positive with respect to the limit cycle) and another for the left coil (negative with respect to the limit cycle). The coil’s step function, either for the right-sided or left-sided coil is given by $s_{r,l}$, taken the equilibrium point as their limit value $\xi$,

In contrast to the work presented in [35] focusing on the swimming modes and gait transition of a robotic fish, the current study, as depicted in Figure 11, introduces a distinctive context. The three oscillatory functions, ${\psi }_{f,r,l}$, are displayed both overlapped and separated, highlighting their unique decoupled step signals. Assuming $\xi =0$ for all $\varphi$ in each case, in Figure 11a, ${\psi }_r={\psi }_l\approx 0$, with $T_{r,l}\ge 6$; in Figure 11b, ${\psi }_f={\psi }_l\approx 0$, with $T_{f,l}\ge 6$; and in Figure 11c, ${\psi }_f={\psi }_r\approx 0$, with $T_{f,r}\ge 6$. For a more comprehensive understanding, Figure 12b presents the electromechanical components of the caudal motion oscillator.

This section introduces the design of the robotic fish mechanism and explores the model of the underactuated physical system to illustrate the fish undulation motions. The conceptualization of the proposed system is inspired by an underactuated structure featuring a links-based caudal spine with passive joints, utilizing helical springs to facilitate undulatory locomotion (see Figure 12a). The robotic fish structure introduces a mechanical oscillator comprising a pair of solenoids activated through coordinated sequences of step signals, as described by (24). Essentially, the electromagnetic coils generate antagonistic attraction/repulsion linear motions, translating into rhythmic oscillations within a mechanized four-bar linkage (depicted in Figure 12b). This linkage takes on the form of a trapezoid, composed of two parallel rigid links and two lateral linear springs functioning as antagonistic artificial muscles. Moreover, beneath the electromagnetic oscillator, there is a ballasting device for either submersion or buoyancy (Figure 12c). The robot’s fixed reference system consists of the X axis, which intersects the lateral sides, and the Y axis aligned with the robot’s longitudinal axis.

In Figure 12a, the electromagnetic oscillator of the robotic avatar responds to opposing coordinated sequences of step signals. The right-sided (R) and left-sided (L) solenoids counteract each other’s oscillations, generating angular moments in the trapezoid linkage (first vertebra). Both solenoids are identical, each comprising a coil and a cylindrical neodymium magnet nucleus. The trapezoid linkage, depicted in Figure 12b, experiences magnetic forces $\pm f_{o{s}_{RL}}$ at the two neodymium magnet attachments situated at a radius of $r_{os}$, resulting in two torques, ${\tau }_{os}$ and ${\tau }_s$, with respect to their respective rotation centers. As input forces $\pm f_{o{s}_{RL}}$ come into play, the linear muscle in its elongated state stores energy. Upon restitution contraction, this stored energy propels the rotation of the link $r_s$, which constitutes the first vertebra of the fish.

Furthermore, the caudal musculoskeletal structure, comprising four links (${\ell }_1,{\ell }_2,{\ell }_3,{\ell }_4$) and three passive joints (${\theta }_1,{\theta }_2,{\theta }_3$), facilitates a sequential rotary motion transmitted from link 1 to link 4. This transmission is accompanied by an incremental storage of energy in each helical spring that is serially connected. Consequently, the last link (link 4) undulates with significantly greater mechanical advantage. In summary, a single electrical pulse in any coil is sufficient to induce a pronounced undulation in the swimming motion of the robot’s skeleton.

As for the ballasting control device situated beneath the floor of the electromechanical oscillator, activation occurs only when either of two possible outputs from the artificial neural network (ANN) is detected: when $y_C$ equals ’sink’ or ’buoyancy.’ However, the fuzzy nature of these inputs results in a gradual slowing down of the fish’s undulation to its minimum speed. Additionally, both actions are independently regulated by a dedicated feedback controller overseeing the ballasting device.

Now, assuming knowledge of the input train of electrical step signals $s_{r,l}$ applied to the coils, let us derive the dynamic model of the biorobot, starting from the electromagnetic oscillator and extending to the motion transmitted to the last caudal link. Thus, as illustrated in Figure 12ab, the force f [N] of the solenoid’s magnetic field oscillator is established on either side (right, denoted as R, or left, denoted as L), in this context, A [m${ }^2$] represents the area of the solenoid’s pole. The symbol ${\mu }_o$ denotes the magnetic permeability of air, expressed as ${\mu }_o=4\pi \times {10}^{-7}$ H/m (henries per meter). Hence, the magnetic field B (measured in Teslas) at one extreme of a solenoid is approximated by: where i represents the coil current [A], N is the number of wire turns in a coil, and l denotes the coil length [m]. Furthermore, a coil’s current is described by the following linear differential equation as a function of time t (in seconds), taking into account a potential difference $\mathrm{v}$ (in volts): here, L represents the coil’s inductance (measured in henries, H) with an initial current condition denoted as $i_0$. Additionally, the coil’s induction model is formulated by:In essence, this study states that both lateral solenoids exhibit linear motion characterized by an oscillator force $f_{os}$. This force is expressed as: and due to the linear impacts of solenoids at both sides R and L (refer to Figure 12a), the first bar of the oscillator mechanism generates a torque, expressed as: it is theorized that the restitution/elongation force along the muscle is denoted as $f_m$ (R or L), with this force being transmitted from the electromechanical oscillator to the antagonistic muscle in the opposite direction (refer to Figure 12b). This implies that the force generated from the linear motion solenoid in the oscillator’s right-sided coil, denoted as $f_{o{s}_R}$, is applied at point R and subsequently reflected towards point L with an opposite direction, represented as $-f_{o{s}_L}$. Similarly, conversely from the oscillator’s left-sided solenoid, the force $f_{o{s}_L}$ is applied at point L and transmitted to point R as $-f_{o{s}_L}$. For $f_{o{s}_L}$ applied in L, Hence, the angles ${\alpha }_{R,L}$ assume significance as the forces acting along the muscles $f_m$ differ, resulting in distinct instant elongations $x_m\left(t\right)$. Consequently, the four-bar trapezoid-shaped oscillator mechanism manifests diverse inner angles, namely ${\theta }_{1,2}$, $\beta$, and ${\gamma }_{1,2}$ as illustrated in Figure 12b.

Thus, prior to deriving an analytical solution for $\alpha$, it is imperative to formulate a theoretical model for the muscle. In this study, a Hill’s model is adopted, as depicted in Figure 12a (on the right side). The model incorporates a serial element SE (overdamped), a contractile element CE (critically damped), and a parallel element PE (critically damped), each representing distinct spring-mass-damper systems.

The generalized model for the antagonistic muscle is conceptualized in terms of the restitution force, and it is expressed as:Therefore, by postulating an equivalent restitution/elongation mass $m_w$ associated with instantaneous weight-force loads w (such as due to hydrodynamic flows), the preceding model is replaced with Newton’s second law of motion, Furthermore, through the independent solution of each element within the system in terms of elongations, the SE model can be expressed as:Here, $s_{1,2}$ represent arbitrary constants representing damping amplitude. The terms ${\lambda }_{1,2}$ denote the root factors, Here, the factors ${\lambda }_{1,2}$ are expressed in relation to the damping coefficient $c_{SE}$ (in kg/s) and the elasticity coefficient $k_{SE}$ (in kg/s²).

Similarly, for the contractile element CE, its elongation is determined by:With amplitude factors $c_{1,2}$ and damping coefficient $c_{CE},$ a similar expression is obtained for the parallel element PE: with amplitude factors $p_{1,2}$ and damping coefficient $c_{PE},$ the next step involves substituting these functional forms into the general muscle model, such that the complete muscle’s force model $f_m$ formulated by Subsequently, simplifying the preceding expression leads to the formulation presented in Proposition 7.

Thus, without loss of generality, considering a muscle model characterized by elongation $x_m$ and a force-based model $f_m$, we proceed to derive the passive angles of the oscillator and the output forces $f_x$ and $f_y$ for ${\ell }_1$.

Under initial conditions, the trapezoid oscillator bars are assumed to have ${\theta }_0=0^o$, aligning the four-bar mechanism with the X axis. As the bars rotate by an angle ${\theta }_1$ due to solenoid impacts at points R or L, the input bar of the oscillator with a radius of $r_{os}$ undergoes an arc displacement $s_1$. Simultaneously, the output bar of shorter radius $r_s$ experiences a displacement rate of $s_2$, such that: the arc displacement at points R or L is given by $s_1=r_{os}{\theta }_{os}$. Consequently, the rotation angle of the input oscillator is expressed as:Therefore, by formulating this relationship in the context of forces and subsequently substituting the newly introduced functions, the resulting expression is Here, $\ddot{y}$ denotes the linear acceleration of either point R or L along the robot’s Y axis. By replacing the solenoid’s mass-force formulation, Hence, the functional expression for $s_1$ takes the form Without loss of generality, the inner angle ${\theta }_1$ of the oscillator mechanism (refer to Figure 12b) is derived as:Initially, when the oscillator bars are aligned with respect to the X axis, an angular displacement denoted by ${\theta }_1$ occurs as a result of the transfer of motion from the solenoid’s tangential linear motion to the input bar. Similarly, in the output bar, the corresponding angular displacement is represented by ${\theta }_2$, here, ${\Delta }_{\theta }$ signifies a minute variation resulting from motion perturbation along the various links of the caudal spine. The selection of the ± operator depends on the robot’s side, whether it is denoted as R or L. As part of the analysis strategy, the four-bar oscillator was geometrically simplified to half a trapezoid for the purpose of streamlining deductions (refer to Figure 12b). Within this reduced mechanism, two triangles emerge. One triangle is defined by the parameters $r_{os}$, ℓ, d, while the other is characterized by $x_m\left(t\right),\ell ,r_s$, where ℓ serves as the hypotenuse and the sides d, $r_{os}$, and $r_s$ remain constant. Consequently, the instantaneous length of the hypotenuse is deduced as follows:Upon determining the value of ℓ, the inner angle ${\gamma }_1$ can be derived as follows:Therefore, by isolating ${\gamma }_1$, Until this point, given the knowledge of ${\gamma }_1$ and ${\theta }_2$, it is feasible to determine the inner complementary angle ${\gamma }_2$ through the following process:Subsequently, the angle formed by the artificial muscle and the output bar can be established according to the following principle: thus, the inner angle $\beta$ is or alternatively, an approximation of the muscle length is this is the mechanism through which the input bar transmits a force $f_{m_1}$, as defined in expression (31), from the tangent $f_{os}$ to the output bar, achieving a mechanical advantage denoted as $f_{m_2}$,

Hence, in accordance with the earlier stipulation in expression (31), Definition 5 delineates the instantaneous angles ${\alpha }_{R,L}$.

It is noteworthy that, owing to the inertial system of the robot, the longitudinal force output component $f_y$ aligns with the input force $f_{os}$ in direction. Consequently, for a right-sided force, we have ${\alpha }_R\doteq {\beta }_R-{\theta }_{1_R}$, where: and Likewise, for the left-sided ${\alpha }_L\doteq {\beta }_L-{\theta }_{1_L}$, as well as

In this scenario, an inverse solution is only applicable for $f_{x_{R,L}}$, with no necessity for determining $f_{y_{R,L}}$. Consequently, the mechanical advantage transferred between the input and output bars can be expressed by a simplified coefficient. Furthermore, through the utilization of the following trigonometric identity, can substitute and streamline the ensuing system of nonlinear equations by solving them simultaneously. Additionally, let ${\theta }_1$ be defined as:Hence, the simultaneous nonlinear system is explicitly presented solely for the force components along the X axis: and Therefore, for the numerical solution of the system, a multidimensional Newton-Raphson approach is employed as outlined in the provided solution: and

Thus, by defining all derivative terms to finalize the system, Therefore, by subsequently organizing and algebraically simplifying, and The objective is to achieve numerical proximity, aiming for ${\beta }_{t+1}\approx {\beta }_t$. Consequently, through this inverse solution, the lateral force components of the first spinal link, denoted by ${\mathbf{f}}_{\mathbf{x}}$, are intended to be estimated because they are perpendicular to the links and produce the angular moments at each passive joint.

Thus, given that the torque of the trapezoid’s second bar is ${\tau }_s=f_s{r}_s$ (see Figure 12b), we establish a torque-angular moment equivalence, denoted as ${\tau }_s\equiv M_1$. Leveraging this equivalence and the prior understanding of the torque ${\tau }_s$ acting on the second bar of the trapezoid, mechanically connected to the first link ${\ell }_1$, we affirm their shared angular moment. Consequently, the general expression for the tangential force $f_k$ applied at the end of each link ${\ell }_k$ is:Yet, considering the angular moment $M_k$ for each helical-spring joint, supporting the mass of the successive links, let’s introduce equivalent inertial moments, starting with $I_{{\epsilon }_1}=I_1+I_2+I_3+I_4$. Subsequently, we define $I_{{\epsilon }_2}=I_2+I_3+I_4$, $I_{{\epsilon }_3}=I_3+I_4$, and finally, $I_4$. Thus, in the continuum of the caudal spine, the transmission of energy to each link is contingent upon the preceding joints, as established by:Each helical spring, connecting pairs of vertebrae, undergoes an input force $f=-kx$, directly proportional to the angular spring deformation indicated by elongation x. Here, k [kg m${ }^2$/s${ }^2$] represents the stiffness coefficient. External forces result in an angular moment, given by $\tau =-k\theta$, where torque serves as an equivalent variable to angular momentum, such that $I\alpha =-k\theta$. Consequently, when expressing the formula as a linear second-order differential equation, we have:Here, ${\ddot{\theta }}_{L_k}$ represents undulatory accelerations arising from external loads or residual motions along the successive caudal links, which are detectable through encoders and IMUs. Assuming an angular frequency ${\omega }^2=k/I$, a period $p=2\pi \sqrt{I/k}$, and moments of inertia expressed as $I_k=r_k^2{m}_k$, the general equation is formulated as follows: where ${\ddot{\theta }}_k$ is replaced by the helical spring expression (69) to derive By algebraically extending, omitting terms, and rearranging for all links in the caudal spine, we arrive at the following matrix-form equation:Hence, in accordance to expression (67), the tangential forces exerted on all the caudal links of the robotic fish are delineated by the following expression:Alternatively, the last expression can be denoted as the following control law:Where $\mathbf{f}={(f_{l1},f_{l2},f_{l3},f_{l4})}^{\top }$, $\mathbf{M}$ represents masses dispersion, and $\ddot{\mathbf{\theta }}={({\ddot{\theta }}_{L_1},{\ddot{\theta }}_{L_2},{\ddot{\theta }}_{L_3},{\ddot{\theta }}_{L_4})}^{\top }$ denotes the vector of angular accelerations for the caudal vertebrae, including external loads. Additionally, $\mathbf{Q}$ stands for the matrix of stiffness coefficients. Therefore, the inverse dynamics control law, presented in a recursive form, is:Finally, for feedback control, both equations are simultaneously employed within a computational recursive scheme, and angular observations are frequently derived from sensors on both joints: encoders and IMUs.