Circulating hot compress control method for improving external application curative effect of Zhuang Yao medicine

Abstract

The invention discloses a circulating hot compress control method for improving the external application curative effect of Zhuang Yao medicine, and relates to the technical field of circulating hot compress. A double-point thermocouple sensor and an NTC thermistor are cooperatively arranged, and time lag and local temperature dynamic characteristics in the skin heat conduction process are accurately recognized; the delay behavior between surface layer heat absorption and deep layer heat diffusion can be distinguished in detail. By collecting and fitting four thermal behavior parameters of skin temperature rising slope delta pT, thermal response inertia time Tr, local heat conduction delay time Tag and maximum temperature difference delta Tu, a complete characteristic curve of a hot compress start-peak value-lag-attenuation process is established, and multi-dimensional comprehensive description of the skin heat sensitivity state is realized. Through the comprehensive evaluation index phi cyc, the system can accurately identify whether the current hot compress is effective or deviates, so that a definite and quantifiable feedback basis is provided for doctors or users, and the formation of a thermal therapy curative effect file and a data record is facilitated.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of circulating hot compress, in particular to a circulating hot compress control method for improving the external application effect of Zhuang and Yao medicine. BACKGROUND

[0002] As an important part of the medical system with Chinese minority characteristics, Zhuang and Yao medicine is often used to treat diseases such as rheumatism, arthralgia and soft tissue injury. Because it emphasizes the penetration of the body surface and the absorption of the meridians, it is often used in combination with hot compress to promote the release of drug efficacy. However, the traditional hot compress method relies on fixed time or fixed temperature control, and does not consider the heat sensitivity changes of skin tissue in the continuous hot compress process, which easily leads to the attenuation of therapeutic effect or the decrease of skin tolerance, and therefore it is urgent to establish a dynamic regulation mechanism coupled with the individual heat response state to ensure the precise matching of the effective release window of the external application medicine and the heat stimulation time.

[0003] In the existing practice of Zhuang and Yao medicine external application heat therapy, hot compress usually adopts a fixed time interval and constant temperature strategy, without distinguishing the heat dissipation speed of the patient's skin, the deep heat accumulation state and the heat sensitivity recovery rhythm, resulting in that most of the hot compress process deviates from the individual true heat response state. Specifically, when the skin has not fully dissipated heat, the hot compress is restarted, which may cause local heat sensitivity numbness or heat retention overload, and vice versa, if the hot compress is delayed for too long, the high efficiency period of drug release will be missed.

[0004] The above technical deficiencies mainly result from the lack of real-time identification capability of the traditional hot compress control system for the cooling state of the skin tissue, and the heat control logic often ignores the dynamic recovery characteristics of the skin after heat stimulation, especially without establishing a logical correlation model between the temperature drop rate, heat response inertia and heat conduction lag. In actual use, this hot compress method which ignores the individual heat response differences is prone to the following abnormal phenomena: the hot compress starts too early, the local heat capacity is not fully released, causing deep tissue heat retention, which may cause microcirculation disorder for a long time; the hot compress is too late, the skin temperature has fallen below the basic value, the drug penetration effect is interrupted, and the therapeutic effect is decreased; the heat rhythm is not matched for a long time, the user's subjective experience is poor, and the skin appears numb or periodic stimulation discomfort. These adverse reactions not only limit the applicability of Zhuang and Yao medicine external application in the treatment of chronic diseases and the promotion of national medicine, but also put forward higher response adaptation requirements for the intelligentization of heat therapy equipment. SUMMARY

[0005] In view of the deficiencies of the prior art, the present application provides a circulating hot compress control method for improving the external application effect of Zhuang and Yao medicine, which solves the problems mentioned in the background art.

[0006] To achieve the above purpose, the present application is implemented by the following technical scheme: a circulating hot compress control method for improving the external application effect of Zhuang and Yao medicine, comprising the following steps:

[0007] S1, when starting the Zhuang Yao medicine hot compress cycle and executing the first round of constant temperature hot compress, collecting local skin thermal response data through a sensor, fitting the raw data set YS, and preprocessing to obtain the hot compress data set RY;

[0008] S2, when the hot compress device completes the hot compress and enters the interval pause period, monitoring the cooling process, collecting monitoring data, fitting the cooling data set QY, and recording the dynamic trend curve of skin temperature drop to determine the local heat retention state;

[0009] S3, feature extraction is performed on the obtained hot compress data RY and cooling data set QY, recovery index IRI and thermal coupling activity index Ψact after the current round of hot compress are calculated, and it is determined whether the current area has recovered to the optimal state for the next heat stimulation;

[0010] S4, according to the obtained recovery index IRI and thermal coupling activity index Ψact, the interval waiting time ΔTds for the next hot compress is calculated;

[0011] S5, after obtaining the interval waiting time ΔTds, the next round of hot compress is started, and the thermal effect indicators under the new round of hot compress are recorded, including the temperature rise slope nΔpT, the new thermal conductance hysteresis time nTag, and the heat stimulation comfort index SRI;

[0012] S6, based on the obtained thermal effect indicators in the complete hot compress cycle, a comprehensive evaluation index Φcyc of cycle thermal control performance is constructed, and the thermal control strategy is evaluated.

[0013] Preferably, S1 includes S11 and S12;

[0014] S11, during the first round of constant temperature stage of Zhuang Yao medicine hot compress, data sampling is performed through a double-point thermocouple sensor and an NTC thermistor, including skin temperature rise slope ΔpT, thermal response inertia time Tr, local thermal conductance hysteresis time Tag, and maximum temperature difference ΔTu;

[0015] Among them, the local thermal conductance hysteresis time Tag is collected by a double-point thermocouple sensor;

[0016] The skin temperature rise slope ΔpT, the thermal response inertia time Tr, and the maximum temperature difference ΔTu are collected by an NTC thermistor;

[0017] The skin temperature rise slope ΔpT is obtained by the following formula:

[0018] ;

[0019] In the formula, Tt1 represents the initial temperature of hot compress start, Tt0 represents the hot compress end temperature, t1 represents the hot compress initial time, and t0 represents the hot compress end time.

[0020] The thermal response inertia time Tr is obtained by the difference between the heat application end time t0 and the time point at which the skin temperature begins to drop;

[0021] The local heat conduction lag time Tag is obtained by the difference between the time at which the maximum temperature rise rate occurs and the heat application initial time t1;

[0022] The maximum temperature difference ΔTu is obtained by the difference between the heat application end temperature Tt0 and the baseline temperature of the skin before heat application;

[0023] The obtained skin temperature rise slope ΔpT, thermal response inertia time Tr, local heat conduction lag time Tag, and maximum temperature difference ΔTu are fitted to obtain an original data set YS;

[0024] S12, noise point elimination and normalization processing are performed on the original data set YS to obtain a heat application data set RY;

[0025] Noise point elimination is performed by using the sliding window mean method and the three times standard deviation method to filter out abnormal values in the original data set YS;

[0026] Normalization processing is performed on the original data set YS by using the Min-Max normalization method to obtain the heat application data set RY;

[0027] The heat application data set RY is obtained by the following formula:

[0028] ;

[0029] In the formula, RYb represents the bth data in the heat application data set RY, YSb represents the bth data in the original data set YS, minYSb represents the valley value of the bth data in the original data set YS, and maxYSb represents the peak value of the bth data in the original data set YS.

[0030] Preferably, S2 includes S21 and S22;

[0031] S21, after the heat application device completes one heating cycle, a natural cooling process is entered, the skin temperature drop process with time is recorded by the NTC thermistor, the skin temperature drop rate ΔsT and the total skin temperature attenuation amplitude Tk are obtained, and a temperature and time trend curve is constructed;

[0032] The skin temperature drop rate ΔsT is obtained by the following formula:

[0033] ;

[0034] In the formula, Tt2 represents the temperature at the end of the cooling sampling, t2 represents the time at the end of the cooling sampling, Tt0 represents the heat application end temperature, and t0 represents the heat application end time.

[0035] The total skin temperature attenuation amplitude Tk is obtained by the following formula:

[0036] ;

[0037] In the formula, N represents the total number of temperature sampling times in the cooling period, Ti+1 represents the skin temperature value at the i+1 sampling time, and Ti represents the skin temperature value at the i sampling time.

[0038] S22, fitting the obtained skin temperature drop rate ΔsT and the total skin temperature attenuation amplitude Tk to obtain a data set W;

[0039] The obtained data set W is normalized to unify the dimension, to obtain a cooling data set QY, and a dynamic trend curve of the skin temperature drop is recorded to determine the local thermal retention state of the skin.

[0040] The cooling data set QY is obtained by the following formula:

[0041] ;

[0042] In the formula, QYg represents the gth data in the cooling data set QY, Wg represents the gth data in the data set W, minWg represents the valley value of the gth data in the data set W, and maxWg represents the peak value of the gth data in the data set W.

[0043] The local thermal retention state of the skin is obtained by matching the following way:

[0044] When 0 < skin temperature drop rate ΔsT < 0.5 and 0 < total skin temperature attenuation amplitude Tk < 0.5, it indicates that the thermal retention of the skin tissue is serious.

[0045] Preferably, S3 comprises S31 and S32.

[0046] S31, feature extraction is performed on the obtained hot compress data RY and the cooling data set QY, including skin temperature drop slope ΔsT, total skin temperature attenuation amplitude Tk, thermal response inertia time Tr, and local thermal conduction lag time Tag, and integration is performed to calculate and obtain thermal coupling activity index Ψact and recovery index IRI.

[0047] The thermal coupling activity index Ψact is obtained by the following formula:

[0048] ;

[0049] In the formula, C represents a non-zero positive number, exp represents an exponential function, and α represents an exponential factor for regulating and controlling the sensitivity of lag.

[0050] The recovery index IRI is obtained by the following formula:

[0051] ;

[0052] In the formula, β represents a curve slope adjustment factor, γ1 represents a thermal response inertia penalty coefficient, γ2 represents a thermal conductance hysteresis penalty coefficient, and e represents a constant.

[0053] Preferably, S32 analyzes the obtained thermal coupling activity index Ψact and combines it with the thermal response inertia time Tr and the local thermal conductance hysteresis time Tag to obtain a recovery state activation value Ψreg.

[0054] The recovery state activation value Ψreg is obtained by the following formula:

[0055] ;

[0056] In the formula, log represents a logarithmic function.

[0057] The obtained recovery index IRI and the recovery state activation value Ψreg are compared with preset index threshold TIR and activation threshold TΨ to determine whether the current region has been recovered to an optimal state for the next thermal stimulation.

[0058] The state determination is obtained by matching in the following manner:

[0059] When the recovery index IRI is greater than or equal to the index threshold TIR and the recovery state activation value Ψreg is greater than or equal to the activation threshold TΨ, it indicates that the current region has been sufficiently recovered and has been recovered to an optimal state for the next thermal stimulation, and the next round of hot compress can be started.

[0060] When the recovery index IRI is less than the index threshold TIR or the recovery state activation value Ψreg is less than the activation threshold TΨ, it indicates that there is a risk of heat retention, and the interval time is extended.

[0061] Preferably, S4 includes S41 and S42.

[0062] S41, the obtained recovery index IRI and the thermal coupling activity index Ψact are combined for joint feature fusion to construct a thermal steady-state regulation strength function Ωc.

[0063] The steady-state regulation strength function Ωc is obtained by the following formula:

[0064] ;

[0065] In the formula, B1 represents a recovery sensitivity factor, and B2 represents a thermal coupling sensitivity factor.

[0066] Preferably, S42, according to the obtained thermal steady-state regulation strength function Ωc, the interval waiting time ΔTds of the next hot compress after the completion of the current hot compress is dynamically calculated.

[0067] The interval waiting time ΔTds of hot compress is obtained by the following formula:

[0068] ;

[0069] In the formula, Tmin represents the minimum allowable interval time, Tmax represents the maximum safety interval upper limit specified by the system, ln represents the natural logarithm function with e as the base, and k represents the curve steepness control factor.

[0070] Preferably, S5 includes S51 and S52.

[0071] S51, after the interval waiting time ΔTds of hot compress ends, a new round of hot compress is started, and the thermal effect indicators under the new round of hot compress are recorded, including the temperature rise slope nΔpT and the new heat flux lag time nTag and the thermal stimulation comfort index SRI.

[0072] The temperature rise slope nΔpT is obtained by the following formula:

[0073] ;

[0074] In the formula, nTt1 represents the initial temperature of the new round of hot compress, nTt0 represents the end temperature of the new round of hot compress, nt1 represents the initial time of the new round of hot compress, and nt0 represents the end time of the new round of hot compress.

[0075] The new heat flux lag time nTag is obtained by the difference between the time when the maximum temperature rise rate occurs and the initial temperature nt1 of the new round of hot compress.

[0076] Preferably, S52, the subjective feeling and objective skin reaction data of the user are synchronously collected during the hot compress process, are fused into the thermal stimulation comfort index SRI, and the state of the next round of hot compress is judged.

[0077] The thermal stimulation comfort index SRI is obtained by the following formula:

[0078] ;

[0079] In the formula, Use represents the user's manual score, which is a percentage ranging from 0 to 1, Ire represents the degree of local skin redness, which is obtained by AI image recognition, Mow represents the degree of microcirculation activity, which is obtained by infrared, near-infrared and skin temperature curve fine fluctuation characteristics, respectively represent the preset weight values of the user's manual score Use, the degree of local skin redness Ire and the degree of microcirculation activity Mow, and ;

[0080] The state of the next round of hot compress is obtained by matching in the following way;

[0081] When 0 < thermal stimulation comfort index SRI < 0.4, it indicates to reduce the intensity of the next round of hot compress and prolong the time interval;

[0082] When 0.4 ≤ thermal stimulation comfort index SRI < 0.7, it indicates to keep the current parameters;

[0083] When 0.7 ≤ thermal stimulation comfort index SRI < 1.0, it indicates to allow temperature to rise and the next round of hot compress to be activated in advance.

[0084] Preferably, S6 includes S61 and S62;

[0085] S61, according to the obtained thermal effect indicators in the complete hot compress cycle, including the temperature rise slope nΔpT, the new heat flux lag time nTag and the thermal stimulation comfort index SRI, a comprehensive evaluation index Φcyc is constructed;

[0086] The comprehensive evaluation index Φcyc is obtained by the following formula:

[0087]

[0088] In the formula, SRI(t) represents the thermal stimulation comfort index at time t, d represents the integral symbol, D1, D2 and D3 respectively represent the preset weight values of the temperature rise slope nΔpT, the new heat flux lag time nTag and the thermal stimulation comfort index SRI, and D1+D2+D3≤1;

[0089] S62, the obtained comprehensive evaluation index Φcyc is compared with the preset comprehensive evaluation threshold Φref, and the state of the current hot compress is judged;

[0090] The state of the hot compress is obtained by matching the following way:

[0091] When the comprehensive evaluation index Φcyc ≥ the comprehensive evaluation threshold Φref, it indicates that the hot compress state is normal and no adjustment is needed

[0092] When the comprehensive evaluation index Φcyc < the comprehensive evaluation threshold Φref, it indicates that the hot compress state is abnormal, and the thermal response inertia penalty coefficient γ1, the heat flux lag penalty coefficient γ2, the recovery sensitive factor B1 and the thermal coupling sensitive factor B2 are adjusted to obtain a new thermal response inertia penalty coefficient nγ1, a new heat flux lag penalty coefficient nγ2, a new recovery sensitive factor nB1 and a new thermal coupling sensitive factor nB2;

[0093] The adjustment formula is as follows:

[0094] nγ1=γ1+(1+η1*(1-Φcyc));

[0095] nγ2=γ2+(1+η1*(1-Φcyc));

[0096] ​nB1 = B1 + (1 + η2 * (1 - Φ cyc ) ) ;

[0097] nB2 = B2 + (1 + η2 * (1 - Φ cyc ) ) ;

[0098] In the formula, η1 represents a recovery factor adjustment step coefficient, and η2 represents a thermal hysteresis factor adjustment step coefficient.

[0099] The application provides a circulation hot compress control method for improving the external application therapeutic effect of Zhuang Yao medicine.

[0100] (1) The double-point thermocouple sensor and NTC thermistor are cooperatively arranged, time lag and local temperature dynamic characteristics in the skin heat conduction process are accurately identified, delay behavior between surface heat absorption and deep heat diffusion can be finely distinguished, and the heat conduction time efficiency identification means cannot be realized by the traditional single-point temperature measurement scheme. Through collection and fitting of four thermal behavior parameters, i.e., skin temperature rising slope ΔpT, thermal response inertia time Tr, local heat conduction lag time Tag and maximum temperature difference ΔTu, a characteristic curve of the hot compress starting-peak-lag-decay process is established, multi-dimensional comprehensive description of the skin heat sensitivity state is realized, and the single index control is significantly superior.

[0101] (2) The data group W is formed and normalized into the cooling data set QY, quantitative and standardized input parameters are provided for the recovery index IRI construction, the adjustment of the subsequent interval time is more scientific and responsive. The traditional thermotherapy strategy only establishes a heat control model for the "heating period", the parameterized modeling of the "cooling period" is realized in the embodiment, the closed loop completion of heating-cooling-judgment in the hot compress cycle is realized for the first time, and the long-term learning type heat control strategy is beneficial to be constructed. The heat retention state is automatically identified through numerical rules, the system can actively prolong the interval or reduce the heating intensity when detecting abnormal cooling behavior, active regulation intervention is realized without increasing the burden of the user, and the system intelligence level is improved.

[0102] (3) The recovery index IRI and the thermal coupling activity index Ψ act are fused, the thermal steady state control strength function Ω c is constructed, and individualized modeling of the hot compress rhythm is realized. Compared with the previous hot compress mode with fixed or experienced interval length, the method can dynamically adjust the next hot compress time according to the actual thermal reaction and recovery capacity of the skin, and the adaptability and safety of the system are significantly enhanced. The logarithmic function and the curve control parameter are used in the embodiment, the thermal steady state control strength function Ω c is nonlinearly mapped, the interval waiting time ΔT ds of the hot compress is generated, the "minimum effective recovery time" and "maximum safe waiting upper limit" can be smoothly connected, and the emergence of potential overheating or invalid window is avoided while the therapeutic effect is ensured.

[0103] (4) Through the comprehensive evaluation index Φcyc, the system can accurately identify whether the current hot compress is effective or whether deviation is generated, provide clear and quantifiable feedback basis for doctors or users, and no longer rely on subjective experience to judge the curative effect, which is beneficial to form a hot compress curative effect file and data record.

[0104] By establishing the thermal steady-state function and the response rule, the scheme provides a structured template for the quantification control and curative effect evaluation of hot compress treatment, and is expected to form an intelligent control standard for the Zhuang Yao medicine hot compress equipment, and support the leap from family auxiliary health care to professional medical auxiliary treatment. BRIEF DESCRIPTION OF DRAWINGS

[0105] Figure 1 A circulation hot compress control method for improving the curative effect of Zhuang Yao medicine external application. DETAILED DESCRIPTION

[0106] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of the present application.

[0107] Embodiment 1

[0108] The present application provides a circulation hot compress control method for improving the curative effect of Zhuang Yao medicine external application, please refer to Figure 1 , comprising the following steps:

[0109] S1, when starting the Zhuang Yao medicine hot compress circulation and performing the first round of constant temperature hot compress, collecting the local skin thermal response data through a sensor, fitting the raw data set YS, and pre-processing to obtain the hot compress data set RY;

[0110] S2, when the hot compress device completes the hot compress and enters the interval pause period, monitoring the cooling process, collecting the monitoring data, fitting the cooling data set QY, recording the dynamic trend curve of the skin temperature drop, and judging the local heat retention state;

[0111] S3, extracting the features of the obtained hot compress data RY and cooling data set QY, calculating the recovery index IRI and the heat coupling activity index Ψact after the current hot compress, and judging whether the current area has recovered to the optimal state for the next heat stimulation;

[0112] S4, according to the obtained recovery index IRI and heat coupling activity index Ψact, calculating the interval waiting time ΔTds for the next hot compress;

[0113] S5, after the acquisition interval waiting time ΔTds, start the next round of hot compress, and record the thermal effect index under the new round of hot compress, including the temperature rise slope nΔpT, the new heat flux lag time nTag and the thermal stimulation comfort index SRI;

[0114] S6, based on the acquired thermal effect index in the complete hot compress cycle, construct the comprehensive evaluation index Φcyc of the cycle thermal control performance, and evaluate the thermal control strategy.

[0115] In this embodiment, by controlling the precise synchronization of the hot compress rhythm and the cooling time, it is ensured that the external application of the medicine is exactly in the strongest interval of skin thermal sensitivity, and the penetration ability of traditional Chinese medicine is enhanced. The hot compress interval time is no longer artificially set, but is generated based on the recovery situation judgment, which can effectively avoid skin tolerance, inflammation or microcirculation disorder caused by heat accumulation.

[0116] Each user has different skin thermal response behavior, and the system can dynamically collect parameters and automatically adjust the control strategy, so that the device is suitable for a wider population, including the elderly, cold and rheumatism and other populations.

[0117] The comfort index is embedded in the regulation link, which can automatically reduce the thermal stimulation intensity or delay the hot compress time according to the user feedback, to prevent discomfort or abandonment caused by "excessive thermal stimulation". The comprehensive evaluation index Φcyc performance function is introduced, which provides a quantitative score for each round of hot compress results, which can not only be used for internal strategy optimization of the device, but also can be used as a reference index for clinicians to judge the efficacy.

[0118] Embodiment 2

[0119] This embodiment is an explanation and description in embodiment 1, please refer to Figure 1 , in particular: S1 includes S11 and S12;

[0120] S11, in the first round of constant temperature phase of Zhuang Yao medicine hot compress, data sampling is carried out through double-point thermocouple sensor and NTC thermistor, including skin temperature rise slope ΔpT, thermal response inertia time Tr, local heat flux lag time Tag and maximum temperature difference ΔTu;

[0121] Among them, the local heat flux lag time Tag is collected by double-point thermocouple sensor;

[0122] The skin temperature rise slope ΔpT, the thermal response inertia time Tr and the maximum temperature difference ΔTu are collected by NTC thermistor;

[0123] The skin temperature rise slope ΔpT is obtained by the following formula:

[0124] ;

[0125] In the formula, Tt1 represents a hot compress starting initial temperature, Tt0 represents a hot compress ending temperature, t1 represents a hot compress initial time, and t0 represents a hot compress ending time.

[0126] The thermal response inertia time Tr is obtained by the difference between the hot compress ending time t0 and the time point at which the skin temperature starts to drop;

[0127] The local thermal conductance lag time Tag is obtained by the difference between the time at which the maximum temperature rise rate occurs and the hot compress initial time t1.

[0128] The maximum temperature difference ΔTu is obtained by the difference between the hot compress ending temperature Tt0 and the baseline temperature of the skin before the hot compress is started.

[0129] The obtained skin temperature rise slope ΔpT, thermal response inertia time Tr, local thermal conductance lag time Tag, and maximum temperature difference ΔTu are fitted to obtain an original data set YS.

[0130] S12, noise point elimination and normalization processing are performed on the original data set YS to obtain a hot compress data set RY.

[0131] The noise point elimination is performed by using a sliding window mean method and a three times standard deviation method to filter out abnormal values in the original data set YS.

[0132] The normalization processing is performed by using a Min-Max normalization method to process the original data set YS to obtain the hot compress data set RY.

[0133] The hot compress data set RY is obtained by the following formula:

[0134]

[0135] In the formula, RYb represents the bth data in the hot compress data set RY, YSb represents the bth data in the original data set YS, minYSb represents the valley value of the bth data in the original data set YS, and maxYSb represents the peak value of the bth data in the original data set YS.

[0136] The embodiment adopts a double-point thermocouple sensor and an NTC thermistor to cooperatively arrange, accurately identifies the time lag and local temperature dynamic characteristics in the skin heat conduction process, can finely distinguish the delay behavior between the surface heat absorption and the deep heat diffusion, and is a heat conductance time effect identification means that cannot be realized by a traditional single-point temperature measurement scheme. Through the collection and fitting of the four thermal behavior parameters of the skin temperature rise slope ΔpT, thermal response inertia time Tr, local thermal conductance lag time Tag, and maximum temperature difference ΔTu, a complete characteristic curve of the hot compress starting-peak-lag-decay process is established, multi-dimensional comprehensive description of the skin heat sensitivity state is realized, and the embodiment is significantly superior to single-index control.

[0137] ​The abnormal points in the original data are removed by the sliding window mean method and the three standard deviation method, combined with Min-Max normalization processing, to construct a standardized hot compress data set RY, which can effectively avoid the misleading of sensor transient fluctuations, environmental interference and other factors on the heat control strategy, and effectively improve the stability and contrast of subsequent modeling. Through high-frequency sampling and multi-parameter synchronous acquisition, this embodiment can clearly capture various physical response processes of the skin after the hot compress is started, including the initial transmission of heat flow, thermal buffering of the organization, peak delay and heat dissipation turning point, etc., to ensure that the subsequent strategy can be based on the real physiological state and not "set assumptions".

[0138] By calculating the thermal response inertia time Tr and the local heat conduction lag time Tag and other parameters, the system can identify the degree of acceptance and recovery speed of the skin to the heat stimulus, and then judge whether it is necessary to extend the interval or adjust the hot compress intensity, get rid of the rigid design of the traditional scheme "time cycle", and truly realize the change according to the person and the adjustment according to the area. The hot compress data set RY after the removal of outliers and normalization processing has good input consistency and contrast reference, and can be directly used to model the dynamics function of the heat response-cooling behavior, improve the accuracy of the recovery index IRI and the heat coupling activity index Ψ act calculation, and ensure that the entire system has continuous closed-loop regulation capability.

[0139] Embodiment 3

[0140] This embodiment is an explanation and description in embodiment 2, please refer to Figure 1 , specifically: S2 includes S21 and S22;

[0141] S21, after the hot compress device completes a heating cycle, enters a natural cooling process, records the skin temperature drop process with time through the NTC thermistor, obtains the skin temperature drop rate ΔsT and the total skin temperature attenuation amplitude Tk, and constructs the temperature and time trend curve;

[0142] The skin temperature drop rate ΔsT is obtained by the following formula:

[0143] ;

[0144] In the formula, Tt2 represents the temperature at the end of the cooling sampling, t2 represents the time at the end of the cooling sampling, Tt0 represents the hot compress end temperature, and t0 represents the hot compress end time;

[0145] The total skin temperature attenuation amplitude Tk is obtained by the following formula:

[0146] ;

[0147] In the formula, N is the total number of temperature sampling times in the cooling period, Ti+1 represents the skin temperature value at the i+1 sampling time, and Ti represents the skin temperature value at the i sampling time.

[0148] S22, fitting the obtained skin temperature drop rate ΔsT and the total skin temperature attenuation amplitude Tk to obtain a data set W;

[0149] The obtained data set W is normalized to unify the dimension, to obtain a cooling data set QY, and a dynamic trend curve of the skin temperature drop is recorded to determine the local thermal retention state of the skin;

[0150] The cooling data set QY is obtained by the following formula:

[0151] ;

[0152] In the formula, QYg represents the gth data in the cooling data set QY, Wg represents the gth data in the data set W, minWg represents the valley value of the gth data in the data set W, and maxWg represents the peak value of the gth data in the data set W;

[0153] The local thermal retention state of the skin is obtained by matching in the following way:

[0154] When 0 < skin temperature drop rate ΔsT < 0.5 and 0 < total skin temperature attenuation amplitude Tk < 0.5, it indicates that the skin tissue is in a serious thermal retention state.

[0155] The embodiment continuously monitors the skin temperature during the "pause period" of the hot compress cycle for the first time, rather than assuming that the skin cooling state is "completed". The cooling process is recorded in real time by the NTC thermistor, the skin temperature drop rate and the total attenuation amplitude are quantified, so as to realize the accurate perception of the skin heat release process. Unlike the traditional thermal control strategy which only focuses on the heating stage, the present scheme constructs a dynamic model of the skin heat dissipation process through the temperature-time curve, which can identify whether the skin is still in a thermal retention state, especially when the skin temperature drop amplitude is small or the speed is slow, and timely capture the "latent heat residual" phenomenon of the tissue which is not fully dissipated.

[0156] The embodiment sets a reasonable two-parameter thermal retention criterion, which is based on the double lower limit judgment of the drop rate and the total attenuation amplitude, to automatically determine whether the current skin has obvious thermal retention, thereby effectively avoiding the thermal control time mismatch problem caused by human subjective intervention or false inference. Through the normalization processing of the multi-source data in the cooling period, the cooling data set QY in a unified scale is formed, which provides a structured input for the subsequent thermal response model, and effectively improves the data compatibility, expandability and control logic consistency of the whole system.

[0157] The embodiment can determine whether the skin has truly completed the heat dissipation process, avoid reheating when the tissue still retains heat, prevent the superposition of heat stimulation from causing microcirculation disorder, heat numbness, redness, or burning pain, and significantly improve the safety of hyperthermia. The cooling rate and decay amplitude, as specific manifestations of the individual skin heat diffusion capacity, can reflect differences in age, skin quality, capillary density, and other factors, enabling the device to adjust the pace based on the "local true response" rather than relying on general settings, thereby enhancing individualized control capabilities.

[0158] By forming the data set W and normalizing it to the cooling data set QY, the embodiment provides quantitative and standardized input parameters for the recovery index IRI construction, making the subsequent interval time adjustment more scientific and responsive. Traditional hyperthermia strategies only establish a thermal control model for the "heating period", and the embodiment parameterizes the "cooling period" modeling for the first time, realizing the closed-loop completion of heating-cooling-judgment in the hot compress cycle, which is conducive to the construction of long-term learning type thermal control strategy. Through numerical rules to automatically identify heat retention state, the system can actively extend the interval or reduce the heating intensity when detecting abnormal cooling behavior, thereby realizing active regulation and intervention without increasing user burden, and improving the system intelligence level.

[0159] Embodiment 4

[0160] This embodiment is an explanation and description in embodiment 3, please refer to Figure 1 , in particular: S3 includes S31 and S32;

[0161] S31, feature extraction is performed on the obtained hot compress data RY and cooling data set QY, including skin temperature drop slope ΔsT, skin temperature total decay amplitude Tk, thermal response inertia time Tr, and local heat conduction lag time Tag, and integration is performed to calculate and obtain thermal coupling activity index Ψact and recovery index IRI;

[0162] The thermal coupling activity index Ψact is obtained by the following formula:

[0163] ;

[0164] In the formula, C represents a non-zero positive number, exp represents an exponential function, and a represents an exponential factor for regulating lag sensitivity;

[0165] The recovery index IRI is obtained by the following formula:

[0166] ;

[0167] In the formula, β represents a curve slope adjustment factor, γ1 represents a thermal response inertia penalty coefficient, γ2 represents a heat conduction lag penalty coefficient, and e represents a constant.

[0168] S32, analyze the obtained thermal coupling activity index Ψact, and combine the thermal response inertia time Tr and the local heat conduction lag time Tag to obtain a recovery state activation value Ψreg;

[0169] The recovery state activation value Ψreg is obtained by the following formula:

[0170] ;

[0171] In the formula, log represents a logarithmic function;

[0172] The obtained recovery index IRI and the recovery state activation value Ψreg are compared with preset index threshold TIR and activation threshold TΨ to determine whether the current area has been recovered to an optimal state for the next heat stimulation;

[0173] The state determination is matched by the following method:

[0174] When the recovery index IRI is greater than or equal to the index threshold TIR, and the recovery state activation value Ψreg is greater than or equal to the activation threshold TΨ, it indicates that the current area has been fully recovered, and is in an optimal state for the next heat stimulation, and the next round of hot compress can be started.

[0175] When the recovery index IRI is less than the index threshold TIR, or the recovery state activation value Ψreg is less than the activation threshold TΨ, it indicates that there is a risk of heat retention, and the interval time is extended.

[0176] In this embodiment, it is first proposed to integrate hot compress data and cooling data in both directions to form two key indicators: thermal coupling activity index Ψact and recovery index IRI, which respectively reflect the response activity of the current skin to heat stimulation and the overall recovery situation, thereby constructing a more accurate heat control state evaluation model, which is significantly better than the traditional single temperature threshold method.

[0177] The thermal coupling activity index Ψact is coupled with the thermal response inertia Tr and the heat conduction lag time Tag to calculate the recovery state activation value Ψreg, which can reflect the dynamic evolution process of the skin from the last round of heat stimulation to the fully recovered state, avoiding the judgment of "whether to reheat" only relying on static indicators. By setting the recovery index threshold TIR and the activation threshold TΨ double judgment standard, the system no longer drives the hot compress rhythm by time, but drives the hot compress decision by the heat response state, ensuring that "only when the skin state is truly ready" the next round of heat stimulation is performed, improving the therapeutic effect while controlling the heat risk.

[0178] By combining the thermal response inertia and the conduction hysteresis characteristics, the embodiment can accurately identify the thermal residual state of the deep skin layer, avoid false triggering of heat application due to the apparent fall of the surface temperature, and effectively prevent the burning sensation, discomfort, and skin damage caused by thermal accumulation. The recovery index IRI and the recovery state activation value Ψreg can be automatically evaluated for different users and different parts of the skin, allowing the heat control strategy to adapt to different physical conditions (such as blood stasis, cold constitution, etc.) and different parts (such as the knee and the waist back) of the cooling recovery rhythm, achieving "decisions based on tissue state rather than fixed time."

[0179] The thermal coupling activity index Ψact reflects the skin thermal conduction state and the activity of thermal dynamic effects, which can be used to determine whether the drug ingredients are in the optimal release interval, achieve "precise coupling" of thermal driving and drug efficacy window, and significantly improve the penetration effect and stability of the Yao medicine.

[0180] Embodiment 5

[0181] This embodiment is an explanation and description in Embodiment 4, please refer to Figure 1 , specifically: S4 includes S41 and S42;

[0182] S41, jointly fusing the obtained recovery index IRI and thermal coupling activity index Ψact, and constructing a thermal steady-state regulation strength function Ωc;

[0183] The steady-state regulation strength function Ωc is obtained by the following formula:

[0184] ;

[0185] In the formula, B1 represents the recovery sensitivity factor, and B2 represents the thermal coupling sensitivity factor.

[0186] S42, according to the obtained thermal steady-state regulation strength function Ωc, dynamically calculating the interval waiting time ΔTds of the next heat application after the completion of the current heat application;

[0187] The interval waiting time ΔTds of the heat application is obtained by the following formula:

[0188] ;

[0189] In the formula, Tmin represents the minimum allowed interval length, Tmax represents the maximum safety interval upper limit specified by the system, ln represents the logarithm function with e as the base, and k represents the curve steepness control factor.

[0190] S5 includes S51 and S52;

[0191] S51, after the interval waiting time ΔTds of the hot compress is ended, a new round of hot compress is started, and the thermal effect indexes under the new round of hot compress are recorded, including the temperature rise slope nΔpT and the new heat flux lag time nTag and the thermal stimulation comfort index SRI;

[0192] The temperature rise slope nΔpT is obtained by the following formula:

[0193] ;

[0194] In the formula, nTt1 represents the initial temperature of the new round of hot compress, nTt0 represents the end temperature of the new round of hot compress, nt1 represents the initial time of the new round of hot compress, and nt0 represents the end time of the new round of hot compress.

[0195] The new heat flux lag time nTag is obtained by the difference between the time when the maximum temperature rise rate appears and the initial temperature nt1 of the new round of hot compress.

[0196] S52, the subjective feeling and the objective skin reaction data of the user are synchronously collected during the hot compress process, are fused into the thermal stimulation comfort index SRI, and the state of the next round of hot compress is judged.

[0197] The thermal stimulation comfort index SRI is obtained by the following formula:

[0198] ;

[0199] In the formula, Use represents the user's manual score, Ire represents the local skin redness degree, Mow represents the microcirculation activity degree, respectively represent the preset weight values of the user's manual score Use, the local skin redness degree Ire and the microcirculation activity degree Mow, and ;

[0200] The state of the next round of hot compress is obtained by matching in the following way:

[0201] When 0 < the thermal stimulation comfort index SRI < 0.4, it means that the next round of hot compress intensity is reduced and the time interval is prolonged.

[0202] When 0.4 ≤ the thermal stimulation comfort index SRI < 0.7, it means that the current parameters are kept.

[0203] When 0.7 ≤ the thermal stimulation comfort index SRI < 1.0, it means that the temperature is allowed to be increased and the next round of hot compress is allowed to be activated in advance.

[0204] In this embodiment, the heat steady-state regulation strength function Ωc is constructed by fusing the characteristics of the recovery index IRI and the heat coupling activity index Ψact, realizing the individualized modeling of the hot compress rhythm. Compared with the previous hot compress mode with fixed or empirical interval length, this method can dynamically adjust the next hot compress timing according to the actual thermal response and recovery ability of the skin, significantly enhancing the system adaptability and safety. In this embodiment, the logarithmic function and the curve control parameter are used to perform nonlinear mapping on the heat steady-state regulation strength function Ωc, generating the interval waiting time ΔTds of hot compress, which can smoothly connect the “minimum effective recovery time” and the “maximum safe waiting upper limit”, ensuring the efficacy while avoiding the occurrence of potential overheating or failure window.

[0205] In S5, the physical indicators such as the heating slope and the heat flux lag are introduced, and the user's subjective score, skin redness and microcirculation activity are collected at the same time, forming a composite heat stimulation comfort index SRI, breaking the previous limitations of relying solely on physical temperature indicators for control, and realizing the intelligent adjustment closed loop of “man-machine-skin state” trinity. According to the SRI score results, three kinds of control instructions are divided: delay cooling, parameter retention and advance heating, forming a clear control response rule chain, ensuring that the system has flexible processing capability for different users, different body feeling states and heat response periods, and continuously improving the individual hot compress experience and the degree of drug efficacy adaptation.

[0206] In this embodiment, the interval waiting time ΔTds of hot compress is regulated by the heat steady-state regulation strength function Ωc, ensuring that the skin state truly completes heat release before each round of heat stimulation starts, effectively avoiding heat overload and efficacy interruption. By collecting user subjective scores and objective responses, the system can no longer rely solely on temperature data as feedback criteria, and can adjust itself according to the actual experience of the user (such as being too hot, tingling, redness), significantly improving the individual adaptability and user stickiness of heat therapy.

[0207] If the interval timing is misaligned with the skin recovery state, it will lead to incomplete drug penetration or premature clearance. In this embodiment, the drug release window and the heat stimulation cycle are precisely aligned through bidirectional adjustment of the recovery index and the comfort index, thereby enhancing the penetration efficiency of the ingredients and prolonging the drug efficacy maintenance time. The heat flux behavior, comfort score and lag response are updated in each round of hot compress, and the system can continuously optimize the control logic accordingly, realizing the strategy evolution and individual learning mechanism, meeting the differentiated needs of long-term users at different stages.

[0208] In this embodiment, by setting the upper and lower limits of the interval time and multiple state judgments, skin damage problems such as red rash and blisters caused by excessive heat sensitivity or insufficient interval can be effectively avoided, which is particularly suitable for the elderly, children or abnormal heat sensitivity groups.

[0209] By establishing the thermal steady-state function and response rules, this scheme provides a structured template for the quantitative control and efficacy evaluation of hot compress treatment, and is expected to form an industry-standard intelligent control standard for Zhuang Yao medicine thermal therapy equipment, supporting its leap from home-assisted healthcare to professional medical auxiliary treatment.

[0210] Embodiment 6

[0211] This embodiment is an explanatory description in Embodiment 5, please refer to Figure 1 , specifically: S6 includes S61 and S62;

[0212] S61, according to the acquired thermal effect indicators in the complete hot compress cycle, including the temperature rise slope nΔpT, the new heat conduction lag time nTag and the thermal stimulation comfort index SRI, the comprehensive evaluation index Φcyc is constructed;

[0213] The comprehensive evaluation index Φcyc is obtained by the following formula:

[0214] ;

[0215] In the formula, SRI(t) represents the thermal stimulation comfort index at time t, d represents the integral symbol, D1, D2 and D3 represent the preset weight values of the temperature rise slope nΔpT, the new heat conduction lag time nTag and the thermal stimulation comfort index SRI respectively, and D1+D2+D3≤1;

[0216] S62, compare the acquired comprehensive evaluation index Φcyc with the preset comprehensive evaluation threshold Φref to judge the state of the current hot compress;

[0217] The state of the hot compress is matched and obtained by the following way:

[0218] When the comprehensive evaluation index Φcyc≥ the comprehensive evaluation threshold Φref, it means that the hot compress state is normal, and there is no need to adjust

[0219] When the comprehensive evaluation index Φcyc< the comprehensive evaluation threshold Φref, it means that the hot compress state is abnormal, and the heat response inertia penalty coefficient γ1, the heat conduction lag penalty coefficient γ2, the recovery sensitivity factor B1 and the thermal coupling sensitivity factor B2 are adjusted to obtain the new heat response inertia penalty coefficient nγ1, the new heat conduction lag penalty coefficient nγ2, the new recovery sensitivity factor nB1 and the new thermal coupling sensitivity factor nB2;

[0220] The adjustment formula is as follows:

[0221] nγ1=γ1+(1+η1*(1-Φcyc));

[0222] nγ2=γ2+(1+η1*(1-Φcyc));

[0223] nB1 = B1 + (1 + η2 * (1 - Φcyc));

[0224] nB2 = B2 + (1 + η2 * (1 - Φcyc));

[0225] In the formula, η1 represents a recovery factor adjustment step coefficient, and η2 represents a thermal hysteresis factor adjustment step coefficient.

[0226] The embodiment first unifies the three heterogeneous indexes of the temperature rise slope nΔpT, the new thermal conductance hysteresis time nTag, and the thermal stimulation comfort index SRI in a unified structure, constructs a comprehensive evaluation index Φcyc, and can not only quantify the comprehensive performance of the current hot compress period, but also breaks the limitation of the traditional thermal control system which takes a single index as a standard. The comfort index is introduced into the comprehensive evaluation index Φcyc and is processed by time integration, so that the user experience is no longer a “static single-point score”, but the fluctuation trend in the whole hot compress period is identified, which reflects the attention of the system to the dynamicity of subjective experience, and makes the evaluation system more human-oriented and time-continuous.

[0227] When the evaluation index is lower than the comprehensive evaluation threshold Φref, the system can automatically dynamically adjust a plurality of core control parameters (including thermal response and recovery factor) according to mathematical formula, form a self-adaptive learning type feedback loop of “evaluation→correction→iterative optimization”, and break the limitation of the traditional control strategy which is static and mechanical. The comprehensive evaluation index Φcyc not only serves as a single-period feedback, but also constitutes a basic trigger for long-term evolution of the system. The evaluation result of each round of thermal control state will act on the penalty coefficient and the sensitive factor of the next period, so that the system can be continuously optimized and gradually fit the skin state, drug effect and subjective feedback of the user.

[0228] Through the comprehensive evaluation index Φcyc, the system can accurately identify whether the current hot compress is effective or deviated, and provide clear and quantifiable feedback basis for doctors or users, so as to avoid relying on subjective experience to judge the curative effect, which is conducive to forming a heat treatment curative effect archive and data record.

[0229] If the system continuously detects that the comprehensive evaluation index Φcyc is low, the penalty coefficient and the sensitive factor can be adjusted in time to avoid the gradual decline of the curative effect without being found, effectively maintain the sustainability of the drug release efficiency and the skin thermal response, and reduce the risk of “cumulative failure”. The continuous feedback and optimization mechanism can make the device gradually form a “user-specific thermal control curve”, and each adjustment is closer to the individual optimal curative effect path, so as to improve the user comfort, curative effect consistency and use trust, and enhance the product stickiness and word-of-mouth effect.

[0230] While embodiments of the application have been shown and described, it is to be understood that the embodiments described are merely exemplary of the principles and application of the present application. Numerous modifications and adaptions can be effected without departing from the spirit and scope of the present application, which is not limited to the exact construction and arrangement described. It is intended, therefore, to cover all modifications and adaptions that fall within the scope of the claims and their equivalents.

Claims

1. A circulating hot compress control method for improving the efficacy of Zhuang and Yao medicine external application, characterized by: The following steps are involved: S1. When starting the Zhuangyao medicine hot compress cycle and performing the first round of constant temperature hot compress, the local skin thermal response data is collected through the sensor, fitted into the original data set YS, and preprocessed to obtain the hot compress data set RY; S2. After the hot compress device completes the hot compress, it enters the interval pause period, monitors the cooling process, collects monitoring data, fits it into the cooling data set QY, and records the dynamic trend curve of the skin temperature drop to determine the local heat retention state; S3. Perform feature extraction on the acquired hot compress data RY and cooling data set QY, calculate the recovery index IRI and thermal coupling activity index Ψact after the current round of hot compress, and determine whether the current area has recovered to the optimal state for the next thermal stimulation; S4. Calculate the waiting time ΔTds for the next hot compress based on the obtained recovery index IRI and thermal coupling activity index Ψact; S5. After obtaining the interval waiting time ΔTds, start the next round of hot compress and record the thermal effect indicators of the new round of hot compress, including the temperature rise slope nΔpT, the new heat conduction hysteresis time nTag and the thermal stimulation comfort index SRI; S6. Based on the thermal effect indicators obtained during the complete hot compress cycle, a comprehensive evaluation index Φcyc of the cycle thermal control performance is constructed to evaluate the thermal control strategy.

2. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 1, characterized in that: S1 includes S11 and S12; S11. During the first constant temperature phase of the Zhuang-Yao medicine hot compress, data is sampled using a dual-point thermocouple sensor and an NTC thermistor, including the skin temperature rise slope ΔpT, thermal response inertia time Tr, local thermal conduction hysteresis time Tag, and maximum temperature difference ΔTu. Among them, the local thermal conduction hysteresis time Tag is collected through a double-point thermocouple sensor; The skin temperature rise slope ΔpT, thermal response inertia time Tr and maximum temperature difference ΔTu are collected through the NTC thermistor; The skin temperature rise slope ΔpT is obtained by the following formula: ； Where, Tt1 represents the initial temperature of hot compress, Tt0 represents the end temperature of hot compress, t1 represents the initial time of hot compress, and t0 represents the end time of hot compress; The thermal response inertia time Tr is obtained by the difference between the end time t0 of the hot compress and the time point when the skin temperature starts to drop; The local heat conduction hysteresis time Tag is obtained by the difference between the time when the maximum heating rate occurs and the initial time t1 of the hot compress; The maximum temperature difference ΔTu is obtained by the difference between the end temperature Tt0 of the hot compress and the baseline temperature of the skin before the hot compress starts; The skin temperature rise slope ΔpT, thermal response inertia time Tr, local thermal conduction hysteresis time Tag and maximum temperature difference ΔTu are fitted to obtain the original data set YS; S12, performing noise point removal and normalization processing on the original data set YS to obtain the hot compress data set RY; Noise point removal is performed by using the sliding window mean method and the triple standard deviation method to filter out the outliers in the original data set YS; Normalization processing is performed by using the Min-Max normalization method to process the original data set YS to obtain the hot compress data set RY; The hot compress dataset RY is obtained using the following formula: ； Where RYb represents the b-th data item in the hot compress dataset RY, YSb represents the b-th data item in the original dataset YS, minYSb represents the valley value of the b-th data item in the original dataset YS, and maxYSb represents the peak value of the b-th data item in the original dataset YS.

3. The cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 1, characterized in that: S2 includes S21 and S22; S21. After the hot compress device completes a heating cycle, it enters a natural cooling process. The skin temperature decreases over time using an NTC thermistor to obtain the skin temperature decrease rate ΔsT and the total skin temperature attenuation amplitude Tk, and construct a temperature and time trend curve. The skin temperature drop rate ΔsT is obtained by the following formula: ； Where, Tt2 represents the temperature at the end of cooling sampling, t2 represents the time at the end of cooling sampling, Tt0 represents the temperature at the end of hot compress, and t0 represents the time at the end of hot compress; The total skin temperature attenuation amplitude Tk is obtained by the following formula: ； Where, N is the total number of temperature sampling times performed during the cooling period, Ti+1 represents the skin temperature value at the i+1th sampling moment, and Ti represents the skin temperature value at the i-th sampling moment; S22, fitting the obtained skin temperature drop rate ΔsT and the total skin temperature attenuation amplitude Tk to obtain a data set W; The acquired data set W is normalized and dimensionally unified to obtain the cooling data set QY. The dynamic trend curve of skin temperature drop is recorded to determine the local heat retention state of the skin. The cooling data set QY is obtained by the following formula: ； Where QYg represents the g-th data item in the cooling data set QY, Wg represents the g-th data item in the data set W, minWg represents the valley value of the g-th data item in the data set W, and maxWg represents the peak value of the g-th data item in the data set W; The local heat retention state of the skin is obtained by matching: When 0 < skin temperature drop rate ΔsT < 0.5, and 0 < total skin temperature attenuation amplitude Tk < 0.5, it indicates that the skin tissue has severe heat retention.

4. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 3, characterized in that: S3 includes S31 and S32; S31. Extract features from the acquired hot compress data RY and cooling data set QY, including the skin temperature drop slope ΔsT, the total skin temperature attenuation amplitude Tk, the thermal response inertia time Tr, and the local thermal conduction hysteresis time Tag, and integrate them to calculate the thermal coupling activity index Ψact and the recovery index IRI. The thermal coupling activity index Ψact is obtained by the following formula: ； Where C represents a non-zero positive number, exp represents an exponential function, and α represents an exponential factor for regulating hysteresis sensitivity; The recovery index IRI is obtained by the following formula: ； Where β represents the curve slope adjustment factor, γ1 represents the thermal response inertia penalty coefficient, γ2 represents the thermal conduction hysteresis penalty coefficient, and e represents a constant.

5. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 4, characterized in that: S32, analyzing the obtained thermal coupling activity index Ψact, and combining it with the thermal response inertia time Tr and the local thermal conduction hysteresis time Tag to obtain a recovery state activation value Ψreg; The restored state activation value Ψreg is obtained by the following formula: ； Where, log represents the logarithmic function; Compare the obtained recovery index IRI and recovery state activation value Ψreg with the preset index threshold TIR and activation threshold TΨ to determine whether the current area has recovered to the optimal state for the next thermal stimulation; Status judgment is obtained by matching in the following ways: When the recovery index IRI ≥ index threshold TIR, and the recovery state activation value Ψreg ≥ activation threshold TΨ, it means that the current area has fully recovered and is back to the optimal state for the next heat stimulation, and the next round of heat compress can be started; When the recovery index IRI is less than the index threshold TIR, or the recovery state activation value Ψreg is less than the activation threshold TΨ, it indicates that there is a risk of thermal retention and the interval time is extended.

6. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 5, characterized in that: S4 includes S41 and S42; S41, performing joint feature fusion on the obtained recovery index IRI and thermal coupling activity index Ψact to construct the thermal steady-state control intensity function Ωc; The steady-state control intensity function Ωc is obtained by the following formula: ； Where B1 represents the recovery sensitivity factor and B2 represents the thermal coupling sensitivity factor.

7. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 6, characterized in that: S42, dynamically calculating the waiting time ΔTds between the next hot compress after the current round of hot compress is completed based on the obtained thermal steady-state control intensity function Ωc; The waiting time ΔTds between hot compresses is obtained by the following formula: ； Where Tmin represents the minimum allowed interval duration, Tmax represents the upper limit of the maximum safety interval specified by the system, ln represents the logarithmic function with e as the base, and k represents the curve steepness control factor.

8. The method for controlling cyclic hot compress for improving the efficacy of external application of Zhuang and Yao medicine according to claim 7, characterized in that: S5 includes S51 and S52; S51. After the waiting time ΔTds between hot compresses ends, a new round of hot compresses is started, and the thermal effect indicators of the new round of hot compresses are recorded, including the temperature rise slope nΔpT, the new heat conduction hysteresis time nTag, and the thermal stimulation comfort index SRI; The heating slope nΔpT is obtained by the following formula: ； Where nTt1 represents the initial temperature for starting a new round of hot compress, nTt0 represents the end temperature for a new round of hot compress, nt1 represents the initial time for a new round of hot compress, and nt0 represents the end time for a new round of hot compress. The new heat conduction hysteresis time nTag is obtained by the difference between the time when the maximum heating rate occurs and the initial temperature nt1 of the new round of hot compress.

9. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 8, characterized in that: S52. During the hot compress process, the user's subjective feelings and objective skin reaction data are simultaneously collected and integrated into a thermal stimulation comfort index (SRI), and the status of the next round of hot compress is judged; The thermal stimulation comfort index SRI is obtained by the following formula: ； In the formula, Use represents the user's manual rating, Ire represents the degree of local skin redness, and Mow represents the degree of microcirculation activity. They represent the preset weight values ​​of the user's manual rating Use, the local skin redness level Ire, and the microcirculation activity level Mow, respectively, and ; The status of the next round of hot compress is obtained by matching the following methods; When 0 < thermal stimulation comfort index SRI < 0.4, it means reducing the intensity of the next round of hot compress and extending the time interval; When 0.4≤SRI<0.7, it means maintaining the current parameters; When 0.7≤SRI<1.0, it means that the temperature can be increased and the next round of hot compress can be activated in advance.

10. A cyclic hot compress control method for improving the efficacy of Zhuang-Yao medicine external application according to claim 9, characterized in that: S6 includes S61 and S62; S61. Construct a comprehensive evaluation index Φcyc based on the thermal effect indicators obtained during the complete hot compress cycle, including the temperature rise slope nΔpT, the new heat conduction hysteresis time nTag, and the thermal stimulation comfort index SRI; The comprehensive evaluation index Φcyc is obtained by the following formula: ； Where SRI(t) represents the thermal comfort index at time t, d represents the integral sign, D1, D2, and D3 represent the preset weight values ​​of the heating slope nΔpT, the new thermal conduction hysteresis time nTag, and the thermal comfort index SRI, respectively, and D1+D2+D3≤1; S62, comparing the obtained comprehensive evaluation index Φcyc with the preset comprehensive evaluation threshold Φref to determine the current state of the hot compress; The status of the hot compress is obtained by matching: When the comprehensive evaluation index Φcyc ≥ comprehensive evaluation threshold Φref, it means that the hot compress state is normal and no adjustment is required. When the comprehensive evaluation index Φcyc is less than the comprehensive evaluation threshold Φref, it indicates that the hot compress state is abnormal. The thermal response inertia penalty coefficient γ1, the thermal conduction hysteresis penalty coefficient γ2, the recovery sensitivity factor B1, and the thermal coupling sensitivity factor B2 are adjusted to obtain new thermal response inertia penalty coefficient nγ1, new thermal conduction hysteresis penalty coefficient nγ2, new recovery sensitivity factor nB1, and new thermal coupling sensitivity factor nB2. The adjustment formula is as follows: nγ1=γ1+(1+η1*(1-Φcyc)); nγ2=γ2+(1+η1*(1-Φcyc)); nB1=B1+(1+η2*(1-Φcyc)); nB2=B2+(1+η2*(1-Φcyc)); Where η1 represents the recovery factor adjustment step coefficient, and η2 represents the thermal hysteresis factor adjustment step coefficient.

Images