CN109634101B - Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal - Google Patents
Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal Download PDFInfo
- Publication number
- CN109634101B CN109634101B CN201910017944.2A CN201910017944A CN109634101B CN 109634101 B CN109634101 B CN 109634101B CN 201910017944 A CN201910017944 A CN 201910017944A CN 109634101 B CN109634101 B CN 109634101B
- Authority
- CN
- China
- Prior art keywords
- photosynthesis
- fluorescence
- chlorophyll fluorescence
- mathematical model
- value
- 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.)
- Active
Links
- 230000029553 photosynthesis Effects 0.000 title claims abstract description 52
- 238000010672 photosynthesis Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 45
- 229930002875 chlorophyll Natural products 0.000 title claims abstract description 33
- 235000019804 chlorophyll Nutrition 0.000 title claims abstract description 33
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 title claims abstract description 33
- 238000013178 mathematical model Methods 0.000 claims abstract description 20
- 238000006243 chemical reaction Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000005516 engineering process Methods 0.000 claims description 16
- 230000009467 reduction Effects 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 238000005286 illumination Methods 0.000 claims description 7
- 230000001105 regulatory effect Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 229910052799 carbon Inorganic materials 0.000 claims description 5
- 230000000875 corresponding effect Effects 0.000 claims description 5
- 230000027756 respiratory electron transport chain Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 5
- 150000001722 carbon compounds Chemical class 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 claims description 4
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000003711 photoprotective effect Effects 0.000 claims description 4
- 235000000346 sugar Nutrition 0.000 claims description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 3
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 239000008103 glucose Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- LJQLQCAXBUHEAZ-UWTATZPHSA-N 3-phospho-D-glyceroyl dihydrogen phosphate Chemical compound OP(=O)(O)OC[C@@H](O)C(=O)OP(O)(O)=O LJQLQCAXBUHEAZ-UWTATZPHSA-N 0.000 claims description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 claims description 2
- ACFIXJIJDZMPPO-NNYOXOHSSA-N NADPH Chemical compound C1=CCC(C(=O)N)=CN1[C@H]1[C@H](O)[C@H](O)[C@@H](COP(O)(=O)OP(O)(=O)OC[C@@H]2[C@H]([C@@H](OP(O)(O)=O)[C@@H](O2)N2C3=NC=NC(N)=C3N=C2)O)O1 ACFIXJIJDZMPPO-NNYOXOHSSA-N 0.000 claims description 2
- 238000005842 biochemical reaction Methods 0.000 claims description 2
- 238000004177 carbon cycle Methods 0.000 claims description 2
- 210000003763 chloroplast Anatomy 0.000 claims description 2
- 230000002596 correlated effect Effects 0.000 claims description 2
- 230000031700 light absorption Effects 0.000 claims description 2
- 229930027945 nicotinamide-adenine dinucleotide Natural products 0.000 claims description 2
- 238000006552 photochemical reaction Methods 0.000 claims description 2
- OSJPPGNTCRNQQC-REOHCLBHSA-N 3-phosphoglyceric acid Chemical compound OC(=O)[C@@H](O)COP(O)(O)=O OSJPPGNTCRNQQC-REOHCLBHSA-N 0.000 claims 1
- 229910019142 PO4 Inorganic materials 0.000 claims 1
- 230000009471 action Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 claims 1
- 239000010452 phosphate Substances 0.000 claims 1
- 238000005070 sampling Methods 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 230000008635 plant growth Effects 0.000 abstract description 10
- 238000011217 control strategy Methods 0.000 abstract description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 3
- 108010060806 Photosystem II Protein Complex Proteins 0.000 abstract description 2
- 241000196324 Embryophyta Species 0.000 description 19
- FKUYMLZIRPABFK-UHFFFAOYSA-N Plastoquinone 9 Natural products CC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCCC(C)=CCC1=CC(=O)C(C)=C(C)C1=O FKUYMLZIRPABFK-UHFFFAOYSA-N 0.000 description 11
- FKUYMLZIRPABFK-IQSNHBBHSA-N plastoquinone-9 Chemical compound CC(C)=CCC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC\C(C)=C\CC1=CC(=O)C(C)=C(C)C1=O FKUYMLZIRPABFK-IQSNHBBHSA-N 0.000 description 11
- 230000000243 photosynthetic effect Effects 0.000 description 7
- YAHZABJORDUQGO-NQXXGFSBSA-N D-ribulose 1,5-bisphosphate Chemical compound OP(=O)(O)OC[C@@H](O)[C@@H](O)C(=O)COP(O)(O)=O YAHZABJORDUQGO-NQXXGFSBSA-N 0.000 description 6
- 150000001720 carbohydrates Chemical class 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 235000014633 carbohydrates Nutrition 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 230000004224 protection Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 150000008163 sugars Chemical class 0.000 description 2
- OSJPPGNTCRNQQC-UWTATZPHSA-N 3-phospho-D-glyceric acid Chemical compound OC(=O)[C@H](O)COP(O)(O)=O OSJPPGNTCRNQQC-UWTATZPHSA-N 0.000 description 1
- 239000002028 Biomass Substances 0.000 description 1
- LXJXRIRHZLFYRP-VKHMYHEASA-N D-glyceraldehyde 3-phosphate Chemical compound O=C[C@H](O)COP(O)(O)=O LXJXRIRHZLFYRP-VKHMYHEASA-N 0.000 description 1
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 1
- 102000004020 Oxygenases Human genes 0.000 description 1
- 108090000417 Oxygenases Proteins 0.000 description 1
- 108010003581 Ribulose-bisphosphate carboxylase Proteins 0.000 description 1
- 241001122767 Theaceae Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000001704 evaporation Methods 0.000 description 1
- 230000008020 evaporation Effects 0.000 description 1
- 238000001506 fluorescence spectroscopy Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000023266 generation of precursor metabolites and energy Effects 0.000 description 1
- 150000002402 hexoses Chemical class 0.000 description 1
- 230000037361 pathway Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000037039 plant physiology Effects 0.000 description 1
- 230000004405 regulation of photosynthesis Effects 0.000 description 1
- 230000008844 regulatory mechanism Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 229910052708 sodium Inorganic materials 0.000 description 1
- 239000011734 sodium Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic 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.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Cultivation Of Plants (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The invention discloses a photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal. The invention relates to a photosynthesis PID closed-loop control method using chlorophyll fluorescence as a feedback signal. It is important to control plant growth in view of limited land and water resources. The current greenhouse control strategy mainly depends on the experience of farmers, does not consider the actual needs of photosynthesis of plants, and therefore cannot achieve the best crop yield and the most effective and reasonable utilization of resources. In this invention, the main process of photosynthesis is represented by an overall mathematical model with chlorophyll fluorescence (ChlF) as output from measurable photosystem ii (psii). The parameters of the model, i.e. the rate constants of the reactions, were further determined by example fitting of experimental data.
Description
Technical Field
The invention relates to the field of plant growth control, in particular to a photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal.
Background
Plants can use solar energy to provide energy and food for organisms. Controlling plant growth is a trend in the future of botany and automatic control, considering limited resources and land. This is one of the 125 global problems promulgated by the "science" journal. Greenhouse control strategies based on plant photosynthesis can reduce waste of resources and increase crop yield. Therefore, it is of great importance to develop plant growth control strategies based on the photosynthetic requirements of plants. According to modern control theory, models covering the main photosynthetic response are crucial for the development of plant growth control strategies. Therefore, it is necessary to design a feasible photosynthesis kinetic model system for achieving the maximized photosynthetic efficiency.
There are three energy pathways for the light energy absorbed by the photosystem ii (psii) of plants, which are involved in photosynthetic chemical reactions, lost as heat energy and dissipated as chlorophyll fluorescence (ChlF). The three pathways are in a trade-off and mutual constraint relation, so that ChlF has very rich photosynthetic information. The process of ChlF loss is complex, but it provides information about the structure and physiology of plant photosynthesis and can be measured by portable instruments. The measurement of ChlF is not invasive to plants and has therefore been widely used as a probe to study photosynthetic energy transfer, physiological processes and intrinsic mechanisms of photosynthesis.
The traditional technology has the following technical problems:
there is a relative lack in the literature of control strategies based on plant physiology to modulate photosynthesis. Although chlorophyll fluorescence can be easily measured, the objective of controlling plant growth is not to obtain the desired fluorescence profile, but to modulate photosynthesis activity, e.g., increase carbohydrate production, etc., by using ChlF as a feedback signal for different purposes. However, a real problem is that in a short time (on the order of seconds) of real-time control, the state quantity of the system, such as the amount of sugars, cannot be measured, and thus the sugars themselves cannot be used as a feedback signal for controlling photosynthesis.
Disclosure of Invention
The invention aims to provide a photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal, and particularly relates to acquisition of chlorophyll fluorescence experimental data, establishment of a photosynthesis model, fitting of the experimental data and a mathematical model, estimation of state quantity by using an extended Kalman filtering technology and regulation of photosynthesis by using a PID control technology.
In order to solve the above technical problems, the present invention provides a photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal, comprising:
acquiring a fluorescence curve;
establishing a mathematical model of the whole process of photosynthesis;
fitting a mathematical model and a fluorescence curve obtained by an experiment, and determining a series of reaction rate constants;
estimating a state variable according to the fluorescence value output by the model and by utilizing an extended Kalman filtering technology;
the photosynthetic activity is regulated using PID control techniques.
In one embodiment, "a fluorescence curve is acquired; "in (1), the chlorophyll fluorometer test is used to obtain the fluorescence curve.
In one embodiment, "a mathematical model of the overall process of photosynthesis is established; "in, establish the photosynthesis mathematical model according to the light reaction electron transfer mechanism, regulation mechanism and dark reaction process of photosynthesis.
In one embodiment, a mathematical model and experimentally derived fluorescence curve are "fitted to determine a series of reaction rate constants; and fitting the experimental data and the established mathematical model by using a Levenberg-Marquardt method to determine the reaction rate constant of the system.
In one embodiment, an extended Kalman filter technique is used to estimate the system's unmeasured state variables from the system's output chlorophyll fluorescence signal.
In one embodiment, the undetectable state variable is the amount of carbohydrate production, electron acceptor QA redox, electron acceptor QB redox, Plastoquinone (PQ).
In one embodiment, "PID control techniques are used to regulate photosynthesis activity. "the plant photosynthesis is regulated using a PID control algorithm using the estimated state variables as feedback signals.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of any of the methods when executing the program.
A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of any of the methods.
A processor for running a program, wherein the program when running performs any of the methods.
The invention has the beneficial effects that:
the method for estimating the photosynthesis state variable by using the extended Kalman filtering technology and adjusting the plant photosynthesis by using the most extensive PID control technology in the current industry has important significance for plant growth control.
Drawings
FIG. 1 is a graph of the fitting effect of experimental data and a mathematical model in the photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal.
FIG. 2 is a flow chart of the extended Kalman filtering technique in the photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal.
FIG. 3 is a comparison graph of the state estimator in the photosynthesis PID control method based on chlorophyll fluorescence as the feedback signal and the results obtained from computer simulation.
FIG. 4 is a block diagram of a closed loop system for regulating plant photosynthesis by the PID control technology in the photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal.
FIG. 5 is an example of a light intensity variation curve (a) and corresponding PQ reduction efficiency (b) obtained by using PID control technique in the photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal according to the present invention.
Detailed Description
The present invention is further described below in conjunction with the following figures and specific examples so that those skilled in the art may better understand the present invention and practice it, but the examples are not intended to limit the present invention.
The ChlF output of the system is measurable, the fluorescence curve can reflect the important characteristics of plant growth, and the ChlF is widely applied to the fields of fault diagnosis, automation, radar systems, missile tracking and the like, and the extended Kalman filtering technology can estimate the state quantity by utilizing the output quantity. The method for estimating the photosynthesis state variable by using the extended Kalman filtering technology provided by the invention adjusts the plant photosynthesis by using the most extensive PID control technology in the industry at present, and has important significance for plant growth control.
Step 1: testing of chlorophyll fluorescence curves using chlorophyll fluorometer
Chlorophyll fluorescence data was obtained from tea leaves collected from university campus (Wuxi, China). The leaves were picked in the morning of 2018, 7 months, and then quickly taken to a laboratory for testing. 12 leaves were used in the experiment, which were floated in water for two hours in order to prevent evaporation of water from the leaves; they were then dark-adapted for at least 15 minutes while consistently held in the clamp before the start of the ChlF measurement. ChlF was measured with FluorPen PSI (Photon Systems Instruments, Czech Reublic) and its OJIP protocol was selected. The intensity of illumination light was set to 750. mu. molphosons m-2s-1(the wavelength of light is in the range of 400-700 nm).
Step 2: and establishing a mathematical model according to the physiological process of plant photosynthesis.
1) Photoreaction stage and electron transfer
Photoreaction and subsequent electron transfer, from light absorption to the formation of NAPDH, also include a photoprotective mechanism when the illumination intensity is relatively strong.
2) Dark reaction stage
The carbon cycle is part of the "dark response" in photosynthesis, with the site of response being the chloroplast stroma. The cycle has three main phases: by CO2Carboxylated RuBP, CO2And reduction of ribulose 1, 5-bisphosphate (RuBP). Most plants incorporate immobilized CO by binding them to the second carbon of 5-carbon RuBP using Rubisco (RuBP carboxylase oxygenase)2A molecule. Thus CO2Is reduced. The six carbon compound thus formed is extremely unstable and immediately decomposes into two molecules of three carbon compound: 3-phosphoglycerate. The latter is then phosphorylated by ATP (formed by ATP synthase) to 1, 3-Diphosphoglycerate (DPGA) and then reduced by NADPH produced by two photoreactions. One of the triosephosphate molecules leaving the cycle is then used to synthesize glucose after a complex series of biochemical reactions. In general, 6 cycles produce one molecule of hexose. In addition, the cycle includes the formed ribulose-1, 5-bisphosphate molecule restarting the cycle. The chemical equation can be expressed as:
TABLE 1 notation of equations (1) - (22), corresponding substance concentrations and initializations
From equations (1) - (22), the following set of differential equations can be established:
and step 3: the experimental chlorophyll fluorescence curve was fitted to the mathematical model using the levenberg-marquardt algorithm.
The Levenberg-Marquardt method is one of the most common methods for system identification, and the fitting result of model simulation and experimental data is shown in the attached drawing 1, so that the model output can be well matched with the experimental data, and the photosynthesis model has good practical significance. Table 2 lists the values of the rate constants used for the model simulations.
Table 2 rate constant values for model simulations
FIG. 1 of the drawings | FIG. 1 of the drawings | FIG. 1 of the drawings | |||
k1u | 0.04 | k14 | 0.09 | k27 | 8.183 |
k2 | 28.34 | k15 | 796.5 | k28 | 2.812 |
k3 | 95.37 | k16 | 27.17 | k29 | 1.878 |
k4 | 1310.2 | k17 | 15.795 | k30 | 8.88 |
k5 | 0.012 | k18 | 38.41 | k31 | 6.13 |
k6 | 144.3 | k19 | 5.514 | PQ0 | 5.25 |
k7 | 0.02 | k20 | 10.53 | PC0 | 1.203 |
k8 | 19299 | k21 | 17.33 | Fd0 | 1.499 |
k9 | 21.31 | k22 | 24.68 | ATP0 | 18.14 |
k10 | 325.7 | k23 | 0.0036 | NADPH0 | 7.5 |
k11 | 481.1 | k24 | 0.227 | PGA0 | 0.024 |
k12 | 590.9 | k25 | 1.23×10-5 | G | 0.75 |
k13 | 8259.1 | k26 | 30.06 |
And 4, step 4: the state quantity is estimated from the chlorophyll fluorescence output of the system using extended kalman filtering techniques.
The Extended Kalman Filter (EKF) flow diagram is shown in fig. 2, and in the present invention, the actual output of the system is chlorophyll fluorescence. In fig. 2, the function f is the established photosynthesis model. The function h is F ═ Gk2x1。xkRepresenting model-based state variables (i.e., predicted values). y iskWhich may include process and measurement noise, corresponds to the actual fluorescence output value. Q and R are the covariance of the process noise and measurement noise possible during the chlorophyll fluorescence measurement (commonly known as gaussian noise). A and H can be calculated from the Jacobian matrix of the mathematical model. P is the chlorophyll fluorescence error estimation covariance matrix, and P0It can be arbitrarily initialized to a matrix of appropriate size because it will converge after a certain number of steps.
Is an estimated value. FIG. 3 shows the estimation of a state variable such as glucose from fluorescence values using EKF. It can be seen that the values obtained by model simulation have a small error with the results obtained by estimation, proving the use of EKF estimation based on measured ChlFThe process of the production amount of saccharides is efficient.
And 5: PID control techniques are used to regulate plant photosynthesis.
Reaction amount v ═ k15x7x8Is related to the reduction efficiency of PQ at a certain time. It can reflect the amount of photoelectrons used for forward photochemical reactions, which is positively correlated with the final sugar production. X can also be estimated using EKF based on measured ChlF7And x8The amount of (c). The dynamics of v are then controlled here using the most common PID (proportional, integral and derivative) closed loop control system. FIG. 4 shows a block diagram of a closed loop system for regulating plant photosynthesis using EKF-based PID control techniques. The controllable input is the light intensity u. The "controller" submodule in the system refers to the PID controller in this operation. By "power supply" sub-module is meant a light energy source such as a Light Emitting Diode (LED) or a high pressure sodium lamp (HPS). "execution unit" refers to a light regulator based on the controller output. The sub-module of the controlled plant refers to the established mathematical model. The "converter" submodule refers to the corresponding photosensor, which is replaced by a gain factor in the simulation. The reduction efficiency (v) of PQ is controlled to a desired reduction efficiency (v ^) using PID control as shown in FIG. 4. For example, we keep the light intensity constant under current greenhouse conditions, the reduction efficiency (v) of PQ will not be kept at a constant level because of the photoprotective mechanisms and other reaction kinetics reactions. However, by setting the PID control to a constant value to achieve the required PQ reduction efficiency (v ^), we can keep v constant using the controller in FIG. 4. Figure 5a shows the controlled post-autotuned lighting dynamics, generating control signals based on the deviation of the set desired PQ recovery efficiency (v ^) from the actual PQ recovery efficiency (v) of the system, automatically adjusting the voltage level of the power supply and gradually reaching the dynamic equilibrium, thus producing the illumination intensity shown in figure 5 (a). As can be seen from FIG. 5(b), the PQ reduction efficiency can be maintained constant in a PID controller in a short time, which is equivalent to artificially regulating the photosynthesis activity of plants, and is significant for the generation of plant photosynthesis efficiency and biomass, and the PQ reduction efficiency can be obtained according to the photosynthesis efficiency requiredThe control system is reasonably designed and the illumination resources are reasonably utilized. Constant levels of lighting intensity are commonly applied in greenhouse production, but too much light can lead to wasted energy, even affecting the photosynthesis mechanism and even damaging the plants due to photoprotection. This example illustrates how the photosynthesis-related quantity can be controlled to the desired dynamics, although it cannot be measured directly. The invention successfully applies the most common extended Kalman filtering technology and PID control technology in industrial application to the regulation of plant photosynthesis, has obvious regulation effect and has important significance for plant growth control.
The above-mentioned embodiments are merely preferred embodiments for fully illustrating the present invention, and the scope of the present invention is not limited thereto. The equivalent substitution or change made by the technical personnel in the technical field on the basis of the invention is all within the protection scope of the invention. The protection scope of the invention is subject to the claims.
Claims (4)
1. A photosynthesis PID control method based on chlorophyll fluorescence as a feedback signal is characterized by comprising the following steps:
step 1, acquiring a fluorescence curve by using a chlorophyll fluorescence instrument;
step 2, establishing a mathematical model of the whole process of photosynthesis;
step 3, fitting a mathematical model and a fluorescence curve obtained by an experiment, and determining a series of reaction rate constants;
step 4, estimating a state variable according to the fluorescence value output by the model and by utilizing an extended Kalman filtering technology;
step 5, regulating the photosynthesis activity by using a PID control technology;
step 2, establishing a mathematical model of the whole process of photosynthesis; the method specifically comprises the following steps:
1) photoreaction stage and electron transfer
Photoreaction and subsequent electron transfer, from light absorption to formation of NAPDH, and photoprotection mechanisms when the illumination intensity is relatively strong;
TABLE 1 notation of equations (1) - (16), corresponding substance concentrations and initializations
2) Dark reaction stage
The carbon cycle is part of the "dark response" in photosynthesis, with the site of response being the chloroplast stroma; the cycle has three main phases: by CO2Carboxylated RuBP, CO2Reduction of (3) and reduction of RuBP; most of themPlants incorporate immobilized CO by binding them to the second carbon of the 5-carbon RuBP using RuBP carboxylase oxygenase2A molecule; thus CO2Is reduced; the six carbon compound thus formed is extremely unstable and immediately decomposes into two molecules of three carbon compound: 3-phosphoglyceric acid; the latter is then phosphorylated by ATP to 1, 3-diphosphoglycerate, and then reduced by NADPH produced by two photoreactions to produce 3-triosephosphate; then, after a series of complex biochemical reactions, one of the triose-3-phosphate molecules is recycled off carbon for the synthesis of glucose; according to the light reaction and dark reaction processes, the following differential equation sets are established:
and step 3: fitting a mathematical model and a fluorescence curve obtained by an experiment to determine a series of reaction rate constants; "in, fit the chlorophyll fluorescence curve and mathematical model using Levenberg-Marquardt algorithm;
the following table lists the values of the rate constants used for the model simulation;
Step 4, estimating the state quantity of chlorophyll fluorescence from the output quantity of the system by using an extended Kalman filtering technology; the method specifically comprises the following steps:
first, settingAnd P0Then according to the formulaComputingAndthen according toCalculating KkAnd then calculate Simultaneous calculationThe actual output of the system is chlorophyll fluorescence; the function f is the established photosynthesis model; the function h is F ═ Gk2x1;xkRepresenting model-based state variables; y iskCorresponding to the actual fluorescence output value, which may include process and measurement noise; q and R are the covariance of the possible process noise and measurement noise during chlorophyll fluorescence measurement; a and H can be calculated from the Jacobian matrix of the mathematical model; p is the chlorophyll fluorescence error estimation covariance matrix, and P0Can be arbitrarily initiatedTo a matrix of appropriate size, since it will converge after a certain number of steps;
step 5, regulating the photosynthesis of the plant by using a PID control technology; "specifically includes
Reaction amount v ═ k15x7x8Relating to the reduction efficiency of PQ at a certain moment; it may reflect the amount of photoelectrons used for the forward photochemical reaction, which is positively correlated with the final sugar production; x can also be estimated using EKF based on measured ChlF7And x8The amount of (c); then the most common PID closed loop control system is adopted to control the dynamics of v; firstly, the difference value between a set value v ^ and the actual output v of a control system is used as the input of a controller, the controller feeds back a control result to an execution unit after generating a control action, so that the execution unit can automatically adjust the power of a power supply to generate different illumination intensities u, the different illumination intensities act on a controlled plant, the blade of the controlled plant scatters out a fluorescence value, meanwhile, the actual output v of the control system is estimated through EKF according to the fluorescence value, and finally, the difference value between the set value and an experimental value is used as the input of a controller of the next sampling time; the controllable input is the light intensity u; the 'controller' submodule in the system refers to a PID controller in the work; the "power supply" submodule refers to a light energy source; "execution unit" refers to a light regulator based on the controller output; the sub-module of the controlled plant refers to an established mathematical model; the "converter" submodule refers to the corresponding photosensor, which is replaced by a gain factor in the simulation.
2. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the steps of the method of claim 1 are performed when the program is executed by the processor.
3. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method as claimed in claim 1.
4. A processor, characterized in that the processor is configured to run a program, wherein the program when running performs the method of claim 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910017944.2A CN109634101B (en) | 2019-01-09 | 2019-01-09 | Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910017944.2A CN109634101B (en) | 2019-01-09 | 2019-01-09 | Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109634101A CN109634101A (en) | 2019-04-16 |
CN109634101B true CN109634101B (en) | 2021-03-30 |
Family
ID=66060377
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910017944.2A Active CN109634101B (en) | 2019-01-09 | 2019-01-09 | Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109634101B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112039967A (en) * | 2020-08-25 | 2020-12-04 | 江南大学 | Open Internet of things platform for feedback control of plant growth |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106130426A (en) * | 2016-07-18 | 2016-11-16 | 南京理工大学 | The permagnetic synchronous motor method for controlling number of revolution of ultrahigh speed without sensor based on EKF |
CN107404787A (en) * | 2017-07-05 | 2017-11-28 | 江苏大学 | A kind of LED light fillings control system and control method based on actual photosynthesis efficiency |
CN108614415A (en) * | 2018-03-27 | 2018-10-02 | 苏州联点数据技术有限公司 | A kind of intelligence control system and control method of plant growth factor |
CN108663928A (en) * | 2017-04-01 | 2018-10-16 | 东北林业大学 | A kind of Greenhouse Intelligent Control System based on fuzzy |
-
2019
- 2019-01-09 CN CN201910017944.2A patent/CN109634101B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106130426A (en) * | 2016-07-18 | 2016-11-16 | 南京理工大学 | The permagnetic synchronous motor method for controlling number of revolution of ultrahigh speed without sensor based on EKF |
CN108663928A (en) * | 2017-04-01 | 2018-10-16 | 东北林业大学 | A kind of Greenhouse Intelligent Control System based on fuzzy |
CN107404787A (en) * | 2017-07-05 | 2017-11-28 | 江苏大学 | A kind of LED light fillings control system and control method based on actual photosynthesis efficiency |
CN108614415A (en) * | 2018-03-27 | 2018-10-02 | 苏州联点数据技术有限公司 | A kind of intelligence control system and control method of plant growth factor |
Non-Patent Citations (1)
Title |
---|
Modelling and simulation of photosystem II chlorophyll fluorescence transition from darkadapted state to light-adapted state;Shaopeng Feng等;《IET Systems Biology》;20181229;第12卷(第6期);第289-293页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109634101A (en) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Automatic carbon dioxide enrichment strategies in the greenhouse: A review | |
Van Henten et al. | Time-scale decomposition of an optimal control problem in greenhouse climate management | |
US20120107921A1 (en) | Model based controls for use with bioreactors | |
CN109634101B (en) | Photosynthesis PID control method based on chlorophyll fluorescence as feedback signal | |
CN111557159A (en) | Water and fertilizer integrated machine control system with crop nutrient demand analysis function and control method thereof | |
US20190204228A1 (en) | Information processing apparatus, information processing method, program, and sensing apparatus | |
Zhang et al. | Robust model-based reinforcement learning for autonomous greenhouse control | |
Yuan et al. | An open Internet of Things (IoT)-based framework for feedback controlof photosynthetic activities | |
Rizkiana et al. | Crop growth prediction model at vegetative phase to support the precision agriculture application in plant factory | |
Outanoute et al. | Synthesis of an Optimal Dynamic Regulator Based on Linear Quadratic Gaussian (LQG) for the Control of the Relative Humidity Under Experimental Greenhouse. | |
Vega et al. | Fuzzy-based automated nutrient solution control for a hydroponic tower system | |
CN117032360A (en) | Plant factory light environment regulation and control system based on photosynthetic simulation and spectrum superposition | |
Gessler et al. | Beyond source and sink control–toward an integrated approach to understand the carbon balance in plants | |
Yu et al. | Crop photosynthetic performance monitoring based on a combined system of measured and modelled chloroplast electron transport rate in greenhouse tomato | |
Ifrim et al. | Multivariable control strategy for the photosynthetic cultures of microalgae | |
Oktarina et al. | Digitized Smart Solar Powered Agriculture Implementation in Palembang, South Sumatra | |
Ödük et al. | The advantages of fuzzy control over traditional control system in greenhouse automation | |
Chakravarty et al. | Quantitative feedback theory-based multi-variable robust control for soil quality improvement in a drip irrigated field | |
Gongxian et al. | Optimizing control of bio-dissimilation process of glycerol to 1, 3-propanediol | |
CN116419455A (en) | Plant environment light supplementing system and method based on sliding mode control | |
Choudhury et al. | Optimized Crop Detection Using IoT and Deep Neural Networks | |
CN117236650B (en) | Water-fertilizer integrated intelligent pump house control method | |
NL2027936B1 (en) | Method of cultivating plants and system therefor | |
Cong et al. | Automated Hydroponic System for Indoor Plant Growth in Urban Area | |
Schaum et al. | Robust adaptive feedforward output-feedback tracking control for microalgae cultures |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |