CN108090693A - The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application - Google Patents
The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application Download PDFInfo
- Publication number
- CN108090693A CN108090693A CN201711494114.6A CN201711494114A CN108090693A CN 108090693 A CN108090693 A CN 108090693A CN 201711494114 A CN201711494114 A CN 201711494114A CN 108090693 A CN108090693 A CN 108090693A
- Authority
- CN
- China
- Prior art keywords
- mrow
- photoresponse
- photosynthetic
- model
- curvature
- 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
Links
- 230000000243 photosynthetic effect Effects 0.000 title claims abstract description 54
- 230000033228 biological regulation Effects 0.000 title claims abstract description 28
- 230000004927 fusion Effects 0.000 title claims abstract description 24
- 238000005286 illumination Methods 0.000 claims abstract description 19
- 238000012937 correction Methods 0.000 claims abstract description 17
- 238000010276 construction Methods 0.000 claims abstract description 10
- 238000000034 method Methods 0.000 claims description 27
- 230000029553 photosynthesis Effects 0.000 claims description 13
- 238000010672 photosynthesis Methods 0.000 claims description 13
- KBPHJBAIARWVSC-RGZFRNHPSA-N lutein Chemical compound C([C@H](O)CC=1C)C(C)(C)C=1\C=C\C(\C)=C\C=C\C(\C)=C\C=C\C=C(/C)\C=C\C=C(/C)\C=C\[C@H]1C(C)=C[C@H](O)CC1(C)C KBPHJBAIARWVSC-RGZFRNHPSA-N 0.000 claims description 10
- 229960005375 lutein Drugs 0.000 claims description 10
- KBPHJBAIARWVSC-XQIHNALSSA-N trans-lutein Natural products CC(=C/C=C/C=C(C)/C=C/C=C(C)/C=C/C1=C(C)CC(O)CC1(C)C)C=CC=C(/C)C=CC2C(=CC(O)CC2(C)C)C KBPHJBAIARWVSC-XQIHNALSSA-N 0.000 claims description 10
- FJHBOVDFOQMZRV-XQIHNALSSA-N xanthophyll Natural products CC(=C/C=C/C=C(C)/C=C/C=C(C)/C=C/C1=C(C)CC(O)CC1(C)C)C=CC=C(/C)C=CC2C=C(C)C(O)CC2(C)C FJHBOVDFOQMZRV-XQIHNALSSA-N 0.000 claims description 10
- 235000008210 xanthophylls Nutrition 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 8
- 230000029058 respiratory gaseous exchange Effects 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000002068 genetic effect Effects 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 abstract description 9
- 230000001105 regulatory effect Effects 0.000 abstract description 5
- 235000000832 Ayote Nutrition 0.000 abstract description 4
- 235000009854 Cucurbita moschata Nutrition 0.000 abstract description 4
- 240000001980 Cucurbita pepo Species 0.000 abstract description 4
- 235000009804 Cucurbita pepo subsp pepo Nutrition 0.000 abstract description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 235000015136 pumpkin Nutrition 0.000 abstract description 4
- 230000007423 decrease Effects 0.000 abstract description 3
- 102000004127 Cytokines Human genes 0.000 abstract 1
- 108090000695 Cytokines Proteins 0.000 abstract 1
- 230000004907 flux Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 244000241257 Cucumis melo Species 0.000 description 4
- 235000009847 Cucumis melo var cantalupensis Nutrition 0.000 description 4
- 108090000790 Enzymes Proteins 0.000 description 3
- 102000004190 Enzymes Human genes 0.000 description 3
- 238000005452 bending Methods 0.000 description 3
- 230000021523 carboxylation Effects 0.000 description 3
- 238000006473 carboxylation reaction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 108010003581 Ribulose-bisphosphate carboxylase Proteins 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000012010 growth Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 206010053759 Growth retardation Diseases 0.000 description 1
- 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 description 1
- 238000009825 accumulation Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 238000005842 biochemical reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 231100000001 growth retardation Toxicity 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 229930027945 nicotinamide-adenine dinucleotide Natural products 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000008635 plant growth Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- A—HUMAN NECESSITIES
- A01—AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
- A01G—HORTICULTURE; CULTIVATION OF VEGETABLES, FLOWERS, RICE, FRUIT, VINES, HOPS OR SEAWEED; FORESTRY; WATERING
- A01G7/00—Botany in general
- A01G7/04—Electric or magnetic or acoustic treatment of plants for promoting growth
- A01G7/045—Electric or magnetic or acoustic treatment of plants for promoting growth with electric lighting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- H—ELECTRICITY
- H05—ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
- H05B—ELECTRIC HEATING; ELECTRIC LIGHT SOURCES NOT OTHERWISE PROVIDED FOR; CIRCUIT ARRANGEMENTS FOR ELECTRIC LIGHT SOURCES, IN GENERAL
- H05B47/00—Circuit arrangements for operating light sources in general, i.e. where the type of light source is not relevant
- H05B47/10—Controlling the light source
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P60/00—Technologies relating to agriculture, livestock or agroalimentary industries
- Y02P60/14—Measures for saving energy, e.g. in green houses
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Biodiversity & Conservation Biology (AREA)
- Forests & Forestry (AREA)
- Primary Health Care (AREA)
- Agronomy & Crop Science (AREA)
- Marine Sciences & Fisheries (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Botany (AREA)
- Ecology (AREA)
- Animal Husbandry (AREA)
- Environmental Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Mining & Mineral Resources (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Cultivation Of Plants (AREA)
- Greenhouses (AREA)
Abstract
The present invention couples photosynthetic desired value model for the problem that the multiple-factor that facilities environment regulator control system lacks fusion regulation and control efficiency, causes production efficiency low, studies the photosynthetic desired value regulation-control model construction method of facility of fusion efficiencies.Pass through intensity of illumination, environment temperature and CO2Concentration obtains the photosynthetic data sample of pumpkin seedling for the multiple-factor coupled nesting experiment of variable, construct photoresponse equilateral hyperbola correction model, upper limit light intensity point is limited using the theoretical light obtained under different specific temperatures of curvature of curve, the pumpkin seedling luminous environment regulation-control model of the fusion efficiencies constraint under being fitted to obtain different temperatures accordingly in this, as goal of regulation and control value.The model is compared with model of the tradition using photosynthetic rate maximum of points as goal of regulation and control value, photosynthetic rate is average only to decline 4.48%, and light requirement then averagely declines 29.91%, it can solve the problems, such as the basic theory of protected crop environment high efficiency regulatory, have directive significance to actual facilities environment cytokine regulatory.
Description
Technical field
The invention belongs to reading intelligent agriculture technical field, more particularly to a kind of photosynthetic desired value of facility of fusion efficiencies constraint
The structure of Optimum Regulation model and application.
Background technology
Industrialized agriculture is to provide optimum environment and condition for crop growth each stage, obtains good quality and high output agricultural product
Modern agriculture business activities, by the end of 2014, China's facility gross area is more than global 85%, up to 410.9 ten thousand hectares, production
800,000,000,000 yuan of value or more, the gross area and total output occupy first, the whole world.But the industrialized agriculture starting in China is being set than later
Applying the countries such as luminous environment intelligent control methods and techniques aspect and the U.S., Britain, Holland, also there is a big difference, so state in recent years
Interior correlation scholar has done numerous studies in this respect.Illumination is to determine the yield of crop as the main factor of photosynthesis is influenced
With the key of quality.Since greenhouse is hidden be subject to obliquity structure, covering material, sun altitude and clean surface degree
The influences such as light, photon flux density are only the 30%-70% in crop field.Photon flux density is difficult to the life for meeting crop in greenhouse
Long demand causes the reduction of various physiological parameters of the facility including plant Net Photosynthetic Rate, will cause growth retardation,
Fruit-setting rate reduces, and can seriously affect the yield and quality of agricultural product, thus builds the luminous environment for being conducive to plant growth, is to realize
The major issue that facility high-efficiency high-quality produces, quickening industry synergy upgrading and transition and upgrade are urgently to be resolved hurrily.
Correlative study shows in addition to intensity of illumination, and it is influence photosynthetic rate main that the factors such as environment temperature, which are also,
Factor, and be at different conditions dynamic change to the influence degree of photosynthetic rate between the different factors.Existing crop is artificial
Light filling technology realizes the regulation and control to facility luminous environment by threshold value mostly, does not account for shadow of the Multi-environment factor to light saturation point
It rings, although it is ineffective to improve photosynthesis high efficiency regulatory;On this basis, this research team proposes multiple-factor fusion
Luminous environment regulates and controls method, by merging multi-environment factor to photosynthetic influence, constructs light saturation point dynamic acquisition method,
The Optimum Regulation on demand of luminous environment is realized to a certain degree.But further analysis crop photoresponse curve can be found that light by
It is cumulative it is strong during plant experienced the process that Rubisco limitations are restricted to from BuBP, wherein rapid increase partial straight lines
Slope is exactly the quantum efficiency of photosynthetic carbon assimilation, represents that photosynthesis is limited by electron transport rate, RuBP carboxylation abilities, electricity
Sub- transfer rate is limited again by photon flux density, just influences the principal element of the part, this partial trace is referred to as light
Restricted part;And curve is slowly converted after reaching saturation state, the effect of enzyme can not be with light reaction in Calvin-Benson cycle
The generation of ATP and NADPH is harmonious in journey, and carboxylation ability of enzyme etc. becomes main limiting factor, and temperature is as influence
An important factor for Rubisco activated enzyme activities etc., also replace light to become the principal element for limiting the stage photosynthetic rate.And light
Saturation point appears in the gradual first transition of curve, and light has not been the principal element for influencing photosynthetic rate variation, but above-mentioned grinds
Main method is studied carefully using light saturation point light filling, does not account for limiting upper limit light intensity point as goal of regulation and control value using light, although can be with
Reach maximum photosynthesis rate, but due to the great decline of the efficiency of light energy utilization, light filling efficiency is caused to significantly reduce.
The content of the invention
The shortcomings that in order to overcome the above-mentioned prior art, it is an object of the invention to provide a kind of facilities of fusion efficiencies constraint
The structure of the Optimum Regulation model of photosynthetic desired value and application limit upper limit light intensity point as light goal of regulation and control value using light, are limited from light
Upper limit light intensity point acquisition methods processed are started with, by analyzing crop photosynthesis mechanism and the Major environment impacts factor, devise with
Intensity of illumination, environment temperature and CO2Concentration is the multiple-factor coupled nesting experiment of experiment variable, utilizes the photosynthetic speed of Li-6400XT
Rate instrument obtains the photosynthetic data sample of pumpkin seedling;Analyze the influence of intensity of illumination, environment temperature to crop photosynthesis rate, structure
Photoresponse equilateral hyperbola correction model is built;Upper limit light intensity is limited using the theoretical light obtained under specific temperature of curvature of curve
Point builds the crop seedling luminous environment regulation-control model under temperature restraint in this, as goal of regulation and control value.
To achieve these goals, the technical solution adopted by the present invention is:
A kind of construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint, including walking as follows
Suddenly:
The photosynthetic data sample of crop seedling is obtained, analyzes the influence of intensity of illumination, environment temperature to crop photosynthesis rate,
Build photoresponse equilateral hyperbola correction model;
Utilize the theoretical curvature maximum obtained under different temperatures on photoresponse curve of curvature of curve, the corresponding light quantity of the value
Sub- flux is the light limitation upper limit light intensity point under relevant temperature, builds the fusion efficiencies constraint under different temperatures on this basis
The photosynthetic desired value regulation-control model of facility.
Using matlab to the photoresponse model regression fit under different temperatures, photochemistry rate, Xanthophyll cycle item, light are satisfied
Photoresponse equilateral hyperbola correction model formula is brought into item and Dark respiration rate obtains photoresponse side under different temperatures
Journey, the photoresponse equilateral hyperbola correction model formula are as follows:
Wherein, PnRepresenting Net Photosynthetic Rate, α represents photochemistry rate, and I represents pharosage, and β represents Xanthophyll cycle item,
γ represents light saturation item, RdRepresent Dark respiration rate.
The process for obtaining point of maximum curvature on photoresponse curve is as follows:
Data are fitted first with the photoresponse equilateral hyperbola correction model, obtain a photoresponse curve,
And it is normalized;
Photoresponse function is defined as sign function, seeks the curvature of curve of photoresponse function after normalization, is calculated using heredity
Method to photoresponse function curvature carry out optimizing, obtain the maximum of points of curvature, on this basis after renormalization curvature maximum
It is light limitation upper limit light intensity point to be worth corresponding illumination spot.
To photoresponse curve normalized under the specific temperature, generation photoresponse function Pn'(I), the normalizing
Change method is linear transformation, Logarithm conversion or arc cotangent are converted.
Curvature acquisition is carried out to photoresponse function under normalized specific temperature, formula is as follows:
Wherein, K represent photoresponse function curvature function, (Pn'(I)) " represent photoresponse function second dervative function,
(Pn'(I)) ' represent photoresponse function first derivative function.
Optimizing is carried out to photoresponse function curvature curve using genetic algorithm, obtains the maximum of points of curvature, it is basic herein
The corresponding illumination spot of the maximum of curvature is light limitation upper limit light intensity point after upper renormalization.
On the basis of obtaining all light limitation upper limit light intensity points, it is because becoming to limit upper limit light intensity as independent variable, light using temperature
Amount, fitting obtain the luminous environment goal of regulation and control value model under different temperatures.
This model can automatically generate the Xanthophyll cycle upper limit under different temperature environments under conditions of temperature is merely entered
Light intensity.Therefore, after this model adds in luminous environment regulator control system, model can complete the real-time automatic calculating of Xanthophyll cycle upper limit light intensity,
The high efficiency regulatory of luminous environment is realized in this, as light goal of regulation and control value.Compared with prior art, the present invention is proposed based on curvature
Theoretical makees object light limitation upper limit light intensity point acquisition methods.It is real according to the photosynthetic rate coupling nested with intensity of illumination of crop temperature
It tests as a result, being fitted photoresponse equilateral hyperbola correction model, introducing curvature concept realizes light by calculating curvature maximum point
The acquisition of upper limit light intensity point is limited, and luminous environment goal of regulation and control value model is constructed with this, compared to using light saturation point as tune
Control desired value point high degree improves regulation and control efficiency.
Description of the drawings
Fig. 1 is the photosynthetic desired value regulation-control model algorithm flow chart of present invention structure fusion efficiencies constraint.
Fig. 2 is curvature maximum optimizing algorithm flow chart of the present invention structure based on genetic algorithm.
Fig. 3 is that the present invention obtains maximum curvature procedure chart.
Fig. 4 is the fusion efficiencies constraint luminous environment Optimum Regulation desired value illustraton of model that the present invention is built.
Specific embodiment
The embodiment that the present invention will be described in detail with reference to the accompanying drawings and examples.
A kind of building process of the photosynthetic desired value regulation-control model of facility of fusion efficiencies constraint of the present invention is as follows:
1st, plan design
This experiment is carried out in May, -2017 in April, 2017 in Xibei Univ. of Agricultural & Forest Science & Technology's agricultural machinery laboratory.Muskmelon Seedlings
It is positioned in the MD1400 incubators of Dutch Sinder companies production and is cultivated, the light source in incubator is by red and blueness
Lamp bead forms, and environment temperature is set to 25 DEG C, and relative humidity is set to 60%, CO2Concentration is set to 400 μm of ol/mol, and the photoperiod is set to
14 hours, incubator have built suitable growing environment for Muskmelon Seedlings.During experiment, crop seedling is carried out normal
Cultivation management.
For trial crops Muskmelon Seedlings, based on the influence of varying environment factor pair photosynthetic rate, U.S. LI-COR is used
The portable photosynthetic speedometers of Li-6800XT of company's production measure, and experiment sets external environment variable to include photon flux
Density, leaf temperature etc. design the nested experiment of each envirment factor under the environmental condition of suitable for crop growth.Wherein photon flux
Density is 0 μm of ol/m2S, 50 μm of ol/m2S, 100 μm of ol/m2S, 200 μm of ol/m2S, 300 μm of ol/m2S, 500 μm of ol/m2S, 700 μ
mol/m2S, 900 μm of ol/m2S, 1100 μm of ol/m2S, 1300 μm of ol/m2S, 1700 μm of ol/m2S, 2000 μm of ol/m212 ladders such as s
Degree, 14 DEG C of temperature setting, 18 DEG C, 22 DEG C, 26 DEG C, 30 DEG C, 34 DEG C, 7 gradients such as 38 DEG C.According to the model for setting each factor above
It encloses and step-length, to reduce the contingency of test data, every group of experiment Nested conditions are repeated three times progress.Based on above experiment side
Case measures photosynthetic rate of the Muskmelon Seedlings under different photon flux densities, leaf temperature, to build photosynthetic rate prediction model
Data are provided to support.
2nd, the structure of photosynthetic desired value regulation-control model
This model obtains the photosynthetic data sample of pumpkin seedling using Li-6400XT photosynthetic rate instrument;It is strong to analyze illumination
Degree, influence of the environment temperature to crop photosynthesis rate, construct photoresponse equilateral hyperbola correction model;It is managed using curvature of curve
Upper limit light intensity point is limited by the light obtained under specific temperature, the fusion efficiencies constraint under building different temperatures on this basis accordingly
The photosynthetic desired value regulation-control model of facility, flow chart is as shown in Figure 1.
2.1st, photoresponse model construction
Charles-Edwards Photosynthesis Models have taken into full account the influence of light respiration, dark respiration and oxygen effect, pass through structure
It builds CO2Concentration from cell inner air gap to mesophyll tissue in photosynthesis position diffusion process with simplify biochemical reaction
The association of system constructs the mathematical model comprising multiple physiological parameters, relation complexity, as shown in formula (1).
Pn 2rd-Pn[αI(rx+rd)+C-Rrd]+αIC-R(αIrx+ C)=0 (1)
In formula:PnRepresenting Net Photosynthetic Rate, α represents photochemistry rate, and I represents pharosage, and β represents Xanthophyll cycle item,
γ represents light saturation item, rdRepresent diffusional resistance;α represents photochemistry rate;rxRepresent carboxylation resistance, C represents environment CO2Concentration,
R represents constant respiratory rate.The model reflects accumulation and the assigning process of photosynthate, is moved to completing crop photosynthesis rate
State emulates, and crop yield and quality-improving is promoted to be of great significance, but there are the problem of physiological parameter is more, model is complicated.For
Simplified model proposes right angle hyperbolic model, on-right angle hyperbolic model, exponential relationship model and modified index model etc. successively
Photosynthetic rate model, model formation is respectively such as (2), (3), (4), (5).
Equilateral hyperbola model:
On-right angle hyperbolic model:
Exponential model:
Equilateral hyperbola correction model:
In formula:K represents convexity, PmaxRepresent maximum photosynthesis rate, RdRepresent Dark respiration rate.Based on experimental data, utilize
Above-mentioned model, which carries out nonlinear regression, can estimate model parameter, and complete value range limitation according to biological significance, so as to
The structure of implementation model.But since the difference of model itself is in processing Xanthophyll cycle type (strong Xanthophyll cycle), (the dim light suppression of light saturation type
System), fitting precision difference is apparent when different types of photoresponse curve for unsaturated type (saturation convergence), and equilateral hyperbola is repaiied
The versatility of its model of positive model is relatively more preferable, therefore the present invention builds photosynthetic respective mode using equilateral hyperbola correction model
Type.
The expression formula of equilateral hyperbola correction model according to formula (5) above is using matlab under different temperatures
Photoresponse model regression fit, fitting result is as shown in table 1, from table root-mean-square error and the goodness of fit analysis it can be found that
This models fitting works well.It is brought into using initial slope, Xanthophyll cycle item, light saturation item and the Dark respiration rate in formula following table
Formula (5) can obtain the photoresponse equation under different temperatures.
1 equilateral hyperbola correction model parameter list of table.
2.2nd, the photosynthetic desired value regulation-control model of facility of fusion efficiencies constraint
Tangent line turns over the limit of angular dimension in curvature of curve expression unit segmental arc, is for describing curved degree
Parameter, curvature is bigger, shows that curve is also bigger in the bending degree of the point, in understanding process, inverse-curvature half of curvature
Footpath may may be understood more readily, and radius of curvature was exactly that inscribed circle is a little done on curve, and the radius of the circle is exactly the point on curve
Radius of curvature, the curvature of the point is smaller, and radius of curvature is bigger, shows that curve is smaller in the bending degree of the point, calculates public
Shown in formula such as formula (6).
Wherein, K represents bent curvature of a curve, and the single order that y ' and y " represent curve respectively is led and led with second order.
And in the such curve being of practical significance of photoresponse curve that the present invention studies, the bending journey of photoresponse curve
Degree represents and is meant that influence of the intensity of illumination to Net Photosynthetic Rate, according to the concept and meaning of above-mentioned curvature, photoresponse curve
Point of maximum curvature just represent the inflection point that light intensity influences Net Photosynthetic Rate, from be affected become influence it is no longer so bright
It is aobvious, that is, light limitation upper limit light intensity point has been corresponded to, so how bent to photoresponse using bent curvature of a curve next will introduce
Line carries out light limitation upper limit light intensity point and obtains.
Its process is fitted data using photoresponse equilateral hyperbola correction model, obtains a photoresponse curve
Afterwards, it is normalized first, because the dimension disunity of photosynthetic rate and intensity of illumination, concentration etc., in order to prevent
Different magnitude of numerical value influences each other, and the phenomenon that " big value eat small value " occurs, the various regressing calculations of progress, matrix operation and
During neural metwork training, data with existing will be normalized first.Common data normalization method is linear turn main
It changes, the conversion of Logarithm conversion and arc cotangent, the present invention uses linear function transformation, shown in linear normalization principle such as formula (7):
Y=(ymax-ymin)*(x-xmin)/(xmax-xmin)+ymin (7)
In formula:xmax、xminThe maximum and minimum value of data, y before respectively normalizingmax、yminBefore respectively normalizing
The upper and lower bound of data, x, y are respectively the data before normalizing and after normalization.
Secondly, curvature estimation is carried out on the basis of normalized photoresponse model is completed.It is asked since curvature needs being asked to use
It leads, only symbolic variable can carry out derivation operations, so photoresponse function is defined as sign function by the present invention middle, utilize
Formula (8) seeks the curvature of photoresponse function after normalization, and formula (8) is as follows:
Wherein, K represent photoresponse function curvature function, (Pn'(I)) " represent photoresponse function second dervative function,
(Pn'(I)) ' represent photoresponse function first derivative function.
Again, using the curvature of photoresponse function under different illumination as fitness value, the curvature based on genetic algorithm is established
Maximum of points optimizing, by repeatedly selecting, intersecting and mutation operation obtains maximum curvature value and corresponding normalization illumination spot,
The corresponding illumination spot of the maximum of curvature is that light limits upper limit light intensity point, flow chart such as Fig. 2 after renormalization on the basis of this
It is shown.
The light limitation upper limit light intensity point under different temperatures is obtained using the above method, wherein the light limitation under the conditions of 14 DEG C
Upper limit light intensity point acquisition process is as shown in Figure 3.According to test data and the above process, the upper limit is limited by independent variable, light of temperature
Light intensity is dependent variable, and fitting obtains luminous environment goal of regulation and control value under different temperatures as shown in figure 3, its model of fit such as formula (9)
It is shown.
PFD=-0.014 × T4+1.37×T3-46.85×T2+752.42×T-3932.77 (9)
The model can obtain the luminous environment goal of regulation and control that fusion efficiencies constrain under arbitrary temp in 14 DEG C of -38 DEG C of temperature ranges
Value, the coefficient of determination 0.995 show that model has fine fitting effect.
3rd, energy-saving effect discussion
In order to further verify the regulating effect of this method, new method is obtained with conventional method in the present invention luminous environment tune
Control desired value is compared, and comparing result is as shown in table 2.
2 two methods photosynthetic rate desired value point of table and photosynthetic rate maximum point contrast table
It can be found that although the decline degree of light filling amount is different under different temperatures, on the whole, under light filling amount from table
The reduction degree of drop degree and photosynthetic rate is much bigger.Maximum intensity of illumination is made under the light restrictive condition obtained using this method
For desired value, compared in the method for traditional light saturation point desired value, photosynthetic rate is average only to reduce 4.48%, and its light requirement
Amount then averagely declines 29.91%, and the protected crop luminous environment regulation-control model for illustrating to be built warm optical coupling using this method is efficiently saved
Can the characteristics of, provide fundamental basis and technical support for the efficient control accurate of protected crop luminous environment.
Claims (7)
- A kind of 1. construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint, which is characterized in that bag Include following steps:Obtain the photosynthetic data sample of crop seedling, the influence of analysis intensity of illumination, environment temperature to crop photosynthesis rate, structure Photoresponse equilateral hyperbola correction model;Using the theoretical curvature maximum obtained under different temperatures on photoresponse curve of curvature of curve, the corresponding light quantum of the value leads to Amount is the light limitation upper limit light intensity point under relevant temperature, builds setting for the fusion efficiencies constraint under different temperatures on this basis Apply photosynthetic desired value regulation-control model.
- 2. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 1, It is characterized in that, using matlab to the photoresponse model regression fit under different temperatures, by photochemistry rate, Xanthophyll cycle item, Light saturation item and Dark respiration rate are brought the light that photoresponse equilateral hyperbola correction model formula is obtained under different temperatures into and are rung Equation is answered, the photoresponse equilateral hyperbola correction model formula is as follows:<mrow> <msub> <mi>P</mi> <mi>n</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&alpha;</mi> <mi>I</mi> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <mi>I</mi> <mo>-</mo> <mi>&beta;</mi> <mi>I</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>I</mi> <mo>+</mo> <mi>&gamma;</mi> <mi>I</mi> </mrow> </mfrac> <mo>-</mo> <msub> <mi>R</mi> <mi>d</mi> </msub> </mrow>Wherein, PnRepresent Net Photosynthetic Rate, α represents photochemistry rate, and I represents pharosage, and β represents Xanthophyll cycle item, and γ is represented Light saturation item, RdRepresent Dark respiration rate.
- 3. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 1, It is characterized in that, the process for obtaining point of maximum curvature on photoresponse curve is as follows:Data are fitted first with the photoresponse equilateral hyperbola correction model, obtain a photoresponse curve, and it is right It is normalized;Photoresponse function is defined as sign function, the curvature of curve of photoresponse function after normalization is sought, utilizes genetic algorithm pair Photoresponse function curvature carry out optimizing, obtain the maximum of points of curvature, on this basis after renormalization curvature maximum pair The illumination spot answered is light limitation upper limit light intensity point.
- 4. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 3, It is characterized in that, to photoresponse curve normalized under specific temperature, generation photoresponse function Pn'(I), the normalization side Method is linear transformation, Logarithm conversion or arc cotangent are converted.
- 5. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 4, It is characterized in that, carrying out curvature acquisition to photoresponse function under normalized specific temperature, formula is as follows:<mrow> <mi>K</mi> <mo>=</mo> <mfrac> <mrow> <mo>|</mo> <msup> <mrow> <mo>(</mo> <msup> <mi>Pn</mi> <mo>&prime;</mo> </msup> <mo>(</mo> <mi>I</mi> <mo>)</mo> <mo>)</mo> </mrow> <mrow> <mo>&prime;</mo> <mo>&prime;</mo> </mrow> </msup> <mo>|</mo> </mrow> <msup> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msup> <mrow> <mo>(</mo> <mrow> <msup> <mi>Pn</mi> <mo>&prime;</mo> </msup> <mrow> <mo>(</mo> <mi>I</mi> <mo>)</mo> </mrow> </mrow> <mo>)</mo> </mrow> <mo>&prime;</mo> </msup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mrow> <mn>3</mn> <mo>/</mo> <mn>2</mn> </mrow> </msup> </mfrac> </mrow>Wherein, K represent photoresponse function curvature function, (Pn'(I)) " represent photoresponse function second dervative function, (Pn' (I)) ' represent photoresponse function first derivative function.
- 6. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 3, It is characterized in that, carrying out optimizing to photoresponse function curvature curve using genetic algorithm, the maximum of points of curvature is obtained, in this base The corresponding illumination spot of the maximum of curvature is light limitation upper limit light intensity point after renormalization on plinth.
- 7. the construction method of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint according to claim 3, It is characterized in that, on the basis of obtaining all light limitation upper limit light intensity points, using temperature as independent variable, light limit upper limit light intensity be because Variable, fitting obtain the luminous environment goal of regulation and control value model under different temperatures.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711494114.6A CN108090693A (en) | 2017-12-31 | 2017-12-31 | The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711494114.6A CN108090693A (en) | 2017-12-31 | 2017-12-31 | The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108090693A true CN108090693A (en) | 2018-05-29 |
Family
ID=62181130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711494114.6A Pending CN108090693A (en) | 2017-12-31 | 2017-12-31 | The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108090693A (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102420A (en) * | 2018-09-05 | 2018-12-28 | 西北农林科技大学 | A kind of acquisition methods of the two-dimentional combined regulating target area based on regulation benefit priority |
CN110786208A (en) * | 2019-11-28 | 2020-02-14 | 贵州光合新植科技有限公司 | Plant planting technology evolution method based on artificial illumination planting technology |
CN112083748A (en) * | 2020-09-18 | 2020-12-15 | 西北农林科技大学 | Facility light environment regulation and control method with priority to light quality |
CN112697702A (en) * | 2020-12-31 | 2021-04-23 | 西北农林科技大学 | CO based on density distribution characteristics2Range finding method |
CN115191258A (en) * | 2022-08-16 | 2022-10-18 | 安徽农业大学 | Facility planting light supplementing system |
CN115226516A (en) * | 2022-07-07 | 2022-10-25 | 西北农林科技大学 | Facility light and carbon dioxide environment cooperative regulation and control method |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105389452A (en) * | 2015-12-31 | 2016-03-09 | 西北农林科技大学 | Cucumber whole-course photosynthetic rate prediction model based on neural network, and establishment method |
CN105678405A (en) * | 2015-12-31 | 2016-06-15 | 西北农林科技大学 | Modeling method for cucumber GA-RBF photosynthetic rate prediction model integrated with stomatal conductance |
WO2016110669A1 (en) * | 2015-01-09 | 2016-07-14 | Tom Robin Caine Boyde | Illumination for horticultural and other applications. |
US20170127622A1 (en) * | 2015-11-10 | 2017-05-11 | Xu Hong | Smart control/iot system for agriculture environment control |
CN107220672A (en) * | 2017-05-31 | 2017-09-29 | 西北农林科技大学 | Acquisition methods between a kind of suitable warm area based on crop demand |
CN107291126A (en) * | 2017-05-31 | 2017-10-24 | 西北农林科技大学 | A kind of facility light supplement control method and system based on crop demand |
CN107329511A (en) * | 2017-05-31 | 2017-11-07 | 西北农林科技大学 | Based on the vegetable aquaculture luminous environment high efficiency regulatory method and system between suitable root warm area |
-
2017
- 2017-12-31 CN CN201711494114.6A patent/CN108090693A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016110669A1 (en) * | 2015-01-09 | 2016-07-14 | Tom Robin Caine Boyde | Illumination for horticultural and other applications. |
US20170127622A1 (en) * | 2015-11-10 | 2017-05-11 | Xu Hong | Smart control/iot system for agriculture environment control |
CN105389452A (en) * | 2015-12-31 | 2016-03-09 | 西北农林科技大学 | Cucumber whole-course photosynthetic rate prediction model based on neural network, and establishment method |
CN105678405A (en) * | 2015-12-31 | 2016-06-15 | 西北农林科技大学 | Modeling method for cucumber GA-RBF photosynthetic rate prediction model integrated with stomatal conductance |
CN107220672A (en) * | 2017-05-31 | 2017-09-29 | 西北农林科技大学 | Acquisition methods between a kind of suitable warm area based on crop demand |
CN107291126A (en) * | 2017-05-31 | 2017-10-24 | 西北农林科技大学 | A kind of facility light supplement control method and system based on crop demand |
CN107329511A (en) * | 2017-05-31 | 2017-11-07 | 西北农林科技大学 | Based on the vegetable aquaculture luminous environment high efficiency regulatory method and system between suitable root warm area |
Non-Patent Citations (4)
Title |
---|
MIROSLAVA RAKOCEVIC等: "Estimating the canopy architecture and photosynthesis of coffea Arabica L. Plants cultivated under long-term elevated air CO2 concentration", 《IEEE》 * |
张海辉等: "黄瓜初花期光合速率主要影响因素分析与模型构建", 《农业机械学报》 * |
胡瑾: "基于作物光合需求的设施光环境调控方法与技术研究", 《中国博士学位论文全文数据库,农业科技辑》 * |
胡瑾等: "基于遗传算法的番茄幼苗光合作用优化调控模型", 《农业工程学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109102420A (en) * | 2018-09-05 | 2018-12-28 | 西北农林科技大学 | A kind of acquisition methods of the two-dimentional combined regulating target area based on regulation benefit priority |
CN110786208A (en) * | 2019-11-28 | 2020-02-14 | 贵州光合新植科技有限公司 | Plant planting technology evolution method based on artificial illumination planting technology |
CN112083748A (en) * | 2020-09-18 | 2020-12-15 | 西北农林科技大学 | Facility light environment regulation and control method with priority to light quality |
CN112083748B (en) * | 2020-09-18 | 2021-06-15 | 西北农林科技大学 | Facility light environment regulation and control method with priority to light quality |
CN112697702A (en) * | 2020-12-31 | 2021-04-23 | 西北农林科技大学 | CO based on density distribution characteristics2Range finding method |
CN115226516A (en) * | 2022-07-07 | 2022-10-25 | 西北农林科技大学 | Facility light and carbon dioxide environment cooperative regulation and control method |
CN115226516B (en) * | 2022-07-07 | 2023-09-26 | 西北农林科技大学 | Cooperative regulation and control method for facility light and carbon dioxide environment |
CN115191258A (en) * | 2022-08-16 | 2022-10-18 | 安徽农业大学 | Facility planting light supplementing system |
CN115191258B (en) * | 2022-08-16 | 2024-04-23 | 安徽农业大学 | Light filling system is planted to facility |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108090693A (en) | The structure of the Optimum Regulation model of the photosynthetic desired value of facility of fusion efficiencies constraint and application | |
McMurtrie et al. | Modifying existing forest growth models to take account of effects of elevated CO2 | |
CN108614601B (en) | Facility light environment regulation and control method fused with random forest algorithm | |
CN104656451B (en) | A kind of closed system envirment factor optimization regulating method based on crop modeling | |
CN108256181B (en) | Construction and application of optimal regulation and control model of facility carbon dioxide target value with integration of efficiency constraint | |
CN107329511B (en) | Based on the vegetable aquaculture luminous environment high efficiency regulatory method and system between suitable root warm area | |
CN101331840A (en) | Method for promoting plant grow using light source and light source products thereof | |
CN109613947A (en) | A kind of embedded facility luminous environment Optimum Regulation system merging illumination frequency and duty ratio | |
CN106842923A (en) | Greenhouse multiple-factor control method for coordinating based on plant physiology and energy optimization | |
CN111915062B (en) | Greenhouse crop water demand regulation and control method with water utilization rate and photosynthesis rate being coordinated | |
CN107390754A (en) | Intelligent plant growth environment adjustment system and method based on Internet of Things cloud platform | |
CN112763427B (en) | Crop growth and fertilization diagnosis simulation method coupled with remote sensing nitrogen information | |
CN108983849A (en) | It is a kind of to utilize compound extreme learning machine ANN Control greenhouse method | |
CN103605385A (en) | CO2 gas fertilizer fine regulation and control method and device used for solar greenhouse | |
CN105389452A (en) | Cucumber whole-course photosynthetic rate prediction model based on neural network, and establishment method | |
Wang et al. | Cotton growth model under drip irrigation with film mulching: A case study of Xinjiang, China | |
CN115310680A (en) | Tomato seedling model modeling and growth prediction method | |
CN114881805A (en) | Forestry carbon sink planning method based on carbon sink density | |
CN114371748B (en) | Greenhouse crop illumination carbon dioxide regulation and control method based on Gaussian curvature | |
CN113625806B (en) | Optimal-efficiency greenhouse plant regulation and control method for carbon neutralization demand | |
CN103444418B (en) | Based on the implementation method of the pipeline system plant factor of plant growth characteristics and rule | |
Wang et al. | Real-time monitoring of environmental information and modeling of the photosynthetic rate of tomato plants under greenhouse conditions | |
Liu et al. | Prediction model of photosynthetic rate based on SOPSO-LSSVM for regulation of greenhouse light environment | |
CN117032360A (en) | Plant factory light environment regulation and control system based on photosynthetic simulation and spectrum superposition | |
CN111667879A (en) | Moso bamboo forest stand structure optimization method based on minimum entropy |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180529 |
|
WD01 | Invention patent application deemed withdrawn after publication |