METHOD AND DEVICE OF ADAPTIVE EMC-EMI RADIO FREQUENCY SIGNAL DATA PROCESSING

Abstract

A method and a device of adaptive EMC-EMI radio frequency signal data processing are provided. The method includes: performing segmentation and preprocessing in response to a radio frequency signal; performing Hilbert-Huang transform on signals after segmentation; calculating EMC power, EMI radio frequency energy and mode, and a radio frequency signal-to-noise mode, and comparing the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode with corresponding thresholds; and adaptively adjusting energy parameters of a radio frequency ablation device, or prompting a user to adjust the energy parameters of the radio frequency ablation device. Instantaneous and dynamic radio frequency plasma is qualitatively and quantitatively detected and classified according to a signal-to-noise mode and pattern recognition of radio frequency emission, and an actual state of the knife head is perceived, to performs adaptive control or prompt the user to perform an adjustment operation.

Claims

1. A method of an adaptive EMC-EMI radio frequency signal data processing, comprising the following steps: S1: responding to, by a processor, a radio frequency signal of a radio frequency ablation device acquired by an EMC-EMI sensor; S2: segmenting, by the processor, the radio frequency signal, wherein a window stage after a segmentation comprises at least one stage of an idling, an energy generation, an ablation and a cutting, and a coagulation and a heat control, and preprocessing signals of each window stage after the segmentation respectively; S3: performing a Hilbert-Huang transform on signals of two window stages of the ablation and the cutting, and the coagulation and the heat control after the segmentation, to obtain time-frequency characteristics of electromagnetic signals of corresponding window stages; S4: calculating an EMC power, an EMI radio frequency energy and mode, and a radio frequency signal-to-noise mode according to the time-frequency characteristics of the electromagnetic signals; and comparing the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode with corresponding thresholds respectively, to obtain knife head states in the two window stages of the ablation and the cutting, and the coagulation and the heat control, wherein the knife head states comprise a touching metal short circuit, an overheating, and an overload; and S5: adjusting, by the processor, energy parameters of the radio frequency ablation device according to the knife head states obtained in step S4, or prompting a user to adjust the energy parameters of the radio frequency ablation device.

2. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 1, wherein the preprocessing in step S2 comprises: step 2.1: removing a high-frequency electrical noise by a low-pass filtering; step 2.2: eliminating a signal drift by a high-pass filtering; and step 2.3: removing a non-physiological pseudo-signal by a band-pass filtering.

3. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 1, wherein the Hilbert-Huang transform in step S3 comprises the following steps: S3.1: performing an empirical mode decomposition (EMD) on the signals of the two window stages of the ablation and the cutting, and the coagulation and the heat control respectively, to obtain limited intrinsic mode function (IMF) component expressions of a corresponding signal s(t);
s(t)=Σ.sub.i=1.sup.kc.sub.i(t)+r.sub.k(t), wherein i represents a number of an IMF component order, k represents the IMF component order, c.sub.i(t) is an i.sup.th order of an IMF component, and r.sub.k(t) is an i.sup.th order of an IMF residual component; S3.2: performing a Hilbert transform on each order of the IMF component c.sub.i(t), to obtain a Hilbert transform expression y.sub.i(t) of the IMF component; $y_i\left(t\right)=\frac{p}{\pi }{\int }_{-\infty }^{+\infty }\frac{c_i\left(t\right)}{t-\tau }d\tau$ wherein p is a Cauchy principal value, and τ represents an amount of change per unit time; S3.3: obtaining an instantaneous phase φ.sub.i(t) and an instantaneous frequency ω.sub.i(t) of each order of the IMF component by using the following formula; ${\phi }_i\left(t\right)=\arctan \left(y_i\left(t\right)/c_i\left(t\right)\right)$ ${\omega }_i\left(t\right)=\frac{d{\phi }_i\left(t\right)}{\mathrm{dt}};$ S3.4: comparing and screening a value of each frequency point in the instantaneous frequency ω.sub.i(t) with an instantaneous frequency threshold λ, if the instantaneous frequency |ω.sub.i|>λ, keeping the IMFi component unchanged, otherwise making the IMFi component be 0; S3.5: distinguishing between a square wave and a harmonic wave according to the value of each frequency point, and recording a spectrum and spectrum characteristics of the IMFi component; and S3.6: reconstructing a signal s′ (t) by using a screened IMF component in step S3.4, to obtain a time-frequency characteristic of an electromagnetic signal of a reconstructed signal.

4. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 3, wherein the obtaining limited intrinsic mode function (IMF) components of the signal s(t) in step S3.1 comprises the following steps: S3.1.1: segmenting an original signal s(t) according to a preset processing window, and taking a maximum value and a minimum value of the original signal s(t) for any processing window; S3.1.2: interpolating the maximum value and the minimum value of the original signal s(t) in the processing window by using a spline interpolation, and fitting upper and lower envelopes; S3.1.3: calculating a mean m.sub.1(t) of the upper and lower envelopes, and obtaining a difference h.sub.1(t) between the original signal s(t) and the mean m.sub.1(t), wherein h.sub.1(t)=s(t)−m.sub.1(t); S3.1.4: calculating a first IMF component c.sub.1(t) and a residual signal r.sub.1(t);
c.sub.1(t)=h.sub.1(t), r.sub.1(t)=s(t)−c.sub.1(t); S3.1.5: determining whether the residual signal r.sub.1(t) meets an empirical mode decomposition (EMD) condition, that is, whether r.sub.1(t) is a monotone function, if the residual signal r.sub.1(t) meets the EMD condition, stopping processing, and marking k=1; otherwise, repeating steps S3.1.1 to S3.1.4 fork times by using h.sub.1(t) as the original signal;
h.sub.k(t)=h.sub.(k−1)(t)−m.sub.k(t), c.sub.k(t)=h.sub.k(t), r.sub.k(t)=h.sub.(k−1)(t)−c.sub.k(t); completing a smooth decomposition processing, and recording k until the residual signal r.sub.k(t) is the monotone function; and S3.1.6: converting s(t) into k intrinsic mode function component expressions: s(t)=Σ.sub.i=1.sup.kc.sub.i(t)+r.sub.k(t), wherein c.sub.i(t) is an i.sup.th IMF component of the signal.

5. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 3, wherein the instantaneous frequency threshold λ in step S3.4 ranges from 0.0094 to 0.01.

6. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 3, wherein the spectrum characteristics in step S3.5 comprise: a frequency characteristic peak value, a peak-to-peak frequency-to-amplitude ratio, and a main frequency component.

7. The method of the adaptive EMC-EMI radio frequency signal data processing according to claim 1, wherein the energy parameters of the radio frequency ablation device in step S5 comprise: a voltage, a current, an impedance, a temperature, a humidity, and a contact force.

8. A device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 1, comprising: a radio frequency signal segmentation unit, configured to respond to the radio frequency signal of the radio frequency ablation device acquired by the EMC-EMI sensor, segment the radio frequency signal, wherein the window stage after the segmentation comprises the at least one stage of the idling, the energy generation, the ablation and the cutting, and the coagulation and the heat control, and preprocess the signals of each window stage after the segmentation respectively; a Hilbert-Huang transformation unit, configured to perform the Hilbert-Huang transform on the signals of the two window stages of the ablation and the cutting, and the coagulation and the heat control after the segmentation, to obtain the time-frequency characteristics of the electromagnetic signals of the corresponding window stages; a threshold comparison unit, configured to calculate the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode according to the time-frequency characteristics of the electromagnetic signals; and compare the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode with the corresponding thresholds respectively, to obtain the knife head states in the two window stages of the ablation and the cutting, and the coagulation and the heat control, wherein the knife head states comprise the touching metal short circuit, the overheating, and the overload; and an output unit, configured to adjust the energy parameters of the radio frequency ablation device according to the knife head states, or prompt the user to adjust the energy parameters of the radio frequency ablation device.

9. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 8, wherein the preprocessing in step S2 comprises: step 2.1: removing a high-frequency electrical noise by a low-pass filtering; step 2.2: eliminating a signal drift by a high-pass filtering; and step 2.3: removing a non-physiological pseudo-signal by a band-pass filtering.

10. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 8, wherein the Hilbert-Huang transform in step S3 comprises the following steps: S3.1: performing an empirical mode decomposition (EMD) on the signals of the two window stages of the ablation and the cutting, and the coagulation and the heat control respectively, to obtain limited intrinsic mode function (IMF) component expressions of a corresponding signal s(t);
s(t)=Σ.sub.i=1.sup.kc.sub.i(t)+r.sub.k(t), wherein I represents a number of an IMF component order, k represents the IMF component order, c.sub.i(t) is an i.sup.th order of an IMF component, and r.sub.k(t) is an i.sup.th order of an IMF residual component; S3.2: performing a Hilbert transform on each order of the IMF component c.sub.i(t), to obtain a Hilbert transform expressiony.sub.i(t) of the IMF component; $y_i\left(t\right)=\frac{p}{\pi }{\int }_{-\infty }^{+\infty }\frac{c_i\left(t\right)}{t-\tau }d\tau$ wherein p is a Cauchy principal value, and τ represents an amount of change per unit time; S3.3: obtaining an instantaneous phase φ.sub.i(t) and an instantaneous frequency ω.sub.i(t) of each order of the IMF component by using the following formula; ${\phi }_i\left(t\right)=\arctan \left(y_i\left(t\right)/c_i\left(t\right)\right)$ ${\omega }_i\left(t\right)=\frac{d{\phi }_i\left(t\right)}{\mathrm{dt}};$ S3.4: comparing and screening a value of each frequency point in the instantaneous frequency ω.sub.i(t) with an instantaneous frequency threshold), λ, if the instantaneous frequency |ω.sub.i|>λ, keeping the IMFi component unchanged, otherwise making the IMFi component be 0; S3.5: distinguishing between a square wave and a harmonic wave according to the value of each frequency point, and recording a spectrum and spectrum characteristics of the IMFi component; and S3.6: reconstructing a signal s′ (t) by using a screened IMF component in step S3.4, to obtain a time-frequency characteristic of an electromagnetic signal of a reconstructed signal.

11. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 10, wherein the obtaining limited intrinsic mode function (IMF) components of the signal s(t) in step S3.1 comprises the following steps: S3.1.1: segmenting an original signal s(t) according to a preset processing window, and taking a maximum value and a minimum value of the original signal s(t) for any processing window; S3.1.2: interpolating the maximum value and the minimum value of the original signal s(t) in the processing window by using a spline interpolation, and fitting upper and lower envelopes; S3.1.3: calculating a mean m.sub.1(t) of the upper and lower envelopes, and obtaining a difference h.sub.1(t) between the original signal s(t) and the mean m.sub.1(t),wherein h.sub.1(t)=s(t)−m.sub.1(t); S3.1.4: calculating a first IMF component c.sub.1(t) and a residual signal r.sub.1(t);
c.sub.1(t)=h.sub.1(t), r.sub.1(t)=s(t)−c.sub.1(t); S3.1.5: determining whether the residual signal r.sub.1(t) meets an empirical mode decomposition (EMD) condition, that is, whether r.sub.1(t) is a monotone function, if the residual signal r.sub.1(t) meets the EMD condition, stopping processing, and marking k=1; otherwise, repeating steps S3.1.1 to S3.1.4 fork times by using h.sub.1(t) as the original signal;
h.sub.k(t)=h.sub.(k−1)(t)−m.sub.k(t), c.sub.k(t)=h.sub.k(t), r.sub.k(t)=h.sub.(k−1)(t)−c.sub.k(t); completing a smooth decomposition processing, and recording k until the residual signal r.sub.k(t) is the monotone function; and S3.1.6: converting s(t) into k intrinsic mode function component expressions: s(t)=Σ.sub.i=1.sup.kc.sub.i(t) r.sub.k(t), wherein c.sub.i(t) is an i.sup.th IMF component of the signal.

12. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 10, wherein the instantaneous frequency threshold λ in step S3.4 ranges from 0.0094 to 0.01.

13. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 10, wherein the spectrum characteristics in step S3.5 comprise: a frequency characteristic peak value, a peak-to-peak frequency-to-amplitude ratio, and a main frequency component.

14. The device for the method of the adaptive EMC-EMI radio frequency signal data processing according to claim 8, wherein the energy parameters of the radio frequency ablation device in step S5 comprise: a voltage, a current, an impedance, a temperature, a humidity, and a contact force.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The foregoing and other objectives, features and advantages of the present invention are more obvious through a more detailed description of exemplary implementations of the present invention with reference to the accompanying drawings. In the exemplary implementations of the present invention, the same reference numerals generally refer to the same parts.

[0049] FIG. 1 is a flow chart of data processing according to the present invention.

[0050] FIG. 2 is a schematic diagram of a waveform comparison between the adaptive method of the present invention when applied and a conventional radio frequency ablation device.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0051] The following describes in detail exemplary implementations of the present invention with reference to the accompanying drawings. Although the accompanying drawings show the exemplary implementations of the present invention, it should be understood that the present invention may be implemented in various manners and should not be limited by the implementations described herein.

[0052] The present invention provides a method of adaptive EMC-EMI radio frequency signal data processing. The method includes the following steps:

[0053] S1. A processor responds to a radio frequency signal of a radio frequency ablation device acquired by an EMC-EMI sensor, including: voltage, current, impedance, temperature, humidity, and contact force.

[0054] S2. The processor segments the radio frequency signal, where a window stage after segmentation includes one or more stages of idling, energy generation, ablation and cutting, and coagulation and heat control, and preprocesses signals of each window stage after segmentation respectively. The preprocessing includes: removing high-frequency electrical noise by low-pass filtering; eliminating signal drift by high-pass filtering; and removing a non-physiological pseudo-signal by band-pass filtering.

[0055] S3. Perform Hilbert-Huang transform on signals of the two window stages of ablation and cutting, and coagulation and heat control after segmentation, to obtain time-frequency characteristics of electromagnetic signals of the corresponding window stages.

[0056] S4. Calculate EMC power, EMI radio frequency energy and mode, and a radio frequency signal-to-noise mode according to the time-frequency characteristics of electromagnetic signals; and compare the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode with corresponding thresholds respectively, to obtain knife head states in the window stages of ablation and cutting, and coagulation and heat control, including touching metal short circuit, overheating and overload.

[0057] S5. The processor adjusts energy parameters of the radio frequency ablation device according to the knife head states obtained in step S4, or prompts a user to adjust the energy parameters of the radio frequency ablation device.

[0058] In this embodiment, according to the method and device of the present invention, the signal can be effectively segmented, the time-frequency characteristics of electromagnetic signals can be accurately recognized, information such as the EMC power, EMI radio frequency energy and mode, and radio frequency signal-to-noise mode can be obtained, and an actual situation of the radio frequency ablation device in each window state can be obtained. According to the actual situation of the device, parameters such as related current and voltage adapted to output are matched, and the user is prompted to implement accurate and safe energy control.

[0059] Further, the Hilbert-Huang transform specifically includes the following steps:

[0060] S3.1. Perform empirical mode decomposition (EMD) on the signals of the window stages of ablation and cutting, and coagulation and heat control respectively, to obtain limited intrinsic mode function (IMF) component expressions of a corresponding signal s(t).

s(t)=Σ.sub.i=1.sup.kc.sub.i(t)+r.sub.k(t),

[0061] where i represents a number of an IMF component order, k represents the IMF component order, c.sub.i(t) is an i.sup.th order of IMF component, and r.sub.k(t) is an i.sup.th order of IMF residual component. The obtaining limited intrinsic mode function (IMF) components of a signal s(t) specifically includes the following steps:

[0062] S3.1.1. Segment an original signal s(t) according to a preset processing window, and take a maximum value and a minimum value of the original signal s(t) for any processing window;

[0063] S3.1.2. Interpolate the maximum value and the minimum value of the original signal s(t) in the processing window by using spline interpolation, and fit upper and lower envelopes.

[0064] S3.1.3. Calculate a mean m.sub.1(t) of the upper and lower envelopes, and obtain a difference h.sub.1(t) between the original signal s(t) and the mean m.sub.1(t); (t)=s(t)−m.sub.1(t);

[0065] S3.1.4. Calculate the first IMF component c.sub.1(t) and a residual signal r.sub.1(t);

c.sub.1(t)=h.sub.1(t), r.sub.1(t)=s(t)−c.sub.1(t);

[0066] S3.1.5. Determine whether the residual signal r.sub.1(t) meets an empirical mode decomposition (EMD) condition, that is, whether r.sub.1(t) is a monotone function, if the residual signal r.sub.1(t) meets the condition, stop processing, and mark k=1;

[0067] otherwise, repeat steps S3.1.1 to S3.1.4 for k times by using h.sub.1(t) as the original signal;

h.sub.k(t)=h.sub.(k−1)(t)−m.sub.k(t), c.sub.k(t)=h.sub.k(t), r.sub.k(t)=h.sub.(k−1)(t)−c.sub.k(t);

[0068] complete smooth decomposition processing, and record k until the residual signal r.sub.k(t) is the monotone function.

[0069] S3.1.6. Convert s(t) into k intrinsic mode function component expressions:

[0070] s(t)=Σ.sub.i=1.sup.kc.sub.i(t)+r.sub.k(t), where c.sub.i(t) is an i.sup.th IMF component of the signal.

[0071] S3.2. Perform Hilbert transform on each order of IMF component c.sub.i(t), to obtain a Hilbert transform expression of the IMF component y.sub.i(t);

[00003]$y_i\left(t\right)=\frac{p}{\pi }{\int }_{-\infty }^{+\infty }\frac{c_i\left(t\right)}{t-\tau }d\tau$

[0072] where p is a Cauchy principal value, and T represents an amount of change per unit time;

[0073] S3.3. Obtain an instantaneous phase φ.sub.i(t) and an instantaneous frequency ω.sub.i(t) of each order of IMF component by using the following formula;

[00004]${\phi }_i\left(t\right)=\arctan \left(y_i\left(t\right)/c_i\left(t\right)\right)$${\omega }_i\left(t\right)=\frac{d{\phi }_i\left(t\right)}{\mathrm{dt}}.$

[0074] S3.4. Compare and screen a value of each frequency point in the instantaneous frequency ω.sub.i(t) with an instantaneous frequency threshold λ, if the instantaneous frequency |ω.sub.i|>λ, keep the IMFi component unchanged, otherwise make the IMFi component be 0, where the instantaneous frequency threshold λ ranges from 0.0094 to 0.01;

[0075] S3.5. Distinguish between a square wave and a harmonic wave according to the value of each frequency point, and record a spectrum and spectrum characteristics of the IMFi component, where the spectrum characteristics include: a frequency characteristic peak value, a peak-to-peak frequency-to-amplitude ratio, and a main frequency component.

[0076] S3.6. Reconstruct a signal s′(t) by using a screened IMF component in step S3.4, to obtain a time-frequency characteristic of an electromagnetic signal of the reconstructed signal.

[0077] In this embodiment, according the first part of the Hilbert-Huang transform, a non-linear and non-stationary signal is decomposed into limited intrinsic mode function (IMF) components, and a mean delay on extremums of two ends of the window signal can be used to achieve a constraint on a signal endpoint and reduce a fitting error near the endpoint. According to the second part of the Hilbert-Huang transform, Hilbert spectrum analysis (HHP) is to perform the Hilbert transform on effective IMFi, use threshold determination before the signal is reconstructed, and separate an effective signal and a decomposition error or noise interference. In data analysis of a plasma knife, signal characterization of few effective IMF components such as electromagnetic conduction and electromagnetic interference are centrally examined. Separated noise interference is eliminated, which effectively improves accuracy of data tracking.

[0078] A device for the method of adaptive EMC-EMI radio frequency signal data processing. The device includes:

[0079] a radio frequency signal segmentation unit, configured to respond to a radio frequency signal of a radio frequency ablation device acquired by an EMC-EMI sensor, segment the radio frequency signal, where a window stage after segmentation includes one or more stages of idling, energy generation, ablation and cutting, and coagulation and heat control, and preprocess signals of each window stage after segmentation respectively;

[0080] a Hilbert-Huang transformation unit, configured to perform Hilbert-Huang transform on signals of the two window stages of ablation and cutting, and coagulation and heat control after segmentation, to obtain time-frequency characteristics of electromagnetic signals of the corresponding window stages;

[0081] a threshold comparison unit, configured to calculate EMC power, EMI radio frequency energy and mode, and a radio frequency signal-to-noise mode according to the time-frequency characteristics of electromagnetic signals; and compare the EMC power, the EMI radio frequency energy and mode, and the radio frequency signal-to-noise mode with corresponding thresholds respectively, to obtain knife head states in the window stages of ablation and cutting, and coagulation and heat control, including touching metal short circuit, overheating and overload; and

[0082] an output unit, configured to adjust energy parameters of the radio frequency ablation device according to the knife head states, or prompt a user to adjust the energy parameters of the radio frequency ablation device.

[0083] During specific implementation,

[0084] as shown in FIG. 1, EMC-EMI signal-noise analysis includes segmentation of a signal waveform in the radio frequency stage, calculation of current/voltage/energy windows, calculation of statistical parameters and indicators (mean, STD, variation and variability, spectral ratio, and the like), multi-calculation indicator combination and information fusion, and the like. A radio frequency ablation process can be divided into four stages from W1 to W4: W1 is an idling stage; W2 is a plasma generation and stabilization stage; W3 is a radio frequency ablation and cutting stage; and W4 is a radio frequency coagulation and heat control stage. Various different methods can be used to calculate, estimate and characterize a radio frequency factor, and to evaluate radio frequency transmission efficiency. Different results can be fused together, which facilitates clinical use to obtain better and more accurate clinical diagnosis.

[0085] In this embodiment, radio frequency of two different types of radio frequency plasma generation mechanisms are used, which are an adaptive radio frequency ablation surgery system and a commonly used radio frequency ablation device. The adaptive radio frequency ablation surgery system uses an adaptive radio frequency plasma generation mechanism based on closed-loop feedback control of wireless EMC-EMI signal acquisition and calculation. The commonly used radio frequency ablation device uses a radio frequency plasma generation mechanism with preset voltage set value and energy level. Comparison is made by monitoring plasma generation and stabilization stages (W2) of the two. As shown in FIG. 2, the upper part in the figure is a conventional radio frequency ablation device; and the lower part in the figure is an adaptive radio frequency ablation surgery system.

[0086] In this embodiment, both radio frequency ablation generators have the same initial radio frequency level set value of 4. After enabled, the commonly used radio frequency ablation device generates a stable 190 V RMS voltage and there are some flashing plasma flames at a plasma knife head in normal saline, which means that when the radio frequency level set value is 4, no sufficient voltage can be transmitted to stably generate the plasma. The adaptive radio frequency ablation surgery system is equipped with a wireless EMC-EMI sensor, converter, calculator, and communication module. An initial phenomenon of the adaptive radio frequency ablation surgery system is similar to that of the commonly used radio frequency ablation device. However, after the radio frequency is activated for twice to three times, plasma of the plasma knife head becomes stable and efficient, and provides feedback control to the clinical user, to increase to 230 V RMS voltage. By monitoring a mode of flashing and ignition of the plasma knife head in the W2 stage, the wireless EMC-EMI module can evaluate and calculate a minimum RMS voltage required to generate stable radio frequency plasma in a particular normal saline environment or clinical situation. From radio frequency output voltage monitoring, it can be seen that an adaptive plasma radio frequency ablation system has a rapid response and adjustment capability for generation and stabilization of the radio frequency plasma.

[0087] In a case that the radio frequency is relatively low, nerves pull muscles due to the stress reaction, resulting in muscular tissue shaking, causing risk in surgery, which is more obvious especially in deep cutting of deep muscles below 9 cm, and about 100 kHz. Through adaptive energy control, the radio frequency is increased to 350 kHz or more, which can effectively reduce neuromuscular shaking by more than 98%, and further reduce the risk in surgery.

[0088] This embodiment merely shows that EMC-EMI signal-noise is used to track/monitor/control voltage transmission stably generated by the plasma. In an actual application, data and waveform of the EMC-EMI signal-noise can provide more detailed information in radio frequency ablation, such as current limitation, ablation duration, and coagulation depth control, and even can predict optimal and minimal energy transmission required for an organ tissues.

[0089] The embodiments of the present invention are described above, and the foregoing descriptions are exemplary but not exhaustive and are not limited to the disclosed embodiments. Without departing from the scope and spirit of the described embodiments, many modifications and variations are apparent to a person of ordinary skill in the technical field.