Extremely weak magnetic determination method and system for traditional Chinese medicine acupuncture and moxibustion

Abstract

The invention relates to the field of signal determination in diagnosis, in particular to an extremely weak magnetic determination method and system for traditional Chinese medicine acupuncture and moxibustion. The method comprises the following steps: during acupuncture treatment, acquiring different signals, amplifying an extremely weak magnetic signal, and displaying the amplified extremely weak magnetic signal in a display; performing low-scale analysis on the signal, and constructing a low-frequency electromagnetic response characteristic factor under a low scale; constructing a biological magnetic interference characteristic value according to the chaos characteristic, the similarity and the low-frequency electromagnetic response characteristic factor between the extremely weak magnetic signal and the other signals; adjusting the Kalman gain according to the biological magnetic interference characteristic value to complete signal denoising; whether acupuncture is effective or not is judged according to the change of the extremely weak magnetic signal after the acupuncture part in the display is denoised. According to the invention, the accuracy of extremely weak magnetic signal measurement in acupuncture and moxibustion composition is improved.

Description

Technical Field

[0001] The present application relates to the field of signal measurement in diagnosis, and in particular to an extremely weak magnetic measurement method and system for traditional Chinese medicine acupuncture. Background Art

[0002] Acupuncture, as an important part of traditional Chinese medicine, regulates the body's qi and blood and treats various diseases by stimulating specific acupuncture points. With the development of science and technology, especially the application of extremely weak magnetic field measurement technology, new tools have been provided for scientific research on the efficacy of acupuncture. This technology can capture the weak magnetic signals generated by the human body during acupuncture, thus providing an objective and quantifiable means for evaluating the efficacy of acupuncture. By measuring the changes in brain magnetic signals after acupuncture at different acupoints, it is possible to optimize treatment plans and improve treatment effects, and also promote a deeper understanding of the mechanism of acupuncture and promote the modernization of traditional Chinese medicine.

[0003] The prospect of extremely weak magnetic field measurement technology is broad. It can not only be used as a new treatment method to provide patients with highly sensitive medical services in a passive, non-destructive and wearable way, but also combine multidisciplinary knowledge such as medicine, physics, and engineering to promote the innovative development of acupuncture. In the extremely weak magnetic determination of acupuncture efficacy, bio-electromagnetic noise is an interference factor that cannot be ignored. The human body will generate weak electromagnetic fields in normal physiological activities, such as heartbeat, nerve impulse transmission, muscle contraction, etc. These electromagnetic fields generated by physiological activities are called bio-electromagnetic noise. They may overlap with the magnetic signals generated by acupuncture efficacy in time and space, thereby mixing with each other, increasing the complexity of magnetic signal processing and greatly affecting the accuracy of extremely weak magnetic signal determination. Summary of the invention

[0004] In order to solve the technical problem of the accuracy of extremely weak magnetic signal measurement, the present application provides an extremely weak magnetic measurement method and system for traditional Chinese medicine acupuncture. The technical solution adopted is as follows: In the first aspect, the present application proposes an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture, the method comprising the following steps: Collect different signals, amplify the extremely weak magnetic signals and display them on the monitor; De-noise the extremely weak magnetic signal based on the interference characteristics of different collected signals on the extremely weak magnetic signal; The effectiveness of acupuncture can be judged by the change of the extremely weak magnetic signal after noise removal at the acupuncture point on the display; The specific steps for denoising extremely weak magnetic signals are as follows: (1) Perform low-scale analysis on the signal. At low scale, construct low-frequency electromagnetic response characteristic factors based on the frequency characteristics of different signals and the signal energy corresponding to the signals; (2) Construct the biomagnetic interference characteristic value based on the difference in the chaotic characteristics of the extremely weak magnetic signal and the rest of the signals at a low scale, the similarity of the signals, the correlation between the rest of the signals, and the low-frequency electromagnetic response characteristic factors; (3) Obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and completing signal denoising based on this.

[0005] In the above scheme, the present application proposes a method and system for extremely weak magnetic measurement for traditional Chinese medicine acupuncture. Aiming at the interference problem of bio-electromagnetic noise on extremely weak magnetic measurement of acupuncture efficacy, based on the method of secondary discrete wavelet transform, information entropy and correlation analysis, the interaction correlation between the magnetic signal generated by acupuncture efficacy and the bio-magnetic signal is analyzed, and a low-frequency electromagnetic response characteristic factor is constructed to characterize the electromagnetic response characteristics and synchronization between different signals; Aiming at the interference characteristic problem of bio-magnetic signals on extremely weak magnetic signals of acupuncture, based on the method of correlation coefficient and least squares linear regression, the correlation strength and trend similarity between bio-magnetic signals are analyzed, and a bio-magnetic interference characteristic value is constructed to characterize the degree of interference of bio-magnetic signals on acupuncture signals; the Kalman filter algorithm is further improved, and weights are added to the calculation of the Kalman gain to achieve the beneficial effect of balancing noise reduction and detail retention in different individuals, thereby solving the challenge of effectively separating signals from bio-electromagnetic noise and ensuring the accuracy of extremely weak magnetic signal measurement in acupuncture efficacy.

[0006] In one embodiment, the collected signals include an extremely weak magnetic signal, a heart rate signal and a blood pressure signal, and all data collected from each signal constitute a sequence, which are respectively recorded as an extremely weak magnetic sequence, a heart rate sequence and a blood pressure sequence.

[0007] In one embodiment, the method of performing low-scale analysis on the signal and constructing a low-frequency electromagnetic response characteristic factor according to the frequency characteristics of different signals and the signal energy corresponding to the signal at the low scale is: The extremely weak magnetic sequence, heart rate sequence and blood pressure sequence are converted into frequency domain signals through wavelet transform, and low-frequency signals are selected at four scales to obtain extremely weak magnetic low-scale signals, heart rate low-scale signals and blood pressure low-scale signals; Perform coherence analysis between low-scale signals to obtain low-frequency coherence eigenvalues ​​of two low-scale signals; The signal energy of each low-scale signal is calculated, and the low-frequency electromagnetic response characteristic factor is obtained according to the low-frequency coherence eigenvalues ​​between all low-scale signals and the signal energy of the low-scale signal.

[0008] In one embodiment, the method of performing coherence analysis between low-scale signals to obtain low-frequency coherence eigenvalues ​​of two low-scale signals is: For any two low-scale signals, obtain the wavelet coherence power spectrum between the two low-scale signals, where each frequency corresponds to a coherence value, set the significance level value, and assign the coherence value of the frequency to 0 when the coherence value is greater than the significance level value among all frequencies, and keep the coherence value unchanged when the coherence value is less than the significance level value; calculate the mean of the re-assigned coherence values ​​corresponding to all frequencies and record it as the low-frequency coherence eigenvalue of the two low-scale signals.

[0009] In one embodiment, the method for obtaining the low-frequency electromagnetic response characteristic factor according to the low-frequency coherence eigenvalues ​​between all low-scale signals and the signal energy of the low-scale signals is: Any low-scale signal is recorded as the target low-scale signal, the ratio of the signal energy of the target low-scale signal to the signal energy of another low-scale signal is taken as the signal energy ratio, the difference between the signal energy ratio and the number 1 is normalized by an inverse proportional function and the product of the corresponding low-frequency coherent eigenvalue is taken as the electromagnetic response feature of the two low-scale signals, and the average of the electromagnetic response features obtained pairwise from all low-scale signals is taken as the low-frequency electromagnetic response characteristic factor.

[0010] In one embodiment, the method for constructing the biomagnetic interference characteristic value according to the difference in chaotic characteristics between the extremely weak magnetic signal and the other signals at a low scale, the similarity of the signals, the correlation between the other signals, and the low-frequency electromagnetic response characteristic factor is: For each low-scale signal, the information entropy of the low-scale signal is calculated, and the low-scale signal of heart rate and the low-scale signal of blood pressure are recorded as the influence signal. For the influence signal and the very weak magnetic signal, a window is constructed based on each data point to obtain the correlation coefficient between the signal slope of the window and the influence signal; The biomagnetic interference characteristic value is constructed based on the difference in signal slope corresponding to the extremely weak magnetic low-scale signal and the signal slope of the influencing signal, the difference in information entropy, the correlation coefficient between the influencing signals corresponding to each window, and the low-frequency electromagnetic response characteristic factor.

[0011] In one embodiment, the method of constructing a window based on each data point and obtaining the signal slope of the window and the correlation coefficient between the influencing signals is: For each impact signal, a 1*N window is constructed with each data point as the center, and the correlation coefficient of the windows corresponding to the data points at the same position between the impact signals is calculated by the Pearson correlation coefficient; For each influencing signal, the data in each window is taken as input, and the least squares linear regression is used to fit it, so as to obtain the straight line fitted by the data in the influencing signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window; For extremely weak magnetic low-scale signals, the data in each window is taken as input and fitted using least squares linear regression to obtain the straight line fitted to the data in the extremely weak magnetic low-scale signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window.

[0012] In one embodiment, the method for constructing the biomagnetic interference characteristic value according to the difference between the signal slope corresponding to the extremely weak magnetic low-scale signal and the signal slope of the influencing signal, the information entropy difference, the correlation coefficient between the influencing signals, and the low-frequency electromagnetic response characteristic factor is: The absolute value of the difference in information entropy between the influencing signal and the extremely weak magnetic low-scale signal is normalized by an inverse proportional function to obtain the information difference, and the ratio of the signal slopes of the windows corresponding to the data points at the same position of the influencing signal and the extremely weak magnetic low-scale signal is recorded as the slope difference; the information difference between the two influencing signals and the extremely weak magnetic low-scale signal and the ratio of each slope difference are recorded as information difference values, and the product of the average of all information difference values, the average of the correlation coefficients between the influencing signals in all windows, and the low-frequency electromagnetic response characteristic factor is taken as the biomagnetic interference characteristic value.

[0013] In one embodiment, the method for obtaining the real-time Kalman gain value by using the initial Kalman gain and the biomagnetic interference characteristic value is: The product of the initial Kalman gain and the biomagnetic interference characteristic value is used as the real-time Kalman gain.

[0014] In a second aspect, the present application also provides an extremely weak magnetic measurement system for traditional Chinese medicine acupuncture, the system comprising: Acupuncture treatment signal filtering module: used to filter the interference characteristics of extremely weak magnetic signals according to the collected different signals, perform low-scale analysis on the signals, and construct low-frequency electromagnetic response characteristic factors according to the frequency characteristics of different signals and the signal energy corresponding to the signals at low scale; The biomagnetic interference characteristic value is constructed according to the difference in the chaotic characteristics between the extremely weak magnetic signal and the other signals at a low scale, the similarity of the signals, the correlation between the other signals and the low-frequency electromagnetic response characteristic factor; Obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and filtering the signal based on the real-time Kalman gain value; Acupuncture treatment effect judgment module: used to judge whether acupuncture is effective through the changes in the extremely weak magnetic signal after filtering at the acupuncture site in the display. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the technical solutions and advantages in the embodiments of the present application or the prior art, the drawings required for use in the embodiments or the prior art descriptions are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0016] Figure 1 A flow chart of an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture provided in one embodiment of the present application; Figure 2 This is the change of the extremely weak magnetic signal for this application. DETAILED DESCRIPTION

[0017] In order to further explain the technical means and effects adopted by the present application to achieve the predetermined invention purpose, the following is a detailed description of the extremely weak magnetic measurement method and system for traditional Chinese medicine acupuncture proposed in accordance with the present application, its specific implementation method, structure, characteristics and effects, in combination with the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" does not necessarily refer to the same embodiment. In addition, specific features, structures or characteristics in one or more embodiments may be combined in any suitable form.

[0018] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs.

[0019] An extremely weak magnetic measurement method and system embodiment for traditional Chinese medicine acupuncture: The specific scheme of the extremely weak magnetic measurement method and system for traditional Chinese medicine acupuncture provided by the present application is described in detail below with reference to the accompanying drawings.

[0020] See also Figure 1 , which shows a flow chart of an extremely weak magnetic measurement method for traditional Chinese medicine acupuncture provided by an embodiment of the present application, the method comprising the following steps: Step S001, during acupuncture treatment, different types of signals at acupuncture positions are obtained, and extremely weak magnetic signals are amplified and displayed on a display.

[0021] Place the ultra-high-sensitivity atomic magnetometer in a stable measurement environment (such as a magnetic shielding room that shields external magnetic field interference). Before starting the measurement, the subjects need to accurately find and mark the acupuncture points. Place the atomic magnetometer near all acupuncture points to capture the extremely weak magnetic signals generated by acupuncture stimulation, realize continuous real-time acquisition of magnetic signal data, and use photoplethysmograph (PPG) to obtain the patient's real-time heart rate and blood pressure data. The collected signal is amplified by a signal amplifier, and the amplified signal is displayed on the display through a signal display to display the extremely weak magnetic signal.

[0022] During the acupuncture process, the signal data will change. The efficacy of acupuncture can be determined by detecting the changes in the extremely weak magnetic signals at the index site. The changes in the signals on the display can provide the physician with necessary feedback information for adjusting the time and depth of acupuncture.

[0023] The extremely weak magnetic signal, heart rate and blood pressure data are uniformly collected at intervals of time t. In this embodiment, the empirical value of t is 0.05s, and the three types of data collected are used as input and are subjected to real-time standardization processing using the Z-Score algorithm to form an extremely weak magnetic sequence, a heart rate sequence and a blood pressure sequence.

[0024] At this point, the signal acquisition is completed.

[0025] Step S002, performing low-scale analysis on the signal, and constructing a low-frequency electromagnetic response characteristic factor according to the frequency characteristics of different signals and the signal energy corresponding to the signal at the low scale.

[0026] In the extremely weak magnetic measurement of acupuncture efficacy, the magnetic field at the acupoints is usually related to acupuncture stimulation, with a weak intensity and limited to a specific area. It may have specific frequency characteristics and obvious time dependence, while the magnetic field generated by the heart and brain is of greater intensity, with a wider frequency range and is relatively stable. When conducting extremely weak magnetic measurement of acupuncture efficacy, the subject's physiological state, including heart rate, blood pressure, etc., may have a significant impact on the measurement results. These physiological parameters are closely related to the electrophysiological activities in the body, and their changes will be reflected in the measurement of bioelectromagnetic fields. The beating of the heart and blood circulation will generate an electromagnetic field. The current generated by each beat propagates in the body to form an electromagnetic field that can be detected by external sensors. Changes in blood pressure may also change the electromagnetic field by affecting blood circulation, thus affecting the display of extremely weak magnetic signals. After being displayed on the monitor, the doctor will make adjustments based on the content of the monitor, which will result in errors. Therefore, it is necessary to suppress the heart rate signal and blood pressure signal.

[0027] Taking extremely weak magnetic sequence, heart rate sequence and blood pressure sequence as input, the three sequences are converted from time domain to frequency domain signals using secondary discrete wavelet transform. After discrete wavelet transform, frequency domain signals at four scales of LL, HL, LH and HH are obtained. Since the electromagnetic field signals generated by the heart and brain are mainly located in the low-frequency region, and these signals may have similar frequency characteristics to the weak magnetic field signals generated by acupuncture stimulation, the LL low-frequency signal after secondary discrete wavelet transform is selected for analysis in the bioelectromagnetic field measurement of acupuncture efficacy.

[0028] The signals of the extremely weak magnetic sequence, the heart rate sequence and the blood pressure sequence at the LL scale after wavelet transformation are extracted and recorded as extremely weak magnetic low-scale signals, heart rate low-scale signals and blood pressure low-scale signals. The signal energy of all low-scale signals is calculated. The calculation method of the signal energy is a well-known technology and will not be described in detail in this embodiment.

[0029] Assume that the maximum frequency range of the extremely weak magnetic low-scale signal, heart rate low-scale signal, and blood pressure low-scale signal is , perform wavelet coherence analysis on any two of the three signals, that is, for any two low-scale signals, obtain the wavelet coherence power spectrum between the two low-scale signals, where each frequency corresponds to a coherence value, set the significance level value, and assign the coherence value of the frequency to 0 when the coherence value is greater than the significance level value among all frequencies, and when the coherence value is less than the significance level value, the coherence value of the frequency remains unchanged. Calculate the mean of the re-assigned coherence values ​​corresponding to all frequencies and record it as the low-frequency coherence eigenvalue of the two low-scale signals.

[0030] Preferably, in one embodiment of the present application, the maximum frequency range is, and the significance level value is 0.5.

[0031] The low-frequency coherence eigenvalue quantifies the interaction between acupuncture stimulation and the subject's physiological state by calculating the average value of the coherence of two low-scale signals within a specific frequency range. The larger the value, the greater the coherence between the two signals, that is, the stronger the synchronization and correlation between the two signals at this frequency.

[0032] The low-frequency electromagnetic response characteristic factor is obtained according to the low-frequency coherence eigenvalues ​​between all low-scale signals and the signal energy of the low-scale signals.

[0033] Any low-scale signal is recorded as the target low-scale signal, the ratio of the signal energy of the target low-scale signal to the signal energy of another low-scale signal is taken as the signal energy ratio, the difference between the signal energy ratio and the number 1 is normalized by an inverse proportional function and the product of the corresponding low-frequency coherent eigenvalue is taken as the electromagnetic response feature of the two low-scale signals, and the average of the electromagnetic response features obtained pairwise from all low-scale signals is taken as the low-frequency electromagnetic response characteristic factor.

[0034] The electromagnetic response characteristic factor is used to quantify the electromagnetic response characteristics between three different low-scale signals, reflecting the interaction correlation between extremely weak magnetic signals and biomagnetic signals. If the signal energy of the two low-scale signals is close, that is, they have similar intensities in low-frequency electromagnetic response, the lower the low-frequency electromagnetic response characteristic factor, the greater the low-frequency coherence characteristic value, which can reflect the synchronization between the electromagnetic field signal generated by acupuncture stimulation and the electromagnetic field signal generated by the heart or brain. The greater the value, the greater the low-frequency electromagnetic response characteristic factor.

[0035] At this point, the low-frequency electromagnetic response characteristic factors are obtained.

[0036] Step S003, constructing a biomagnetic interference characteristic value according to the difference in chaos characteristics between the extremely weak magnetic signal and the remaining signals at a low scale, the similarity of the signals, the correlation between the remaining signals, and the low-frequency electromagnetic response characteristic factor.

[0037] For all signals, all data collected from the beginning of collection to the current moment is the signal corresponding to the current moment.

[0038] For each low-scale signal, the information entropy of the low-scale signal is calculated, and the low-scale signal of heart rate and the low-scale signal of blood pressure are recorded as influencing signals. For each influencing signal, a 1*N window is constructed with each data point as the center. In this embodiment, N is 11. The correlation coefficient of the corresponding windows of the data points at the same position between the influencing signals is calculated using the Pearson correlation coefficient.

[0039] It should be noted that the present implementation uses the Pearson correlation coefficient to calculate the correlation coefficient corresponding to two windows, and the implementer may use other methods of calculating the correlation coefficient, such as the Spearman rank correlation coefficient.

[0040] For each influencing signal, the data in each window is taken as input and fitted using least squares linear regression to obtain the straight line fitted to the data in the influencing signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window.

[0041] For extremely weak magnetic low-scale signals, the data in each window is also used as input and fitted using least squares linear regression to obtain the straight line fitted to the data in the extremely weak magnetic low-scale signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window.

[0042] In this way, the signal slope of the window corresponding to each data point of the extremely weak magnetic low-scale signal is obtained.

[0043] The biomagnetic interference characteristic value is constructed based on the difference in signal slope corresponding to the extremely weak magnetic low-scale signal and the signal slope of the influencing signal, the difference in information entropy, the correlation coefficient between the influencing signals corresponding to each window, and the low-frequency electromagnetic response characteristic factor.

[0044] The absolute value of the difference in information entropy between the influencing signal and the extremely weak magnetic low-scale signal is normalized by an inverse proportional function to obtain the information difference, and the ratio of the signal slopes of the windows corresponding to the data points at the same position of the influencing signal and the extremely weak magnetic low-scale signal is recorded as the slope difference; the information difference between the two influencing signals and the extremely weak magnetic low-scale signal and the ratio of each slope difference are recorded as information difference values, and the product of the average of all information difference values, the average of the correlation coefficients between the influencing signals in all windows, and the low-frequency electromagnetic response characteristic factor is taken as the biomagnetic interference characteristic value.

[0045] The biomagnetic interference characteristic value reflects the interference characteristics of the biomagnetic signal on the extremely weak magnetic signal of acupuncture. The larger the value, the more significant the interference of the biomagnetic signal on the acupuncture signal, which indicates that the acupuncture effect is more affected by the electromagnetic fields of the heart and brain. The correlation coefficient between the influencing signals is used as a weight, which characterizes the correlation strength between the influencing signals when the electromagnetic response characteristics are considered. The low-frequency electromagnetic response characteristic factor quantifies the electromagnetic response characteristics between the three different signals, reflecting the interaction correlation between the extremely weak magnetic signal and the biomagnetic signal. The larger the value, the stronger the correlation between the signals and the larger the biomagnetic interference characteristic value. The information difference gives a higher weight to the signal combination with a smaller information entropy difference. If the information entropy of the influencing signal is close to the information entropy of the extremely weak magnetic signal, it means that the uncertainty of the influencing signal and the extremely weak magnetic signal are similar. Its value will be close to 1, and the contribution to the biomagnetic interference characteristic value will be greater. On the contrary, if the information entropy difference is large, this value will be close to 0, reducing the contribution to the biomagnetic interference characteristic value. The signal slope measures the similarity of the trends of the two signals. When its value is larger, the trends of the two signals are not similar, and the biomagnetic interference characteristic value is smaller.

[0046] At this point, the biomagnetic interference characteristic value is obtained.

[0047] Step S004, obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and filtering the signal based on the real-time Kalman gain value.

[0048] In order to eliminate the influence of the electromagnetic fields of the heart and brain on the extremely weak magnetic signals of acupuncture, the Kalman filter algorithm is used to denoise the acquired extremely weak magnetic sequences. The Kalman gain K in the Kalman filter algorithm controls the balance between noise reduction and detail retention. A larger K value can remove noise more effectively, but may lead to the loss of useful signals; a smaller K value helps to retain signal details, but the denoising effect may be weakened. Since extremely weak magnetic signals are often interfered by the strong electromagnetic fields generated by organs such as the heart and brain, in order to improve the signal-to-noise ratio of the signal and accurately extract the acupuncture signal, the method of adjusting the K value is adopted to adapt to the noise level and signal characteristics of different individuals, ensuring that the weak magnetic signals generated by acupuncture are retained to the maximum extent while removing the electromagnetic interference of the heart and brain.

[0049] The initial empirical value of the Kalman gain of the algorithm is calculated as a fixed formula, which is obtained by the ratio of the prediction error covariance and the observation noise covariance. At this time, the extremely weak magnetic sequence is used as the state input of the Kalman filter and the object of Kalman filter denoising. The extremely weak magnetic sequence, heart rate sequence and blood pressure sequence are used as the reference input of the Kalman filter to improve the calculation of the Kalman gain, calibrate and adjust the performance of the filter, calculate the real-time biomagnetic interference characteristic value through the reference input, and make real-time adjustments to the calculation of the K value.

[0050] The product of the normalized biomagnetic interference characteristic value and the initial Kalman gain is used as the real-time Kalman gain.

[0051] The extremely weak magnetic signals are denoised using the Kalman filter algorithm.

[0052] Through real-time improvement of the Kalman gain, the algorithm can adaptively adjust the Kalman gain according to the strength of interference of the biomagnetic signal on the extremely weak magnetic signal of acupuncture. The stronger the interference of the biomagnetic signal on the extremely weak magnetic signal of acupuncture, the larger the real-time Kalman gain of the algorithm, which makes the algorithm more capable of denoising the state input to ensure the elimination of the noise influence caused by the biomagnetic signal. The weaker the interference of the biomagnetic signal on the extremely weak magnetic signal of acupuncture, the smaller the real-time Kalman gain of the algorithm, which allows the algorithm to retain more details when denoising the state input to ensure the accuracy of the extremely weak magnetic signal.

[0053] At this point, the real-time Kalman gain is obtained and the denoising of the extremely weak magnetic signal is completed.

[0054] Step S005, judging whether acupuncture is effective by the change of the extremely weak magnetic signal after denoising at the acupuncture site in the display.

[0055] Through the above steps, the extremely weak magnetic signal after denoising is displayed on the monitor. When the doctor is performing acupuncture, he adjusts the depth of acupuncture and judges the efficacy of acupuncture by observing the peak change of the extremely weak magnetic signal. As acupuncture proceeds, the extremely weak magnetic signal changes. The higher the peak value, the stronger the magnetic field and the better the efficacy. The extremely weak magnetic signal changes as shown in the figure. Figure 2 The depth with the highest peak value is stored in the memory for reference in the next acupuncture.

[0056] Based on the same inventive concept as the above method, the embodiment of the present invention also provides an extremely weak magnetic field measurement system for traditional Chinese medicine acupuncture, including the following modules: Acupuncture treatment signal sensing module: used to obtain different types of signals at acupuncture positions during acupuncture treatment, amplify extremely weak magnetic signals and display them on the display; Acupuncture treatment signal filtering module: used to filter the interference characteristics of extremely weak magnetic signals according to the different signals collected, perform low-scale analysis on the signals, and construct low-frequency electromagnetic response characteristic factors according to the frequency characteristics of different signals and the signal energy corresponding to the signals at low scale; The biomagnetic interference characteristic value is constructed according to the difference in chaotic characteristics between the extremely weak magnetic signal and the rest of the signals at low scale, the similarity of the signals, the correlation between the rest of the signals and the low-frequency electromagnetic response characteristic factor; Obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and filtering the signal based on the real-time Kalman gain value; Acupuncture treatment effect judgment module: used to judge whether acupuncture is effective through the changes in the extremely weak magnetic signal after filtering at the acupuncture site in the display.

[0057] It should be noted that the above-described embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the scope of the technical solutions of the embodiments of the present application, and should all be included in the protection scope of the present application.

[0058] The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referenced to each other, and each embodiment focuses on the differences from other embodiments.

Claims

1. A very weak magnetic measurement method for traditional Chinese medicine acupuncture, characterized in that: The method comprises the following steps: During acupuncture treatment, different types of signals are obtained at the acupuncture positions, and the extremely weak magnetic signals are amplified and displayed on the monitor; Signal filtering is performed on the interference characteristics of extremely weak magnetic signals according to the different signals collected; The effectiveness of acupuncture can be judged by the changes in the extremely weak magnetic signals filtered at the acupuncture points on the display; The signal filtering process is as follows: (1) Perform low-scale analysis on the signal. At low scale, construct low-frequency electromagnetic response characteristic factors based on the frequency characteristics of different signals and the signal energy corresponding to the signals; (2) Construct the biomagnetic interference characteristic value based on the difference in the chaotic characteristics of the extremely weak magnetic signal and the rest of the signals at a low scale, the similarity of the signals, the correlation between the rest of the signals, and the low-frequency electromagnetic response characteristic factors; (3) Obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and filtering the signal based on the real-time Kalman gain value.

2. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 1, characterized in that: The collected signals include extremely weak magnetic signals, heart rate signals and blood pressure signals. All data collected from each signal constitute a sequence, which are recorded as extremely weak magnetic sequence, heart rate sequence and blood pressure sequence respectively.

3. A very weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 2, characterized in that: The method of performing low-scale analysis on the signal and constructing a low-frequency electromagnetic response characteristic factor according to the frequency characteristics of different signals and the signal energy corresponding to the signal at a low scale is: The extremely weak magnetic sequence, heart rate sequence and blood pressure sequence are converted into frequency domain signals through wavelet transform, and low-frequency signals are selected at four scales to obtain extremely weak magnetic low-scale signals, heart rate low-scale signals and blood pressure low-scale signals; Perform coherence analysis between low-scale signals to obtain low-frequency coherence eigenvalues ​​of two low-scale signals; The signal energy of each low-scale signal is calculated, and the low-frequency electromagnetic response characteristic factor is obtained according to the low-frequency coherence eigenvalues ​​between all low-scale signals and the signal energy of the low-scale signal.

4. A very weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 3, characterized in that: The method for performing coherence analysis between low-scale signals to obtain low-frequency coherence eigenvalues ​​of two low-scale signals is: For any two low-scale signals, obtain the wavelet coherence power spectrum between the two low-scale signals, where each frequency corresponds to a coherence value, set the significance level value, and among all frequencies, when the coherence value is greater than the significance level value, assign the coherence value of the frequency to 0, and when the coherence value is less than the significance level value, the coherence value of the frequency remains unchanged; The average of the re-assigned coherence values ​​corresponding to all frequencies is calculated and recorded as the low-frequency coherence eigenvalue of the two low-scale signals.

5. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 3, characterized in that: The method for obtaining the low-frequency electromagnetic response characteristic factor according to the low-frequency coherent eigenvalues ​​between all low-scale signals and the signal energy of the low-scale signals is: Any low-scale signal is recorded as the target low-scale signal, the ratio of the signal energy of the target low-scale signal to the signal energy of another low-scale signal is taken as the signal energy ratio, the difference between the signal energy ratio and the number 1 is normalized by an inverse proportional function and the product of the corresponding low-frequency coherent eigenvalue is taken as the electromagnetic response feature of the two low-scale signals, and the average of the electromagnetic response features obtained pairwise from all low-scale signals is taken as the low-frequency electromagnetic response characteristic factor.

6. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 3, characterized in that: The method for constructing the biomagnetic interference characteristic value according to the difference in chaotic characteristics between the extremely weak magnetic signal and the other signals at a low scale, the similarity of the signals, the correlation between the other signals and the low-frequency electromagnetic response characteristic factor is: For each low-scale signal, the information entropy of the low-scale signal is calculated, and the low-scale signal of heart rate and the low-scale signal of blood pressure are recorded as the influence signal. For the influence signal and the very weak magnetic signal, a window is constructed based on each data point to obtain the correlation coefficient between the signal slope of the window and the influence signal; The biomagnetic interference characteristic value is constructed based on the difference in signal slope corresponding to the extremely weak magnetic low-scale signal and the signal slope of the influencing signal, the difference in information entropy, the correlation coefficient between the influencing signals corresponding to each window, and the low-frequency electromagnetic response characteristic factor.

7. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 6, characterized in that: The method of constructing a window based on each data point and obtaining the signal slope of the window and the correlation coefficient between the influencing signals is: For each impact signal, a 1*N window is constructed with each data point as the center, and the correlation coefficient of the windows corresponding to the data points at the same position between the impact signals is calculated by the Pearson correlation coefficient; For each influencing signal, the data in each window is taken as input, and the least squares linear regression is used to fit it, so as to obtain the straight line fitted by the data in the influencing signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window; For extremely weak magnetic low-scale signals, the data in each window is taken as input and fitted using least squares linear regression to obtain the straight line fitted to the data in the extremely weak magnetic low-scale signal window, and the slope of the straight line is obtained through the fitted straight line, thereby obtaining the signal slope corresponding to each window.

8. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 6, characterized in that: The method for constructing the biomagnetic interference characteristic value according to the difference between the signal slope corresponding to the extremely weak magnetic low-scale signal and the signal slope of the influencing signal, the information entropy difference, the correlation coefficient between the influencing signals, and the low-frequency electromagnetic response characteristic factor is: The absolute value of the difference in information entropy between the influencing signal and the extremely weak magnetic low-scale signal is normalized by an inverse proportional function to obtain the information difference, and the ratio of the signal slopes of the windows corresponding to the data points at the same position of the influencing signal and the extremely weak magnetic low-scale signal is recorded as the slope difference; the information difference between the two influencing signals and the extremely weak magnetic low-scale signal and the ratio of each slope difference are recorded as information difference values, and the product of the average of all information difference values, the average of the correlation coefficients between the influencing signals in all windows, and the low-frequency electromagnetic response characteristic factor is taken as the biomagnetic interference characteristic value.

9. The extremely weak magnetic measurement method for traditional Chinese medicine acupuncture as claimed in claim 1, characterized in that: The method for obtaining the real-time Kalman gain value by using the initial Kalman gain and the biomagnetic interference characteristic value is: The product of the initial Kalman gain and the biomagnetic interference characteristic value is used as the real-time Kalman gain.

10. An extremely weak magnetic measurement system for traditional Chinese medicine acupuncture, characterized in that: The system comprises: Acupuncture treatment signal sensing module: used to obtain different types of signals at acupuncture positions during acupuncture treatment, amplify extremely weak magnetic signals and display them on the display; Acupuncture treatment signal filtering module: used to filter the interference characteristics of extremely weak magnetic signals according to the different signals collected, perform low-scale analysis on the signals, and construct low-frequency electromagnetic response characteristic factors according to the frequency characteristics of different signals and the signal energy corresponding to the signals at low scale; The biomagnetic interference characteristic value is constructed according to the difference in chaotic characteristics between the extremely weak magnetic signal and the rest of the signals at low scale, the similarity of the signals, the correlation between the rest of the signals and the low-frequency electromagnetic response characteristic factor; Obtaining an initial Kalman gain value, using the initial Kalman gain and the biomagnetic interference characteristic value to obtain a real-time Kalman gain value, and filtering the signal based on the real-time Kalman gain value; Acupuncture treatment effect judgment module: used to judge whether acupuncture is effective through the changes in the extremely weak magnetic signal after filtering at the acupuncture site in the display.

Images