Harmonic Noise-Tolerant ZNN for Dynamic Matrix Pseudoinversion and Its Application to Robot Manipulator.

As we know, harmonic noises widely exist in industrial fields and have a crucial impact on the computational accuracy of the zeroing neural network (ZNN) model. For tackling this issue, by combining the dynamics of harmonic signals, two harmonic noise-tolerant ZNN (HNTZNN) models are designed for the dynamic matrix pseudoinversion. In the design of HNTZNN models, an adaptive compensation term is adopted to eliminate the influence of harmonic noises, and a Li activation function is introduced to further improve the convergence rate. The convergence and robustness to harmonic noises of the proposed HNTZNN models are proved through theoretical analyses. Besides, compared with the ZNN model without adaptive compensation term, the HNTZNN models are more effective for tacking the problem of dynamic matrix pseudoinverse under harmonic noises environments. Moreover, HNTZNN models are further applied to the kinematic control of a four-link planar robot manipulator under harmonic noises. In general, the experimental results verify the effectiveness, superiority, and broad application prospect of the models.

1. Introduction

The matrix pseudoinverse (i.e., Moore Penrose generalized inverse) is the generalized inverse of a singular matrix or non-square matrix. Similar to the inverse of the matrix, the matrix pseudoinverse is an important subject in the fields of science and engineering. It is considered a powerful formula and design technology for image reduction (Juang and Wu, 2010), signal processing (Van der Veen et al., 1997), linear classifier (Skurichina and Duin, 2002), and associative memory (Zhang et al., 2004). Due to its importance, numerous efforts have been made in the calculation of matrix pseudoinverse in the past few decades. For instance, Perković and Stanimirović (2011) developed an iterative algorithm to estimate the Moore Penrose generalized inverse. Zhou et al. (2002) presented a sequential recursive formula to calculate the pseudo inverse of the matrix based on the famous Greville formula. Courrieu (2008) proposed an algorithm based on full rank Cholesky decomposition for fast calculation of Moore Penrose inverse matrix. However, most of these serial processing algorithms need to be executed in a single sampling period, and when the system order becomes large, these algorithms may fail. Therefore, these studies are neither feasible nor valid for real-time applications (Zhang and Yi, 2011). While the real-time computing of dynamic matrix pseudoinverse problems often exists in industrial applications. For example, the inverse kinematics control problem of a redundant robot manipulator needs to calculate the pseudoinverse of the dynamic matrix in real time (Liao et al., 2016).

Due to the nature of high-speed parallel, distributed processing and the convenience of hardware implementation, neural networks have been used to replace the numerical algorithm of serial processing, and widely used in the fields of scientific and engineering (Jin et al., 2017; Xu et al., 2019; Jafari, 2020; Liu et al., 2020). Scholars have deeply studied many recurrent neural networks and parallel computing schemes for computing scientific problems including matrix pseudoinverse in real time (Li et al., 2021; Kornilova et al., 2022). In particular, a special kind of recurrent neural network termed zeroing neural network (ZNN), has achieved remarkable success in solving dynamic problems (Zhang et al., 2002; Ding et al., 2017; Xiao et al., 2017; Jin, 2021). Ding et al. (2017) investigated a new neural network based on ZNN to obtain the online solution of complex-valued systems of linear equation in a complex domain with higher precision and higher convergence rate. Xiao et al. (2017) presented a ZNN to solve effectively the unified QP problem and applied it to the coordinated path tracking of dual robot manipulators. Specifically, ZNN successfully solved the lag error generated by traditional recurrent neural networks which uses the indefinite error monitoring function. In addition, the traditional ZNN does not have obvious advantages in real-time calculation of large-scale problems due to its slow convergence speed. Thus, based on the essential nonlinear method of ZNN design, Liao and Zhang (2014b) proposed two finite time convergent ZNN models to solve the real time dynamic matrix pseudoinverse effectively.

In industrial applications, harmonic noises exist widely (Du et al., 2017; Karsli and Dondurur, 2018). In addition, any type of signal can be expressed as the superposition of a sinusoidal signal based on the Fourier transform. However, most of the current studies are aimed at constant noises or linear noises, and there is a lack of an adaptive mechanism to suppress harmonic noises. Furthermore, high frequency and large amplitude harmonic noises will seriously affect the calculation accuracy of the ZNN model. Therefore, by incorporating the dynamics of harmonic signals, we design and use the adaptive compensation term to learn harmonic noises and compensate for the influence of harmonic noises adaptively. Moreover, a well-defined activation function is helpful to accelerate the convergence rate of the ZNN model (Xiao et al., 2018). Therefore, on the basis of harmonic noise adaptation, we employ the Li activation function (Li et al., 2013) to speed up the convergence rate of ZNN and finally design two harmonic noise-tolerant ZNN (HNTZNN) models for solving the dynamic matrix pseudoinverse problem in the harmonic-noise environment. Furthermore, these HNTZNN models are applied to the kinematic control of a four-link planar robot manipulator.

The remainder of this article is organized into six sections. Section 2 introduces the problem formulation and preliminaries. HNTZNN models are developed and studied in Section 3. In Section 4, the convergence and noise suppression ability of HNTZNN are proved theoretically. In Section 5, numerical examples are used to prove the efficacy and superiority of HNTZNN models for calculating the dynamic matrix pseudoinverse. Section 6 shows the application of two ZNN models on a four-link planar robot manipulator. Finally, the conclusion of this article is given in Section 7. The main contributions of this article are as follows.

In this article, two novel harmonic noise-tolerant ZNN models with a fast convergence rate are first proposed and investigated for computing dynamic matrix pseudoinversion by combining the dynamic properties of harmonic signals and the Li activation function.

Theoretical analyses are conducted, which deduce the excellent convergence and robustness of these HNTZNN models in coping with single-harmonic and multiple-harmonic noises.

The experiment results are illustrated, which further substantiate the efficacy and superiority of the proposed HNTZNN models for the dynamic matrix pseudoinversion under low-frequency, high-frequency, periodic, and aperiodic harmonic noises.

The proposed HNTZNN models are applied to the kinematic control of a four-link planar robot manipulator in the presence of harmonic noise, thereby depicting the application prospect of the proposed models.

2. Problem Formulation and Preliminaries

Definition 1. (Wang, B(t) ∈ ℝ, its pseudoinverse Y(t) ∈ ℝ satisfies the following Penrose equations:

where (·)T denotes the transpose of a matrix, Y(t) is called the dynamic pseudoinverse of B(t), which is often denoted by B†(t).

It is worth noting that the dynamic pseudoinverse B†(t) always exists and is unique (Liao and Zhang, 2014a). In particular, if the matrix B(t) is a full rank matrix at any time instant t, i.e., rank[B(t)] = min{m, n} ∀t ∈ [0, ∞), the dynamic pseudoinverse of B(t) is obtained by the following lemma.

Lemma 1. (Wang, B(t) ∈ ℝ, if rank[B(t)] = min{m, n} ∀t ∈ [0, ∞), the unique dynamic pseudoinverse B†(t) can be expressed as:

When m > n, it represents the left pseudoinverse of the matrix B(t). When m ≤ n, it represents the right pseudoinverse of the matrix B(t). In the case of m ≤ n, the procedure of obtaining the dynamic pseudoinverse B†(t) is similar to that of m > n, and thus this paper only considers the case of the left pseudoinverse of the matrix B(t). Therefore, the dynamic pseudoinverse problem studied in this study is in the form of

Our objective in this study is to find the Y(t) of the problem (1).

3. Harmonic Noise-Tolerant ZNN Models

In this section, the improved ZNN models against harmonic noises for dynamic pseudoinverse are detailed.

3.1. ZNN Model Against Single-Harmonic Noise

In order to monitor the calculation process of dynamic matrix pseudoinverse (1) in real time, the error function is defined as

To force E(t) to converge to zero, we choose

where τ > 0 ∈ ℝ is a positive design parameter, and Ψ(·) denotes the activation function of the neural network, which is a monotonically increasing odd function. Then,

In this article, we choose the Li activation function (Li et al., 2013), and its expression is as follows:

where e represents the ijth element of E(t), parameter λ ∈ (0, 1), and the function Lipλ(·) is defined as follows:

In view of the noise exiting, formula (2) can be rewritten as

where O(t) ∈ ℝ represents the matrix-form harmonic noises and each element is defined as follows:

where i ∈ {1, ⋯ , m}, j ∈ {1, ⋯ , n}. It is assumed that the frequency f in (7) is known, the amplitude Γ and phase φ are unknown. To adaptively learn the f-frequency harmonics noise, an additional quantity −X(t) ∈ ℝ is introduced into equation (6):

The purpose of this study is to design and use −X(t) to adaptively learn O(t) and compensate for the impact of O(t). Differentiating (7) twice with respect to time t, the following result is obtained:

Thus, the dynamics of o(t) can be formulated as follows:

Rewrite formula (9) into matrix form:

where N(t) ∈ ℝ with n(t) as its element. The unknown amplitude and phase information of a harmonic signal can be adaptively eliminated from (10) when the frequency f is known.

Therefore, combing formulas (3), (8), and (10), the following single-harmonic noise-tolerant ZNN (HNTZNN) model is obtained:

where λ ≥ 1 ∈ ℝ.

3.2. ZNN Model Against Multiple-Harmonic Noise

In this section, the proposed HNTZNN model (11) is extended to multiple-harmonic noises, and a neural network model against multiple-harmonic noises is obtained. Let Õ(t) ∈ ℝ be the matrix-form multiple-harmonic noises. Each element of such matrix is defined as

Similarly, in (12), it is assumed that the frequency f (with k ∈ 1, ⋯ , l) is known, the amplitude Γ and phase φ are unknown.

Differentiating (12) twice with respect to time t, we have the following results:

The dynamics of õj(t) are obtained as follows:

which is rewritten in the matrix form as:

where and are the matrices composed of n(t) and o(t), respectively.

Combined with the dynamic analysis in Formula (13), a multiple-harmonic noise-tolerant ZNN model is as follows:

When l = 1, (14) can be reduced to (11).

The corresponding theoretical analysis and proof are carried out in the next section.

4. Theoretical Analysis

Three theorems are proposed to verify the convergence and robustness performance of the proposed HNTZNN model.

Theorem 1. For any smooth dynamic full rank matrix B(t) ∈ ℝ (m > n), the state matrix Y(t) of HNTZNN model (11) globally converges to the dynamic theoretical pseudoinverse B†(t) of (1) in the free-noise situation.

Proof: In the free-noise situation, (11) can be simplified as

where E(t) = BT(t)B(t)Y(t) − BT(t), (15) is a compact matrix form of the following set of mn-decoupled equations:

Then, we could define the Lyapunov function candidate (Xiang et al., 2019) to analyze the ijth subsystem (16) as

which guarantees v(t) ≥ 0, i.e., v(t) > 0 for any or or , and v(t) = 0 if and only if , , . The time derivative of (17) could be obtained as

Substituting (5) into the above equation, we further have

Thus, when e(t) ≠ 0, , when e(t) = 0, . According to Lyapunov stability theory (Khalil, 2001), it can be concluded that e(t) globally to 0 for any i ∈ {1, ⋯ , m}, j ∈ {1, ⋯ , n}, which means that

In view of E(t) = BT(t)B(t)Y(t) − BT(t) and the nonsingularity of BT(t)B(t), we can deduce Y(t) → (BT(t)B(t))−1BT(t) as t → ∞. Therefore, the state matrix Y(t) of HNTZNN model (11) globally converges to the dynamic theoretical pseudoinverse B†(t) of (1). The proof is, thus, completed.

Theorem 2. For any smooth dynamic full rank matrix B(t) ∈ ℝ (m > n), the state matrix Y(t) of HNTZNN model (11) globally converges to the dynamic theoretical pseudoinverse B†(t) of (1) when the single-harmonic noise is considered.

Proof: Let us consider a single harmonic noise O(t) = [o(t)] = [Γsin(2πft + φ)] with known frequency, and the unknown amplitude and phase. Model (11) is simplified and rewritten as

Defining Δ(t) = X(t) − O(t) and Φ(t) = C(t) − N(t) yields the following reformulation of (18):

Noting that Equations (15) and (19) have exactly the same form, therefore, one can also obtain

for any i ∈ {1, ⋯ , m}, j ∈ {1, ⋯ , n}. In addition, according to (16), we have

According to the dynamic description in (9), Equation (20) can generate a harmonic signal x(t) that can adaptively compensate for the influence of noise over time. Then, Equation (18) is a compact matrix form of the following set of mn-decoupled equations:

Because of , the harmonic signal x(t) is generated based on the analysis of (20), which meets the following results:

In particular, although only the frequency of o(t) is known, the second and third dynamics in (21) can adaptively predict the unknown parameters of o(t), i.e., amplitude Γ and phase φ. The resultant signal x(t) can adaptively compensate the impact of o(t).

In summary, the state matrix Y(t) of the HNTZNN model (11) globally converges to the dynamic theoretical pseudoinverse B†(t) of (1) when the single-harmonic noise is considered. The proof is, thus, completed.

Theorem 3. For any smooth dynamic full rank matrix B(t) ∈ ℝ (m > n), the state matrix Y(t) of the HNTZNN model (14) globally converges to the dynamic theoretical pseudoinverse B†(t) of (1) when the multiple-harmonic noise is considered, i.e., , i ∈ {1, ⋯ , m}, j ∈ {1, ⋯ , n}.

Proof: This theorem can be obtained through the proof of the above two theorems and the superposition principle. Therefore, it is omitted. The proof is, thus, completed.

Remarks: Generally, for any periodic noises which can be decomposed into a series of harmonics through Fourier transform, the proposed model is effective. The limitation of the current work may lie in the suppression of random noises.

5. Numerical Verifications

Previously, we analyzed the convergence and noise tolerance of models (11) and (14). In this section, we verify the efficacy and the superiority of the proposed models for solving the following dynamic matrix pseudoinverse:

To check the correctness of the presented HNTZNN models, the theoretical dynamic matrix pseudoinverse of (22) is directly given as

In order to verify the superiority of the proposed model, the ZNN-1 model (Liao and Zhang, 2014a) for dynamic matrix pseudoinverse is directly given in this section:

The following experiments are carried out to solve the pseudoinverse of the dynamic matrix (22) in the case of four kinds of harmonic noise.

Single-harmonic noise with a low-frequency: Figures 1, 2 present the experimental result synthesized by HNTZNN model (11) and ZNN-1 model (24) for the dynamic matrix (21) pseudoinverse under the single-harmonic noise with a near-zero frequency, i.e.,

Figure 1

Results synthesized by the HNTZNN model (11) using τ = 10 and λ = 1 for computing the pseudoinverse of the dynamic matrix (21) under the single-harmonic noise with frequency f = 0.001 Hz. (A) State trajectory. (B) Residual error.

Figure 2

Residual errors of ZNN-1 model (24) for computing the pseudoinverse of the dynamic matrix (21) under the single-harmonic noise with different frequencies. (A) f = 0.001 Hz. (B) f = 50 Hz.

where frequency f = 0.002π/(2π) = 0.001Hz. In Figure 1A, the red dashed-dotted line denotes the theoretical dynamic solution, and the blue solid line denotes the neural-state solution. In addition, we set τ = 10 and λ = 1. As shown in Figure 1A, for a random initial value Y(0) ∈ ℝ2×3, state matrix Y(t) ∈ ℝ2×3 of the proposed HNTZNN model (11) converges to the theoretical pseudoinverse (23) accurately and rapidly in a short period of time. Figure 1B shows that the residual error ‖E(t)‖F of the HNTZNN model (11) can converge to zero within 1 s. However, with the same parameters, the residual error ‖E(t)‖F of the ZNN-1 model (24) can not converge to zero displayed in Figure 2A. These experimental results verify the efficacy and superiority of model (11) for solving the dynamic matrix pseudoinverse under the single-harmonic noise.

Single-harmonic noise with high-frequency: HNTZNN model (11) and ZNN-1 model (24) are tested for the following harmonic noise:

where frequency f = 50 Hz, which is widely encountered in power systems (Karsli and Dondurur, 2018). The corresponding results are shown in Figures 2, 3. As shown in Figure 3A, the state matrix Y(t) ∈ ℝ2×3 of the proposed HNTZNN model (11) still can converge to the dynamic theoretical pseudoinverse (23) accurately and rapidly. Figure 3B shows that the residual error ‖E(t)‖F of HNTZNN model (11) can converge to zero rapidly. On the contrary, Figure 2B shows that the residual error ‖E(t)‖F of ZNN-1 model (24) has a large value under harmonic noise and cannot converge to zero.

Figure 3

Results synthesized by the HNTZNN model (11) using τ = 10 and λ = 1 for computing the pseudoinverse of the dynamic matrix (21) under the single-harmonic noise with frequency f = 50 Hz. (A) State trajectory. (B) Residual error.

For further investigation, let λ = 1 and vary the value of τ (i.e., let τ = 50 and τ = 100). The experimental results are shown in Figure 4. It can be seen from Figures 3B, 4 that the larger value of τ, the faster convergence rate of HNTZNN model (11).

Figure 4

Residual errors of HNTZNN model (11) using different τ for computing the pseudoinverse of the dynamic matrix (21) under the single-harmonic noise with frequency f = 50 Hz. (A) τ = 50, λ = 1. (B) τ = 100, λ = 1.

Multiple-harmonic noise with a periodic characteristic: HNTZNN model (14) is tested for the following harmonic noise:

where frequencies are denoted by f1 = 1, f2 = 2, f3 = 5 Hz. As shown in Figure 5, the state matrix Y(t) of proposed HNTZNN model (14) converges to the dynamic theoretical pseudoinverse (23), and the residual error ‖E(t)‖F of HNTZNN model (11) is convergent to zero.

Figure 5

Results synthesized by the HNTZNN model (14) using τ = 10 and λ = 1 for computing the pseudoinverse of the dynamic matrix (21) under the multiple-harmonic noise with a periodic characteristic. (A) State trajectory. (B) Residual error.

Multiple-harmonic noise with an aperiodic characteristic: HNTZNN model (14) is tested for the following harmonic noise:

where frequencies are denoted by f1 = 0Hz, f2 = 2.5/πHz, f3 = 5Hz. The corresponding results are demonstrated in Figure 6. As shown in Figure 6, the state matrix Y(t) of proposed HNTZNN model (14) converges to the dynamic theoretical pseudoinverse (23), and the residual error ‖E(t)‖F of HNTZNN model (11) is convergent to zero.

Figure 6

Results synthesized by the HNTZNN model (14) using τ = 10 and λ = 1 for computing the pseudoinverse of the dynamic matrix (21) under the multiple-harmonic noise with an aperiodic characteristic. (A) State trajectory. (B) Residual error.

In summary, the experimental results verify the effectiveness and superiority of HNTZNN models (11) and (14) for solving the dynamic matrix pseudoinverse in the presence of harmonic noises.

6. Robot Manipulator Application

In this section, HNTZNN models (11) and (14) are applied to the kinematic control of a four-link planar robot manipulator under the presence of harmonic noises.

6.1. ZNN-Combined Kinematic Control

The geometry of the four-link planar robot manipulator is shown in Guo et al. (2018a). The Jacobian matrix J(θ(t)) ∈ ℝ2×4 of this robot manipulator is:

where θ(t) ∈ ℝ4 denotes the joint-angle vector, l (with i = 1, 2, 3, 4) denotes the length of the ith link, , and . The redundancy-resolution problem of such a robot manipulator, which corresponds to the kinematic control, is described as follows. Given the desired end-effector path , we need to obtain θ(t) in real time t. In particular, the following acceleration-level redundancy-resolution problem is formulated by adopting the solutions presented in previous work (Siciliano et al., 2009; Guo et al., 2018b):

where is the joint-acceleration vector, denotes actual end-effector path and with as the joint-velocity vector and as the time derivative of the Jacobian matrix J(t). In addition, k1 and k2 > 0 are the feedback gains, ψ(·) is a differentiable nonlinear mapping (Siciliano et al., 2009; Guo et al., 2018b), and ṙd(t) and are first-order and second-order time derivatives of rd(t), respectively. By effectively solving the acceleration-level kinematic equation in (25), the purpose of kinematic control for robot manipulator is achieved.

According to Siciliano et al. (2009), the following pseudoinverse-type (P-type) scheme, which is a typical solution to (25), is presented for the four-link planar robot manipulator:

The dynamic matrix inversion is integrated into (26) to achieve the kinematic control of the four-link planar robot manipulator. By defining B(t) = J(t), combining (11) with (26) yields the following the new ZNN-combined kinematic control method that is tolerant to single-harmonic noise:

Furthermore, by combining (14) with (26), the new ZNN-combined kinematic control method for multiple-harmonic noise can be obtained as follows:

To provide a better understanding, we show Figure 7 to depict the block diagram of the ZNN-combined kinematic control method in (28) [including (27) as a special case] for the four-link planar robot manipulator. The robot application in this study is clearly limited by the harmonic noise associated with the ZNN model instead of that of the robot manipulator.

Figure 7

Block diagram of the new ZNN-combined kinematic control method for the four-link planar root manipulator.

6.2. Simulation Verification

In this subsection, the simulation results are presented to show the effectiveness of the presented ZNN-combined kinematic control methods in (27) and (28) for the four-link planar robot manipulator. In the simulations, the lengths of the robot's links are l1 = l2 = l3 = l4 = 1 m, and the initial joint state is set to θ(0) = [π/9, π/12, π/12, π/12]T rad.

Tricuspid path-tracking example: In this example, the presented ZNN-combined kinematic control methods in (27) and (28) are simulated for the four-link planar robot manipulator with its end-effector tracking a tricuspid path.

★ Single-harmonic noise: The ZNN-combined kinematic control method in (27) is investigated for the single-harmonic noise of , in which the harmonic frequency is f = 2π/(2π) = 1 Hz. For comparison, the following noise-polluted kinematic control method based on the original ZNN design in (2) is also simulated:

Figure 8 presents the simulation results synthesized by (29) using τ = 10. The simulated end-effector trajectory of the robot manipulator does not match the desired tricuspid path, which means that (29) has failed to achieve the kinematic control of the robot manipulator due to the existence of harmonic noise. As shown in Figures 9A,B, the simulated end-effector trajectory is very close to the desired tricuspid path. Figure 9C shows that the maximal end-effector positioning error is <2 × 10−7 m, suggesting that the kinematic control is achieved successfully via (27) even if harmonic noise exists. These comparative results substantiate the effectiveness and noise-suppression capability of the ZNN-combined kinematic control method in (27) for the four-link planar robot manipulator.

Figure 8

End-effector of the four-link planar robot manipulator tracking the tricuspid path synthesized by the noise-polluted kinematic control method in (29) using τ = 10, where single-harmonic noise is considered. (A) Simulated motion trajectories. (B) Path and trajectory.

Figure 9

End-effector of the four-link planar robot manipulator tracking the tricuspid path synthesized by the ZNN-combined kinematic control method in (27) using τ = 10 and λ = 1, where single-harmonic noise is considered. (A) Simulated motion trajectories. (B) Path and trajectory. (C) End-effector positioning errors.

★ Multiple-harmonic noise: The ZNN-combined kinematic control method in (28) is investigated for the multiple-harmonic noise:

in which the harmonic frequencies are denoted by f1 = 1, f2 = 5, and f3 = 10 Hz. The corresponding simulation results synthesized by (28) are presented in Figure 10. The simulated end-effector trajectory is close to the desired tricuspid path, with the maximal positioning error being <5 × 10−6 m. Owing to the robustness of (28), this method is tolerant of multiple-harmonic noise, and the robot manipulator can track the desired tricuspid path successfully. The effectiveness of the ZNN-combined kinematic control method in (28) for the four-link planar robot manipulator is, thus, demonstrated.

Figure 10

End-effector of the four-link planar robot manipulator tracking the tricuspid path synthesized by the ZNN-combined kinematic control method in (28) using τ = 10 and λ = 1, where multiple-harmonic noise is considered. (A) Simulated motion trajectories. (B) Path and trajectory. (C) End-effector positioning errors.

Circular path-tracking example: In this example, the presented ZNN-combined kinematic control methods in (27) and (28) are simulated for the four-link planar robot manipulator with its end-effector tracking a circular path. Similarly, the following two harmonic noises are considered in the investigation of (27) and (28):

The frequency of the single-harmonic noise approximates the bandwidth of the robot manipulator. The corresponding simulation results are illustrated in Figures 11, 12.

Figure 11

End-effector of the four-link planar robot manipulator tracking the circular path synthesized by the ZNN-combined kinematic control method in (28) using τ = 10 and λ = 1, where single-harmonic noise is considered. (A) Simulated motion trajectories. (B) Path and trajectory. (C) End-effector positioning errors.

Figure 12

End-effector of the four-link planar robot manipulator tracking the circular path synthesized by the ZNN-combined kinematic control method in (28) using τ = 10 and λ = 1, where multiple-harmonic noise is considered. (A) Simulated motion trajectories. (B) Path and trajectory. (C) End-effector positioning errors.

As shown in Figure 11, the simulated end-effector trajectory of the robot manipulator matches the desired circular path, and the maximal positioning error is less than 4 × 10−7 m. These results substantiate the successful application of the ZNN-combined kinematic control method in (27) to the robot manipulator. The results presented in Figure 11 also show that (27) has achieved robustness against single-harmonic noise, even if the noise frequency is nearly the same as the robot's bandwidth. Furthermore, the results presented in Figure 12 suggest that the path-tracking task has been successfully completed even with the presence of multiple-harmonic noise. The effectiveness of the ZNN-combined kinematic control method in (28) for the robot manipulator is further demonstrated.

In summary, the simulation results (i.e., Figures 8–12) prove that the presented ZNN-combined kinematic control methods in (27) and (28) are more effective for the robot manipulator under the presence of harmonic noise compared with the noise-polluted the kinematic control method in (29). More importantly, the results indicate the application prospect of the proposed HNTZNN models in (11) and (14) via dynamic matrix pseudoinverse.

7. Conclusion

In this article, by combining the harmonic signal dynamics and the Li activation function, two harmonic noises tolerant ZNN models with fast convergence have been proposed. The theoretical analysis demonstrates that the proposed HNTZNN models possess satisfactory convergence and robustness properties. Furthermore, an example has been given to verify its effectiveness in solving the dynamic matrix pseudoinverse under different harmonic noises. Besides, such two HNTZNN models have been applied to the kinematic control of the four-link planar robot manipulator. In summary, the corresponding experimental results substantiate the effectiveness, superiority, and application prospect of the proposed HNTZNN models even with the existence of harmonic noises. Future study may lie in studying the proposed ZNN models activated by different activation functions. Another future research direction is to extend the proposed ZNN models to other application fields.

Data Availability Statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author Contributions

BL and DG presented the scheme. BL and YW designed experiments and wrote the manuscript. YW and JL carried out experiments. JL and YH analyzed the experimental results. All authors contributed to the article and approved the submitted version.

Funding

This study was supported in part by the National Natural Science Foundation of China (62066015 and 61962023), the Hunan Natural Science Foundation of China (2020JJ4511), the Research Foundation of Education Bureau of Hunan Province, China (20A396), and the Scientific Research Foundation of Jishou University (Jdy20063).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.