A neurology nursing hidden danger monitoring system and a monitoring method thereof

Abstract

The application relates to the technical field of medical monitoring, and discloses a neurology nursing hidden danger monitoring system and a monitoring method thereof, which comprises the following steps: separating static posture pressure values and tremor fluctuation pressure values through sliding average; extracting original mass center space sliding amounts and thermal drift space false image sliding amounts respectively, and calculating a dynamic contact thermal resistance modulation rate by using tremor energy; extracting a real space sliding amount after difference correction; further, calculating a tremor body energy consumption danger index by performing attenuation integral on the tremor energy, and calculating an edge sliding bed falling danger index by using the real space sliding amount and combining a safety buffer constant; and finally, multiplying the two indexes, combining an adaptive background, and generating a final hidden danger alarm state quantity.

Description

Technical Field

[0001] This invention relates to the field of medical monitoring technology, specifically to a neurological nursing risk monitoring system and its monitoring method. Background Technology

[0002] In neurology wards, patients often exhibit involuntary, high-frequency movements such as pathological tremors or spasms due to their illness. Existing patient fall risk monitoring solutions largely rely on pressure sensor arrays installed on the bed surface, triggering fall warnings by calculating the shift in the center of gravity due to pressure. However, such solutions face several intractable problems in practical applications. First, the patient's continuous tremors generate high-frequency fluctuating pressure, which can easily mask the patient's slow, static sliding tendency, preventing the system from accurately extracting the true posture reference and leading to misjudgments. Second, when the patient's local body temperature is transferred to the sensor through the bed surface, the elastic modulus of the sensor material changes due to thermo-pressure coupling effects, generating thermodynamic drift pressure artifacts unrelated to the actual force. This thermal drift causes a false shift in the calculated center of gravity. Furthermore, existing technologies typically assume a relatively stable contact state between the patient and the bed surface when performing environmental compensation. However, in a neurology setting, the patient's continuous tremors cause drastic changes in the contact area between the body and the bed surface. This dynamic change in air gap renders the thermal compensation mechanism ineffective in static environments, resulting in the inability to effectively correct for thermal artifacts. Finally, prolonged uncontrollable tremors can also lead to the risk of physical exhaustion in patients, and existing systems often overlook this hidden danger, failing to simultaneously assess the dual risks of falling out of bed and physical exhaustion.

[0003] In summary, there is a need for a monitoring scheme that can separate tremor noise, correct for thermal drift artifacts, and comprehensively assess the risks of falls and physical exhaustion in order to improve the safety early warning capabilities of neurological nursing. Summary of the Invention

[0004] This invention provides a monitoring system and method for potential risks in neurological nursing care, which helps to solve the problems mentioned in the background art.

[0005] This invention provides the following technical solution: a method for monitoring potential risks in neurological nursing, comprising:

[0006] The moving average of the original pressure acquisition values ​​within the physiological characteristic time window is calculated to generate the static posture pressure value; the static posture pressure value is calculated by subtracting the static posture pressure value from the original pressure acquisition value.

[0007] The coordinates of the spatial pressure centroid are calculated by weighting the static attitude pressure values; the original centroid spatial slip is calculated based on the difference between the coordinates of the spatial pressure centroids at adjacent time points.

[0008] The thermodynamic drift pressure artifact value is calculated by combining the difference between the local temperature acquisition value and the environmental reference temperature value with the thermo-pressure coupling coefficient; the thermal drift space artifact slip amount is calculated by combining the polar coordinate range-weighted value with the static attitude pressure value.

[0009] The global tremor energy value is generated by summing the squares of the tremor fluctuation pressure values ​​over the entire domain; the dynamic contact thermal resistance modulation rate is calculated using the ratio of the static attitude pressure value to the total pressure including fluctuations; the corrected thermal drift artifact slip is generated by multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate.

[0010] The difference between the original centroid spatial slip and the differential variation between the thermal drift artifact slip at adjacent time points is calculated to extract the true spatial slip.

[0011] Based on an integral mechanism with time decay weighting, the global tremor energy values ​​within the physiological characteristic time window are accumulated to calculate the tremor physical energy consumption risk index.

[0012] The edge slippage and bed fall risk index is calculated by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed and the small constant of spatial safety buffer.

[0013] The tremor physical exertion risk index is multiplied by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. The difference between the comprehensive hazard assessment value and the sum of squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window is used to generate the final hazard alarm state quantity through a mathematical symbol function.

[0014] Optionally, the step of calculating the moving average of the original pressure acquisition values ​​within a physiological characteristic time window to generate a static posture pressure value; and calculating the tremor fluctuation pressure value by subtracting the static posture pressure value from the original pressure acquisition value, includes:

[0015] During the system initialization phase, a grid-like array of piezoelectric pressure sensors and a grid of temperature sensors are laid out on the hospital bed.

[0016] A two-dimensional plane coordinate system parallel to the bed surface is established with the geometric center of the bed surface as the origin. The horizontal coordinate and vertical coordinate of the sensor array corresponding to each sensor in the grid-distributed piezoelectric pressure sensor array and temperature sensor array laid on the bed are set, and the pressure sensor and temperature sensor are guaranteed to correspond one-to-one on the same spatial coordinate node.

[0017] Obtain the raw pressure data directly measured at each node;

[0018] Set the length of the time window for physiological characteristics;

[0019] The original pressure acquisition values ​​at each spatial coordinate node at the current time point and at historical time points within the physiological characteristic time window are summed and then divided by the physiological characteristic time window length to generate the static posture pressure value of each spatial coordinate node at the current time point.

[0020] The tremor pressure value of each spatial coordinate node at the current time point is calculated by subtracting the corresponding static attitude pressure value from the original pressure acquisition value of each spatial coordinate node at the current time point.

[0021] Optionally, the step of calculating the spatial pressure centroid coordinates by weighting the static attitude pressure values ​​and calculating the original centroid spatial slip based on the difference in spatial pressure centroid coordinates at adjacent time points includes:

[0022] The lateral product is obtained by multiplying the lateral coordinate of each spatial coordinate node with its corresponding static attitude pressure value.

[0023] The global horizontal sum is obtained by summing the horizontal products of all spatial coordinate nodes in the entire domain.

[0024] The global pressure sum is obtained by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain.

[0025] Divide the global horizontal sum by the global pressure sum to calculate the horizontal coordinate of the spatial pressure centroid at the current time point;

[0026] The longitudinal product is obtained by multiplying the longitudinal coordinate of each spatial coordinate node by its corresponding static attitude pressure value.

[0027] The global vertical sum is obtained by summing the vertical products of all spatial coordinate nodes in the entire domain.

[0028] Divide the global longitudinal sum by the global pressure sum to calculate the spatial pressure centroid ordinate at the current time point;

[0029] Calculate the difference between the abscissa of the spatial pressure centroid at the current time point and the abscissa of the spatial pressure centroid at the previous time point, and square it to obtain the squared value of the change in the abscissa.

[0030] Calculate the difference between the ordinate of the spatial pressure centroid at the current time point and the ordinate of the spatial pressure centroid at the previous time point, and square it to obtain the squared change value of the ordinate.

[0031] Add the squared change value of the horizontal axis to the squared change value of the vertical axis to obtain the sum of changes;

[0032] The original centroid spatial slip at the current time point is obtained by taking the square root of the sum of the variations.

[0033] Optionally, the step of calculating the thermodynamic drift pressure artifact value based on the difference between the local temperature acquisition value and the environmental reference temperature value, combined with the thermo-pressure coupling coefficient; and calculating the thermal drift spatial artifact slip amount using polar coordinate range-weighted averaging combined with static attitude pressure value, includes:

[0034] The local temperature data of the corresponding nodes is synchronously acquired by a grid-like temperature sensor array pre-installed on the hospital bed.

[0035] Obtain the ambient reference temperature value and the thermo-pressure coupling coefficient;

[0036] Calculate the difference between the local temperature acquisition value and the environmental reference temperature value at each spatial coordinate node;

[0037] Multiplying this difference by the thermo-pressure coupling coefficient yields the illusory thermodynamic drift pressure value for each spatial coordinate node at the current time point;

[0038] Squaring the horizontal coordinate of each spatial coordinate node and squaring the vertical coordinate;

[0039] The square root of the sum of the squares of the horizontal and vertical coordinates is used to calculate the polar coordinate distance.

[0040] The node weighting value is obtained by multiplying the polar coordinate distance of each spatial coordinate node by its corresponding thermodynamic drift pressure illusory value.

[0041] The global weighted sum is obtained by summing the node weights of all spatial coordinate nodes in the entire domain.

[0042] Divide the global weighted sum by the sum of the static attitude pressure values ​​of all spatial coordinate nodes in the global domain to obtain the thermal drift spatial illusion slip at the current time point.

[0043] Optionally, the step of summing the squares of the tremor fluctuation pressure values ​​globally to generate a global tremor energy value; calculating the dynamic contact thermal resistance modulation rate using the ratio of the static attitude pressure value to the total pressure including fluctuations; and multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate to generate a corrected thermal drift artifact slip, including:

[0044] Squaring the vibrational pressure value at each spatial coordinate node;

[0045] The global tremor energy value at the current time point is generated by summing the squares of the tremor fluctuation pressure values ​​at all spatial coordinate nodes across the entire domain.

[0046] The sum of the absolute values ​​of the static attitude pressure value and the corresponding tremor fluctuation pressure value of each spatial coordinate node in the whole domain is accumulated and summed in the whole domain to obtain the total pressure of the whole domain including fluctuation.

[0047] The dynamic contact thermal resistance modulation rate at the current time point is calculated by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain and dividing the sum by the total pressure including fluctuations in the entire domain.

[0048] Multiply the current thermal drift spatial artifact slip by the current dynamic contact thermal resistance modulation rate to generate the corrected thermal drift artifact slip at the current time.

[0049] Optionally, the step of calculating the difference between the original centroid spatial slip and the differential variation of the thermal drift artifact slip corrected at adjacent time points to extract the true spatial slip includes:

[0050] Calculate the difference between the corrected thermal drift artifact slip at the current time point and the corrected thermal drift artifact slip at the previous time point;

[0051] Take the absolute value of this difference as the amount of change in the difference;

[0052] The true spatial slip difference is obtained by subtracting the differential change from the original centroid spatial slip at the current time point.

[0053] The absolute value of the real spatial slip difference is taken to obtain the real spatial slip amount at the current time point.

[0054] Optionally, the calculation of the tremor energy depletion risk index based on the time-decay-weighted integral mechanism involves accumulating the global tremor energy values ​​within a physiological characteristic time window, including:

[0055] In the current time point and the historical time points within the length of the physiological characteristic time window, for each time point, calculate the difference between the integral vertex corresponding to that time point and the current time point;

[0056] Add the physiological characteristic time window length to the difference and divide by the physiological characteristic time window length to obtain the time decay weight;

[0057] The weighted energy value is obtained by multiplying the global tremor energy value at each time point with its corresponding time decay weight;

[0058] The weighted energy values ​​at all time points within the specified physiological characteristic time window are summed to calculate the tremor energy consumption risk index at the current time point.

[0059] Optionally, the step of calculating the edge slippage and fall-from-bed risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed and the small constant of spatial safety buffer includes:

[0060] Set the physical distance limit from the geometric center to the edge of the hospital bed and the small constant of spatial safety buffer;

[0061] Squaring the horizontal axis of the spatial pressure centroid at the current time point and squaring the vertical axis of the spatial pressure centroid at the current time point;

[0062] The square root of the sum of the squares of the abscissa and ordinate of the centroid of spatial pressure is used to calculate the actual distance from the centroid of spatial pressure to the geometric center at the current time point.

[0063] The remaining spatial distance is calculated by subtracting the actual distance from the physical distance limit from the geometric center to the edge of the bed.

[0064] Compare the remaining spatial distance with the magnitude of the spatial safety buffer constant;

[0065] Extract the maximum value as the denominator of the safety boundary distance;

[0066] The edge slippage and bed fall risk index at the current time point is calculated by dividing the actual spatial slippage at the current time point by the denominator of the safety boundary distance.

[0067] Optionally, the step of multiplying the tremor energy consumption hazard index with the edge slippage and bed fall hazard index to generate a comprehensive hazard evaluation value; and using the difference between the comprehensive hazard evaluation value and the sum of squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window, a final hazard alarm state quantity is generated through a mathematical symbol function, including:

[0068] Multiply the tremor physical exertion risk index at the current time point by the edge slippage and bed fall risk index at the current time point to generate the comprehensive hazard assessment value at the current time point;

[0069] The static attitude pressure value of each spatial coordinate node at the end of the first complete physiological characteristic time window is squared;

[0070] The base sum is obtained by summing the squared static attitude pressure values ​​of all spatial coordinate nodes in the entire domain.

[0071] Multiply the sum of the baselines by the length of the physiological characteristic time window to obtain the baseline value;

[0072] The alarm judgment difference is obtained by subtracting the baseline value from the comprehensive hazard assessment value at the current time point;

[0073] The sign direction value of the alarm judgment difference is extracted using a mathematical symbol function;

[0074] Add one to the symbol direction value and then multiply it by one-half to calculate and generate the final hidden danger alarm status value at the current time point;

[0075] When the final hazard alarm status value equals 1, the nursing station will activate the audible and visual alarm.

[0076] Optionally, a system for implementing the neurological nursing risk monitoring method includes:

[0077] Separation module: used to calculate the tremor fluctuation pressure value using the raw pressure acquisition values ​​within the physiological characteristic time window;

[0078] Slip tracking module: used to calculate the coordinates of the spatial pressure centroid using the weighted average of static attitude pressure values, and to calculate the original centroid spatial slip amount based on the coordinate difference;

[0079] The illusory mapping module is used to calculate the illusory value of thermodynamic drift pressure and the illusory slip of thermal drift space based on the local temperature acquisition value, the ambient reference temperature value, and the thermo-pressure coupling coefficient.

[0080] Thermal resistance correction module: used to calculate the global vibration energy value and dynamic contact thermal resistance modulation rate based on the vibration fluctuation pressure value, and generate the corrected thermal drift artifact slip amount;

[0081] The true slip extraction module is used to calculate the difference between the original centroid spatial slip and the differential change between the thermal drift artifact slip at adjacent time points, and extract the true spatial slip.

[0082] Energy Consumption Assessment Module: Used to integrate and accumulate global tremor energy values ​​within a physiological characteristic time window to calculate the tremor energy consumption risk index;

[0083] The bed fall warning module is used to calculate the edge slippage bed fall risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance from the spatial pressure centroid coordinate to the bed distance limit and the small constant of spatial safety buffer.

[0084] The comprehensive alarm module is used to multiply the tremor physical energy consumption hazard index and the edge slippage and bed fall hazard index to obtain a comprehensive hazard evaluation value, and combine it with the static posture pressure value to generate the final hazard alarm status quantity.

[0085] The present invention has the following beneficial effects:

[0086] 1. By calculating the moving average of the raw pressure data collected within a physiological characteristic time window to separate static posture pressure values ​​from tremor fluctuation pressure values, this solution addresses the problem in neurological nursing environments where high-frequency fluctuations caused by pathological nerve tremors mask the slow static slippage trend of the patient's body. In this specific monitoring environment, the patient's continuous tremor and spasm not only interfere with conventional spatial pressure centroid coordinate tracking, but the local body temperature transfer to the sensor also causes a thermo-pressure coupling effect, producing a thermodynamic drift pressure artifact. Furthermore, continuous tremor leads to drastic dynamic changes in the contact area between the patient's body and the bed surface, making conventional static thermal compensation methods ineffective. Therefore, this solution uses the dynamic contact thermal resistance modulation rate calculated from the tremor fluctuation pressure value to generate a corrected thermal drift artifact slippage amount, and calculates the difference between the original centroid spatial slippage amount and the differential change amount of the corrected thermal drift artifact slippage amount at adjacent time points to extract the true spatial slippage amount, thereby eliminating interference caused by thermal drift and dynamic air gap changes. Furthermore, since prolonged tremors themselves can cause harm, this scheme calculates the edge slippage and fall risk index using the actual spatial slippage amount and the small constant of spatial safety buffer. It also performs attenuation-weighted integration on the global tremor energy value to calculate the tremor energy consumption risk index. Finally, the two are multiplied together and combined with baseline data to generate a comprehensive hazard assessment value and a final hazard alarm status. This scheme, by eliminating the dynamic interference of environmental thermal drift and physiological tremors, achieves a comprehensive and objective assessment of the risk of bed slippage and the hazard of tremor energy consumption.

[0087] 2. By establishing a two-dimensional planar coordinate system with the geometric center of the bed surface as the origin, and ensuring that the piezoelectric pressure sensor array and the temperature sensor array correspond one-to-one at the same spatial coordinate nodes, a unified reference benchmark is provided for the spatial registration of subsequent multi-source physical signals. Furthermore, the original pressure acquisition values ​​at each spatial coordinate node at the current time point and at historical time points within the set physiological characteristic time window are accumulated, and the moving average is calculated to generate the static posture pressure value of each spatial coordinate node at the current time point. Then, the static posture pressure value is subtracted from the original pressure acquisition value to obtain the tremor fluctuation pressure value. This process preserves the patient's long-term posture characteristics while removing the pressure disturbance caused by high-frequency pathological nerve tremors, thereby separating the complex superimposed pressure signal into a static component representing slow displacement and a dynamic component representing spastic tremor. This eliminates the masking interference of high-frequency tremor noise on the true posture distribution and provides a structurally clear basic parameter for subsequent assessment of slippage trend and tremor energy consumption.

[0088] 3. By multiplying the horizontal and vertical coordinates of each spatial coordinate node by its corresponding static posture pressure value and summing them across the entire domain, and then dividing by the sum of the static posture pressure values ​​of all spatial coordinate nodes across the entire domain, the horizontal and vertical coordinates of the spatial pressure centroid at the current time point are calculated. This coordinate weighting process uses pure posture pressure data filtered of tremor disturbances for centroid positioning, objectively reflecting the actual physical spatial distribution of the patient on the bed. Subsequently, the difference between the horizontal and vertical coordinates of the spatial pressure centroid at the current time point and the previous time point is calculated, and the original centroid spatial slip is calculated using the Euclidean distance algorithm of square root of the square. This transforms the static spatial position state into a dynamic relative displacement state, quantifies the instantaneous spatial change rate of the patient moving towards the edge of the bed, and enables the monitoring system to capture the slow bed slip trajectory masked by high-frequency tremors, providing kinematic reference for subsequent safety boundary approach determination.

[0089] 4. By acquiring local temperature data and environmental baseline temperature data, the difference between the two is calculated and multiplied by the thermo-pressure coupling coefficient to obtain the illusory thermodynamic drift pressure value at the current time point. This process maps the sensor modulus change caused by the patient's local body temperature conduction into a precise physical force error, clarifying the interference of non-gravity environmental factors on the sensor acquisition signal. Subsequently, the square root of the sum of the squares of the horizontal and vertical coordinates is taken to obtain the polar coordinate distance. The polar coordinate distance is multiplied by the illusory thermodynamic drift pressure value to obtain the node weighted value. The global weighted sum is divided by the global static attitude pressure sum to calculate the thermal drift spatial illusory slip. This process projects the local thermal deformation error to the overall spatial displacement dimension, quantifies the illusory spatial offset introduced by temperature deviation in global centroid tracking, and provides a quantitative correction benchmark for subsequently eliminating environmental temperature interference in the motion trajectory.

[0090] 5. By summing the squared tremor fluctuation pressure values ​​of all spatial coordinate nodes in the global domain to generate the global tremor energy value at the current time point, the severity of pathological tremor spasm in neurological patients was quantitatively assessed. Furthermore, the dynamic contact thermal resistance modulation rate was calculated by dividing the sum of the global static posture pressure values ​​by the sum of the absolute values ​​of the global static posture pressure values ​​and the tremor fluctuation pressure values. This ratio accurately represents the dynamic loss of area due to the continuous separation and contact between the patient's body and the bed surface caused by persistent tremor, reflecting the fluctuation state of the system's contact thermal resistance. Finally, the thermal drift spatial artifact slip at the current time point was multiplied by this dynamic contact thermal resistance modulation rate to generate the corrected thermal drift artifact slip. This step uses the dimensionless area modulation rate to compensate for the basic thermal drift error, eliminating the nonlinear interference of dynamic contact area changes on thermodynamic conduction and solving the concurrent problem of the failure of the conventional static environment thermal compensation mechanism caused by persistent tremor.

[0091] 6. By calculating the difference between the corrected thermal drift illusory slip at the current time point and the corrected thermal drift illusory slip at the previous time point, and taking the absolute value of this difference to extract the differential variation, this process transforms the thermal drift amplitude, which originally represented the absolute spatial position error, into a differential error component representing the relative position change. This ensures that the environmental disturbance variable is consistent with the spatial displacement increment representing the sliding velocity in the kinematic dimension. Subsequently, the above differential variation is subtracted from the original centroid spatial slip at the current time point, and the absolute value of the calculation result is taken to obtain the true spatial slip at the current time point. This operation removes the thermal illusory component modulated by tremor nonlinearity from the original motion trajectory, which is a mixture of physiological displacement and physical thermal deformation, and eliminates the differential mismatch contradiction caused by the direct subtraction of relative displacement and absolute position error, thus extracting the purified true physical spatial displacement vector.

[0092] 7. By calculating the difference between the integral cursor corresponding to each historical time point and the current time point within a set physiological characteristic time window, and adding the time window length to this difference and dividing by the time window length, a time decay weight that gradually decreases over time is constructed. Subsequently, the global tremor energy value corresponding to each time point is multiplied by its corresponding time decay weight to obtain a weighted energy value. All weighted energy values ​​within the time window are summed to calculate the tremor physical exertion risk index at the current time point. This time-constrained integral mechanism, when quantifying the physical function consumed by a patient's prolonged uncontrollable tremor spasm, gives higher weight to recent tremor movements, so that the destructive effect of long-term tremor gradually and smoothly decays. It objectively reflects the potential risk of physical exhaustion caused by the continuous accumulation of tremor and provides a physiological state assessment dimension independent of spatial displacement monitoring.

[0093] 8. The actual distance from the centroid to the geometric center of the bed is obtained by calculating the square root of the sum of the squares of the horizontal and vertical coordinates of the centroid at the current time point. The remaining spatial distance is obtained by subtracting this actual distance from the set physical distance limit from the geometric center of the bed to the edge. Furthermore, the remaining spatial distance is compared with the size of the set spatial safety buffer constant, and the maximum value is extracted as the denominator of the safety boundary distance. The actual spatial slippage is divided by the denominator of the safety boundary distance to calculate the edge slippage and bed fall risk index. This process correlates the patient's actual slippage rate with the physical buffer margin approaching the edge of the bed, and uses the maximum value function to truncate the remaining distance. This avoids the situation where the denominator is zero or negative when the patient's centroid just moves to the physical edge or goes beyond the edge, preventing the physical meaning of the bed fall risk index from being reversed and ensuring the arithmetic stability of the monitoring system in the extreme danger boundary area.

[0094] 9. By multiplying the tremor physical exertion risk index at the current time point with the edge slippage and bed fall risk index to generate a comprehensive hazard evaluation value, the superposition effect of the patient's physical exhaustion and bed slippage concurrent states is measured. At the same time, the static posture pressure values ​​of each spatial coordinate node at the moment when the system fills the first physiological characteristic time window are summed by squares and multiplied by the time window length to generate a dimension-matched baseline value. The alarm judgment difference is obtained by subtracting this baseline value from the comprehensive hazard evaluation value. Finally, the sign direction value is extracted using a mathematical symbol function to generate the final hazard alarm state quantity. This operation strictly anchors the adaptive safety baseline to the moment when the system completes the initial historical cache, eliminating the calculation deviation problem caused by the lack of data in the initial sampling stage, and constructing the judgment boundary by using basic physical parameters to realize the step-like state warning of the comprehensive hazard state. Attached Figure Description

[0095] Figure 1 This is a schematic diagram of the basic process of the present invention.

[0096] Figure 2 This is a schematic diagram of the pressure separation and centroid slip tracking process of the present invention.

[0097] Figure 3 This is a schematic diagram of the thermal drift artifact correction and real slip extraction process of the present invention.

[0098] Figure 4 This is a schematic diagram of the vibration-induced energy consumption and bed fall early warning and integrated alarm process of the present invention. Detailed Implementation

[0099] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0100] Example 1, refer to Figure 1 A method for monitoring potential risks in neurological nursing care, comprising:

[0101] The moving average of the original pressure acquisition values ​​within the physiological characteristic time window is calculated to generate the static posture pressure value; the static posture pressure value is calculated by subtracting the static posture pressure value from the original pressure acquisition value.

[0102] The coordinates of the spatial pressure centroid are calculated by weighting the static attitude pressure values; the original centroid spatial slip is calculated based on the difference between the coordinates of the spatial pressure centroids at adjacent time points.

[0103] The thermodynamic drift pressure artifact value is calculated by combining the difference between the local temperature acquisition value and the environmental reference temperature value with the thermo-pressure coupling coefficient; the thermal drift space artifact slip amount is calculated by combining the polar coordinate range-weighted value with the static attitude pressure value.

[0104] The global tremor energy value is generated by summing the squares of the tremor fluctuation pressure values ​​over the entire domain; the dynamic contact thermal resistance modulation rate is calculated using the ratio of the static attitude pressure value to the total pressure including fluctuations; the corrected thermal drift artifact slip is generated by multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate.

[0105] The difference between the original centroid spatial slip and the differential variation between the thermal drift artifact slip at adjacent time points is calculated to extract the true spatial slip.

[0106] Based on an integral mechanism with time decay weighting, the global tremor energy values ​​within the physiological characteristic time window are accumulated to calculate the tremor physical energy consumption risk index.

[0107] The edge slippage and bed fall risk index is calculated by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed and the small constant of spatial safety buffer.

[0108] The tremor physical exertion risk index is multiplied by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. The difference between the comprehensive hazard assessment value and the sum of squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window is used to generate the final hazard alarm state quantity through a mathematical symbol function.

[0109] The calculation of the moving average of the original pressure acquisition values ​​within a physiological characteristic time window generates a static posture pressure value; the calculation of the tremor fluctuation pressure value is achieved by subtracting the static posture pressure value from the original pressure acquisition value, including:

[0110] During the system initialization phase, a grid-like array of piezoelectric pressure sensors and a grid of temperature sensors are laid out on the hospital bed.

[0111] A two-dimensional plane coordinate system parallel to the bed surface is established with the geometric center of the bed surface as the origin. The horizontal coordinate and vertical coordinate of the sensor array corresponding to each sensor in the grid-distributed piezoelectric pressure sensor array and temperature sensor array laid on the bed are set, and the pressure sensor and temperature sensor are guaranteed to correspond one-to-one on the same spatial coordinate node.

[0112] Obtain the raw pressure data directly measured at each node;

[0113] Set the length of the time window for physiological characteristics;

[0114] The original pressure acquisition values ​​at each spatial coordinate node at the current time point and at historical time points within the physiological characteristic time window are summed and then divided by the physiological characteristic time window length to generate the static posture pressure value of each spatial coordinate node at the current time point.

[0115] The tremor pressure value of each spatial coordinate node at the current time point is calculated by subtracting the corresponding static attitude pressure value from the original pressure acquisition value of each spatial coordinate node at the current time point.

[0116] The calculation of spatial pressure centroid coordinates using coordinate weighting based on static attitude pressure values, and the calculation of the original centroid spatial slip based on the difference in spatial pressure centroid coordinates at adjacent time points, includes:

[0117] The lateral product is obtained by multiplying the lateral coordinate of each spatial coordinate node with its corresponding static attitude pressure value.

[0118] The global horizontal sum is obtained by summing the horizontal products of all spatial coordinate nodes in the entire domain.

[0119] The global pressure sum is obtained by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain.

[0120] Divide the global horizontal sum by the global pressure sum to calculate the horizontal coordinate of the spatial pressure centroid at the current time point;

[0121] The longitudinal product is obtained by multiplying the longitudinal coordinate of each spatial coordinate node by its corresponding static attitude pressure value.

[0122] The global vertical sum is obtained by summing the vertical products of all spatial coordinate nodes in the entire domain.

[0123] Divide the global longitudinal sum by the global pressure sum to calculate the spatial pressure centroid ordinate at the current time point;

[0124] Calculate the difference between the abscissa of the spatial pressure centroid at the current time point and the abscissa of the spatial pressure centroid at the previous time point, and square it to obtain the squared value of the change in the abscissa.

[0125] Calculate the difference between the ordinate of the spatial pressure centroid at the current time point and the ordinate of the spatial pressure centroid at the previous time point, and square it to obtain the squared change value of the ordinate.

[0126] Add the squared change value of the horizontal axis to the squared change value of the vertical axis to obtain the sum of changes;

[0127] The original centroid spatial slip at the current time point is obtained by taking the square root of the sum of the variations.

[0128] The process involves calculating the illusory thermodynamic drift pressure value based on the difference between the local temperature acquisition value and the environmental reference temperature value, combined with the thermo-pressure coupling coefficient; and calculating the illusory thermal drift space slip amount using polar coordinate range-weighted averaging combined with static attitude pressure values, including:

[0129] The local temperature data of the corresponding nodes is synchronously acquired by a grid-like temperature sensor array pre-installed on the hospital bed.

[0130] Obtain the ambient reference temperature value and the thermo-pressure coupling coefficient;

[0131] Calculate the difference between the local temperature acquisition value and the environmental reference temperature value at each spatial coordinate node;

[0132] Multiplying this difference by the thermo-pressure coupling coefficient yields the illusory thermodynamic drift pressure value for each spatial coordinate node at the current time point;

[0133] Squaring the horizontal coordinate of each spatial coordinate node and squaring the vertical coordinate;

[0134] The square root of the sum of the squares of the horizontal and vertical coordinates is used to calculate the polar coordinate distance.

[0135] The node weighting value is obtained by multiplying the polar coordinate distance of each spatial coordinate node by its corresponding thermodynamic drift pressure illusory value.

[0136] The global weighted sum is obtained by summing the node weights of all spatial coordinate nodes in the entire domain.

[0137] Divide the global weighted sum by the sum of the static attitude pressure values ​​of all spatial coordinate nodes in the global domain to obtain the thermal drift spatial illusion slip at the current time point.

[0138] The global tremor energy value is generated by summing the squares of the tremor fluctuation pressure values ​​over the entire domain; the dynamic contact thermal resistance modulation rate is calculated using the ratio of the static attitude pressure value to the total pressure including fluctuations; and the corrected thermal drift artifact slip is generated by multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate, including:

[0139] Squaring the vibrational pressure value at each spatial coordinate node;

[0140] The global tremor energy value at the current time point is generated by summing the squares of the tremor fluctuation pressure values ​​at all spatial coordinate nodes across the entire domain.

[0141] The sum of the absolute values ​​of the static attitude pressure value and the corresponding tremor fluctuation pressure value of each spatial coordinate node in the whole domain is accumulated and summed in the whole domain to obtain the total pressure of the whole domain including fluctuation.

[0142] The dynamic contact thermal resistance modulation rate at the current time point is calculated by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain and dividing the sum by the total pressure including fluctuations in the entire domain.

[0143] Multiply the current thermal drift spatial artifact slip by the current dynamic contact thermal resistance modulation rate to generate the corrected thermal drift artifact slip at the current time.

[0144] The method of calculating the difference between the original centroid spatial slip and the differential variation of the thermal drift artifact slip at adjacent time points to extract the true spatial slip includes:

[0145] Calculate the difference between the corrected thermal drift artifact slip at the current time point and the corrected thermal drift artifact slip at the previous time point;

[0146] Take the absolute value of this difference as the amount of change in the difference;

[0147] The true spatial slip difference is obtained by subtracting the differential change from the original centroid spatial slip at the current time point.

[0148] The absolute value of the real spatial slip difference is taken to obtain the real spatial slip amount at the current time point.

[0149] The integral mechanism with time decay weighting accumulates the global tremor energy values ​​within the physiological characteristic time window to calculate the tremor physical exertion risk index, including:

[0150] In the current time point and the historical time points within the length of the physiological characteristic time window, for each time point, calculate the difference between the integral vertex corresponding to that time point and the current time point;

[0151] Add the physiological characteristic time window length to the difference and divide by the physiological characteristic time window length to obtain the time decay weight;

[0152] The weighted energy value is obtained by multiplying the global tremor energy value at each time point with its corresponding time decay weight;

[0153] The weighted energy values ​​at all time points within the specified physiological characteristic time window are summed to calculate the tremor energy consumption risk index at the current time point.

[0154] The method of calculating the edge slippage and fall-from-bed risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed, and the small constant of the spatial safety buffer, includes:

[0155] Set the physical distance limit from the geometric center to the edge of the hospital bed and the small constant of spatial safety buffer;

[0156] Squaring the horizontal axis of the spatial pressure centroid at the current time point and squaring the vertical axis of the spatial pressure centroid at the current time point;

[0157] The square root of the sum of the squares of the abscissa and ordinate of the centroid of spatial pressure is used to calculate the actual distance from the centroid of spatial pressure to the geometric center at the current time point.

[0158] The remaining spatial distance is calculated by subtracting the actual distance from the physical distance limit from the geometric center to the edge of the bed.

[0159] Compare the remaining spatial distance with the magnitude of the spatial safety buffer constant;

[0160] Extract the maximum value as the denominator of the safety boundary distance;

[0161] The edge slippage and bed fall risk index at the current time point is calculated by dividing the actual spatial slippage at the current time point by the denominator of the safety boundary distance.

[0162] The process involves multiplying the tremor energy depletion risk index by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. Then, using the difference between the comprehensive hazard assessment value and the sum of the squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window, a final hazard alarm state quantity is generated through a mathematical symbolic function, including:

[0163] Multiply the tremor physical exertion risk index at the current time point by the edge slippage and bed fall risk index at the current time point to generate the comprehensive hazard assessment value at the current time point;

[0164] The static attitude pressure value of each spatial coordinate node at the end of the first complete physiological characteristic time window is squared;

[0165] The base sum is obtained by summing the squared static attitude pressure values ​​of all spatial coordinate nodes in the entire domain.

[0166] Multiply the sum of the baselines by the length of the physiological characteristic time window to obtain the baseline value;

[0167] The alarm judgment difference is obtained by subtracting the baseline value from the comprehensive hazard assessment value at the current time point;

[0168] The sign direction value of the alarm judgment difference is extracted using a mathematical symbol function;

[0169] Add one to the symbol direction value and then multiply it by one-half to calculate and generate the final hidden danger alarm status value at the current time point;

[0170] When the final hazard alarm status value equals 1, the nursing station will activate the audible and visual alarm.

[0171] Example 2: A method for monitoring potential risks in neurological nursing care, comprising:

[0172] Reference Figure 2 It illustrates the continuous process of pressure signal separation, centroid localization, and generation of the original centroid spatial slip. The calculation of the sliding mean of the original pressure acquisition values ​​within a physiological characteristic time window generates a static posture pressure value; the calculation of the tremor fluctuation pressure value is performed by subtracting the static posture pressure value from the original pressure acquisition value.

[0173] This step aims to address the problem of tremor noise masking the true static posture. Since neural tremor is a high-frequency disturbance, a stable posture reference is obtained by calculating the moving average within a physiological characteristic time window.

[0174] During the system initialization phase, a grid-like array of piezoelectric pressure sensors and a grid of temperature sensors are laid out on the hospital bed.

[0175] A two-dimensional plane coordinate system parallel to the bed surface is established with the geometric center of the bed surface as the origin. The horizontal coordinate and vertical coordinate of the sensor array corresponding to each sensor in the grid-distributed piezoelectric pressure sensor array and temperature sensor array laid on the bed are set, and the pressure sensor and temperature sensor are guaranteed to correspond one-to-one on the same spatial coordinate node.

[0176] Obtain the raw pressure data directly measured at each node;

[0177] Setting the length of the physiological characteristic time window ,in This represents the number of consecutive sampling points, including the current time point. Its value is based on the product of the system sampling frequency and the typical tremor cycle time in neurology. The larger the value, the smoother the filtering of high-frequency noise, but the more severe the response delay to real bed slippage.

[0178] First, in order to filter high-frequency neural impulses through mean integral to obtain a long-term posture reference, a formula for calculating static posture stress is given:

[0179]

[0180] in, Represents the coordinates in the sensor array At the current point in time The static attitude pressure value, Indicates the length of the time window for physiological characteristics. Represents the coordinates in the sensor array At a point in time The original pressure values ​​obtained;

[0181] Subsequently, in order to isolate purely pathological wave characteristics for downstream energy calculations, a formula for extracting thrill wave pressure is given by subtracting the static baseline from the current acquired values:

[0182]

[0183] in, Represents the coordinates in the sensor array At the current point in time The vibrational pressure value, Represents the coordinates in the sensor array At the current point in time The raw pressure values ​​obtained, Represents the coordinates in the sensor array At the current point in time The static attitude pressure value.

[0184] By establishing a two-dimensional planar coordinate system with the geometric center of the bed surface as the origin, and ensuring that the piezoelectric pressure sensor array and the temperature sensor array correspond one-to-one at the same spatial coordinate nodes, a unified reference benchmark is provided for the spatial registration of subsequent multi-source physical signals. Furthermore, the original pressure acquisition values ​​at each spatial coordinate node at the current time point and at historical time points within the set physiological characteristic time window are accumulated, and the moving average is calculated to generate the static posture pressure value of each spatial coordinate node at the current time point. Then, the static posture pressure value is subtracted from the original pressure acquisition value to obtain the tremor fluctuation pressure value. This process preserves the patient's long-term posture characteristics while removing the pressure disturbance caused by high-frequency pathological nerve tremors, thereby separating the complex superimposed pressure signal into a static component representing slow displacement and a dynamic component representing spastic tremor. This eliminates the masking interference of high-frequency tremor noise on the true posture distribution and provides a structurally clear basic parameter for subsequent assessment of slippage trend and tremor energy consumption.

[0185] The calculation of spatial pressure centroid coordinates using coordinate weighting based on static attitude pressure values, and the calculation of the original centroid spatial slip based on the difference in spatial pressure centroid coordinates at adjacent time points, includes:

[0186] The purpose of this step is to quantitatively track the slow bed slippage motion masked by the elimination of vibrations. By obtaining the static attitude pressure, the absolute physical center of gravity of the entire bed surface is calculated, and the displacement vector is extracted by comparing adjacent time points.

[0187] First, in order to locate the patient's spatial center of gravity in the absence of high-frequency noise interference, pressure-weighted normalization using planar coordinates is required, and the formula for calculating the abscissa of the center of gravity is given:

[0188]

[0189] in, Indicates the current time point The horizontal axis of the centroid of spatial pressure. Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At the current point in time The static attitude pressure value;

[0190] Subsequently, the formula for calculating the centroid's ordinate is given:

[0191]

[0192] in, Indicates the current time point The spatial pressure centroid ordinate, Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At the current point in time The static attitude pressure value;

[0193] Next, in order to extract the instantaneous rate of change of the patient sliding towards the edge of the bed due to muscle weakness, the Euclidean spatial difference of the center of mass between two consecutive time points is calculated, and the formula for calculating the original center of mass spatial slip is given:

[0194]

[0195] in, Indicates the current time point The original centroid space slip, Indicates the current time point The horizontal axis of the centroid of spatial pressure. Indicates the current time point The spatial pressure centroid ordinate, Indicates the previous time point The horizontal axis of the centroid of spatial pressure. Indicates the previous time point The ordinate of the centroid of spatial pressure.

[0196] By multiplying the x and y coordinates of each spatial coordinate node by its corresponding static posture pressure value and summing them across the entire domain, and then dividing by the sum of the static posture pressure values ​​of all spatial coordinate nodes across the entire domain, the x and y coordinates of the spatial pressure centroid at the current time point are calculated. This coordinate weighting process uses pure posture pressure data filtered of tremor disturbances for centroid positioning, objectively reflecting the actual physical spatial distribution of the patient on the bed. Subsequently, the difference between the x and y coordinates of the spatial pressure centroid at the current time point and the previous time point is calculated, and the original centroid spatial slip is calculated using the Euclidean distance algorithm of square root of the square. This transforms the static spatial position state into a dynamic relative displacement state, quantifies the instantaneous spatial change rate of the patient moving towards the edge of the bed, and enables the monitoring system to capture the slow bed slip trajectory masked by high-frequency tremors, providing kinematic reference for subsequent safety boundary approach determination.

[0197] Reference Figure 3 It illustrates the process of temperature acquisition, thermal drift artifact mapping, thermal resistance modulation correction, and extraction of real space slip. The process involves calculating the thermodynamic drift pressure artifact value based on the difference between the local temperature acquisition value and the environmental reference temperature value, combined with the thermo-pressure coupling coefficient; and calculating the thermal drift space artifact slip using polar coordinate range-domain weighted summation combined with static attitude pressure values, including:

[0198] The purpose of this step is to address the problem that the patient's local body temperature changes the elastic modulus of the sensor, thereby creating a false displacement illusion.

[0199] The local temperature data of the corresponding nodes is synchronously acquired by a grid-like temperature sensor array pre-installed on the hospital bed.

[0200] The ambient reference temperature value was obtained using an independent thermometer in the ward. ;

[0201] Obtaining the thermo-compression coupling coefficient Its value is obtained through sensor factory calibration;

[0202] First, to clarify the direct physical deformation error caused by temperature deviation, a formula for calculating the illusory value of thermodynamic drift pressure is given using the thermo-pressure coupling characteristic mapping:

[0203]

[0204] in, Represents the coordinates in the sensor array At the current point in time The illusory value of thermodynamic drift pressure, Indicates the thermo-compression coupling coefficient. Represents the coordinates in the sensor array At the current point in time The acquired local temperature data Indicates the ambient reference temperature value;

[0205] Furthermore, in order to determine the centroid shift illusion caused by this pressure error in global space, a formula for calculating the spatial illusory slip of thermal drift is given by distance weighting in polar coordinates:

[0206]

[0207] in, Indicates the current time point The thermal drift space illusion slip amount, Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At the current point in time The illusory value of thermodynamic drift pressure, Represents the coordinates in the sensor array At the current point in time The static attitude pressure value.

[0208] By acquiring local temperature data and environmental baseline temperature data, calculating the difference between the two and multiplying it by the thermo-pressure coupling coefficient, the illusory thermodynamic drift pressure value at the current time point is obtained. This process maps the sensor modulus change caused by the patient's local body temperature conduction into a precise physical force error, clarifying the interference of non-gravitational environmental factors on the sensor acquisition signal. Subsequently, the square root of the sum of the squares of the horizontal and vertical coordinates is taken to obtain the polar coordinate distance. The polar coordinate distance is multiplied by the illusory thermodynamic drift pressure value to obtain the node weighting value. The global weighted sum is divided by the global static attitude pressure accumulation to calculate the thermal drift spatial illusory slip. This process projects the local thermal deformation error to the overall spatial displacement dimension, quantifies the illusory spatial offset introduced by temperature deviation in global centroid tracking, and provides a quantitative correction benchmark for subsequently eliminating environmental temperature interference in the motion trajectory.

[0209] The global tremor energy value is generated by summing the squares of the tremor fluctuation pressure values ​​over the entire domain; the dynamic contact thermal resistance modulation rate is calculated using the ratio of the static attitude pressure value to the total pressure including fluctuations; and the corrected thermal drift artifact slip is generated by multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate, including:

[0210] This step is specifically designed to address the problem that persistent tremors in neurological patients cause continuous changes in the air gap between the body and the bed surface, leading to the direct failure of thermal drift compensation in a static environment.

[0211] First, to quantitatively assess the severity of air gap separation caused by neurotic tremor, the squares of local fluctuation pressures are summed globally, giving a formula for calculating the global tremor energy value:

[0212]

[0213] in, Indicates the current time point The global tremor energy value, Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At the current point in time The vibrational pressure value;

[0214] Next, to accurately reflect the proportion of dynamic contact area loss caused by vibration, a formula for calculating the dynamic contact thermal resistance modulation rate is given by using the relative ratio of static effective support pressure to the absolute total pressure including fluctuations:

[0215]

[0216] in, Indicates the current time point Dynamic contact thermal resistance modulation rate, Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At the current point in time The static attitude pressure value, Represents the coordinates in the sensor array At the current point in time The vibrational pressure value;

[0217] Finally, in order to eliminate the nonlinear interference of dynamic contact area changes on thermal drift, the modulation rate is applied to the basic thermal drift amount, and the calculation formula for the corrected thermal drift spurious slip is given:

[0218]

[0219] in, Indicates the current time point Corrected thermal drift artifact slip amount, Indicates the current time point The thermal drift space illusion slip amount, Indicates the current time point Dynamic contact thermal resistance modulation rate.

[0220] By summing the squared tremor fluctuation pressure values ​​of all spatial coordinate nodes in the global domain to generate the global tremor energy value at the current time point, the severity of pathological tremor spasm in neurological patients was quantitatively assessed. Furthermore, the dynamic contact thermal resistance modulation rate was calculated by dividing the sum of the global static posture pressure values ​​by the sum of the absolute values ​​of the global static posture pressure values ​​and the tremor fluctuation pressure values. This ratio accurately represents the dynamic loss of area due to the continuous separation and contact between the patient's body and the bed surface caused by persistent tremor, reflecting the fluctuation state of the system's contact thermal resistance. Finally, the thermal drift spatial artifact slip at the current time point was multiplied by this dynamic contact thermal resistance modulation rate to generate the corrected thermal drift artifact slip. This step uses the dimensionless area modulation rate to compensate for the basic thermal drift error, eliminating the nonlinear interference of dynamic contact area changes on thermodynamic conduction and resolving the concurrent problem of the failure of conventional static environmental thermal compensation mechanisms caused by persistent tremor.

[0221] The method of calculating the difference between the original centroid spatial slip and the differential variation of the thermal drift artifact slip at adjacent time points to extract the true spatial slip includes:

[0222] The purpose of this step is to completely eliminate the difference mismatch caused by the direct subtraction between the kinematic relative displacement and the absolute spatial error of thermal drift through the negative feedback loop of the data; by introducing the difference change at the previous time point, the highly purified real spatial displacement is extracted from the mixed error trajectory, providing an objective physical basis for judging whether the patient has fallen out of bed.

[0223] To obtain the unbiased, true pre-fall motion vector, the differential variation of the thermal drift spurious slip amount after correction at adjacent time points is subtracted from the original slip amount, and the formula for calculating the true spatial slip amount is given:

[0224]

[0225] in, Indicates the current time point The actual spatial slip, Indicates the current time point The original centroid space slip, Indicates the current time point Corrected thermal drift artifact slip amount, Indicates a point in time The corrected thermal drift artifact slip amount.

[0226] By calculating the difference between the corrected thermal drift spurious slip at the current time point and the corrected thermal drift spurious slip at the previous time point, and taking the absolute value of this difference to extract the differential variation, this process transforms the thermal drift amplitude, which originally represented the absolute spatial position error, into a differential error component representing the relative position change. This ensures that the environmental disturbance variable remains consistent with the spatial displacement increment representing the sliding velocity in the kinematic dimension. Subsequently, the above differential variation is subtracted from the original centroid spatial slip at the current time point, and the absolute value of the calculation result is taken to obtain the true spatial slip at the current time point. This operation removes the thermal spurious component modulated by tremor nonlinearity from the original motion trajectory, which is a mixture of physiological displacement and physical thermal deformation, and eliminates the differential mismatch contradiction caused by directly subtracting the relative displacement from the absolute position error, thus extracting the purified true physical spatial displacement vector.

[0227] Reference Figure 4 It illustrates the process of tremor energy integration, fall risk index calculation, and final hazard warning state generation. The integration mechanism, based on time decay weighting, accumulates global tremor energy values ​​within a physiological characteristic time window to calculate the tremor physical exertion risk index, including:

[0228] The purpose of this step is to monitor another type of hidden danger unique to neurology, namely the continuous energy consumption abnormalities and nursing risks that may be caused by prolonged uncontrollable tremor spasm.

[0229] To highlight the cumulative destructive effect of persistent tremor over time, a structured attenuation weight based on the relative lag time of window sliding is introduced to gradually reduce the destructiveness of long-term tremor. A formula for calculating the tremor physical exertion risk index is given:

[0230]

[0231] in, Indicates the current time point The risk index of physical exertion from tremors. Indicates the length of the time window for physiological characteristics. Indicates a point in time The global tremor energy value, Represents an integral cursor. Indicates the current time point.

[0232] By calculating the difference between the integral cursor corresponding to each historical time point and the current time point within a set physiological characteristic time window, and adding the time window length to this difference and dividing by the time window length, a time decay weight that gradually decreases over time is constructed. Subsequently, the global tremor energy value corresponding to each time point is multiplied by its corresponding time decay weight to obtain a weighted energy value. All weighted energy values ​​within the time window are summed to calculate the tremor physical exertion risk index at the current time point. This time-constrained integral mechanism, when quantifying the physical function consumed by a patient's prolonged uncontrollable tremor spasm, assigns a higher weight to recent tremor movements, so that the destructive effect of long-term tremor gradually and smoothly decays. It objectively reflects the potential risk of physical exhaustion caused by the continuous accumulation of tremor and provides a physiological state assessment dimension independent of spatial displacement monitoring.

[0233] The method of calculating the edge slippage and fall-from-bed risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed, and the small constant of the spatial safety buffer, includes:

[0234] The purpose of this step is to use a highly accurate real physical displacement vector to determine its approximation to the safety boundary of the bed edge.

[0235] Set the physical distance limit from the geometric center to the edge of the hospital bed. Its value is half the actual width of the medical bed surface;

[0236] Set space safety buffer small constant Its value is based on the sensor's minimum physical resolution, in order to prevent the denominator from collapsing when the centroid approaches the edge;

[0237] To prevent the denominator from becoming negative and reversing the physical meaning of the danger index when the center of mass crosses the physical edge, a formula for calculating the danger index of edge slippage and bed fall is given by introducing a maximum value function to truncate the remaining spatial distance:

[0238]

[0239] in, Indicates the current time point Edge slippage and fall risk index Indicates the current time point The actual spatial slip, This represents the limit of the physical distance from the geometric center of the hospital bed to its edge. Indicates the current time point The horizontal axis of the centroid of spatial pressure. Indicates the current time point The spatial pressure centroid ordinate, It provides a small constant for space safety buffering.

[0240] The actual distance from the centroid to the geometric center of the bed is obtained by calculating the square root of the sum of the squares of the horizontal and vertical coordinates of the centroid at the current time point. The remaining spatial distance is obtained by subtracting this actual distance from the set physical distance limit from the geometric center of the bed to the edge. Furthermore, the remaining spatial distance is compared with the size of the set spatial safety buffer constant, and the maximum value is extracted as the denominator of the safety boundary distance. The actual spatial slippage is divided by the denominator of the safety boundary distance to calculate the edge slippage and bed fall risk index. This process correlates the patient's actual slippage rate with the physical buffer margin approaching the edge of the bed, and uses the maximum value function to truncate the remaining distance. This avoids the situation where the denominator is zero or negative when the patient's centroid just moves to the physical edge or goes beyond the edge, preventing the physical meaning of the bed fall risk index from being reversed and ensuring the arithmetic stability of the monitoring system in the extreme danger boundary area.

[0241] The process involves multiplying the tremor energy depletion risk index by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. Then, using the difference between the comprehensive hazard assessment value and the sum of the squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window, a final hazard alarm state quantity is generated through a mathematical symbolic function, including:

[0242] The purpose of this step is to cross-couple displacement hazard with energy exhaustion hazard to complete the final state machine transition determination.

[0243] First, to measure the amplifying effect of concurrent dual hazards, the hazard indices of the two dimensions are multiplicatively coupled, and a formula for calculating the comprehensive hazard evaluation value is given:

[0244]

[0245] in, Indicates the current time point The comprehensive hazard assessment value, Indicates the current time point The risk index of physical exertion from tremors. Indicates the current time point Edge slippage and fall risk index;

[0246] Finally, the final state, which requires no manual optimization, is output through a step function, and the calculation formula for the final hazard alarm state quantity is given:

[0247]

[0248] in, Indicates the current time point The final hazard alarm status quantity, Indicates the current time point The comprehensive hazard assessment value, Indicates the length of the time window for physiological characteristics. Represents the lateral coordinates of the sensor array. Represents the longitudinal coordinates of the sensor array. Represents the coordinates in the sensor array At a point in time The static attitude pressure value;

[0249] When the final hidden danger alarm status quantity When the value equals 1, the nursing station will activate an audible and visual alarm;

[0250] The final hazard alarm status is used to prompt the nursing station to conduct safety inspections or assist in nursing procedures, and is not used to directly output disease diagnosis conclusions.

[0251] By multiplying the tremor physical exertion risk index at the current time point with the edge slippage and bed fall risk index to generate a comprehensive hazard evaluation value, the superposition effect of the dual concurrent states of patient physical exhaustion and bed slippage is measured. At the same time, the static posture pressure values ​​of each spatial coordinate node at the moment when the system fills the first physiological characteristic time window are summed by squares and multiplied by the time window length to generate a dimension-matched baseline value. The alarm judgment difference is obtained by subtracting this baseline value from the comprehensive hazard evaluation value. Finally, the sign direction value is extracted using a mathematical symbol function to generate the final hazard alarm state quantity. This operation strictly anchors the adaptive safety baseline to the moment when the system completes the initial historical cache, eliminating the calculation deviation problem caused by the lack of data in the initial sampling stage, and constructs the judgment boundary by using basic physical parameters to realize the step-like state warning of the comprehensive hazard state.

[0252] Example 3: A system for implementing the aforementioned neurological nursing risk monitoring method, comprising:

[0253] Separation module: used to calculate the tremor fluctuation pressure value using the raw pressure acquisition values ​​within the physiological characteristic time window;

[0254] Slip tracking module: used to calculate the coordinates of the spatial pressure centroid using the weighted average of static attitude pressure values, and to calculate the original centroid spatial slip amount based on the coordinate difference;

[0255] The illusory mapping module is used to calculate the illusory value of thermodynamic drift pressure and the illusory slip of thermal drift space based on the local temperature acquisition value, the ambient reference temperature value, and the thermo-pressure coupling coefficient.

[0256] Thermal resistance correction module: used to calculate the global vibration energy value and dynamic contact thermal resistance modulation rate based on the vibration fluctuation pressure value, and generate the corrected thermal drift artifact slip amount;

[0257] The true slip extraction module is used to calculate the difference between the original centroid spatial slip and the differential change between the thermal drift artifact slip at adjacent time points, and extract the true spatial slip.

[0258] Energy Consumption Assessment Module: Used to integrate and accumulate global tremor energy values ​​within a physiological characteristic time window to calculate the tremor energy consumption risk index;

[0259] The bed fall warning module is used to calculate the edge slippage bed fall risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance from the spatial pressure centroid coordinate to the bed distance limit and the small constant of spatial safety buffer.

[0260] The comprehensive alarm module is used to multiply the tremor physical energy consumption hazard index and the edge slippage and bed fall hazard index to obtain a comprehensive hazard evaluation value, and combine it with the static posture pressure value to generate the final hazard alarm status quantity.

[0261] It should be noted that, in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article, or apparatus.

[0262] The above description is only a preferred embodiment of the present invention. It should be noted that for those skilled in the art, several improvements and modifications can be made without departing from the technical principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for monitoring potential nursing risks in neurology, characterized in that, include: The moving average of the original pressure acquisition values ​​within the physiological characteristic time window is calculated to generate the static posture pressure value; the static posture pressure value is calculated by subtracting the static posture pressure value from the original pressure acquisition value. The coordinates of the spatial pressure centroid are calculated by weighting the static attitude pressure values; the original centroid spatial slip is calculated based on the difference between the coordinates of the spatial pressure centroids at adjacent time points. The thermodynamic drift pressure artifact value is calculated by combining the difference between the local temperature acquisition value and the environmental reference temperature value with the thermo-pressure coupling coefficient; the thermal drift space artifact slip amount is calculated by combining the polar coordinate range-weighted value with the static attitude pressure value. The global tremor energy value is generated by summing the squares of the tremor fluctuation pressure values ​​over the entire domain. The dynamic contact thermal resistance modulation rate is calculated by using the ratio of the static attitude pressure value to the total pressure including fluctuations; the thermal drift space illusory slip is multiplied by the dynamic contact thermal resistance modulation rate to generate the corrected thermal drift illusory slip. The difference between the original centroid spatial slip and the differential variation between the thermal drift artifact slip at adjacent time points is calculated to extract the true spatial slip. Based on an integral mechanism with time decay weighting, the global tremor energy values ​​within the physiological characteristic time window are accumulated to calculate the tremor physical energy consumption risk index. The edge slippage and bed fall risk index is calculated by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed and the small constant of spatial safety buffer. The tremor physical exertion risk index is multiplied by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. The difference between the comprehensive hazard assessment value and the sum of squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window is used to generate the final hazard alarm state quantity through a mathematical symbol function.

2. The method for monitoring potential risks in neurological nursing care according to claim 1, characterized in that, The calculation of the moving average of the original pressure acquisition values ​​within a physiological characteristic time window generates a static posture pressure value; the calculation of the tremor fluctuation pressure value is achieved by subtracting the static posture pressure value from the original pressure acquisition value, including: During the system initialization phase, a grid-like array of piezoelectric pressure sensors and a grid of temperature sensors are laid out on the hospital bed. A two-dimensional plane coordinate system parallel to the bed surface is established with the geometric center of the bed surface as the origin. The horizontal coordinate and vertical coordinate of the sensor array corresponding to each sensor in the grid-distributed piezoelectric pressure sensor array and temperature sensor array laid on the bed are set, and the pressure sensor and temperature sensor are guaranteed to correspond one-to-one on the same spatial coordinate node. Obtain the raw pressure data directly measured at each node; Set the length of the time window for physiological characteristics; The original pressure acquisition values ​​at each spatial coordinate node at the current time point and at historical time points within the physiological characteristic time window are summed and then divided by the physiological characteristic time window length to generate the static posture pressure value of each spatial coordinate node at the current time point. The tremor pressure value of each spatial coordinate node at the current time point is calculated by subtracting the corresponding static attitude pressure value from the original pressure acquisition value of each spatial coordinate node at the current time point.

3. The method for monitoring potential risks in neurological nursing care according to claim 2, characterized in that, The calculation of spatial pressure centroid coordinates using coordinate weighting based on static attitude pressure values, and the calculation of the original centroid spatial slip based on the difference in spatial pressure centroid coordinates at adjacent time points, includes: The lateral product is obtained by multiplying the lateral coordinate of each spatial coordinate node with its corresponding static attitude pressure value. The global horizontal sum is obtained by summing the horizontal products of all spatial coordinate nodes in the entire domain. The global pressure sum is obtained by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain. Divide the global horizontal sum by the global pressure sum to calculate the horizontal coordinate of the spatial pressure centroid at the current time point; The longitudinal product is obtained by multiplying the longitudinal coordinate of each spatial coordinate node by its corresponding static attitude pressure value. The global vertical sum is obtained by summing the vertical products of all spatial coordinate nodes in the entire domain. Divide the global longitudinal sum by the global pressure sum to calculate the spatial pressure centroid ordinate at the current time point; Calculate the difference between the abscissa of the spatial pressure centroid at the current time point and the abscissa of the spatial pressure centroid at the previous time point, and square it to obtain the squared value of the change in the abscissa. Calculate the difference between the ordinate of the spatial pressure centroid at the current time point and the ordinate of the spatial pressure centroid at the previous time point, and square it to obtain the squared change value of the ordinate. Add the squared change value of the horizontal axis to the squared change value of the vertical axis to obtain the sum of changes; The original centroid spatial slip at the current time point is obtained by taking the square root of the sum of the variations.

4. The method for monitoring potential risks in neurological nursing care according to claim 3, characterized in that, The process involves calculating the illusory thermodynamic drift pressure value based on the difference between the local temperature acquisition value and the environmental reference temperature value, combined with the thermo-pressure coupling coefficient; and calculating the illusory thermal drift space slip amount using polar coordinate range-weighted averaging combined with static attitude pressure values, including: The local temperature data of the corresponding nodes is synchronously acquired by a grid-like temperature sensor array pre-installed on the hospital bed. Obtain the ambient reference temperature value and the thermo-pressure coupling coefficient; Calculate the difference between the local temperature acquisition value and the environmental reference temperature value at each spatial coordinate node; Multiplying this difference by the thermo-pressure coupling coefficient yields the illusory thermodynamic drift pressure value for each spatial coordinate node at the current time point; Squaring the horizontal coordinate of each spatial coordinate node and squaring the vertical coordinate; The square root of the sum of the squares of the horizontal and vertical coordinates is used to calculate the polar coordinate distance. The node weighting value is obtained by multiplying the polar coordinate distance of each spatial coordinate node by its corresponding thermodynamic drift pressure illusory value. The global weighted sum is obtained by summing the node weights of all spatial coordinate nodes in the entire domain. Divide the global weighted sum by the sum of the static attitude pressure values ​​of all spatial coordinate nodes in the global domain to obtain the thermal drift spatial illusion slip at the current time point.

5. The method for monitoring potential risks in neurological nursing care according to claim 4, characterized in that, The squares of the tremor fluctuation pressure values ​​are summed over the entire domain to generate a global tremor energy value; the dynamic contact thermal resistance modulation rate is calculated using the ratio of the static attitude pressure value to the total pressure including fluctuations. Multiplying the thermal drift space artifact slip by the dynamic contact thermal resistance modulation rate generates the corrected thermal drift artifact slip, which includes: Squaring the vibrational pressure value at each spatial coordinate node; The global tremor energy value at the current time point is generated by summing the squares of the tremor fluctuation pressure values ​​at all spatial coordinate nodes across the entire domain. The sum of the absolute values ​​of the static attitude pressure value and the corresponding tremor fluctuation pressure value of each spatial coordinate node in the whole domain is accumulated and summed in the whole domain to obtain the total pressure of the whole domain including fluctuation. The dynamic contact thermal resistance modulation rate at the current time point is calculated by summing the static attitude pressure values ​​of all spatial coordinate nodes in the entire domain and dividing the sum by the total pressure including fluctuations in the entire domain. Multiply the current thermal drift spatial artifact slip by the current dynamic contact thermal resistance modulation rate to generate the corrected thermal drift artifact slip at the current time.

6. The method for monitoring potential risks in neurological nursing care according to claim 5, characterized in that, The method of calculating the difference between the original centroid spatial slip and the differential variation of the thermal drift artifact slip at adjacent time points to extract the true spatial slip includes: Calculate the difference between the corrected thermal drift artifact slip at the current time point and the corrected thermal drift artifact slip at the previous time point; Take the absolute value of this difference as the amount of change in the difference; The true spatial slip difference is obtained by subtracting the differential change from the original centroid spatial slip at the current time point. The absolute value of the real spatial slip difference is taken to obtain the real spatial slip amount at the current time point.

7. The method for monitoring potential risks in neurological nursing care according to claim 6, characterized in that, The integral mechanism with time decay weighting accumulates the global tremor energy values ​​within the physiological characteristic time window to calculate the tremor physical exertion risk index, including: In the current time point and the historical time points within the length of the physiological characteristic time window, for each time point, calculate the difference between the integral vertex corresponding to that time point and the current time point; Add the physiological characteristic time window length to the difference and divide by the physiological characteristic time window length to obtain the time decay weight; The weighted energy value is obtained by multiplying the global tremor energy value at each time point with its corresponding time decay weight; The weighted energy values ​​at all time points within the specified physiological characteristic time window are summed to calculate the tremor energy consumption risk index at the current time point.

8. A method for monitoring potential risks in neurological nursing care according to claim 7, characterized in that, The method of calculating the edge slippage and fall-from-bed risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance between the spatial pressure centroid coordinates and the limit value of the distance to the bed, and the small constant of the spatial safety buffer, includes: Set the physical distance limit from the geometric center to the edge of the hospital bed and the small constant of spatial safety buffer; Squaring the horizontal axis of the spatial pressure centroid at the current time point and squaring the vertical axis of the spatial pressure centroid at the current time point; The square root of the sum of the squares of the abscissa and ordinate of the centroid of spatial pressure is used to calculate the actual distance from the centroid of spatial pressure to the geometric center at the current time point. The remaining spatial distance is calculated by subtracting the actual distance from the physical distance limit from the geometric center to the edge of the bed. Compare the remaining spatial distance with the magnitude of the spatial safety buffer constant; Extract the maximum value as the denominator of the safety boundary distance; The edge slippage and bed fall risk index at the current time point is calculated by dividing the actual spatial slippage at the current time point by the denominator of the safety boundary distance.

9. A method for monitoring potential nursing risks in neurology according to claim 8, characterized in that, The process involves multiplying the tremor energy depletion risk index by the edge slippage and bed fall risk index to generate a comprehensive hazard assessment value. Then, using the difference between the comprehensive hazard assessment value and the sum of the squares of the static posture pressure values ​​at the end of the first complete physiological characteristic time window, a final hazard alarm state quantity is generated through a mathematical symbolic function, including: Multiply the tremor physical exertion risk index at the current time point by the edge slippage and bed fall risk index at the current time point to generate the comprehensive hazard assessment value at the current time point; The static attitude pressure value of each spatial coordinate node at the end of the first complete physiological characteristic time window is squared; The base sum is obtained by summing the squared static attitude pressure values ​​of all spatial coordinate nodes in the entire domain. Multiply the sum of the baselines by the length of the physiological characteristic time window to obtain the baseline value; The alarm judgment difference is obtained by subtracting the baseline value from the comprehensive hazard assessment value at the current time point; The sign direction value of the alarm judgment difference is extracted using a mathematical symbol function; Add one to the symbol direction value and then multiply it by one-half to calculate and generate the final hidden danger alarm status value at the current time point; When the final hazard alarm status value equals 1, the nursing station will activate the audible and visual alarm.

10. A system employing the neurological nursing risk monitoring method of claim 1, characterized in that, include: Separation module: used to calculate the tremor fluctuation pressure value using the raw pressure acquisition values ​​within the physiological characteristic time window; Slip tracking module: used to calculate the coordinates of the spatial pressure centroid using the weighted average of static attitude pressure values, and to calculate the original centroid spatial slip amount based on the coordinate difference; The illusory mapping module is used to calculate the illusory value of thermodynamic drift pressure and the illusory slip of thermal drift space based on the local temperature acquisition value, the ambient reference temperature value, and the thermo-pressure coupling coefficient. Thermal resistance correction module: used to calculate the global vibration energy value and dynamic contact thermal resistance modulation rate based on the vibration fluctuation pressure value, and generate the corrected thermal drift artifact slip amount; The true slip extraction module is used to calculate the difference between the original centroid spatial slip and the differential change between the thermal drift artifact slip at adjacent time points, and extract the true spatial slip. Energy Consumption Assessment Module: Used to integrate and accumulate global tremor energy values ​​within a physiological characteristic time window to calculate the tremor energy consumption risk index; The bed fall warning module is used to calculate the edge slippage bed fall risk index by dividing the actual spatial slippage by the maximum value of the remaining spatial distance from the spatial pressure centroid coordinate to the bed distance limit and the small constant of spatial safety buffer. The comprehensive alarm module is used to multiply the tremor physical energy consumption hazard index and the edge slippage and bed fall hazard index to obtain a comprehensive hazard evaluation value, and combine it with the static posture pressure value to generate the final hazard alarm status quantity.

Images