CN115509115A - Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method - Google Patents

Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method Download PDF

Info

Publication number
CN115509115A
CN115509115A CN202210824562.2A CN202210824562A CN115509115A CN 115509115 A CN115509115 A CN 115509115A CN 202210824562 A CN202210824562 A CN 202210824562A CN 115509115 A CN115509115 A CN 115509115A
Authority
CN
China
Prior art keywords
fuzzy
temperature
tissue
pid
pid control
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210824562.2A
Other languages
Chinese (zh)
Inventor
汤云东
陈鸣
苏航
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fuzhou University
Original Assignee
Fuzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fuzhou University filed Critical Fuzhou University
Priority to CN202210824562.2A priority Critical patent/CN115509115A/en
Publication of CN115509115A publication Critical patent/CN115509115A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N2/00Magnetotherapy
    • A61N2/02Magnetotherapy using magnetic fields produced by coils, including single turn loops or electromagnets

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Temperature (AREA)

Abstract

The invention relates to a fuzzy self-adaptive PID control-based magnetic nanoparticle heat production optimization method, which comprises the following steps of: step S1: constructing a geometric model of the biological tissue; step S2: constructing a Pennes biological heat transfer model based on a geometric model of the biological tissue, and predicting the temperature distribution in the biological tissue; and step S3: controlling the heating power of the alternating magnetic field through a PID algorithm; and step S4: controlling the heating power of the alternating magnetic field based on the PID, judging whether the highest temperature of the area is converged to a set value, if not, setting the parameters of the PID algorithm by using a fuzzy controller, and if so, jumping to the step S5; step S5: and using the set parameters in a PID algorithm and outputting an optimized temperature control curve. The invention realizes the control of the calorific value of the magnetic nanoparticles under the action of the alternating magnetic field through the fuzzy self-adaptive PID, thereby effectively controlling the temperature distribution of the tissue area.

Description

Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method
Technical Field
The invention relates to the technical field of modeling of magnetic nanoparticles in an alternating magnetic field, in particular to a fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method.
Background
The emerging technology is continuously developed along with the development of science and technology, and in recent years, the magnetic nano thermotherapy becomes a potential tissue thermal ablation technology with the advantages of high safety, few side effects, high targeting and the like. Under the action of the alternating magnetic field, the magnetic nano particles can convert magnetic field energy into heat energy, so that the temperature of a target area is increased to 42-46 ℃, and the purpose of local tissue ablation is achieved.
In magnetic nano hyperthermia, it is desirable to maintain the temperature of the target area in the range of 42-46 ℃ so that the temperature in the tissue can kill specific cells without damaging normal tissue. Therefore, the heat generation of the magnetic nanoparticles needs to be controlled and optimized in the magnetic nano hyperthermia, so that the accuracy and stability of the treatment temperature in the magnetic nano hyperthermia are improved.
Disclosure of Invention
In view of the above, the present invention aims to provide a method for optimizing heat generation of magnetic nanoparticles based on fuzzy adaptive PID control, which realizes controlling the heat generation value of magnetic nanoparticles under the action of an alternating magnetic field by fuzzy adaptive PID, and further effectively controls the temperature distribution of a tissue region.
In order to achieve the purpose, the invention adopts the following technical scheme:
a magnetic nanoparticle heat generation optimization method based on fuzzy self-adaptive PID control comprises the following steps:
step S1: constructing a geometric model of the biological tissue;
step S2: constructing a Pennes biological heat transfer model based on a geometric model of the biological tissue, and predicting the temperature distribution in the biological tissue;
and step S3: controlling the heating power of the alternating magnetic field through a PID algorithm;
and step S4: controlling the heating power of the alternating magnetic field based on the PID, judging whether the highest temperature of the area is converged to a set value, if not, setting the parameters of the PID algorithm by using a fuzzy controller, and if so, jumping to the step S5;
step S5: and using the set parameters in a PID algorithm and outputting an optimized temperature control curve.
Further, the geometric model of the biological tissue is composed of a circle and an ellipse, wherein the radius is R 1 The circle of (a) represents a first tissue region, and the ellipse having a major axis a and a minor axis b represents a second tissue region, wherein the first tissue region is contained in the second tissue region.
Further, the Pennes biological heat transfer model is as follows:
Figure RE-GDA0003952578630000021
wherein ρ, c, T and k represent the density, specific heat capacity, absolute temperature and heat conductivity coefficient of the tissue, ω b 、ρ b 、c b 、T b Respectively representing blood perfusion rate, blood density, blood specific heat capacity and blood temperature, t represents heating time of the magnetic nanoparticles under the action of an alternating magnetic field, and symbols
Figure RE-GDA0003952578630000022
Representing Hamiltonian, Q m Represents the metabolic heat per unit volume, alpha represents the correction factor for power dissipation, and P represents power dissipation.
Further, the predicting the temperature distribution in the biological tissue specifically includes: setting attribute parameters of the tissue region, and solving the Pennes biological heat transfer model by adopting a finite element method, wherein the tissue region parameters comprise the density, specific heat capacity and heat conductivity coefficient of the tissue.
Further, the control equation of the PID algorithm is:
error(t)=y d (t)-y(t)
Figure RE-GDA0003952578630000031
wherein error (t) represents the set temperature y d The difference between (t) and the actual output temperature y (t), u (t) represents the PID control equation, k p Denotes the proportionality coefficient, k i Denotes the integration time constant, k d Representing the differential time constant.
Further, the step S4 specifically includes the following steps:
step S41: initial parameter k for a given PID control p 、k i 、k d
Step S42: derivation of error and error
Figure RE-GDA0003952578630000032
Step S43: determining input membership functions, derivatives of error and error
Figure RE-GDA0003952578630000033
Fuzzification is carried out;
step S44: determining fuzzy rule for parameter k p 、k i 、k d Fuzzy setting is carried out;
step S45: determining an output membership function, and deblurring the output fuzzy quantity to obtain an accurate quantity;
step S46: will k p 、k i 、k d The initial value of the PID controller and the accurate value obtained in the step S45 are linearly combined to form a new controlled variable of the PID equation, the new controlled variable is brought back to the PID control equation, and whether the temperature is converged to a set value or not is judged; if the temperature does not converge to the set point, the above steps need to be repeated until the temperature converges to the set point.
Further, the fuzzification is to describe the error and the magnitude of the derivative of the error through a language form, and the fuzzy subset is as follows: { NB, NM, NS, ZO, PS, PM, PB }, representing { big negative, middle negative, small negative, zero, small positive, middle positive, big positive }, respectively.
Further, the defuzzification adopts a gravity center method to solve fuzzy output, and the expression is as follows:
Figure RE-GDA0003952578630000041
wherein,z 0 representing the exact value, mu, of the output of the fuzzy controller after the output is deblurred c (z i ) Denotes z i Of membership value of z i Are values in the fuzzy control theoretic domain.
Compared with the prior art, the invention has the following beneficial effects:
the invention realizes the control of the calorific value of the magnetic nanoparticles under the action of the alternating magnetic field through the fuzzy self-adaptive PID, and further effectively controls the temperature distribution of the tissue area, thereby effectively improving the effect of tissue thermal ablation and improving the control precision.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of a geometric model constructed in an embodiment of the present invention
FIG. 3 is a schematic diagram of temperature distribution prediction without fuzzy adaptive PID control according to an embodiment of the present invention
FIG. 4 is a temperature profile without fuzzy adaptive PID control in accordance with an embodiment of the invention;
FIG. 5 is a temperature profile of fuzzy adaptive PID control in one embodiment of the invention.
Detailed Description
The invention is further explained below with reference to the drawings and the embodiments.
Referring to fig. 1, the invention provides a fuzzy adaptive PID control-based magnetic nanoparticle heat generation optimization method, which includes the following steps:
step S1: constructing a geometric model of the biological tissue;
step S2: constructing a Pennes biological heat transfer model based on a geometric model of the biological tissue, and predicting the temperature distribution in the biological tissue;
and step S3: controlling the heating power of the alternating magnetic field through a PID algorithm;
and step S4: controlling the heating power of the alternating magnetic field based on the PID, judging whether the highest temperature of the area is converged to a set value, if not, setting the parameters of the PID algorithm by using a fuzzy controller, and if so, jumping to the step S5;
step S5: and using the set parameters in a PID algorithm and outputting an optimized temperature control curve.
Preferably, in this embodiment, the geometric model described in step S1 is shown in fig. 2, and the model is composed of a circle and an ellipse, wherein the radius is R 1 A circle of =5mm represents a first tissue area, and an ellipse having a major axis of a =20mm and a minor axis of b =15mm represents a second tissue area, wherein the first tissue area is included in the second tissue area.
Preferably, in this embodiment, in step S2, a Pennes biological heat transfer model is constructed to predict the temperature distribution in the biological tissue, where the Pennes biological heat transfer model is:
Figure RE-GDA0003952578630000051
wherein ρ, c, T and k represent the density, specific heat capacity, absolute temperature and heat conductivity coefficient of the tissue, ω b 、ρ b 、c b 、T b Respectively representing blood perfusion rate, blood density, blood specific heat capacity and blood temperature, t represents heating time of the magnetic nanoparticles under the action of an alternating magnetic field, and symbols
Figure RE-GDA0003952578630000052
Representing Hamiltonian, Q m Represents the metabolic heat per unit volume, alpha represents the correction factor for power dissipation, and P represents power dissipation.
Preferably, in this embodiment, the setting of the property parameters of the tissue region and the solving of the biological heat transfer equation by using the finite element method are performed, and the temperature distribution in the tissue is as shown in fig. 3, wherein the tissue region parameters include the density, specific heat capacity and thermal conductivity of the tissue. The parameter data corresponding to the first tissue region are: 1060Kg · m -3 、3540J·Kg -1 ·K -1 、0.52W·m -1 · K -1 The parameter data corresponding to the second tissue region is: 1064Kg m -3 、4500J·Kg -1 · K -1 、0.59W·m -1 ·K -1
Preferably, in this embodiment, in the step S3, the alternating magnetic field heating power is controlled by a PID algorithm, and a control equation of the PID algorithm is error (t) = y d (t)-y(t)
Figure RE-GDA0003952578630000061
Wherein error (t) represents the set temperature y d The difference between (t) and the actual output temperature y (t), u (t) represents the PID control equation, k p Denotes the proportionality coefficient, k i Denotes the integration time constant, k d Representing the differential time constant.
Preferably, in this embodiment, step S41: initial parameter k for a given PID control p = 1、k i =0.05、k d =0.05;
Step S42: derivation of error and error
Figure RE-GDA0003952578630000062
Step S43: determining input membership functions, derivatives of error and error
Figure RE-GDA0003952578630000063
Fuzzification is carried out;
step S44: determining fuzzy rule for parameter k p 、k i 、k d Fuzzy setting is carried out;
step S45: determining an output membership function, and deblurring the output fuzzy quantity to obtain an accurate quantity;
step S46: will k p 、k i 、k d The initial value of the PID controller and the accurate value obtained in step S45 are linearly combined to form a new controlled variable of the PID equation, and the new controlled variable is brought back to the PID control equation to judge whether the temperature is converged to the set value. If the temperature does not converge to 46 deg.C, the above steps need to be repeated until the temperature converges to the set point. FIG. 4 is a temperature profile with center point not being fuzzy adaptive PID controlled, FIG. 5 is a center pointPoint-fuzzy adaptive PID controlled temperature profile.
Preferably, in this embodiment, the step S4 of blurring is to describe the error and the magnitude of the derivative of the error by a language, and the subset of blurring is: { NB, NM, NS, ZO, PS, PM, PB }, representing { big negative, middle negative, small negative, zero, small positive, middle positive, big positive }, respectively.
Preferably, in this embodiment, the defuzzification in step S4 uses a gravity center method to solve the fuzzy output, and the expression is:
Figure RE-GDA0003952578630000071
wherein z is 0 Representing the exact value, mu, of the output of the fuzzy controller after the output is deblurred c (z i ) Denotes z i Of membership value of z i Are values in the fuzzy control theoretic domain.
Preferably, in this embodiment, in the step S5, the setting k in the step S4 is adjusted p 、 k i 、k d And the temperature control curve is used in a PID algorithm and the optimized temperature control curve is output.
The above description is only a preferred embodiment of the present invention, and all equivalent changes and modifications made in accordance with the claims of the present invention should be covered by the present invention.

Claims (8)

1. A magnetic nanoparticle heat generation optimization method based on fuzzy self-adaptive PID control is characterized by comprising the following steps:
step S1: constructing a geometric model of the biological tissue;
step S2: constructing a Pennes biological heat transfer model based on a geometric model of the biological tissue, and predicting the temperature distribution in the biological tissue;
and step S3: controlling the heating power of the alternating magnetic field through a PID algorithm;
and step S4: controlling the heating power of the alternating magnetic field based on the PID, judging whether the highest temperature of the area is converged to a set value, if not, setting the parameters of the PID algorithm by using a fuzzy controller, and if so, jumping to the step S5;
step S5: and using the set parameters in a PID algorithm and outputting an optimized temperature control curve.
2. The method for optimizing the generation of heat of magnetic nanoparticles based on fuzzy adaptive PID control as claimed in claim 1, wherein the geometric model of the biological tissue is composed of a circle and an ellipse, wherein the radius is R 1 The circle of (b) represents a first tissue area, and the ellipse with major axis a and minor axis b represents a second tissue area, wherein the first tissue area is contained in the second tissue area.
3. The magnetic nanoparticle heat generation optimization method based on the fuzzy adaptive PID control according to claim 1, wherein the Pennes biological heat transfer model is:
Figure FDA0003745931220000011
wherein ρ, c, T and k represent the density, specific heat capacity, absolute temperature and heat conductivity coefficient of the tissue, ω b 、ρ b 、c b 、T b Respectively representing blood perfusion rate, blood density, blood specific heat capacity and blood temperature, t represents heating time of the magnetic nanoparticles under the action of an alternating magnetic field, and symbol
Figure FDA0003745931220000024
Representing Hamiltonian, Q m Represents the metabolic heat per unit volume, alpha represents the correction factor for power dissipation, and P represents power dissipation.
4. The fuzzy adaptive PID control based magnetic nanoparticle heat generation optimization method according to claim 3, wherein the predicting of temperature distribution in biological tissue is specifically: setting attribute parameters of the tissue region, and solving the Pennes biological heat transfer model by adopting a finite element method, wherein the tissue region parameters comprise the density, specific heat capacity and heat conductivity coefficient of the tissue.
5. The fuzzy adaptive PID control based magnetic nanoparticle heat generation optimization method according to claim 1, wherein the control equation of the PID algorithm is:
error(t)=y d (t)-y(t)
Figure FDA0003745931220000021
wherein error (t) represents the set temperature y d The difference between (t) and the actual output temperature y (t), u (t) represents the PID control equation, k p Denotes the proportionality coefficient, k i Denotes the integration time constant, k d Representing the differential time constant.
6. The fuzzy adaptive PID control based magnetic nanoparticle heat generation optimization method according to claim 1, wherein the step S4 specifically comprises the following steps:
step S41: initial parameter k for a given PID control p 、k i 、k d
Step S42: derivation of error and error
Figure FDA0003745931220000022
Step S43: determining input membership functions, derivatives of error and error
Figure FDA0003745931220000023
Fuzzification is carried out;
step S44: determining fuzzy rule for parameter k p 、k i 、k d Fuzzy setting is carried out;
step S45: determining an output membership function, and deblurring the output fuzzy quantity to obtain an accurate quantity;
step S46: will k p 、k i 、k d The initial value of the PID controller and the accurate value obtained in the step S45 are linearly combined to form a new controlled variable of the PID equation, the new controlled variable is brought back to the PID control equation, and whether the temperature is converged to a set value or not is judged; if the temperature does not converge to the set point, the above steps need to be repeated until the temperature converges to the set point.
7. The fuzzy adaptive PID control based magnetic nanoparticle heat generation optimization method of claim 6, wherein the fuzzification is to describe the error and the derivative of the error in terms of language, and the fuzzy subset is: { NB, NM, NS, ZO, PS, PM, PB }, representing { big negative, middle negative, small negative, zero, small positive, middle positive, big positive }, respectively.
8. The fuzzy adaptive PID control based magnetic nanoparticle heat generation optimization method of claim 6, wherein the defuzzification uses a gravity center method to solve fuzzy output, and the expression is as follows:
Figure FDA0003745931220000031
wherein z is 0 Representing the exact value, mu, of the output of the fuzzy controller after the output has been deblurred c (z i ) Denotes z i Degree of membership of z i Are values in the fuzzy control theoretic domain.
CN202210824562.2A 2022-07-14 2022-07-14 Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method Pending CN115509115A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210824562.2A CN115509115A (en) 2022-07-14 2022-07-14 Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210824562.2A CN115509115A (en) 2022-07-14 2022-07-14 Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method

Publications (1)

Publication Number Publication Date
CN115509115A true CN115509115A (en) 2022-12-23

Family

ID=84502682

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210824562.2A Pending CN115509115A (en) 2022-07-14 2022-07-14 Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method

Country Status (1)

Country Link
CN (1) CN115509115A (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100204840A1 (en) * 2007-05-25 2010-08-12 Carrier Corporation Modified Fuzzy Control for Chiller Electronic Expansion Valve
CN105629720A (en) * 2015-12-30 2016-06-01 常州联德电子有限公司 Wide-band oxygen sensor temperature control method
CN108775596A (en) * 2018-05-25 2018-11-09 中国科学院深圳先进技术研究院 A kind of heat energy recycle device
CN111258213A (en) * 2020-03-09 2020-06-09 深圳市锐同技术有限公司 Fuzzy self-tuning PID-based temperature control method
CN112947091A (en) * 2021-03-26 2021-06-11 福州大学 PID control-based method for optimizing heat production of magnetic nanoparticles in biological tissues

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100204840A1 (en) * 2007-05-25 2010-08-12 Carrier Corporation Modified Fuzzy Control for Chiller Electronic Expansion Valve
CN105629720A (en) * 2015-12-30 2016-06-01 常州联德电子有限公司 Wide-band oxygen sensor temperature control method
CN108775596A (en) * 2018-05-25 2018-11-09 中国科学院深圳先进技术研究院 A kind of heat energy recycle device
CN111258213A (en) * 2020-03-09 2020-06-09 深圳市锐同技术有限公司 Fuzzy self-tuning PID-based temperature control method
CN112947091A (en) * 2021-03-26 2021-06-11 福州大学 PID control-based method for optimizing heat production of magnetic nanoparticles in biological tissues

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈岩;杜晓明;: "模糊PID控制在温室环境中的应用", 农机化研究, no. 08, 1 August 2010 (2010-08-01), pages 179 - 183 *

Similar Documents

Publication Publication Date Title
CN108134114B (en) Proton exchange membrane fuel cell temperature control method
WO2021196997A1 (en) Method for predicting and controlling water level of series of water conveyance canals on basis of fuzzy neural network
Jiang et al. Design of an intelligent temperature control system based on the fuzzy self-tuning PID
CN110888317A (en) PID controller parameter intelligent optimization method
Salahshour et al. Designing multi-layer quantum neural network controller for chaos control of rod-type plasma torch system using improved particle swarm optimization
CN113552797A (en) Heating furnace temperature control method and system based on improved particle swarm optimization
KR101388069B1 (en) Method for controlling constant pressure using fuzzy pid controller
Murshed et al. Estimation and control of solid oxide fuel cell system
CN112035949B (en) Real-time fuzzy energy management method combining Q reinforcement learning
CN102778844A (en) Induction heating closed loop simulation method based on finite element model and system identification
Tavoosi An experimental study on inverse adaptive neural fuzzy control for nonlinear systems
Xiao Parameter tuning of PID controller for beer filling machine liquid level control based on improved genetic algorithm
CN115509115A (en) Fuzzy self-adaptive PID control-based magnetic nanoparticle heat generation optimization method
Castillo et al. Hierarchical Genetic Algorithms for topology optimization in fuzzy control systems
CN107045290A (en) Reaction-regeneration system optimal control method based on MQPSO DMPC
Whitelam How to train your demon to do fast information erasure without heat production
CN113886990A (en) Reheating fuzzy PID temperature adjusting method for steel billet produced by strip steel hot continuous rolling
Rüttgers et al. Machine-learning-based control of perturbed and heated channel flows
CN117117992B (en) Output power adjusting method, device, computer equipment and storage medium
Zhang et al. An efficient optimization algorithm for extreme value of nonlinear function based on the SAGA and BP algorithm
Li et al. Optimal trajectory exploration large-scale deep reinforcement learning tuned optimal controller for proton exchange membrane fuel cell
Wang et al. Using different Heuristic strategies and an adaptive Neuro-Fuzzy inference system for multi-objective optimization of Hybrid Nanofluid to provide an efficient thermal behavior
CN107461977A (en) A kind of intelligent temperature control method of semiconductor refrigeration temperature control case
Wu et al. Stable fuzzy neural tracking control of a class of unknown nonlinear systems based on fuzzy hierarchy error approach
Cizmeci et al. Improved electro search algorithm with intelligent controller control system: ESPID algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination