Research on the recognition model of exercise fatigue based on the fusion of sEMG and ECG signals

Summary This study significantly enhances the accuracy of exercise state identification in wearable devices through improved denoising techniques for sEMG and ECG signals. By adopting an optimized Variational Mode Decomposition (VMD) method, combined with the Improved Sparrow Search Algorithm and Second Generation Wavelet Transform (ISSA-VMD-SWT), and introducing chaos mapping to strengthen the algorithm’s initial population, this approach effectively reduces noise while preserving key fatigue-related features. In tests conducted on data from 32 participants, the method achieved accuracy rates of 93.25%, 95.16%, and 93.05% for identifying “Easy,” “Transition,” and “Tired” exercise states, respectively, showing significant advantages over traditional denoising techniques. These results indicate that the denoising technology developed in this study represents a significant technological advancement for the application of ECG and sEMG fatigue identification technologies in wearable health monitoring devices.


INTRODUCTION
Muscle fatigue is defined as the decrease in the efficiency of action potential propagation between the sarcolemma and muscle fibers, resulting from the accumulation of lactate, hydrogen ions, and inorganic phosphates in the blood during exercise, leading to a decrease in muscle pH, maximum voluntary contraction capacity, and maximum power output of the muscle. 1 Accurate identification of muscle fatigue status is of great significance in preventing sports injuries and has important implications in the fields of kinematics and rehabilitation medicine.Non-invasive detection techniques can be employed for precise identification of muscle fatigue, offering broad application prospects. 2hese techniques primarily include surface electromyography (sEMG), 3 sonomyography (SMG), near-infrared spectroscopy (NIRS) analysis, 4 and Mechanomyography (MMG). 5In recent years, there has been increasing attention from researchers on the combined use of sEMG and electrocardiogram (ECG) for the identification of muscle fatigue status.These two signals respectively represent the superficial EMG of the skin and the sustained depolarization of myocardial cells, both containing rich information on human neural function.The combination of these two methods can provide a highly accurate and robust approach for identifying muscle fatigue status.
When applying these methods, the first challenge is to address the denoising issue in weak bioelectrical signals of the human body.Noise in sEMG and ECG signals is considered as one of the crucial factors that affect the recognition performance.Chatterjee et al. emphasized that denoising is a key preprocessing step. 6The sources of noise mainly include baseline drift, power line interference, muscle artifacts, and channel noise.Several researchers have conducted studies on denoising weak bioelectrical signals, employing methods such as empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD), improved empirical mode decomposition, and wavelet-based digital filtering.In comparison to these methods, the more recent denoising technique, variational mode decomposition (VMD), can overcome drawbacks like mode mixing and endpoint effects, making it one of the recognized effective decomposition methods for weak bioelectrical signals.Xiao et al. proposed a method that combines wavelet soft thresholding with VMD for denoising sEMG signals, where the key parameters of VMD were selected using an empirical approach.Comparative results demonstrated that this method outperformed Empirical Mode Decomposition and wavelet-based methods, highlighting the importance of parameter selection in VMD for achieving superior model performance. 7However, the selection of VMD parameters relies heavily on subjective judgment, making it less user-friendly for non-experts or sports practitioners.Therefore, many researchers have devoted significant efforts to introducing adaptive mechanisms for parameter selection in VMD.
The denoising efficacy of VMD heavily relies on the choice of penalty factors a and decomposition levels k.Swarm intelligence optimization algorithms, like Beetle Antennae Search, Gray Wolf, and Cuckoo Search, have been employed to optimize VMD parameters.Wang et al. proposed an improved VMD method based on the Beetle Antennae Search algorithm, where the kurtosis of intrinsic mode functions was used ll OPEN ACCESS as the fitness function during the search process. 8Gai et al. optimized VMD parameters using a hybrid Gray Wolf algorithm, significantly improving the optimization speed. 9Yan et al. optimized VMD parameters using the Cuckoo Search algorithm. 10Numerous research findings have shown that swarm intelligence parameter optimization algorithms can indeed introduce adaptive mechanisms for parameter selection in VMD.However, swarm intelligence algorithms have drawbacks such as a tendency to get stuck in local optima and slow convergence speed, leading to different optimization results and a lack of unified denoising effects when using these algorithms to optimize the penalty factor a and k.Prior to this study, there have been few research efforts on the improvement of swarm intelligence optimization algorithms in the field of VMD for denoising muscle signals.
In this paper, the fast swarm intelligent optimization algorithm is applied to the parameter optimization of denoising algorithm, which proposes an improved denoising method called Improved Sparrow Search algorithm to optimize VMD combined with Second Generation Wavelet Transform (ISSA-VMD-SWT).Through the comparison of denoising effects on simulated noisy ECG and sEMG signals, analysis of denoising performance on actual collected signals, and the analysis of the improvement in fatigue recognition performance by combining the denoising algorithm with sEMG and ECG, it is demonstrated that ISSA-VMD-SWT can effectively handle both local and global characteristics of the signals, making it one of the suitable methods for preprocessing muscle signals.The overall arrangement of this paper is as follows: the second part presents the materials and methods, including data collection and a detailed description of ISSA-VMD-SWT.The third part compares the denoising effects on simulated noisy signals and real collected signals, as well as the analysis of the improvement in muscle fatigue recognition using the proposed denoising method.The fourth part presents the conclusions.This paper will provide a reliable and efficient denoising method for the joint use of sEMG and ECG, supporting the development and promotion of wearable devices in the field of sports science in the future.

RESULT ECG signal characteristics analysis under different fatigue states
Due to the significant variations in ECG signals caused by individual and environmental differences, Heart Rate Variability (HRV) analysis is employed in this study to analyze the ECG signals.HRV provides a comprehensive set of evaluation metrics derived from the ECG.It can serve as a reference for assessing mental and physiological stress. 11This paper utilizes several time-domain and frequency-domain HRV indices.The time-domain indices include the mean of normal to normal intervals (NN.mean) and the proportion of NN50 divided by the total number of NNs (PNN50).The frequency-domain indices include the total power (TP) and the High frequency from 0.15 Hz to 0.4 Hz (HF).The characteristic parameters of HRV are shown in Table 1.
HRV exhibits different patterns under three different states.(1) For the parameter NN.mean, the average value is higher in the Easy and Transition states, but slightly lower in the Tired state.However, the variance increases, indicating greater fluctuations in NN.mean during the Tired state.(2) rMSSD is used to evaluate the degree of heart rate variation between adjacent heartbeats, reflecting the rapid or short-term changes in HRV.The calculated values for rMSSD show slight differences but with minimal changes.The values are generally close to each other.(3) PNN50 represents the percentage of adjacent RR intervals differing by more than 50 ms.In the Easy state, PNN50 has a higher value, while in the Transition and Tired states, the value is lower.The low variance also indicates less variability in this parameter.This suggests a higher stability in the autonomic nervous system regulation of the participants.(4) total spectral power gradually increases from the Easy to Tired state.This indicates that overall HRV increases with the level of muscle fatigue.(5) HF and LF values increase sequentially with increasing muscle fatigue.This suggests a greater contribution of HF and LF to HRV. (6) VLF also increases with muscle fatigue, indicating an imbalance in sympathetic and parasympathetic nervous system activity and an impact on autonomic nervous system regulation.(7) the LF/HF ratio indicates the balance between sympathetic and parasympathetic activity.When LF/HF is lower, it indicates a relatively balanced state between sympathetic and parasympathetic activity.As LF/HF increases, sympathetic activity becomes relatively enhanced while The "mean" represents the average value, "var" denotes the variance.
parasympathetic activity gradually weakens.These findings illustrate the diverse patterns of HRV parameters under different states of muscle fatigue.
Based on the comprehensive data analysis, the development of muscle fatigue is a relatively complex process with various changes.The sympathetic and parasympathetic nervous systems have a significant impact on ECG signals, and multiple HRV parameters reflect the individual's sympathetic and parasympathetic responses to muscle fatigue.However, one obvious problem is that HRV parameters vary from person to person, making it impossible to determine muscle fatigue from just one or a few parameters.

Analysis of sEMG signal characteristics under different fatigue conditions
We collected sEMG signals from the rectus femoris (RF), vastus lateralis (VL), vastus medialis (VM), and gastrocnemius (GA) muscles through experiments, and conducted an analysis on the collected data.The sEMG signals exhibited complexity, and fatigue cannot be analyzed solely based on the raw signals.Therefore, we selected commonly used sEMG signal features for analysis, including integrated electromyography (IEMG), root-mean-square (RMS), mean power frequency (MPF), and median frequency (MF).The calculation methods for these four indices can be found in ref. 12.
Figure 1 shows the sEMG signal features of RF in different fatigue states (Easy, Transition, Tired).In Figure 1A, sEMG IEMG can easily distinguish Tired and Easy states.In the Tired state, sEMG IEMG has higher values and larger variance, indicating greater fluctuations.In contrast, in the Easy state, sEMG IEMG has smaller values and less variability.Therefore, sEMG IEMG exhibits good discriminative ability and reflects the overall muscle activation level, where higher values correspond to stronger muscle stimulation.In Figure 1B, sEMG RMS is mainly used to quantify the amplitude characteristics and reflects the level of muscle contraction force.The results show that in the Tired state, sEMG RMS amplitude increases because more neural activation is required to maintain the same force output during fatigue.However, the separation of RMS between Tired and Easy states is not significant, indicating that this parameter cannot accurately reflect the degree of muscle fatigue.In Figure 1C, sEMG MPF represents the median frequency, which plays a role in reflecting muscle fatigue level and fiber type.According to the results, the MPF decreases when the body is in a fatigued state.This is because as fatigue progresses, the activation pattern of muscle fibers changes, leading to an increase in low-frequency information.In the Easy stage, the MPF value is relatively high, indicating a more normal activation pattern of muscle fibers, and the signal tends to contain higher-frequency information.The variation in MPF between individuals may differ depending on the different exercise states.In Figure 1D, sEMG MF shows a clear pattern with changes in exercise states.As the exercise continues, sEMG MF gradually decreases.This is because muscle contractions become less coordinated and the activation pattern changes.The signal shifts toward lower frequencies, resulting in a gradual decrease in sEMG MF values.
By analyzing the HRV parameters of ECG and the parameters of sEMG, such as IEMG, RMS, MPF, and MF, the body experiences exerciseinduced fatigue during physical activity.This fatigue is mainly characterized by a decrease in the stability of the sympathetic and parasympathetic nervous systems and an increase in the low-frequency component of muscle signals.These findings provide a basis for a deeper understanding of ECG and sEMG signals and the development of suitable denoising methods.

Comparison of different optimization algorithms
To evaluate the optimization effect of ISSA on VMD parameters and the advantages of the SWT in denoising, this study used sEMG signals collected from a healthy volunteer in the NinaPro electromyography database as simulated data.The data had minimal noise and served as clean signals for experimentation.Both ECG and sEMG signals contain a significant number of low-frequency components during the Tired stage, which should be preserved while removing noise in signal processing.In the simulated signals, 5 dB and 10 dB noise signals were added.Common algorithms, including genetic algorithm (GA), particle swarm optimization (POS), and conventional sparrow search algorithm (SSA), were also used as comparisons.The parameter settings for the optimization algorithms are as follows: population size of 50, iteration count of 50, optimization range for the number of modal components [1, 2000], optimization range for the penalty factor [2, 20], safety thresholds of 0.8 for SSA and ISSA, and discoverers and informers set to 10.For GA, the crossover probability and mutation probability were set to 0.85 and 0.1, respectively.For the POS algorithm, the learning factors C 1 and C 2 were set to 1.5, and the initial and final inertia weights were set to 0.8 and 0.4, respectively.The sym8 wavelet was chosen as the basis function for the SWT, and the decomposition level was set to 3. The original sEMG and ECG signals, as well as the signals with added 5 dB and 10 dB noise, are shown in Figure 2.
The noisy signals were input into four optimization algorithms.These algorithms were combined with VMD for decomposition, resulting in GA-VMD, PSO-VMD, SSA-VMD, and ISSA-VMD models.To reduce randomness during the optimization process, each model was run five times, and the average fitness value was obtained.The fitness variation curves of the different optimization algorithms at different iteration counts are shown in the Figure 3.The results show that at the first iteration, the initial values of the optimal Q f for GA-VMD, PSO-VMD, SSA-VMD, and ISSA-VMD models were slightly different, with values of 1.85, 1.88, 2.00, and 1.81.ISSA-VMD had the lowest initial value among them.These differences in initial values are related to the optimization strategies for the initial population in each algorithm.ISSA uses an improved Tent chaotic mapping to generate sequences, and then maps the sequences to different individuals to improve the distribution of the initial population.Furthermore, the Q f values of the four methods varied between 2.01 and 0.15, and the convergence was achieved in less than 30 iterations for all models.This indicates that all four models are effective and can be used to optimize the parameters of VMD.The Q f of VMD in the four models decreases as the iteration count increases.In the initial stage of the models, Q f decreases rapidly.Specifically, in the GA-VMD model, Q f decreases to 1.15 after 6 iterations and reaches convergence at 0.2 after 26 iterations.In the PSO-VMD model, Q f is 0.89 after 11 iterations and reaches 0.23 after 28 iterations, achieving convergence.In the SSA-VMD model, the initial Q f is relatively high at 2.01.After 9 iterations, Q f is 1.01, and after 20 iterations, Q f is 0.47.Finally, Q f is 0.35 after 26 iterations.In the ISSA-VMD model, Q f rapidly decreases to 0.5 after 14 iterations, and the model converges and stabilizes after 19 iterations with a final Q f value of 0.15.The ISSA-VMD model is the fastest in terms of convergence speed among the four models and exhibits the best Q f value in terms of parameter optimization.During the experiment, we found that the convergence speed and final Q f of the four models are not unique, which may be related to the grid search step size adopted by each model during the optimization process.The search mechanism of the ISSA-VMD model allows it to easily escape local optima and further optimize in the search space.

Comparison of denoising effect of denoised signal
The above comparison pertains to the parameter optimization aspect of the ISSA method.Now, we will compare the denoising performance of the established model.We will apply the ISSA-VMD method, as well as two other methods, namely the SWT with soft thresholding and the SWT with hard thresholding, to noisy ECG and sEMG signals.
Figure 4 shows the denoising effects of five methods: SWT soft thresholding, SWT hard thresholding, ISSA-VMD, ISSA-VMD with SWT soft thresholding, and ISSA-VMD with SWT hard thresholding, on sEMG signals with 5dB and 10dB noise.Ten models were  2.
According to the evaluation metrics of the different models after denoising, when using the SWT with soft thresholding or hard thresholding alone, the denoising results are not satisfactory.The soft thresholding method shows certain advantages over the hard thresholding method, as it achieves higher SNR values and lower MSE values.However, when comparing the PE values, the second-generation wavelet hard thresholding method performs better, with smaller peak errors.This is because the soft thresholding method, while removing noise and improving the signal-to-noise ratio, also reduces the peak values, resulting in increased peak errors between the reconstructed signal and the original signal.
When ISSA-VMD is used alone, the denoising effect is not ideal.Compared to the soft and hard thresholding methods of the second-generation wavelet, ISSA-VMD achieves better SNR, PE, and NCC, indicating that ISSA-VMD can improve the signal-to-noise ratio of the reconstructed signal while minimizing its impact on the waveform.However, the MSE value of ISSA-VMD is slightly higher than that of the soft and hard thresholding of the second-generation wavelet, indicating that using ISSA-VMD alone may not achieve a clean removal of noise details in the signal.The combination of ISSA-VMD and the second-generation wavelet in the soft and hard thresholding model shows a better overall performance.In terms of NCC and SNR, the difference between ISSA-VMD and the soft and hard thresholding is not significant, indicating that the waveform and signal-to-noise ratio after the SWT remain largely unchanged.In terms of PE, there is a noticeable decrease after the soft and hard thresholding, indicating effective correction of waveform distortion.In terms of MSE, ISSA-VMD-SWT-Hard achieves good results.Overall, when processing sEMG signals, using ISSA-VMD first helps protect the original signal waveform and suppress distortion in the reconstructed signal, and then using hard thresholding emphasizes the protection of signal details, leading to better noise suppression and higher similarity to the original signal.However, we also observed that the difference between soft and hard thresholding is very small, and it is difficult to determine the absolute advantage through a single experiment.Nevertheless, by combining the advantages of different methods, employing a strategy that combines multiple methods can achieve better results in denoising tasks.
When using the SWT soft and hard thresholding methods alone, the denoising performance is relatively poor, regardless of whether the noise level is 5 dB or 10 dB.The NCC is low, indicating low similarity between the reconstructed signal and the original signal.The MSE is high, indicating a large mean square error between the reconstructed signal and the original signal.The PE is also high, indicating a large peak error.This result may be because the SWT method may not handle signal details well when processing ECG signals, leading to a decrease in the quality of the reconstructed signal.
The denoising performance of the ISSA-VMD method alone is also not satisfactory.Although the NCC value is high, indicating a high similarity between the reconstructed signal and the original signal, the MSE value is still large, indicating a lower quality of the reconstructed signal.The PE value is also relatively high, indicating a large peak error.This suggests that the ISSA-VMD method may not effectively remove noise when processing ECG signals, leading to a decrease in the quality of the reconstructed signal.However, when combining ISSA-VMD with the second-generation wavelet soft and hard thresholding methods, namely the ISSA-VMD-SWT-Soft and ISSA-VMD-SWT-Hard models, the denoising performance improves.These models demonstrate better performance in terms of NCC, MSE, and PE metrics.The NCC value is relatively high, the MSE value is relatively low, and the PE value is relatively small, indicating a high similarity between the reconstructed signal and the original signal, good quality of the reconstructed signal, and small peak error.In summary, for denoising ECG signals, using the SWT soft and hard thresholding methods or the ISSA-VMD method alone yields unsatisfactory results.However, combining ISSA-VMD with the SWT soft and hard thresholding methods, especially with the soft thresholding, can achieve better denoising performance and improve the quality and similarity of the reconstructed signal.The evaluation metrics for each model are shown in Table 4.
In conclusion, for denoising sEMG signals, combining ISSA-VMD with the SWT and the hard thresholding method yields good results.For denoising ECG signals, ISSA-VMD combining with SWT the soft thresholding method achieves satisfactory results.It is important to note that these comparative results are based on experiments conducted with simulated noisy data.When applied to real signals, the effectiveness may be influenced by the characteristics of the actual signals and the types of noise present.

Comparison of real sEMG and ECG noise reduction
To validate the superiority of the ISSA-VMD and SWT algorithms in real sEMG and ECG signals, we collected sEMG and ECG signals from subjects and evaluated their performance using various metrics, as shown in Figure 6.Detailed parameters and standard deviations for each model are presented in Table S1.Concerning sEMG signal denoising, the method that combines ISSA-VMD with the SWT hard threshold technique outperforms other methods, exhibiting the highest NCC (0.988) and the lowest MSE (0.001), along with the highest SNR (24.060) and the lowest PE (0.020).It demonstrates excellent denoising performance, handling signal details better while effectively preserving useful signals.For ECG signal denoising, the method combining ISSA-VMD with SWT soft threshold technique shows better performance, with higher NCC (0.388) and lower MSE (0.527), along with higher SNR (14.622) and lower PE (0.273).These results align with denoising conclusions drawn from simulated noisy signals.

Overall performance analysis
Based on the previous research, the workflow for fatigue classification was adjusted using the combined sEMG and ECG signals.For sEMG, The ISSA-VMD-SWT method is combined with hard threshold technology, while for ECG, the ISSA-VMD-SWT method is combined with the soft thresholding technique for denoising.The classification model established using these methods is named ISSA-VMD-SWT.To compare the performance, we also established three other models: SWT-Soft, SWT-Hard, and ISSA-VMD.The accuracy of these models is shown in Figure 7.
Based on the statistical results, the models are ranked in descending order of average accuracy as ISSA-VMD-SWT, ISSA-VMD, SWT-Soft, and SWT-Hard.The larger the distribution on the Y axis of the violin plot, the higher the variance in model accuracy, indicating significant The range of recognition rates for the three states further narrowed, indicating improved adaptability of the model after denoising.However, the improvement was not significant, which may be due to the differences in the fatigue information used in the two experiments.The previous study mainly relied on self-assessment using the Rating of Perceived Exertion (RPE) scale, which may have had issues with unclear transition states.In contrast, this study employed the more scientific Ultima GX system.Additionally, a comparison was made between the average accuracies of different genders.The average accuracy for males was 82.9% and for females was 82.91%.The two-factor ANOVA for gender and fatigue status and denoising methods are shown in Tables S2 and S3 No significant difference was found in the evaluation of exercise fatigue between genders, indicating that this denoising method is applicable to both male and female athletes.

Confusion matrix evaluation
Different models established by various denoising methods correspond to confusion matrices shown in Figure 8. Precision, Recall, and F1 Score were calculated based on statistical principles, and detailed results are in Table S4.Precision, Recall, and F1 Score are crucial metrics in assessing classification algorithm quality, especially considering F1 Score, which integrates Precision and Recall, serving as a comprehensive evaluation metric.Regarding the SWT-Soft model, it exhibits higher Precision and Recall in the Easy category, indicating its strong predictive ability for this class.However, Precision and Recall are slightly lower for the Transition and Tired categories.SWT-Hard generally demonstrates slightly lower Precision, Recall, and F1 Score in the Easy, Transition, and Tired categories compared to SWT-Soft, particularly evident in the Transition and Tired categories.This suggests potential errors or biases in SWT-Hard when predicting these categories, resulting in slightly inferior performance compared to SWT-Soft.The ISSA-VMD-SWT model showcases relatively higher Precision, Recall, and F1 Score across the Easy, Transition, and Tired categories.These results indicate that these models perform well across multiple categories, especially achieving a high Precision level in the Easy category, signifying high predictive accuracy in this class.These outcomes highlight performance disparities among different models across various categories, indicating that certain models might excel in specific fatigue types.Further exploration might be necessary to understand the reasons for differences among fatigue states.

Variance analysis
Conducting two-way ANOVA on the data, where the first factor is set as the denoising method and the second factor as the fatigue state.The interaction plot for the two factors is illustrated as Figure S1.The results of two-factor ANOVA are shown in Table 3.At a significance level of p < 0.05, there exists a significant difference in the overall means between the Denoising Method and State of Fatigue.The interaction between these two factors has been investigated.Statistical analysis of differences among various groups has been carried out to deeply comprehend the impact of different factors on the study outcomes.
The denoising method significantly impacts the recognition rate (Fð3; 378Þ = 2224:87, p < 0:0001).Simultaneously, the fatigue state also significantly influences the recognition rate (Fð2; 378Þ = 20:91, p < 0:0001).There is an interaction effect between these two factors (Fð6;378Þ = 2:43, p = 0:0321).Among them, 'ISSA-VMD-SWT' exhibits a relatively higher recognition rate.The Tired state exhibits a higher recognition rate compared to "Easy" and "Transition"; this might be attributed to a sudden increase in oxygen consumption and significant muscle contractions.Figure 1 vividly displays the differences in the Tired state.Follow-up post-hoc analyses were conducted, and the results are shown in Tables S5 and S6.The q values and probabilities in the tables indicate significant differences between various denoising methods and fatigue states in many cases.Confidence interval data display considerable differences among different denoising methods in the Transition state, while the variation in confidence intervals is smaller in the Tired state.Based on the grouping description of the interaction, it can be concluded that ISSA-VMD-SWT has relatively minor statistical differences in the Transition, Easy, and Tired states.In essence, ISSA-VMD-SWT can effectively overcome muscle conditions under different fatigue states, demonstrating stronger generalization capabilities.

ISSA-VMD-SWT compared with other methods
Select two more advanced methods of human weak current signal denoising, Combines the traditional wavelet threshold denoising with wavelet digital filter threshold denoising (WDFTD) and empirical wavelet transform (EWT) and improved interval thresholding (EWT-IIT). 14ased on the fusion sEMG and ECG experimental data and ISSA-VMD-SWT compared processing signal parameters, as well as the recognition accuracy of muscle fatigue state.The differences among the three algorithms are shown in Table 5. NCC and MSE indicators show that the ISA-VMD-SWT algorithm performs well in preserving the original signal characteristics.SNR and PE display success This method successfully reduces signal noise and has little change in signal amplitude.Table S2 shows the difference in Precision, Recall and F1 Score caused by

Conclusion
This study introduces a novel fatigue recognition model by fusing sEMG and ECG signals to address the robustness against noises interfering weak electrical signals.The innovative method integrates ISSA-VMD and SWT denoising techniques.The optimization speed is accelerated by introducing the Tent chaotic map.The soft or hard threshold of SWT can effectively retain the fatigue-related features and reduce the noise level.The comprehensive evaluations involving both simulated noisy signals and real-world data demonstrate substantial advantages.Different from the GA, PSO, and SSA, the ISSA exhibits a faster convergence rate and lower fitness value.The searching mechanism of the ISSA-VMD model makes it easy to escape from local optima.The combination of ISSA-VMD with SWT soft or hard threshold achieves superior values in terms of NCC, MSE, SNR, and PE.The ISSA-VMD-SWT obtains significantly higher recognition rates for three fatigue levels.
These results indicate that the proposed model possesses higher recognition accuracy and better generalization performance.Specifically, the recognition rates for Easy, Transition, and Tired levels are 93.25%,95.16%, and 93.05%, respectively.The reduced inter-subject variance of fatigue levels highlights the enhanced generalization ability of the model.Furthermore, gender shows slight influence on fatigue classification, indicating that there is no significant difference in fatigue estimation between males and females.In conclusion, the signal denoising method based on ISSA-VMD-SWT soft and hard threshold exhibits superior performance in muscle fatigue recognition models using sEMG and ECG signals.This method holds great potential for practical applications in the fields of exercise physiology, engineering, and rehabilitation medicine.Future research of our team will focus on improving feature extraction and classification models for weak biological electrical signals.The eventual goal is to integrate the entire technology, involving signal acquisition, processing, and recognition, into embedded devices such as portable wearable devices.This will facilitate the advancement of sports science toward informatization.

Limitations of the study
Despite the promising results, this study has several limitations.First, the method's effectiveness is primarily validated in controlled environments with a limited number of participants, which may not fully represent the diverse conditions encountered in real-world applications.Second, while the ISSA-VMD-SWT approach shows significant improvements in noise reduction and fatigue recognition, its computational complexity is higher compared to simpler models, which may impact its real-time application in wearable devices.Third, the study focuses on the general population without considering specific groups that may have different physiological responses, such as athletes or individuals with health conditions.Lastly, the influence of gender on fatigue classification, although slight, suggests potential variability in the model's performance across different demographic groups, which was not extensively explored.
Future studies will focus on broadening the model's applicability through diversifying participant demographics, enhancing computational efficiency for real-time applications, tailoring models for specific populations, conducting in-depth demographic variability analysis, 82.16  7.12 10.72 to testing, participants familiarized themselves with the equipment and procedures.After a 3-min warm-up on a cycle ergometer, each participant performed incremental exercise at room temperature (26 G 3 C), starting at an initial workload of 100 W and increasing by 25 W/m.Participants maintained a cycling speed of 70-75 r/m throughout the entire testing process.The test was terminated when participants were unable to maintain 70 r/m.Each participant underwent the same test for three consecutive days, with a 24-h interval between tests.After the experiment, we analyzed the slope of the VCO 2 and VO 2 kinetics curves 17 and calculated the anaerobic threshold (AT) using the V-slope method.The sEMG and ECG corresponding to the time of AT were categorized as "Easy," "Transition," and "Tired."

Signal preprocessing
The original sEMG and ECG signal contains noise interference, which needs to be preprocessed.Firstly, the original ECG and sEMG signals were filtered by 0$100 Hz and 0-500 Hz low-pass filters to remove high-frequency interference.Secondly, 49.5-50 Hz adaptive notch filters were used to filter the power frequency and harmonic interference in the signal.Finally, empirical mode decomposition(EMD) and discrete wavelet transform(DWT) domains were used to reduce the noise11.which reduce the noise from the initial IMFs instead of discarding them completely thus yielding a relatively cleaner ECG signal.

Noise reduction algorithm based on ISSA-VMD-SWT Variational mode decomposition (VMD)
VMD as a method with the characteristics of Wiener filtering and adaptive signal decomposition, effectively addresses the issues of mode mixing and high-frequency signal loss.It performs well in handling nonlinear signals. 18The specific decomposition process of the signal is as follows.
(1) The initialization of parameters in VMD is based on the optimization model and constraints of the variational problem.The optimization model and constraints for the variational problem are as follows: (2) By introducing an augmented Lagrange function, the constrained variational problem can be transformed into an unconstrained variational problem as follows: (3) The augmented Lagrange function in Equation 2 can be solved using the alternating direction method of multipliers (ADMM).This iterative algorithm aims to find the optimal solution of the model by updating the mode components and Lagrange multipliers alternately.The steps of the ADMM algorithm for solving Equation 2 can be outlined as follows: 3) 4) Where u n+1 k ðwÞ and w n+1 k are the wiener filter and frequency center of each component.It can be seen from the VMD decomposition process that the empirical selection of parameters k and a in practical applications may not achieve the best decomposition performance or even appear mode aliasing.Therefore, how to self-apply the selection of VMD parameters is the key to accurate decomposition.

Improved sparrow search algorithm (ISSA)
In order to achieve better decomposition and denoising effects, this study utilizes intelligent optimization algorithms to optimize the key parameters, k and a, of VMD within a reasonable range.The SSA is derived from the foraging and anti-predation processes of sparrows, possesses strong global optimization ability and local search capability. 19However, it suffers from issues such as uneven initial population distribution, slow convergence speed, and susceptibility to local optima in later stages. 20In this study, an improved Tent chaotic mapping is introduced to optimize the initial population, thereby improving the optimization speed.Simultaneously, the number of discoverers, joiners, and scouts is optimized to enhance the global search ability in the early stages and the search capability in the later stages, thus avoiding being trapped in local optima.The specific steps are as follows: (1) Based on the improved Tent chaotic mapping, a sequence y i is generated, and this sequence is mapped to different individual sparrows to improve the population distribution.

> > > <
> > > : 5) 6) Where y i represents the chaotic sequence; i = 1;2;:::;N; y i is a random number between 0 and 1. N represents the population size, T max is the maximum number of iterations; X i represents the initial population; b min and b max denote the lower and upper limits of the search space.
(2) Adaptive weights are introduced to update the finder position as: > > : 7) Where X n i is the position of the i sparrow of the N-th generation, X n best is the best position for individual sparrows, n indicates the number of iterations, g 1 and g 2 are random numbers of (0,1) and random numbers that follow a normal distribution, respectively, L is a row of multidimensional all-one matrix, R 2 is an alarm value [0,1].ST is a security threshold [0.5,1].When R 2 ＜ST, the discoverer continued to search for food in situ, guiding the population to get better fitness, When R 2 R ST, the finder should immediately move to a safe area to prevent predators.
(3) During the foraging process, if the discoverer finds a better food source, the joiner will immediately move toward it, and its position is updated as follows: (Equation 8) Where X n worst is the worst position for individual sparrows, When i＞N=2, the i-th entrant will leave the original feeding ground and move to other food sources.When i % N=2, the i-th entrant stays in place and forages.
(4) During the foraging process, when danger is detected, the warner will leave its current location and move toward safer food sources or approach other sparrows.The position of the warner is updated as follows: X n best + g 3 3 X n i À X n bot ; C i sC g X n i + g 4 3 X n i À X n best ðC i À C w Þ+ε ; C i = C g (Equation 9) Where g 3 and g 4 represent random numbers following a normal distribution and [-1,1] uniform distribution.ε is a very small value, typically on the order of 10 À 50 .C i is the fitness value of the current sparrow individual, and C w represents the fitness value of the global worst individual.
To improve the algorithm's initial global search capability and avoid falling into local optima later on, the optimization of the number of explorers and sentinels is implemented.The numbers of explorers and sentinels are dynamically adjusted during the algorithm's execution.fatigue recognition.Referring to our previous research, 20 we will use a feature fusion model based on Support Vector Machine (SVM) to perform fatigue recognition by combining preprocessed and denoised ECG and sEMG signals.By comparing the experimental results using different denoising methods, we can evaluate the impact of each method on the recognition rate.A higher recognition rate indicates a better improvement in the assessment of muscular weak electrical signal fatigue using the denoising method, enabling more accurate identification of fatigue states.At the same time, we also need to consider factors such as algorithm complexity, real-time performance, and applicability to different signal types to comprehensively evaluate the advantages and disadvantages of various methods and select the most suitable method for muscular weak electrical signal fatigue assessment tasks.The formula for calculating the recognition rate is as follows: Recognition rate = number of samples correctly identified total number of test samples 3 100% (Equation 21)

QUANTIFICATION AND STATISTICAL ANALYSIS
Statistical analyses were performed using custom python.Details of all statistical analyses can be found above in the relevant subsections of the method details section.Sample number (n) and statistical methods used to assess differences between groups are indicated in the relevant subsections of the result section.
Research on the recognition model of exercise fatigue based on the fusion of sEMG and ECG signals

Figure 4 .
Figure 4. Different algorithm for noise reduction of sEMG signals with 5dB and 10dB noise 5dB indicates the denoised 5dB, 10dB indicates the denoised 10dB, Soft indicates the soft threshold, and Hard indicates the hard threshold.

Figure 7 .
Figure 7.The recognition rate of sports fatigue of different models The values marked in the figure are the first quartile value, the average value and the third quartile value of the model respectively.

u k ðtÞ ! e À jw k t 2 X k u k ðtÞ = f ðtÞ (Equation 1 )
where u k ðtÞ represents the k-th mode component obtained from the decomposition of the original signal f ðtÞ, w k represents the central frequency of each component.

Table 1 .
The statistical analysis of HRV feature parameters under different levels of muscle fatigue

Table 4 .
Denoise indicates value of the denoised ECG signal

Table 3 .
The results of two-factor analysis of variance different denoising methods.Compared with the three algorithms, ISA-VMD-SWT has certain advantages in different parameters of the three states.

Table 5 .
Comparison of denoising effects of various algorithmsintegrating technology into wearable devices, and refining feature extraction and classification techniques.These efforts aim to solidify the model's utility in exercise physiology, rehabilitation, and healthcare monitoring, pushing forward the practical implementation of advanced fatigue recognition in wearable technology.