CN114237031A - Space ozone concentration control method - Google Patents

Space ozone concentration control method Download PDF

Info

Publication number
CN114237031A
CN114237031A CN202111483201.8A CN202111483201A CN114237031A CN 114237031 A CN114237031 A CN 114237031A CN 202111483201 A CN202111483201 A CN 202111483201A CN 114237031 A CN114237031 A CN 114237031A
Authority
CN
China
Prior art keywords
concentration
value
output
deviation
representing
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.)
Granted
Application number
CN202111483201.8A
Other languages
Chinese (zh)
Other versions
CN114237031B (en
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.)
Shanghai Kangjiu Disinfection Technology Co ltd
Original Assignee
Shanghai Kangjiu Disinfection Technology Co ltd
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 Shanghai Kangjiu Disinfection Technology Co ltd filed Critical Shanghai Kangjiu Disinfection Technology Co ltd
Priority to CN202111483201.8A priority Critical patent/CN114237031B/en
Publication of CN114237031A publication Critical patent/CN114237031A/en
Application granted granted Critical
Publication of CN114237031B publication Critical patent/CN114237031B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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)
  • Feedback Control In General (AREA)

Abstract

The invention discloses a method for controlling ozone concentration in a space, and belongs to the technical field of ozone control algorithms. The invention comprises the following steps: measuring static loss, measuring control margin, adjusting PID output upper calibration value, adjusting static loss, controlling actual output of PID algorithm and fuzzy control PID parameter; the control algorithm takes a self-adaptive algorithm as a core, allows the control margin to freely fluctuate to a certain degree, and improves the adaptive degree and the reaction time of the control algorithm to the environment; by using double closed-loop control as an execution algorithm, on one hand, the ozone fast response to the environmental ozone concentration is fast achieved, and on the other hand, the control degree of the output concentration is enhanced; through calculation and adjustment of static loss, the adaptation degree of a control algorithm to the environment is improved, and the accuracy degree of control of the ozone concentration is provided; controlling inaccurate parameters through a fuzzy control algorithm, and enhancing the robustness of the algorithm; through the arrangement of dead zone detection, the accuracy of the control algorithm on concentration control is improved.

Description

Space ozone concentration control method
Technical Field
The invention relates to the technical field of ozone control algorithms, in particular to a method for controlling the concentration of ozone in a space.
Background
Ozone is a strong oxidizing gas, widely used in medicine, agriculture, catering industry, sterilization and disinfection and the like, and has special chemical properties: in the air, the ozone decomposition speed is high when the ozone concentration is low, the ozone decomposition speed is slow along with the increase of the ozone concentration, and a plurality of ozone generation modes are provided, wherein a plurality of problems exist in the process of preparing ozone by corona discharge;
the ozone decomposition speed is related to the outlet concentration of an ozone generator, the outlet concentration of the ozone generator is closely related to the cleanliness of a gas source used by a machine, the cleanliness comprises oil content, air dew point and the like, the outlet concentration of the ozone generator is related to the cooling water temperature of the machine, the gas source temperature, the grid voltage of a power grid and the temperature of a space where the machine is located, the outlet concentration of the machine is a time variable, all related quantities cannot be collected and calculated when a controller is designed, and even if all related quantities are calculated, a listed equation is a transcendental equation, so that a result is difficult to solve;
the existing algorithm only collects the ozone concentration in the space when controlling the ozone concentration, which requires that the reaction speed of the control parameters is high, the sampled ozone concentration can quickly react when changing, and the output parameters are quickly updated, which results in that the PID execution time requires a short control period; however, the concentration of ozone is a large hysteresis quantity, and the parameter change can be sampled only after the adjustment of the pre-stage actuator, so that the condition frequently occurs, when the adjuster finds that the output is over-adjusted, the serious overshoot actually occurs, and when the underadjustment is sampled, the serious undershoot occurs; on the other hand, the half-life period of ozone in the air is inversely related to the current space concentration and is in a nonlinear relation; therefore, according to the traditional control method, the controlled variable can vibrate up and down in a large range and cannot meet the requirement;
the ozone concentration in the space follows a plurality of parameters, such as: changes such as environment temperature, environment humidity, environment cleanliness, ozone machine air outlet temperature, working pressure, air source temperature, air source humidity and cooling water temperature are correlated with one another, a specific model is difficult to establish, stable parameters are reluctantly adjusted by an existing control algorithm, and the control parameters can be seriously distorted along with time accumulation, so that control is out of control.
Disclosure of Invention
The present invention is directed to a method for controlling ozone concentration in a space, which solves the above-mentioned problems of the prior art.
In order to solve the technical problems, the invention provides the following technical scheme: the method for controlling the ozone concentration in the space comprises the following steps:
s100: measuring static loss;
s200: calculating a control margin;
s300: adjusting the PID output upper calibration value by controlling the PID output value calibration coefficient;
s400: static loss is adjusted up;
s500: fuzzy control PID parameters;
s600: controlling the actual output of the PID algorithm.
The static loss measurement includes: opening an air inlet valve, closing the output of a high-voltage power supply, resetting the ozone concentration in the environment, opening the output of the ozone high-voltage power supply when the ozone concentration is less than a measurement low value, setting the output of the ozone high-voltage power supply as a minimum output value, obtaining static loss by measuring the running trend of the ozone concentration, and outputting after averaging multiple times of static loss measurement;
ozone can act with the dust particles in the space, acts with the pipeline and objects in the environment, so that loss is caused by loss, air flow in the space or gas leakage, the size of the space is in direct proportion to the loss, the loss of the environment space is static loss, the static loss is almost invariable when the spatial arrangement is determined, when the input quantity of the ozone is static loss, the maximum value of the ozone concentration rise is not generated after the ozone is introduced into the environment, and the environment concentration does not rise obviously within 24 hours; when the input amount of ozone exceeds the static loss, the concentration of the environmental ozone has a stable rising trend;
the specific calculation method of the static loss measurement comprises the following steps:
s101: calculating the concentration difference per unit time:
Figure BDA0003396264920000021
wherein, cMeasureErrRepresents the difference in concentration, cPreRepresenting the current sample concentration, cLastRepresenting the concentration of the last sample per unit time;
s102: when the concentration difference is smaller than the lower limit of the concentration difference, if the lower counter value is smaller than the maximum lower counter value, the lower counter value is increased by one;
nMeasureLow＝nMeasureLow+1，nMeasureLow<NMeasureLowerand c isMeasureErr<CErrLower
Wherein n isMeasureLowRepresents the Down counter value, NMeasureLowerDenotes the maximum value of the count-down, CErrLowerDenotes the lower limit of the concentration difference, NMeasureLowerAnd CErrLowerIs a fixed threshold;
when the concentration difference is smaller than the lower limit of the concentration difference, if the lower counter value is larger than or equal to the lower counting maximum value, the current output concentration is judged to be too low, the trend counter is cleared, and the output of the high-voltage power supply is increased:
Figure BDA0003396264920000031
nMeasureLow≥NMeasureLowerand c isMeasureErr<CErrLower
Wherein n isMeasureUpRepresenting the up counter value, nMeasureZeroIndicating a non-changing counter value, nMeasureLowRepresents the down counter value, oPreIndicating the high voltage power supply output, OUpperThe step length of output power adjustment is represented, and P represents the control precision of the maximum power;
s103: when the concentration difference is greater than or equal to the lower concentration difference limit and the concentration difference is less than or equal to the upper concentration difference limit, if the value of the non-change counter is less than the maximum value of the non-change count, the value of the non-change counter is equal to the value of the non-change counter plus the value of an upper counter, the value of the upper counter is cleared, and the value of the non-change counter plus one;
Figure BDA0003396264920000032
nMeasureZero<NMeasureZeroand CErrLower≤cMeasureErr≤CErrUpper
Wherein N isMeasureZeroRepresenting the maximum value of the count without change, CErrUpperDenotes the upper limit of the concentration difference, NMeasureZeroAnd CErrUpperIs a fixed threshold;
when the concentration difference is more than or equal to the lower limit of the concentration difference, and the concentration difference is less than or equal to the upper limit of the concentration difference, if the value of the non-change counter is more than or equal to the maximum value of the non-change count, the current output concentration is judged to be too low, the trend counter is cleared, and the output of the high-voltage power supply is increased:
Figure BDA0003396264920000033
nMeasureZero≥NMeasureZeroand CErrLower≤cMeasureErr≤CErrUpper；
S104: when the concentration difference is greater than the concentration difference upper limit, if the up counter value is less than the up counter value maximum, then the up counter value is incremented by one:
nMeasureUp＝nMeasureUp+1，nMeasureUp<NMeasureUpperand c isMeasureErr>CErrUpper
Wherein N isMeasureUpperRepresents the upper limit of the count, CErrUpperDenotes the upper limit of the concentration difference, NMeasureUpperAnd CErrUpperIs a fixed threshold;
when the concentration difference is greater than the concentration difference upper limit, if the upper counter value is greater than or equal to the upper counting maximum value, the trend counter is reset, the static loss measured value of the time is output, the single static loss control bit is set, the next measurement period is waited, and the set single static loss control bit represents that the static loss of the time is completed:
Figure BDA0003396264920000041
nMeasureUp≥NMeasureUpper，cMeasureErr>CErrUpper
wherein n isMeasureZeroIndicating a non-changing counter value, nMeasureLowRepresents the Down counter value,/ZeroRepresents the static loss, oPreThe output of a high-voltage power supply is represented, the lossSingle represents a single static loss control bit, and the lossSingle is Boolean quantity;
s104: measuring the static loss for multiple times, calculating an average value, outputting the final static loss, and averaging the static loss measured for multiple times to obtain a stable static loss value so as to improve the accuracy of static measurement;
when the single static loss control position I is in a state that the counter value of the measuring times is smaller than the maximum value of the measuring count, the single static loss control position is set to zero, the next static loss measurement is continued, the counter value of the measuring is increased, and the static loss value is accumulated:
Figure BDA0003396264920000042
nZero<NZeroMax
wherein n isZeroRepresenting the value of the meter,/ZeroSumRepresents the total static loss, NZeroMaxRepresenting the maximum value of the measuring counter, NZeroMaxIs a fixed threshold;
when the single static loss control position is one, if the value of the measurement time counter is greater than or equal to the maximum value of the measurement count, the measurement count value is cleared, and the static loss is output:
Figure BDA0003396264920000043
nZero≥NZeroMax。
the trend counters include an up counter, a down counter, and a no change counter.
The accuracy of ozone output and the speed of ozone output need to be controlled, the PID output value is required to change in an interval, static loss is set as a static parameter of an output lower limit, the static loss can fluctuate due to factors such as air change, and a control algorithm adjusts the static loss according to actual output feedback;
the upper limit of the region cannot be directly obtained, an initial value is preset as the upper limit of the region, and the initial value is automatically fluctuated within a certain range according to the reaction speed and the control precision of the system in the actual control process;
setting a control margin as a fluctuation range of an upper limit of a PID output value, and when the environmental concentration is lower than a set value and the PID output value reaches the upper limit, the environmental concentration cannot be increased to the set value within a preset time, the control upper limit is too small and needs to be increased; when the environmental concentration is lower than a set value, the PID output value is reduced to a lower limit, and the concentration can not be reduced to be lower than the set value within a preset time, the control upper limit is over large, and the control upper limit needs to be reduced; meanwhile, when the system is started initially, the concentration of the system is slowly accumulated, the adjustment time limit value needs to be increased, and the normal adjustment time limit value is recovered when the system is operated normally;
the specific calculation mode of the control margin comprises the following steps:
s201: controlling and adjusting a time limit value according to the system starting time;
when the system is initially started, the adjustment time limit is increased:
Figure BDA0003396264920000051
FirstScal＝0
wherein n isUpMaxRepresents the maximum value of the up-regulation count, nLowMaxRepresenting the maximum value of the Down count, NUpFirstLimitDenotes the initial Up-regulation Limit, NLowFirstLimitDenotes the initial turndown limit, NUpFirstLimitAnd NLowFirstLimitThe first Scal is a fixed threshold value and represents the first adjustment of the control bit, the first Scal is a Boolean quantity, and the default value is 0;
and (3) restoring and adjusting the time limit value during normal operation:
Figure BDA0003396264920000052
FirstScal＝1
wherein N isUpLimitDenotes the upper regulation limit, NLowLimitIndicating the turndown limit, NUpLimitAnd NLowLimitIs a fixed threshold;
s202: calculating the deviation on control:
Figure BDA0003396264920000053
wherein, cUpErrTo representDeviation in control, cSetIndicating a set concentration, cPerRepresenting the concentration of the sample;
s203: when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is greater than the upper deviation limit value, the lower deviation detection is closed, when the upper deviation detection control position is zero, the up-regulation dead zone detection is started, and the up-regulation dead zone count value is cleared:
Figure BDA0003396264920000054
cUpErr>CUpErrMax，oPID>OPIDUpand c isPre<cSet
Wherein LowerMeasure represents lower deviation detection control bit, Uppermeasure represents upper deviation detection control bit, LowerMeasure and Uppermeasure are Boolean, nUpIndicates an up dead zone count value, CUpErrMaxRepresents the maximum value of the upper deviation, oPIDRepresents PID output, OPIDUpRepresents the upper limit of PID output, CUpErrMaxAnd OPIDUpIs a fixed threshold;
when the upper deviation detection control position is zero, the dead zone detection is not carried out at the moment; when the upper deviation detection control position is set, the dead zone detection is performed at the moment, the dead zone count value is effective, and the dead zone count value zero clearing operation is not continued;
when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is less than or equal to the upper deviation limit value, closing the upper deviation detection, and clearing the up-regulation dead zone count value:
Figure BDA0003396264920000061
cUpErr≤CUpErrMax，oPID>OPIDUpand c isPre<cSet；
S204: calculating the deviation under control, and calculating a formula:
Figure BDA0003396264920000062
wherein, cLowErrIndicating a deviation under control;
s205: when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is greater than the lower deviation limit value, the upper deviation detection is closed, and when the lower deviation detection control position is zero, the lower dead zone detection is started, and the lower dead zone count value is cleared:
Figure BDA0003396264920000063
cLowErr>CLowErrMax，oPID≤OPIDLowand c isPre>cSet
Wherein n isLowIndicating a down dead zone count value, CLowErrMaxRepresents the maximum value of the lower deviation, OPIDLowRepresents the lower limit of PID output, CLowErrMaxAnd OPIDLowIs a fixed threshold;
when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is less than or equal to the lower deviation limit value, the lower deviation detection is closed, and the down-regulation dead zone count value is cleared:
Figure BDA0003396264920000064
cLowErr≤CLowErrMax，oPID≤OPIDLowand c isPre>cSet；
S206: when the sampling concentration is greater than the set concentration, if the PID output is greater than the output lower limit, and the PID output is less than or equal to the output upper limit, closing the upper limit and the lower limit for detection:
Figure BDA0003396264920000071
OPIDLow<oPID≤OPIDUpand c isPre>cSet；
S207: when the upper limit detection is allowed, the upper limit detection is counted, after the maximum value is counted, an upper deviation detection control bit is set to zero, and a first adjustment control bit is set:
nUp＝nUp+1，nUp<nUpMaxand UpperMeasure ═ 1
Figure BDA0003396264920000072
nUp≥nUpMaxAnd UpperMeasure ═ 1
Wherein n isUpIndicates the up dead zone count value, nUpMaxThe maximum value of the up-regulation count, wherein Upperenable represents an up-regulation output upper limit allowable control bit;
s208: when the lower limit detection is allowed, counting the lower limit detection, and setting a down control bit after the maximum value is counted:
nLow＝nLow+1，nLow<nLowMaxand LowerMeasure is 1
Figure BDA0003396264920000073
nLow≥NLowMaxAnd LowerMeasure is 1
Wherein, LowerEnable represents a lower output upper limit permission control bit.
And setting a dead zone detection mechanism, carrying out time delay detection to prevent data from shaking, if a dead zone detection area is started, controlling the upper deviation to continuously deviate upwards or controlling the lower deviation to continuously deviate downwards, indicating the continuous upper deviation or lower deviation of the concentration, judging that the control margin is too large or too small if no fluctuation occurs, and adjusting the PID to output an upper calibration value.
And adjusting the PID output upper calibration value by controlling the PID output value calibration coefficient, wherein the specific algorithm of the PID output upper calibration value comprises the following steps:
s301: when the first adjustment output upper limit allowable control position is located, resetting the upper adjustment output upper limit allowable control position, and performing up-adjustment on the PID output value calibration coefficient:
Figure BDA0003396264920000074
UpperEnable＝1
wherein v isScalIndicating PID output calibration, VUpperPrecisionIndicating PID output calibration up-regulation step, VPWMMAXIndicating a PID output calibration upper limit value;
s302: resetting the lower-regulation output upper limit allowable control bit when the lower-regulation output upper limit allowable control position is one, and calibrating the PID output value to be regulated downwards if the PID output value is greater than the output lower limit; and if the PID output value is less than or equal to the output lower limit and the static loss is greater than the static loss minimum value, the static loss is adjusted downwards:
Figure BDA0003396264920000081
LowerEnable＝1
wherein o isPIDRepresents PID output, OPIDLowIndicating the lower limit of PID output, VLowerPrecisionIndicating PID output calibration down-regulation step, VScalLowIndicating the lower limit of the PID output calibration,/ZeroRepresents the static loss, LLowerPrecisionRepresenting the static loss down-regulation step, LZeroLowerRepresenting the static loss minimum.
When the concentration rises too fast, the current concentration exceeds a set value and is difficult to drop below a lower limit in preset time, so that the output lower limit is too large, namely the static loss is too large, and the static loss is adjusted downwards; when the output power exceeds the set output power and the concentration cannot reach the set value within the preset time, the static loss is over small, and the static loss is adjusted upwards at the moment.
The specific algorithm for static loss up-regulation comprises the following steps:
s401: when the sampling concentration is greater than the set concentration, starting concentration maximum value measurement:
MeasureMaxEnable＝1,cPre>cSet
wherein MeasureMeasureMaxEnable represents a maximum value measurement control bit, MeasureMeasureMaxEnable is a Boolean quantity, cPreRepresenting the concentration of the sample, cSetRepresents the set concentration;
s402: when the measurement of the maximum concentration value is started, if the maximum value of single sampling is less than the sampling concentration, the maximum value is updated to be equal to the sampling concentration:
cMeasureMax＝cPre,cMeasureMax<cPre，MeasureMaxEnable＝1
wherein, cMeasureMaxRepresenting the maximum of a single sample
S403: when the sampling concentration is less than the set concentration and the concentration maximum value is measured and opened, the maximum value is closed for measurement, the concentration maximum value is output, the single sampling maximum value is emptied, and the single sampling minimum value is opened for measurement:
Figure BDA0003396264920000082
measurmaxenable 1 and cPre<cSet
Wherein, cPreMaxDenotes the maximum concentration, cMeasureMaxRepresenting the maximum value of single sampling, wherein MeasureMinEnable represents the minimum measurement control bit and is Boolean;
s404: when the measurement of the concentration minimum value is started, if the single sampling minimum value is greater than the sampling concentration, the updated minimum value is equal to the sampling concentration:
cMeasureMin＝cPre,cMeasureMin>cPreand measureinEnable 1
Wherein, cMeasureMinIndicating the minimum value of single sampling;
s405: when the measurement of concentration minimum is opened, if the sampling concentration is greater than the set concentration, the measurement of maximum and minimum is closed, the single sampling minimum is output, and the single sampling minimum is set to be equal to 1000:
Figure BDA0003396264920000091
measureinenable 1 and cPre>cSet
Wherein, cPreMinRepresents the minimum concentration;
s406: when the concentration maximum value and the concentration minimum value are obtained, the root mean square of the concentration extreme value difference is calculated:
Figure BDA0003396264920000092
nRMS<NRMSMax
Figure BDA0003396264920000093
nRMS≥NRMSMax
wherein n isRMSRepresents the number of extreme samples, NRMSMaxRepresents the upper limit of the extremum sample, NRMSMaxTo fix the threshold value, cRMSSumRepresents the sum of extrema samples, cRMSRepresenting an extreme root mean square;
s407: and comparing the extreme value difference root-mean-square with a set difference value, and if the extreme value difference root-mean-square is larger than a limit value, adjusting the static loss value upwards.
Figure BDA0003396264920000094
Wherein lZeroRepresents the running static loss, LUpperPrecisionRepresents the static loss up-regulation step, CRMSRepresenting an extreme root mean square setting, LZeroUpperRepresents the maximum allowable value of static loss, LUpperPrecision、CRMSAnd LZeroUpperThe threshold is fixed.
The fuzzy control PID parameter specific algorithm comprises the following steps:
s501: sampling control deviation and deviation change rate:
Figure BDA0003396264920000101
wherein, cErrRepresents the difference in concentration, cPreRepresenting the concentration of the sample, cSetIndicates the set concentration, eEIndicating concentration calibrationDeviation, EErrIndicating the degree of deviation of the calibration, eEcIndicating the rate of change of the deviation of the calibration of the concentration, cErrLastIndicating the limit of the calibration deviation, EEcRepresenting the limit of the rate of change of the calibration deviation, cErrShowing the historical concentration difference, EErr、EEcIs a fixed threshold;
s502: converting the control interval of the concentration calibration deviation and the change rate of the concentration calibration deviation to-3 to 3, dispersing the continuous number into eight numbers, and then calibrating the data to an integer of 1 to 7:
Figure BDA0003396264920000102
Figure BDA0003396264920000103
Figure BDA0003396264920000104
Figure BDA0003396264920000105
wherein e isFuzzyRepresenting the fuzzified output value, eltRepresenting a left value formatting output value, ertRepresenting the right value formatted output value, x, y, z representing the blurred input value, e1Representing the first column of the fuzzy matrix, e2Representing the second index of the fuzzy matrix, eC1Representing the first row index of the fuzzy matrix, eC2Representing the second row index of the fuzzy matrix, eE1Representing the first row and column elements of the matrix M, eE2Representing the first row and the second column of elements of the matrix M, eEc1Representing the first row and column elements of the matrix N, eEc2Representing the second row and the first column of elements of the matrix M;
s503: according to the influence strength of the concentration calibration deviation and the change rate of the concentration calibration deviation on proportion, integration and differentiation, three 7 multiplied by 7 fuzzy matrixes P are made7×7、I7×7、D7×7Calculating the row and column values of variable coordinates of the fuzzy matrix according to the concentration parameters, and deriving three new fuzzy matrices A2×2、B2×2、C2×2：
Figure BDA0003396264920000111
Wherein
Figure BDA0003396264920000112
Figure BDA0003396264920000113
Wherein
Figure BDA0003396264920000114
Figure BDA0003396264920000115
Wherein
Figure BDA0003396264920000116
Wherein, P7×7A fuzzy matrix representing the ratio coefficients of the deviation of the calibration of the concentration and the rate of change of the deviation of the calibration of the concentration, I7×7Fuzzy matrices representing the deviation of the calibration of the concentration and the rate of change of the deviation of the calibration of the concentration versus the integral coefficient, D7×7A fuzzy matrix representing the deviation of the concentration calibration and the rate of change of the deviation of the concentration calibration versus the differential coefficient, A2×2Fuzzy matrix representing the ratio of the matrix P to the concentration calibration deviation and the rate of change of the concentration calibration deviation, B2×2A fuzzy matrix representing the ratio of the matrix I to the concentration calibration deviation and the rate of change of the concentration calibration deviation, C2×2Fuzzy matrix, p, representing the ratio of the matrix D to the concentration calibration deviation and the rate of change of the concentration calibration deviationi,j、ii,j、di,jRespectively representing elements, a, of a matrix P, I, Di,j、bi,j、ci,jRespectively represent elements of matrix A, B, C;
s504: constructing a concentration calibration deviation fuzzification matrix M and a concentration calibration deviation change rate fuzzification matrix N:
M1×2＝[eE1 eE2]
N1×2＝[eEc1 eEc2]T
wherein M is1×2Indicating the concentration calibration deviation matrix, N1×2Representing a concentration calibration deviation change rate matrix;
s505: determining a fuzzy PID parameter (k)p,tI,td)：
Figure BDA0003396264920000117
Figure BDA0003396264920000118
Figure BDA0003396264920000119
Wherein k isPDenotes the proportionality coefficient, tIRepresenting the integral coefficient, tDDenotes the differential coefficient, KPRepresenting a scale factor setting parameter, TIRepresenting an integral coefficient setting parameter, TDRepresenting a differential coefficient setting parameter, KP、TIAnd TDIs a fixed threshold.
The actual output of the control PID algorithm is calibrated and limited through static loss and PID output:
Figure BDA0003396264920000121
wherein o isPIDSumRepresenting calculated values of the PID algorithm, i.e. output values, k, without nominal clippingPDenotes the proportionality coefficient, PWRepresenting the proportional weight, tIRepresenting the integral coefficient, tDRepresenting a differential coefficient, s representing the Laplace operator, DDelayRepresenting a differential delay systemNumber, DWRepresenting a differential weight, cSetIndicates a set value, cPreRepresenting the current sample concentration, oPIDRepresents the PID algorithm output, oOutRepresenting the actual power control value, l, delivered to the power supplyZeroRepresenting static losses.
The maximum output power is set to ten thousand, the step length is set to ten thousandth, and the precision of the control algorithm is improved.
Compared with the prior art, the invention has the following beneficial effects: the control algorithm only collects the ozone concentration in the space, avoids controlling the ozone concentration through calculating the concentration at the outlet of a machine, reduces the calculated amount, improves the calculation cost and the reaction time of the control algorithm, and avoids data deviation caused by slow reaction of the control algorithm on the ozone concentration of the environment due to calculation; the control algorithm takes a self-adaptive algorithm as a core, allows the control margin to freely fluctuate to a certain degree, and improves the adaptive degree and the reaction time of the control algorithm to the environment; by using double closed-loop control as an execution algorithm, on one hand, the ozone fast response to the environmental ozone concentration is fast achieved, and on the other hand, the control degree of the output concentration is enhanced; through calculation and adjustment of static loss, the adaptation degree of a control algorithm to the environment is improved, and the accuracy degree of control of the ozone concentration is provided; controlling inaccurate parameters through a fuzzy control algorithm, and enhancing the robustness of the algorithm; through the arrangement of dead zone detection, excessive reaction of a control algorithm due to data fluctuation is avoided, and the accuracy of the control algorithm on concentration control is improved; the control algorithm has high precision and adjustable precision.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a schematic flow chart of a method for controlling ozone concentration in a space according to the present invention;
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Referring to fig. 1, the present invention provides a technical solution:
the first embodiment is as follows: the method for controlling the ozone concentration in the space comprises the following steps:
s100: measuring static loss;
s200: calculating a control margin;
s300: adjusting the PID output upper calibration value by controlling the PID output value calibration coefficient;
s400: static loss is adjusted up;
s500: fuzzy control PID parameters;
s600: controlling the actual output of the PID algorithm.
The static loss measurement includes: opening an air inlet valve, closing the output of a high-voltage power supply, resetting the ozone concentration in the environment, opening the output of the ozone high-voltage power supply when the ozone concentration is less than a measurement low value, setting the output of the ozone high-voltage power supply as a minimum output value, obtaining static loss by measuring the running trend of the ozone concentration, and outputting after averaging multiple times of static loss measurement;
ozone can act with the dust particles in the space, acts with the pipeline and objects in the environment, so that loss is caused by loss, air flow in the space or gas leakage, the size of the space is in direct proportion to the loss, the loss of the environment space is static loss, the static loss is almost invariable when the spatial arrangement is determined, when the input quantity of the ozone is static loss, the maximum value of the ozone concentration rise is not generated after the ozone is introduced into the environment, and the environment concentration does not rise obviously within 24 hours; when the input amount of ozone exceeds the static loss, the concentration of the environmental ozone has a stable rising trend;
the specific calculation method of the static loss measurement comprises the following steps:
s101: calculating the concentration difference per unit time:
Figure BDA0003396264920000131
wherein, cMeasureErrRepresents the difference in concentration, cPreRepresenting the current sample concentration, cLastRepresenting the concentration of the last sample per unit time;
s102: when the concentration difference is smaller than the lower limit of the concentration difference, if the lower counter value is smaller than the maximum lower counter value, the lower counter value is increased by one;
nMeasureLow＝nMeasureLow+1，nMeasureLow<NMeasureLowerand c isMeasureErr<CErrLower
Wherein n isMeasureLowRepresents the Down counter value, NMeasureLowerDenotes the maximum value of the count-down, CErrLowerDenotes the lower limit of the concentration difference, NMeasureLowerAnd CErrLowerIs a fixed threshold;
when the concentration difference is smaller than the lower limit of the concentration difference, if the lower counter value is larger than or equal to the lower counting maximum value, the current output concentration is judged to be too low, the trend counter is cleared, and the output of the high-voltage power supply is increased:
Figure BDA0003396264920000141
nMeasureLow≥NMeasureLowerand c isMeasureErr<CErrLower
Wherein n isMeasureUpRepresenting the up counter value, nMeasureZeroIndicating a non-changing counter value, nMeasureLowRepresents the down counter value, oPreIndicating the high voltage power supply output, OUpperThe step length of output power adjustment is represented, and P represents the control precision of the maximum power;
s103: when the concentration difference is greater than or equal to the lower concentration difference limit and the concentration difference is less than or equal to the upper concentration difference limit, if the value of the non-change counter is less than the maximum value of the non-change count, the value of the non-change counter is equal to the value of the non-change counter plus the value of an upper counter, the value of the upper counter is cleared, and the value of the non-change counter plus one;
Figure BDA0003396264920000142
nMeasureZero<NMeasureZeroand CErrLower≤cMeasureErr≤CErrUpper
Wherein N isMeasureZeroRepresenting the maximum value of the count without change, CErrUpperDenotes the upper limit of the concentration difference, NMeasureZeroAnd CErrUpperIs a fixed threshold;
when the concentration difference is more than or equal to the lower limit of the concentration difference, and the concentration difference is less than or equal to the upper limit of the concentration difference, if the value of the non-change counter is more than or equal to the maximum value of the non-change count, the current output concentration is judged to be too low, the trend counter is cleared, and the output of the high-voltage power supply is increased:
Figure BDA0003396264920000143
nMeasureZero≥NMeasureZeroand CErrLower≤cMeasureErr≤CErrUpper；
S104: when the concentration difference is greater than the concentration difference upper limit, if the up counter value is less than the up counter value maximum, then the up counter value is incremented by one:
nMeasureUp＝nMeasureUp+1，nMeasureUp<NMeasureUpperand c isMeasureErr>CErrUpper
Wherein N isMeasureUpperRepresents the upper limit of the count, CErrUpperDenotes the upper limit of the concentration difference, NMeasureUpperAnd CErrUpperIs a fixed threshold;
when the concentration difference is greater than the concentration difference upper limit, if the upper counter value is greater than or equal to the upper counting maximum value, the trend counter is reset, the static loss measured value of the time is output, the single static loss control bit is set, the next measurement period is waited, and the set single static loss control bit represents that the static loss of the time is completed:
Figure BDA0003396264920000151
nMeasureUp≥NMeasureUpper，cMeasureErr>CErrUpper
wherein n isMeasureZeroIndicating a non-changing counter value, nMeasureLowRepresents the Down counter value,/ZeroRepresents the static loss, oPreThe output of a high-voltage power supply is represented, the lossSingle represents a single static loss control bit, and the lossSingle is Boolean quantity;
s104: measuring the static loss for multiple times, calculating an average value, outputting the final static loss, and averaging the static loss measured for multiple times to obtain a stable static loss value so as to improve the accuracy of static measurement;
when the single static loss control position I is in a state that the counter value of the measuring times is smaller than the maximum value of the measuring count, the single static loss control position is set to zero, the next static loss measurement is continued, the counter value of the measuring is increased, and the static loss value is accumulated:
Figure BDA0003396264920000152
nZero<NZeroMax
wherein n isZeroRepresenting the value of the meter,/ZeroSumRepresents the total static loss, NZeroMaxRepresenting the maximum value of the measuring counter, NZeroMaxIs a fixed threshold;
when the single static loss control position is one, if the value of the measurement time counter is greater than or equal to the maximum value of the measurement count, the measurement count value is cleared, and the static loss is output:
Figure BDA0003396264920000153
nZero≥NZeroMax。
the trend counters include an up counter, a down counter, and a no change counter.
The accuracy of ozone output and the speed of ozone output need to be controlled, the PID output value is required to change in an interval, static loss is set as a static parameter of an output lower limit, the static loss can fluctuate due to factors such as air change, and a control algorithm adjusts the static loss according to actual output feedback;
the upper limit of the area can not be directly obtained, an initial value is preset as the upper limit of the area, and the initial value automatically fluctuates within a certain range according to the reaction speed and the control precision of the system in the actual control process;
setting a control margin as a fluctuation range of an upper limit of a PID output value, and when the environmental concentration is lower than a set value and the PID output value reaches the upper limit, the environmental concentration cannot be increased to the set value within a preset time, the control upper limit is too small and needs to be increased; when the environmental concentration is lower than a set value, the PID output value is reduced to a lower limit, and the concentration can not be reduced to be lower than the set value within a preset time, the control upper limit is over large, and the control upper limit needs to be reduced; meanwhile, when the system is started initially, the concentration of the system is slowly accumulated, the adjustment time limit value needs to be increased, and the normal adjustment time limit value is recovered when the system is operated normally;
the specific calculation mode of the control margin comprises the following steps:
s201: controlling and adjusting a time limit value according to the system starting time;
when the system is initially started, the adjustment time limit is increased:
Figure BDA0003396264920000161
FirstScal＝0
wherein n isUpMaxRepresents the maximum value of the up-regulation count, nLowMaxRepresenting the maximum value of the Down count, NUpFirstLimitDenotes the initial Up-regulation Limit, NLowFirstLimitDenotes the initial turndown limit, NUpFirstLimitAnd NLowFirstLimitThe first Scal is a fixed threshold value and represents the first adjustment of the control bit, the first Scal is a Boolean quantity, and the default value is 0;
and (3) restoring and adjusting the time limit value during normal operation:
Figure BDA0003396264920000162
FirstScal＝1
wherein N isUpLimitDenotes the upper regulation limit, NLowLimitIndicating the turndown limit, NUpLimitAnd NLowLimitIs a fixed threshold;
s202: calculating the deviation on control:
Figure BDA0003396264920000163
wherein, cUpErrIndicating deviation in control, cSetIndicating a set concentration, cPerRepresenting the concentration of the sample;
s203: when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is greater than the upper deviation limit value, the lower deviation detection is closed, when the upper deviation detection control position is zero, the up-regulation dead zone detection is started, and the up-regulation dead zone count value is cleared:
Figure BDA0003396264920000164
cUpErr>CUpErrMax，oPID>OPIDUpand c isPre<cSet
Wherein LowerMeasure represents lower deviation detection control bit, Uppermeasure represents upper deviation detection control bit, LowerMeasure and Uppermeasure are Boolean, nUpIndicates an up dead zone count value, CUpErrMaxRepresents the maximum value of the upper deviation, oPIDRepresents PID output, OPIDUpRepresents the upper limit of PID output, CUpErrMaxAnd OPIDUpIs a fixed threshold;
when the upper deviation detection control position is zero, the dead zone detection is not carried out at the moment; when the upper deviation detection control position is set, the dead zone detection is performed at the moment, the dead zone count value is effective, and the dead zone count value zero clearing operation is not continued;
when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is less than or equal to the upper deviation limit value, closing the upper deviation detection, and clearing the up-regulation dead zone count value:
Figure BDA0003396264920000171
cUpErr≤CUpErrMax，oPID>OPIDUpand c isPre<cSet；
S204: calculating the deviation under control, and calculating a formula:
Figure BDA0003396264920000172
wherein, cLowErrIndicating a deviation under control;
s205: when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is greater than the lower deviation limit value, the upper deviation detection is closed, and when the lower deviation detection control position is zero, the lower dead zone detection is started, and the lower dead zone count value is cleared:
Figure BDA0003396264920000173
cLowErr>CLowErrMax，oPID≤OPIDLowand c isPre>cSet
Wherein n isLowIndicating a down dead zone count value, CLowErrMaxRepresents the maximum value of the lower deviation, OPIDLowRepresents the lower limit of PID output, CLowErrMaxAnd OPIDLowIs a fixed threshold;
when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is less than or equal to the lower deviation limit value, the lower deviation detection is closed, and the down-regulation dead zone count value is cleared:
Figure BDA0003396264920000174
cLowErr≤CLowErrMax，oPID≤OPIDLowand c isPre>cSet；
S206: when the sampling concentration is greater than the set concentration, if the PID output is greater than the output lower limit, and the PID output is less than or equal to the output upper limit, closing the upper limit and the lower limit for detection:
Figure BDA0003396264920000181
OPIDLow<oPID≤OPIDUpand c isPre>cSet；
S207: when the upper limit detection is allowed, the upper limit detection is counted, after the maximum value is counted, an upper deviation detection control bit is set to zero, and a first adjustment control bit is set:
nUp＝nUp+1，nUp<nUpMaxand UpperMeasure ═ 1
Figure BDA0003396264920000182
nUp≥nUpMaxAnd UpperMeasure ═ 1
Wherein n isUpIndicates the up dead zone count value, nUpMaxThe maximum value of the up-regulation count, wherein Upperenable represents an up-regulation output upper limit allowable control bit;
s208: when the lower limit detection is allowed, counting the lower limit detection, and setting a down control bit after the maximum value is counted:
nLow＝nLow+1，nLow<nLowMaxand LowerMeasure is 1
Figure BDA0003396264920000183
nLow≥NLowMaxAnd LowerMeasure is 1
Wherein, LowerEnable represents a lower output upper limit permission control bit.
And setting a dead zone detection mechanism, carrying out time delay detection to prevent data from shaking, if a dead zone detection area is started, controlling the upper deviation to continuously deviate upwards or controlling the lower deviation to continuously deviate downwards, indicating the continuous upper deviation or lower deviation of the concentration, judging that the control margin is too large or too small if no fluctuation occurs, and adjusting the PID to output an upper calibration value.
The PID output calibration value is adjusted by controlling the PID output value calibration coefficient, and the specific algorithm of the PID output calibration value comprises the following steps:
s301: when the first adjustment output upper limit allowable control position is located, resetting the upper adjustment output upper limit allowable control position, and performing up-adjustment on the PID output value calibration coefficient:
Figure BDA0003396264920000184
UpperEnable＝1
wherein v isScalIndicating PID output calibration, VUpperPrecisionIndicating PID output calibration up-regulation step, VPWMMAXIndicating a PID output calibration upper limit value;
s302: resetting the lower-regulation output upper limit allowable control bit when the lower-regulation output upper limit allowable control position is one, and calibrating the PID output value to be regulated downwards if the PID output value is greater than the output lower limit; and if the PID output value is less than or equal to the output lower limit and the static loss is greater than the static loss minimum value, the static loss is adjusted downwards:
Figure BDA0003396264920000191
LowerEnable＝1
wherein o isPIDRepresents PID output, OPIDLowIndicating the lower limit of PID output, VLowerPrecisionIndicating PID output calibration down-regulation step, VScalLowIndicating the lower limit of the PID output calibration,/ZeroRepresents the static loss, LLowerPrecisionRepresenting the static loss down-regulation step, LZeroLowerRepresenting the static loss minimum.
When the concentration rises too fast, the current concentration exceeds a set value and is difficult to drop below a lower limit in preset time, so that the output lower limit is too large, namely the static loss is too large, and the static loss is adjusted downwards; when the output power exceeds the set output power and the concentration cannot reach the set value within the preset time, the static loss is over small, and the static loss is adjusted upwards at the moment.
The specific algorithm for static loss up-regulation comprises the following steps:
s401: when the sampling concentration is greater than the set concentration, starting concentration maximum value measurement:
MeasureMaxEnable＝1,cPre>cSet
wherein MeasureMeasureMaxEnable represents a maximum value measurement control bit, MeasureMeasureMeasureMaxEnable is a Boolean quantity, cPreRepresenting the concentration of the sample, cSetRepresents the set concentration;
s402: when the measurement of the maximum concentration value is started, if the maximum value of single sampling is less than the sampling concentration, the maximum value is updated to be equal to the sampling concentration:
cMeasureMax＝cPre,cMeasureMax<cPre，MeasureMaxEnable＝1
wherein, cMeasureMaxRepresenting the maximum of a single sample
S403: when the sampling concentration is less than the set concentration and the concentration maximum value is measured and opened, the maximum value is closed for measurement, the concentration maximum value is output, the single sampling maximum value is emptied, and the single sampling minimum value is opened for measurement:
Figure BDA0003396264920000192
measurmaxenable 1 and cPre<cSet
Wherein, cPreMaxDenotes the maximum concentration, cMeasureMaxRepresenting the maximum value of single sampling, wherein MeasureMinEnable represents the minimum measurement control bit and is Boolean;
s404: when the measurement of the concentration minimum value is started, if the single sampling minimum value is greater than the sampling concentration, the updated minimum value is equal to the sampling concentration:
cMeasureMin＝cPre,cMeasureMin>cPreand measureinEnable 1
Wherein, cMeasureMinIndicating the minimum value of single sampling;
s405: when the measurement of concentration minimum is opened, if the sampling concentration is greater than the set concentration, the measurement of maximum and minimum is closed, the single sampling minimum is output, and the single sampling minimum is set to be equal to 1000:
Figure BDA0003396264920000201
measureinenable 1 and cPre>cSet
Wherein, cPreMinRepresents the minimum concentration;
s406: when the concentration maximum value and the concentration minimum value are obtained, the root mean square of the concentration extreme value difference is calculated:
Figure BDA0003396264920000202
nRMS<NRMSMax
Figure BDA0003396264920000203
nRMS≥NRMSMax
wherein n isRMSRepresents the number of extreme samples, NRMSMaxRepresents the upper limit of the extremum sample, NRMSMaxTo fix the threshold value, cRMSSumRepresents the sum of extrema samples, cRMSRepresenting an extreme root mean square;
s407: and comparing the extreme value difference root-mean-square with a set difference value, and if the extreme value difference root-mean-square is larger than a limit value, adjusting the static loss value upwards.
Figure BDA0003396264920000204
Wherein lZeroRepresents the running static loss, LUpperPrecisionRepresents the static loss up-regulation step, CRMSRepresenting an extreme root mean square setting, LZeroUpperRepresents the maximum allowable value of static loss, LUpperPrecision、CRMSAnd LZeroUpperThe threshold is fixed.
The fuzzy control PID parameter specific algorithm comprises the following steps:
s501: sampling control deviation and deviation change rate:
Figure BDA0003396264920000205
wherein, cErrRepresents the difference in concentration, cPreRepresenting the concentration of the sample, cSetIndicates the set concentration, eEIndicating deviation of calibration of concentration, EErrIndicating the degree of deviation of the calibration, eEcIndicating the rate of change of the deviation of the calibration of the concentration, cErrLastIndicating the limit of the calibration deviation, EEcRepresenting the limit of the rate of change of the calibration deviation, cErrShowing the historical concentration difference, EErr、EEcIs a fixed threshold;
s502: converting the control interval of the concentration calibration deviation and the change rate of the concentration calibration deviation to-3 to 3, dispersing the continuous number into eight numbers, and then calibrating the data to an integer of 1 to 7:
Figure BDA0003396264920000211
Figure BDA0003396264920000212
Figure BDA0003396264920000213
Figure BDA0003396264920000214
wherein e isFuzzyThe output value of the fuzzification is represented,eltrepresenting a left value formatting output value, ertRepresenting the right value formatted output value, x, y, z representing the blurred input value, e1Representing the first column of the fuzzy matrix, e2Representing the second index of the fuzzy matrix, eC1Representing the first row index of the fuzzy matrix, eC2Representing the second row index of the fuzzy matrix, eE1Representing the first row and column elements of the matrix M, eE2Representing the first row and the second column of elements of the matrix M, eEc1Representing the first row and column elements of the matrix N, eEc2Representing the second row and the first column of elements of the matrix M;
s503: according to the influence strength of the concentration calibration deviation and the change rate of the concentration calibration deviation on proportion, integration and differentiation, three 7 multiplied by 7 fuzzy matrixes P are made7×7、I7×7、D7×7Calculating the row and column values of variable coordinates of the fuzzy matrix according to the concentration parameters, and deriving three new fuzzy matrices A2×2、B2×2、C2×2：
Figure BDA0003396264920000215
Wherein
Figure BDA0003396264920000216
Figure BDA0003396264920000221
Wherein
Figure BDA0003396264920000222
Figure BDA0003396264920000223
Wherein
Figure BDA0003396264920000224
Wherein, P7×7A fuzzy matrix representing the ratio coefficients of the deviation of the calibration of the concentration and the rate of change of the deviation of the calibration of the concentration, I7×7Indicating concentration calibration deviation and concentrationFuzzy matrices of degree-scale deviation change rate versus integral coefficient, D7×7A fuzzy matrix representing the deviation of the concentration calibration and the rate of change of the deviation of the concentration calibration versus the differential coefficient, A2×2Fuzzy matrix representing the ratio of the matrix P to the concentration calibration deviation and the rate of change of the concentration calibration deviation, B2×2A fuzzy matrix representing the ratio of the matrix I to the concentration calibration deviation and the rate of change of the concentration calibration deviation, C2×2Fuzzy matrix, p, representing the ratio of the matrix D to the concentration calibration deviation and the rate of change of the concentration calibration deviationi,j、ii,j、di,jRespectively representing elements, a, of a matrix P, I, Di,j、bi,j、ci,jRespectively represent elements of matrix A, B, C;
s504: constructing a concentration calibration deviation fuzzification matrix M and a concentration calibration deviation change rate fuzzification matrix N:
M1×2＝[eE1 eE2]
N1×2＝[eEc1 eEc2]T
wherein M is1×2Indicating the concentration calibration deviation matrix, N1×2Representing a concentration calibration deviation change rate matrix;
s505: determining a fuzzy PID parameter (k)p,tI,td)：
Figure BDA0003396264920000225
Figure BDA0003396264920000226
Figure BDA0003396264920000227
Wherein k isPDenotes the proportionality coefficient, tIRepresenting the integral coefficient, tDDenotes the differential coefficient, KPRepresenting a scale factor setting parameter, TIRepresenting an integral coefficient setting parameter, TDRepresents a differential systemNumber setting parameter, KP、TIAnd TDIs a fixed threshold.
And controlling the actual output of the PID algorithm to calibrate amplitude limiting through static loss and PID output:
Figure BDA0003396264920000228
wherein o isPIDSumRepresenting calculated values of the PID algorithm, i.e. output values, k, without nominal clippingPDenotes the proportionality coefficient, PWRepresenting the proportional weight, tIRepresenting the integral coefficient, tDRepresenting a differential coefficient, s representing the Laplace operator, DDelayRepresenting the differential delay coefficient, DWRepresenting a differential weight, cSetIndicates a set value, cPreRepresenting the current sample concentration, oPIDRepresents the PID algorithm output, oOutRepresenting the actual power control value, l, delivered to the power supplyZeroRepresenting static losses.
The maximum output power precision is adjustable, the satisfactory control precision can be obtained by combining a control algorithm, the maximum output power is set to ten thousand, the step length is set to one ten thousand, and the control precision is less than 10 ppb.
Compared with the prior art, the invention has the following beneficial effects: the control algorithm only collects the ozone concentration in the space, avoids controlling the ozone concentration through calculating the concentration at the outlet of a machine, reduces the calculated amount, improves the calculation cost and the reaction time of the control algorithm, and avoids data deviation caused by slow reaction of the control algorithm on the ozone concentration of the environment due to calculation; the control algorithm takes a self-adaptive algorithm as a core, allows the control margin to freely fluctuate to a certain degree, and improves the adaptive degree and the reaction time of the control algorithm to the environment; by using double closed-loop control as an execution algorithm, on one hand, the ozone fast response to the environmental ozone concentration is fast achieved, and on the other hand, the control degree of the output concentration is enhanced; through calculation and adjustment of static loss, the adaptation degree of a control algorithm to the environment is improved, and the accuracy degree of control of the ozone concentration is provided; controlling inaccurate parameters through a fuzzy control algorithm, and enhancing the robustness of the algorithm; through the arrangement of dead zone detection, excessive reaction of a control algorithm due to data fluctuation is avoided, and the accuracy of the control algorithm on concentration control is improved; the control algorithm has high precision and adjustable precision.
It is noted that, herein, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions. Also, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it will be apparent to those skilled in the art that changes may be made in the embodiments and/or equivalents thereof without departing from the spirit and scope of the invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. The method for controlling the ozone concentration in the space is characterized in that: the control method comprises the following steps:
s100: measuring static loss;
s200: calculating a control margin;
s300: adjusting the PID output upper calibration value by controlling the PID output value calibration coefficient;
s400: static loss is adjusted up;
s500: fuzzy control PID parameters;
s600: controlling the actual output of the PID algorithm.
2. The method of claim 1, wherein: the static loss measurement includes: opening an air inlet valve, closing the output of a high-voltage power supply, resetting the ozone concentration in the environment, opening the output of the ozone high-voltage power supply when the ozone concentration is less than a measurement low value, setting the output of the ozone high-voltage power supply as a minimum output value, obtaining static loss by measuring the running trend of the ozone concentration, and outputting after averaging multiple times of static loss measurement;
the specific calculation method of the static loss measurement comprises the following steps:
s101: obtaining the concentration difference of unit time;
s102: when the concentration difference is smaller than the lower limit of the concentration difference, if the lower counter value is smaller than the maximum lower counter value, the lower counter value is increased by one;
when the concentration difference is smaller than the lower limit of the concentration difference, if the value of the lower counter is larger than or equal to the maximum value of the lower count, the current output concentration is judged to be too low, the trend counter is cleared, and the output of the high-voltage power supply is increased;
s103: when the concentration difference is greater than or equal to the lower concentration difference limit and the concentration difference is less than or equal to the upper concentration difference limit, if the non-change counter value is less than the maximum non-change count value, the non-change counter value continues to accumulate and inherits the upper counter value, and the upper counter value is cleared;
when the concentration difference is greater than or equal to the lower concentration difference limit and the concentration difference is less than or equal to the upper concentration difference limit, if the value of the non-change counter is greater than or equal to the maximum non-change count value, judging that the current output concentration is too low, resetting the trend counter, and increasing the output of the high-voltage power supply;
s104: when the concentration difference is larger than the concentration difference upper limit, if the upper counter value is smaller than the maximum value of the upper counter value, the upper counter value is increased by one;
when the concentration difference is greater than the concentration difference upper limit, if the upper counter value is greater than or equal to the upper counting maximum value, resetting the trend counter, outputting the static loss measured value of the time, and setting a single static loss control bit;
s104: and measuring the static loss for multiple times, calculating an average value, and outputting the final static loss.
3. The method of claim 2, wherein: the specific calculation mode of the control margin comprises the following steps:
s201: when the system is initially started, increasing the adjustment time limit value, and recovering the adjustment time limit value when the system is normally operated;
s202: calculating the deviation on control:
Figure FDA0003396264910000021
wherein, cUpErrIndicating deviation in control, cSetIndicating a set concentration, cPerRepresenting the concentration of the sample;
s203: when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is greater than the upper deviation limit value, closing the lower deviation detection, when the upper deviation detection control position is zero, starting the up dead zone detection, and resetting the up dead zone count value;
Figure FDA0003396264910000022
cUpErr>CUpErrMax，oPID>OPIDUpand c isPre<cSet
Wherein LowerMeasure represents lower deviation detection control bit, Uppermeasure represents upper deviation detection control bit, LowerMeasure and Uppermeasure are Boolean, nUpIndicates an up dead zone count value, CUpErrMaxRepresents the maximum value of the upper deviation, oPIDRepresents PID output, OPIDUpRepresents the upper limit of PID output, CUpErrMaxAnd OPIDUpIs a fixed threshold;
when the PID output is greater than the output upper limit and the sampling concentration is less than the set concentration, if the output upper deviation is less than or equal to the upper deviation limit value, closing the upper deviation detection, and clearing the up-regulation dead zone count value:
Figure FDA0003396264910000023
cUpErr≤CUpErrMax，oPID>OPIDUpand c isPre<cSet；
S204: calculating the deviation under control, and calculating a formula:
Figure FDA0003396264910000024
wherein, cLowErrIndicating a deviation under control;
s205: when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is greater than the lower deviation limit value, the upper deviation detection is closed, and when the lower deviation detection control position is zero, the lower dead zone detection is started, and the lower dead zone count value is cleared:
Figure FDA0003396264910000031
cLowErr>CLowErrMax，oPID≤OPIDLowand c isPre>cSet
Wherein n isLowIndicating a down dead zone count value, CLowErrMaxRepresents the maximum value of the lower deviation, OPIDLowRepresents the lower limit of PID output, CLowErrMaxAnd OPIDLowIs a fixed threshold;
when the PID output is less than or equal to the output lower limit and the sampling concentration is greater than the set concentration, if the output lower deviation is less than or equal to the lower deviation limit value, the lower deviation detection is closed, and the down-regulation dead zone count value is cleared:
Figure FDA0003396264910000032
cLowErr≤CLowErrMax，oPID≤OPIDLowand c isPre>cSet；
S206: when the sampling concentration is greater than the set concentration, if the PID output is greater than the output lower limit, and the PID output is less than or equal to the output upper limit, closing the upper limit and the lower limit for detection:
Figure FDA0003396264910000033
OPIDLow<oPID≤OPIDUpand c isPre>cSet；
S207: when the upper limit detection is allowed, the upper limit detection is counted, after the maximum value is counted, an upper deviation detection control bit is set to zero, and a first adjustment control bit is set:
nUp＝nUp+1，nUp<nUpMaxand UpperMeasure ═ 1
Figure FDA0003396264910000034
nUp≥nUpMaxAnd UpperMeasure ═ 1
Wherein n isUpIndicates the up dead zone count value, nUpMaxThe maximum value of the up-regulation count, wherein Upperenable represents an up-regulation output upper limit allowable control bit;
s208: when the lower limit detection is allowed, counting the lower limit detection, and setting a down control bit after the maximum value is counted:
nLow＝nLow+1，nLow<nLowMaxand LowerMeasure is 1
Figure FDA0003396264910000041
nLow≥NLowMaxAnd LowerMeasure is 1
Wherein, LowerEnable represents a lower output upper limit permission control bit.
4. The method of claim 3, wherein: the specific algorithm of the calibration value on the PID output comprises the following steps:
s301: when the first adjusting output upper limit allowable control position is located, resetting the upper adjusting output upper limit allowable control position, and adjusting the PID output value calibration coefficient;
s302: resetting the lower-regulation output upper limit allowable control bit when the lower-regulation output upper limit allowable control position is one, and calibrating the PID output value to be regulated downwards if the PID output value is greater than the output lower limit; and if the PID output value is less than or equal to the output lower limit and the static loss is greater than the static loss minimum value, the static loss is adjusted downwards.
5. The method of claim 4, wherein: the specific algorithm for static loss up-regulation comprises the following steps:
s401: when the sampling concentration is greater than the set concentration, starting concentration maximum value measurement;
s402: when the measurement of the maximum concentration value is started, if the maximum value of single sampling is less than the sampling concentration, the maximum value is updated to be equal to the sampling concentration;
s403: when the sampling concentration is less than the set concentration and the concentration maximum value measurement is started, closing the maximum value measurement, outputting the concentration maximum value, emptying the single sampling maximum value, and starting the single sampling minimum value measurement;
s404: when the measurement of the minimum concentration value is started, if the single sampling minimum value is greater than the sampling concentration, the updated minimum value is equal to the sampling concentration;
s405: when the measurement of the minimum concentration value is started, if the sampling concentration is greater than the set concentration, the measurement of the maximum value and the minimum value is closed, and the single sampling minimum value is output;
s406: when the concentration maximum value and the concentration minimum value are obtained, the root mean square of the concentration extreme value difference is obtained;
s407: and comparing the extreme value difference root-mean-square with a set difference value, and if the extreme value difference root-mean-square is larger than a limit value, adjusting the static loss value upwards.
6. The method of claim 5, wherein: the fuzzy control PID parameter specific algorithm comprises the following steps:
s501: sampling control deviation and deviation change rate:
Figure FDA0003396264910000051
wherein, cErrRepresents the difference in concentration, cPreRepresenting the concentration of the sample, cSetIndicates the set concentration, eEIndicating deviation of calibration of concentration, EErrIndicating the degree of deviation of the calibration, eEcIndicating the rate of change of the deviation of the calibration of the concentration, cErrLastIndicating the limit of the calibration deviation, EEcRepresenting the limit of the rate of change of the calibration deviation, cErrShowing the historical concentration difference, EErr、EEcIs a fixed threshold;
s502: converting the control interval of the concentration calibration deviation and the change rate of the concentration calibration deviation to-3 to 3, dispersing the continuous number into eight numbers, and then calibrating the data to an integer of 1 to 7:
Figure FDA0003396264910000052
Figure FDA0003396264910000053
Figure FDA0003396264910000054
Figure FDA0003396264910000055
wherein e isFuzzyRepresenting the fuzzified output value, eltRepresenting a left value formatting output value, ertRepresenting the right value formatted output value, x, y, z representing the blurred input value, e1Representing the first column of the fuzzy matrix, e2Representing the second index of the fuzzy matrix, eC1Representing the first row index of the fuzzy matrix, eC2Representing the second row index of the fuzzy matrix, eE1Representing the first row and column elements of the matrix M, eE2Representing the first row and the second column of elements of the matrix M, eEc1Representing the first row and column elements of the matrix N, eEc2Representing the second row and the first column of elements of the matrix M;
s503: according to the influence strength of the concentration calibration deviation and the change rate of the concentration calibration deviation on proportion, integration and differentiation, three 7 multiplied by 7 fuzzy matrixes P are made7×7、I7×7、D7×7Calculating the row and column values of variable coordinates of the fuzzy matrix according to the concentration parameters, and deriving three new fuzzy matrices A2×2、B2×2、C2×2：
Figure FDA0003396264910000061
Wherein
Figure FDA0003396264910000062
Figure FDA0003396264910000063
Wherein
Figure FDA0003396264910000064
Figure FDA0003396264910000065
Wherein
Figure FDA0003396264910000066
Wherein, P7×7A fuzzy matrix representing the ratio coefficients of the deviation of the calibration of the concentration and the rate of change of the deviation of the calibration of the concentration, I7×7Fuzzy matrices representing the deviation of the calibration of the concentration and the rate of change of the deviation of the calibration of the concentration versus the integral coefficient, D7×7A fuzzy matrix representing the deviation of the concentration calibration and the rate of change of the deviation of the concentration calibration versus the differential coefficient, A2×2Ambiguity representing the matrix P about the concentration calibration deviation and the rate of change of the concentration calibration deviationMatrix, B2×2A fuzzy matrix representing the ratio of the matrix I to the concentration calibration deviation and the rate of change of the concentration calibration deviation, C2×2Fuzzy matrix, p, representing the ratio of the matrix D to the concentration calibration deviation and the rate of change of the concentration calibration deviationi,j、ii,j、di,jRespectively representing elements, a, of a matrix P, I, Di,j、bi,j、ci,jRespectively represent elements of matrix A, B, C;
s504: constructing a concentration calibration deviation fuzzification matrix M and a concentration calibration deviation change rate fuzzification matrix N:
M1×2＝[eE1 eE2]
N1×2＝[eEc1 eEc2]T
wherein M is1×2Indicating the concentration calibration deviation matrix, N1×2Representing a concentration calibration deviation change rate matrix;
s505: determining a fuzzy PID parameter (k)p,tI,td)：
Figure FDA0003396264910000067
Figure FDA0003396264910000068
Figure FDA0003396264910000069
Wherein, KPRepresenting a scale factor setting parameter, TIRepresenting an integral coefficient setting parameter, TDRepresenting a differential coefficient setting parameter, KP、TIAnd TDIs a fixed threshold.
7. The method of claim 6, wherein: the actual output of the control PID algorithm is calibrated and limited through static loss and PID output:
Figure FDA0003396264910000071
wherein o isPIDSumRepresenting calculated values of the PID algorithm, i.e. output values, k, without nominal clippingPDenotes the proportionality coefficient, PWRepresenting the proportional weight, tIRepresenting the integral coefficient, tDRepresenting a differential coefficient, s representing the Laplace operator, DDelayRepresenting the differential delay coefficient, DWRepresenting a differential weight, cSetIndicates a set value, cPreRepresenting the current sample concentration, oPIDRepresents the PID algorithm output, oOutRepresenting the actual power control value, l, delivered to the power supplyZeroRepresenting static losses and P representing the control accuracy of the maximum power.
8. The method of controlling ozone concentration in a space according to any one of claims 1 to 7, wherein: the space ozone concentration control method takes PID double closed-loop control as an execution algorithm, the inner loop control quantity of the execution algorithm is a PWM value of a power supply, the inner loop process quantity is power supply power, the outer loop control quantity of the execution algorithm is power supply power, and the outer loop process quantity is air ozone concentration.
CN202111483201.8A 2021-12-07 2021-12-07 Space ozone concentration control method Active CN114237031B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111483201.8A CN114237031B (en) 2021-12-07 2021-12-07 Space ozone concentration control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111483201.8A CN114237031B (en) 2021-12-07 2021-12-07 Space ozone concentration control method

Publications (2)

Publication Number Publication Date
CN114237031A true CN114237031A (en) 2022-03-25
CN114237031B CN114237031B (en) 2024-05-10

Family

ID=80753595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111483201.8A Active CN114237031B (en) 2021-12-07 2021-12-07 Space ozone concentration control method

Country Status (1)

Country Link
CN (1) CN114237031B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6060020A (en) * 1998-04-09 2000-05-09 S.P.M. Recovery Technologies Ltd Method and apparatus for treating objects with ozone
CN1393746A (en) * 2001-07-02 2003-01-29 李晓枫 Self-adaptive fuzzy controller for PID parameters
JP2003146624A (en) * 2001-11-08 2003-05-21 Mitsubishi Electric Corp Apparatus for controlling ozone concentration
CN1544912A (en) * 2003-11-28 2004-11-10 华中科技大学 Concentration changing and sending adjusting device
CN102759506A (en) * 2011-04-27 2012-10-31 中国科学院大连化学物理研究所 Method and device for measuring absolute concentration of singlet-state oxygen
CN109839963A (en) * 2019-01-23 2019-06-04 北京雪迪龙科技股份有限公司 A kind of ultraviolet irradiation formula ozone generator concentration control method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6060020A (en) * 1998-04-09 2000-05-09 S.P.M. Recovery Technologies Ltd Method and apparatus for treating objects with ozone
CN1393746A (en) * 2001-07-02 2003-01-29 李晓枫 Self-adaptive fuzzy controller for PID parameters
JP2003146624A (en) * 2001-11-08 2003-05-21 Mitsubishi Electric Corp Apparatus for controlling ozone concentration
CN1544912A (en) * 2003-11-28 2004-11-10 华中科技大学 Concentration changing and sending adjusting device
CN102759506A (en) * 2011-04-27 2012-10-31 中国科学院大连化学物理研究所 Method and device for measuring absolute concentration of singlet-state oxygen
CN109839963A (en) * 2019-01-23 2019-06-04 北京雪迪龙科技股份有限公司 A kind of ultraviolet irradiation formula ozone generator concentration control method

Also Published As

Publication number Publication date
CN114237031B (en) 2024-05-10

Similar Documents

Publication Publication Date Title
US6554198B1 (en) Slope predictive control and digital PID control
US20090312876A1 (en) Mass flow controller
EP1196833B1 (en) System and method for a variable gain proportional-integral (pi) controller
CN109630446B (en) Gas estimation device and vacuum exhaust device
US5511006A (en) Method and apparatus for determining an air quality level
US20040256473A1 (en) Slope predictive control and digital PID control for a variable temperature control system
CN101836173A (en) Method and device for controlling pressure of vacuum container
CN111579013B (en) Gas mass flow controller and flow calibration method thereof
JP2010079827A (en) Mass flow controller
CN114237031A (en) Space ozone concentration control method
WO2016178739A1 (en) Nonlinear control of mass flow controller devices using sliding mode
JP2007034550A (en) Mass flow controller
JP2005538462A (en) Method and apparatus for determining hysteresis
US8056579B2 (en) Mass flow controller
JP2009015514A (en) Pid control device
JP7473397B2 (en) Mass flow controller and hunting suppression method
CN112748660A (en) Memory, heating furnace outlet temperature control method, device and equipment
CN113739385B (en) Air conditioner control method, device and system and storage medium
CN110531715B (en) Control method and system for LPCVD (low pressure chemical vapor deposition) process production environment
JP7089336B2 (en) Control device
KR101668483B1 (en) Mass flow controller
JP7089335B2 (en) Control device
JP2024048549A (en) Flow Control Device
JPS61198302A (en) Regulator
JP2024067716A (en) Flow Control System

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