CN116679026A - Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method - Google Patents
Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method Download PDFInfo
- Publication number
- CN116679026A CN116679026A CN202310772649.4A CN202310772649A CN116679026A CN 116679026 A CN116679026 A CN 116679026A CN 202310772649 A CN202310772649 A CN 202310772649A CN 116679026 A CN116679026 A CN 116679026A
- Authority
- CN
- China
- Prior art keywords
- unbiased
- window size
- impulse response
- finite impulse
- dissolved oxygen
- 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
Links
- 239000010865 sewage Substances 0.000 title claims abstract description 109
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 title claims abstract description 94
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 94
- 239000001301 oxygen Substances 0.000 title claims abstract description 94
- 230000004044 response Effects 0.000 title claims abstract description 93
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000001914 filtration Methods 0.000 title claims abstract description 18
- 238000011282 treatment Methods 0.000 claims abstract description 42
- 238000013178 mathematical model Methods 0.000 claims abstract description 15
- 230000007246 mechanism Effects 0.000 claims abstract description 12
- 238000012544 monitoring process Methods 0.000 claims abstract description 9
- 239000011159 matrix material Substances 0.000 claims description 29
- 230000008569 process Effects 0.000 claims description 17
- 238000005259 measurement Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 244000005700 microbiome Species 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 6
- 239000010802 sludge Substances 0.000 claims description 5
- 239000000758 substrate Substances 0.000 claims description 4
- 238000011166 aliquoting Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000003102 growth factor Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 239000005416 organic matter Substances 0.000 claims description 3
- 238000004062 sedimentation Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims 9
- 238000004590 computer program Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 239000000126 substance Substances 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000003860 storage Methods 0.000 description 4
- 230000006866 deterioration Effects 0.000 description 3
- 239000003344 environmental pollutant Substances 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000011221 initial treatment Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 231100000719 pollutant Toxicity 0.000 description 2
- 238000004065 wastewater treatment Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 241001148470 aerobic bacillus Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006065 biodegradation reaction Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000013526 transfer learning Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
- Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
- Y02W10/00—Technologies for wastewater treatment
- Y02W10/10—Biological treatment of water, waste water, or sewage
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Chemical & Material Sciences (AREA)
- Computing Systems (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Medicinal Chemistry (AREA)
- Pathology (AREA)
- Activated Sludge Processes (AREA)
Abstract
The application relates to the technical field of sewage treatment, and discloses a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method, which comprises the following steps: establishing a mathematical model of a sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration; estimating a state value of a system by using an unbiased finite impulse response filter in a window interval, acquiring estimated values of system states of the unbiased finite impulse response filter under different window sizes, and selecting a window size corresponding to the window with highest estimated value precision as an optimal window size according to a Markov distance; and estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter under a window interval corresponding to the optimal window size. The application can realize real-time estimation of the concentration of the dissolved oxygen in the sewage, reduce the estimation error of the filter and increase the estimation precision.
Description
Technical Field
The application relates to the technical field of sewage treatment, in particular to a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method.
Background
Sewage treatment means that harmful substances in sewage are removed through a series of physical, chemical and biological treatment processes, so that the harmful substances reach the emission standards regulated by the national and local governments, and the influence on the environment and human health is reduced. The sewage treatment process generally comprises a plurality of stages of pretreatment, primary treatment, secondary treatment and tertiary treatment, wherein the pretreatment mainly comprises the removal of large-particle substances and sediments, the primary treatment mainly comprises the removal of suspended substances and biochemical oxygen demand through precipitation and filtration, the secondary treatment mainly comprises the removal of nutrient substances such as nitrogen and phosphorus through biodegradation, and the tertiary treatment mainly comprises the removal of refractory substances and trace pollutants through chemical treatment. The sewage treatment is an important environmental protection work, and can effectively reduce the pollution to the environment and the threat to human health.
The sewage treatment has the characteristics of large scale, multiple pollutant types, large water quality fluctuation, large treatment difficulty, high operation cost and the like in industry, and the characteristics of the sewage treatment increase the treatment difficulty and the cost to a certain extent, so that advanced technology and management means are needed to ensure the balance of treatment effect and economic benefit. By 2022, 3 months, there are ten thousands of wastewater treatment plants taking a pollution discharge license nationwide, wherein 98.6% of the wastewater treatment plants have published limit information on the total amount of self-pollutant discharge or the concentration of discharge. According to the sewage treatment capacity, the sewage treatment plant of 1-5 ten thousand tons/day is the most, and the ratio is 33.2%; a sewage treatment plant of 5000-10000 tons/day accounts for 11.7%; the sewage treatment plant of 1000-5000 tons/day accounts for 25.8%; a sewage treatment plant of 1000 tons/day or less accounts for 16.4%; the sewage treatment plant with the sewage treatment plant weight of less than 1 ten thousand tons/day accounts for 52.9 percent.
The sewage treatment process has the characteristics of nonlinearity, time variation, randomness and large time lag, the modeling and null value process of the sewage treatment process is still not easy at present, and in a specific sewage treatment field environment, the sensor can cause errors between the calculated state value and the actual value due to the interference effect. The quality of the water quality is directly determined by the quality of the control of the concentration of the dissolved oxygen in the sewage treatment process. Too little or too much dissolved oxygen concentration in the activated sludge basin can lead to deterioration of the sludge environment: the insufficient concentration of dissolved oxygen can cause the reduction of the growth rate of aerobic bacteria, thereby reducing the quality of effluent; conversely, if the concentration of dissolved oxygen is too high, the flocculant is destroyed, resulting in deterioration of the settling property of suspended solids and energy waste.
The finite impulse response is a popular research direction in recent years, and particularly in the unbiased finite impulse response filter, the body and shadow of which can be seen in many related fields, such as the fields of track tracing, transfer learning and the like, and the sensor can be used for measuring various parameters of water quality, such as PH value, temperature, dissolved oxygen, conductivity and the like, and the data can be input into the filter for predicting the future state of the water quality and controlling the future state of the water quality. However, the concentration of dissolved oxygen in sewage is a dynamic process, and real-time monitoring and data acquisition are required to ensure that a filter can effectively process dissolved oxygen signals in time, but real-time monitoring and data acquisition cannot be realized by simply using an unbiased finite impulse response filter. In addition, the prior art lacks consideration and evaluation on stability, instantaneity and the like of the filter, and cannot ensure the application effect and stability of the filter, so that the estimation on the concentration of the dissolved oxygen in the sewage is inaccurate.
Disclosure of Invention
Therefore, the application aims to solve the technical problems of overcoming the defects in the prior art, and providing the self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method which can realize real-time estimation of the sewage dissolved oxygen concentration, reduce filter estimation errors and increase estimation accuracy.
In order to solve the technical problems, the application provides a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method, which comprises the following steps:
establishing a mathematical model of a sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration;
estimating a state value of the system by using an unbiased finite impulse response filter in a window interval by combining a mathematical model of the system and an observation model of the sensor;
acquiring estimated values of system states of unbiased finite impulse response filters under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size;
and estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter under a window interval corresponding to the optimal window size.
In one embodiment of the application, the mathematical model of the sewage dissolved oxygen concentration system according to the sewage treatment mechanism comprises the following steps:
the state space equation of the concentration of the dissolved oxygen is established according to the treatment mechanism of the sewage treatment by the activated sludge process and is as follows:
x k ＝Ax k-1 +w k ，
wherein ：
k is the time index, x k Is the system state at time k, x k ＝[x 1,k ,x 2,k ,x 3,k ] T； wherein ,x1,k Is the mass concentration of microorganisms, x 2,k Is the mass concentration of the substrate, x 3,k T represents transposition for mass concentration of dissolved oxygen;
w k system noise term at time k, w k -N (0, Q), where Q is the covariance matrix of the system noise term;
a is the state transition matrix of the system, wherein ,uH Is the maximum growth rate of microorganisms, k d Is the endogeneous hysteresis parameter, C is the concentration factor of the secondary sedimentation tank, Q w Is the flow of the sewage, V is the volume of the reactor, Y NH For the observed growth factor, Q in For inflow, f is the organic matter of interestFactor with oxygen demand, f x Is the water pump factor, delta is the impulse coefficient set for dissolved oxygen.
In one embodiment of the present application, the method for modeling the observation of the sensor for dissolved oxygen concentration of sewage includes:
the method for establishing the observation model of the sensor comprises the following steps:
y k ＝Cx k +v k ，
wherein ,yk For the observation value of the dissolved oxygen of the sewage measured by a sensor at the moment k, C is the observation matrix of the system, x k Is the system state at time k, v k Noise term, v, which is an observation of time k k N (0, R), where R is the covariance matrix of the observed noise term.
In one embodiment of the present application, the estimating the state value of the system using the unbiased finite impulse response filter in the window section includes:
batch processing of unbiased finite impulse response filters using standard least squares in window intervals [ m, k]The internal estimated value to obtain the estimated value of the system state at k timeThe method comprises the following steps:
wherein ,Hm,k Representing a mapping matrix, Y m,k Representing an extended observation vector, T representing a transpose; m represents the initial time, m=k-n+1, N being the window size;
transforming the system state estimation value to obtain a recursive Kalman system state estimation valueThe method comprises the following steps:
wherein ,Ak Representing the state transition matrix of the system, G k Representing generalized noise power gain matrix, C k Representing an observation matrix of the system;
prediction prior state estimationThe method comprises the following steps:
wherein ,is a priori state estimate, l is an auxiliary variable, iterating from l=m+k to l=k ending to obtain an estimate at the time scale K;
iteratively updating the prior state estimation to obtain the posterior state estimationThe method comprises the following steps:
wherein ,Kl Is the offset correction gain and,G l is the generalized noise power gain,/->z l Is the updated measurement residual,/>
In one embodiment of the present application, the obtaining the estimated value of the system state of the unbiased finite impulse response filter under different window sizes, selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance, includes:
randomly selecting an estimated value of a system state of the unbiased finite impulse response filter in a preset range of the window size, and primarily evaluating the basic set by using a mahalanobis distance to obtain an initial optimal window size;
expanding the initial optimal window size, re-evaluating the estimated value of the system state corresponding to the expanded initial optimal window size by using the mahalanobis distance, and selecting the window size corresponding to the highest estimated value precision at the moment as the final optimal window size.
In one embodiment of the present application, the estimating values of the system states of the unbiased finite impulse response filter are randomly selected in groups within a preset range of window sizes to obtain a basic set, and the basic set is first estimated by using a mahalanobis distance to obtain an initial optimal window size, specifically:
setting the maximum value N of the window size N preset range max And a minimum value N min Interval [ N ] min ,N max ]Aliquoting into c+1 parts to obtain a set wherein ,/>The method comprises the steps of equally dividing a set into c+1 intervals, wherein i is { 1..c }, and c is a preset natural number;
at the collectionIn each of the c+1 intervals of (2) takes a value +.>i.e { 1..c+1 }, run +.>The corresponding unbiased finite impulse response filter generates respective state estimation values, and calculates each state estimation valueThe mahalanobis distance value is:
wherein ,yk Is the measurement vector of the data set,is an estimated vector, C k Is a measurement matrix, R k Is the measurement noise covariance;
and selecting the estimated value with the minimum Markov distance value as an initial optimal estimated value, and taking the window size corresponding to the initial optimal estimated value as the initial optimal window size.
In one embodiment of the present application, the expanding the initial optimal window size, re-evaluating the estimated value of the system state corresponding to the expanded initial optimal window size by using the mahalanobis distance, and selecting the window size corresponding to the highest estimated value precision at this time as the final optimal window size, which specifically includes:
at an initial optimal window sizeThe basic filter bank is obtained by expanding the basic filter bank to the left and the right sides respectively, and the basic filter bank is integrated as follows:
wherein ,respectively correspond to->d is the expansion distance of window size;
operationCorresponding unbiased finite pulseGenerating respective state estimation values by response filters, and calculating corresponding Markov distances to obtain +.>According to->Is selected to the final optimal window size.
In one embodiment of the application, the methodWhen the final optimal window size is selected, three conditions are selected, specifically:
first case: if it isThe final optimal window size
Second case: if it isAn unbiased finite impulse response filter with a smaller window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->The optimal window size at the current moment is the optimal window size, the expansion is stopped, and the final optimal window size is +.>
If it isThen the "add window smaller unbiased finite impulse response filter is repeated with a window size ofContinuing to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance;
third case: if it isAn unbiased finite impulse response filter with a larger window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->The optimal window size at the current moment is the optimal window size, the expansion is stopped, and the final optimal window size is +.>
If it isThen the "add window larger unbiased finite impulse response filter is repeated with a window size ofThe mahalanobis distance is continued to be used to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance.
In one embodiment of the present application, the following conditions are satisfied when the unbiased finite impulse response filter with smaller add window and the unbiased finite impulse response filter with larger add window are used simultaneously:
at time scale k, the window size of the unbiased finite impulse response filter must be smaller than k;
the mahalanobis distance calculated by adding a state estimation of an unbiased finite impulse response filter with a smaller or larger window needs to be smaller than the mahalanobis distance calculated from the previous state estimation;
the newly added window size must be within the initial preset range.
The application also provides a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation system, which is characterized by comprising:
the modeling module is used for establishing a mathematical model of the sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration;
the estimated value acquisition module is used for estimating the state value of the system by using an unbiased finite impulse response filter in a window interval by combining the mathematical model of the system and the observation model of the sensor;
the optimal window size acquisition module is used for acquiring estimated values of the system states of the unbiased finite impulse response filter under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size;
and the real-time estimation module is used for estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter in a window interval corresponding to the optimal window size.
Compared with the prior art, the technical scheme of the application has the following advantages:
the application combines an unbiased finite impulse response filter and utilizes a sensor to obtain data, introduces the difference between the Marsdian distance estimated value and the measured value, updates the filter window on line, realizes the real-time estimation of the concentration of the dissolved oxygen in the sewage, reduces the estimation error of the filter and increases the estimation precision.
Drawings
In order that the application may be more readily understood, a more particular description of the application will be rendered by reference to specific embodiments thereof which are illustrated in the appended drawings, in which:
fig. 1 is a flow chart of the present application.
FIG. 2 is a schematic diagram of a process for window size primary screening in accordance with the present application.
FIG. 3 is a schematic diagram of a window size expansion reevaluation process in accordance with the present application.
FIG. 4 is a graph showing the effect of the method of the present application on the concentration of dissolved oxygen in sewage.
Fig. 5 is a graph showing a comparison of Root Mean Square Error (RMSE) of the concentration of dissolved oxygen in sewage estimated by using the method of the present application, the kalman method, and the fixed window UFIR filtering in the example of the present application.
Detailed Description
The present application will be further described with reference to the accompanying drawings and specific examples, which are not intended to be limiting, so that those skilled in the art will better understand the application and practice it.
Referring to FIG. 1, the application discloses a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method. In order to obtain more accurate dissolved oxygen concentration of sewage and solve the problem of environmental deterioration, the application introduces an unbiased finite impulse response filter, and the filter needs to obtain an optimal estimation window firstly in the running process so as to drive an estimated value to approach to the optimal under the meaning of minimum mean square error. However, for the sewage oxygen concentration, when the characteristics of the sewage oxygen concentration are interfered by the outside or the modeling is inaccurate, the filter capable of adjusting the window on line is more advantageous than the filter with the fixed window, and can be better adapted to the actual situation of the complex and intricate, so that the accuracy of the sewage dissolved oxygen concentration estimation value can be improved by introducing the mahalanobis distance measurement estimation precision, and performing the two stages of primary screening and expansion and re-estimation, so that the optimal window at each moment can be selected and estimated on the premise of acceptable calculated quantity, the estimation of the sewage dissolved oxygen concentration at each moment can be obtained, and the accuracy of the sewage dissolved oxygen concentration estimation value can be improved. The method specifically comprises the following steps:
s1: and establishing a mathematical model of a sewage dissolved oxygen concentration system according to a sewage treatment mechanism.
The state space equation of the concentration of the dissolved oxygen is established according to the treatment mechanism of the sewage treatment by the activated sludge process and is as follows:
x k ＝Ax k-1 +w k ，
wherein ：
k is the time index, x k Is the system state at time k, x k ＝[x 1,k ,x 2,k ,x 3,k ] T； wherein ,x1,k Is the mass concentration of microorganisms, x 2,k Is the mass concentration of the substrate, x 3,k T represents transposition for mass concentration of dissolved oxygen;
w k system noise term at time k, w k Obeying a gaussian distribution with zero mean value, w k -N (0, Q), where Q is the covariance matrix of the system noise term;
a is the state transition matrix of the system, wherein ,uH Is the maximum growth rate of microorganisms, k d Is the endogeneous hysteresis parameter, C is the concentration factor of the secondary sedimentation tank, Q w Is the flow of the sewage, V is the volume of the reactor, Y NH For the observed growth factor, Q in As inflow amountF is a factor linking organic matter and oxygen demand, f x Is the water pump factor, delta is the impulse coefficient set for dissolved oxygen.
S2: and monitoring the concentration of the dissolved oxygen in the sewage by using a sensor, and establishing an observation model of the sensor of the concentration of the dissolved oxygen in the sewage.
The method for establishing the observation model of the sensor comprises the following steps:
y k ＝Cx k +v k ，
wherein ,yk For the observation value of the dissolved oxygen in the sewage measured by the sensor at the time k, C is the observation matrix of the system, and in this embodiment, c= [0, 1]；x k Is the system state at time k, v k Noise term, v, which is an observation of time k k Obeying a gaussian distribution with zero mean value v k N (0, R), where R is the covariance matrix of the observed noise term.
S3: an Unbiased Finite Impulse Response (UFIR) filter is used to estimate a state value of the system within a window interval in combination with a mathematical model of the system and an observation model of the sensor.
The UFIR filter uses N measured values from the initial time m=k-N+1 to k to estimate, and the value of N is XX; no initial state x is required in the iteration interval m And an initial error covariance matrix p m Nor does it require a noise covariance matrix Q i and Ri Wherein i is e [ m, k]. UFIR filter can forget the measured value outside the window before m time, and determine the state estimated value at m+K-1 time in a batch processing modeWhere K is the dimension of the system state. The derivation process of the iterative form of the UFIR filter specifically comprises the following steps:
s3-1: batch processing of unbiased finite impulse response filters using standard least squares in window intervals [ m, k]The internal estimated value to obtain the estimated value of the system state at k timeThe method comprises the following steps:
wherein ,Hm,k Representing a mapping matrix, Y m,k Representing an extended observation vector, T representing a transpose; m represents the initial time, m=k-n+1, N being the window size;
s3-2: transforming the system state estimation value to obtain a recursive Kalman system state estimation valueThe method comprises the following steps:
wherein ,Ak Representing the state transition matrix of the system, G k Representing generalized noise power gain matrix, C k Representing the observation matrix of the system.
S3-3: each recursion in the UFIR filter has two phases: prediction and updating. The UFIR algorithm does not require noise statistics, so the prediction stage only requires prediction of prior state estimationThe method comprises the following steps:
wherein ,is a priori state estimate, l is an auxiliary variable, iterating from l=m+k to l=k ending to obtain an estimate at the time scale K.
S3-4: state estimation is improved by combining current state observation, prior state estimation is updated in an iterative mode, and posterior state estimation is obtainedThe method comprises the following steps:
wherein ,Kl Is the offset correction gain and,G l is generalized noise power gain (generalizednoise power gain, GNPG), and is>z l Is the updated measurement residual,/>
S4: and searching an optimal filtering window in real time. And acquiring estimated values of the system state of the unbiased finite impulse response filter under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size. The adjustment parameter required for this process is the window size N, which must be optimal, i.e. N opt It is guaranteed that the MSE is minimized. Based on this, a UFIR (MUFIR) filter based on mahalanobis distance window adaptation is proposed, the basic idea is: using a filter bank, wherein the filter bank comprises UFIR filters with different window sizes, aiming at obtaining a plurality of estimates from the UFIR filters with different window sizes, measuring the estimation accuracy by using the Markov distance, and selecting the window size of the UFIR filter with the smallest Markov distance, i.e. the highest estimation accuracy, as N of the current moment opt I.e.The method comprises the following steps:
s4-1: window size primary screening: and randomly selecting an estimated value of the system state of the unbiased finite impulse response filter in a preset range of the window size, and primarily evaluating the basic set by using the mahalanobis distance to obtain the initial optimal window size. As shown in fig. 2, specifically:
s4-1-1: firstly, design parameters are put forward, and the maximum value N of a preset range of window size N is set max And a minimum value N min To minimize dead space, N is typically set min Set to minimum allowable value (same as system order), for N max It has an impact on the computational burden of the system and therefore needs to be selected according to the particular system. In the range of preset N [ N min ,N max ]Thereafter, section [ N ] min ,N max ]Aliquoting into c+1 parts to obtain a set wherein ,/>For equally dividing the set into c+1 intervals as much as possible, i e { 1..c }; c is a preset natural number, and the value of c has a certain influence on the calculation load and the estimation precision, and in the embodiment, the value of c is c epsilon [2,5 ]]。
S4-1-2: at the collectionIn each of the c+1 intervals of (2) takes a value +.>There is->i.e. { 1..c+1 }. Run->The corresponding unbiased finite impulse response filter generates the respective state estimate +.>The mahalanobis distance value of each state estimation value is calculated as follows:
wherein ,yk Is the measurement vector of the data set,is an estimated vector, C k Is a measurement matrix, R k Is the measurement noise covariance;
s4-1-3: selecting the estimated value with the minimum Markov distance value as an initial optimal estimated value, and taking the window size corresponding to the initial optimal estimated value as an initial optimal window size, namely the initial optimal window size
S4-2: window size expansion reevaluation: expanding the initial optimal window size, re-evaluating the estimated value of the system state corresponding to the expanded initial optimal window size by using the mahalanobis distance, and selecting the window size corresponding to the highest estimated value precision at the moment as the final optimal window size. As shown in fig. 3, in particular
S4-2-1: in order to expand, a basic window needs to be obtained first, and the basic window is obtained at the end of the operation of the processAt an initial optimal window size->The basic filter bank is obtained by expanding the basic filter bank to the left and the right sides respectively, and the basic filter bank is integrated as follows:
wherein ,respectively correspond to->d is the expansion distance of window size; d is a design parameter, typically set to 1 or 2.
S4-2-2: operationThe corresponding unbiased finite impulse response filter generates respective state estimation values, and the corresponding mahalanobis distance is calculated to obtain +.>
S4-2-3: according toIs selected to the final optimal window size. The selection is carried out in three cases, specifically:
first case: if it isThe final optimal window size
Second case: if it isAn unbiased finite impulse response filter with a smaller window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->The optimal window size at the current moment is the optimal window size, the expansion is stopped, and the final optimal window size is +.>
If it isThen the "add window smaller unbiased finite impulse response filter is repeated with a window size ofContinuing to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance;
third case: if it isAn unbiased finite impulse response filter with a larger window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->The optimal window size at the current moment is the optimal window size, the expansion is stopped, and the final optimal window size is +.>
If it isThen the "add window larger unbiased finite impulse response filter is repeated with a window size ofThe mahalanobis distance is continued to be used to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance.
In short, if a new UFIR filter is added, the calculated D k Calculated D compared with the original UFIR filter k Small, adding UFIR filters with new window sizes to the filter bank; otherwise, stopping expansion, and selecting the D with the smallest current moment k The corresponding window is N of the current moment opt The corresponding UFIR filter is the current time filter, and the corresponding state estimate is the current time optimal estimate.
In adding an unbiased finite impulse response filter with a smaller window and adding an unbiased finite impulse response filter with a larger window, the following conditions need to be satisfied at the same time:
(1) At time scale k, the window size of the unbiased finite impulse response filter must be smaller than k;
(2) The mahalanobis distance calculated by adding a state estimation of an unbiased finite impulse response filter with a smaller or larger window needs to be smaller than the mahalanobis distance calculated from the previous state estimation;
(3) Newly added window size N new Must be within an initial preset range, i.e., N new ∈[N min ,N max ]。
When any one of the conditions is not satisfied, stopping expansion to obtain N opt 。
S5: and finally estimating online. And estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter under a window interval corresponding to the optimal window size. The unbiased finite impulse response filter capable of adaptively changing the window on line can be obtained by searching the optimal window size, and the filter can estimate the state value of the dissolved oxygen concentration of the sewage in a complex sewage environment, so that on-line evaluation is realized.
The application also discloses a self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation system, which comprises a modeling module, an estimated value acquisition module, an optimal window size acquisition module and a real-time estimation module. The modeling module is used for establishing a mathematical model of the sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration. The estimated value acquisition module is used for estimating the state value of the system by using an unbiased finite impulse response filter in a window interval in combination with the mathematical model of the system and the observation model of the sensor. The optimal window size acquisition module is used for acquiring estimated values of the system states of the unbiased finite impulse response filter under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size. The real-time estimation module is used for estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter under the window interval corresponding to the optimal window size.
The application also discloses a computer readable storage medium, on which a computer program is stored, which when being executed by a processor, realizes the self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method.
The application also discloses a device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor realizes the self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method when executing the computer program.
Compared with infinite impulse response (such as Kalman filtering) and an unbiased finite impulse response filtering method using a fixed window, the method combines an unbiased finite impulse response filter and utilizes a sensor to obtain data, and the filter window is updated online by introducing the difference between the Mars distance estimated value and the measured value, so that the real-time estimation of the concentration of the dissolved oxygen in the sewage is realized, the estimation error of the filter is reduced, the estimation precision is increased, and the method has very important significance for accurately estimating the concentration of the dissolved oxygen.
In order to further prove the beneficial effects of the application, simulation experiments are carried out by taking the following parameter settings as examples in the embodiment:
the parameters of the sensor are respectively as follows:the system has 3 states, namely the first state: mass concentration of microorganisms mg/L, second state: mass concentration of substrate mg/L and third state: the concentration of dissolved oxygen in sewage is mg/L. The initial values for the 3 states were all set to 0mg/L.
The method of the application is used for tracking the concentration state of the dissolved oxygen in the sewage, and the tracking effect is shown in figure 4. MUFIR in fig. 4 shows the results using the present application, and TURE shows the true value. It can be seen that the estimated value of the present application substantially coincides with the actual value.
The method, the Kalman method and the fixed window UFIR filter are respectively used for estimating the concentration state of the dissolved oxygen in the sewage, and the Root Mean Square Error (RMSE) between each method and the actual value is calculated, and the result is shown in figure 5. In fig. 5 KF shows the result of using the kalman method, MUFIR shows the result of using the method of the present application, UFIR shows the result of using the fixed window UFIR filtering. It can be seen that the root mean square error of the application is minimum, the difference between the root mean square error and the actual value is minimum, and the estimation accuracy is effectively improved.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
It is apparent that the above examples are given by way of illustration only and are not limiting of the embodiments. Other variations and modifications of the present application will be apparent to those of ordinary skill in the art in light of the foregoing description. It is not necessary here nor is it exhaustive of all embodiments. And obvious variations or modifications thereof are contemplated as falling within the scope of the present application.
Claims (10)
1. The method for estimating the concentration of the dissolved oxygen in the sewage by self-adaptive unbiased finite impulse response filtering is characterized by comprising the following steps of:
establishing a mathematical model of a sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration;
estimating a state value of the system by using an unbiased finite impulse response filter in a window interval by combining a mathematical model of the system and an observation model of the sensor;
acquiring estimated values of system states of unbiased finite impulse response filters under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size;
and estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter under a window interval corresponding to the optimal window size.
2. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 1, characterized in that: the mathematical model for establishing a sewage dissolved oxygen concentration system according to a sewage treatment mechanism comprises the following steps:
the state space equation of the concentration of the dissolved oxygen is established according to the treatment mechanism of the sewage treatment by the activated sludge process and is as follows:
wherein ：
k is the time index, x k Is the system state at time k, x k ＝[x 1,k ,x 2,k ,x 3,k ] T； wherein ,x1,k Is the mass concentration of microorganisms, x 2,k Is the mass concentration of the substrate, x 3,k T represents transposition for mass concentration of dissolved oxygen;
w k as a system noise term at time k,w k -N (0, Q), where Q is the covariance matrix of the system noise term;
a is the state transition matrix of the system, wherein ,uH Is the maximum growth rate of microorganisms, k d Is the endogeneous hysteresis parameter, C is the concentration factor of the secondary sedimentation tank, Q w Is the flow of the sewage, V is the volume of the reactor, Y NH For the observed growth factor, Q in For inflow, f is a factor linking organic matter to oxygen demand, f x Is the water pump factor, delta is the impulse coefficient set for dissolved oxygen.
3. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 2, characterized in that: the observation model for establishing the sensor of the dissolved oxygen concentration of the sewage comprises the following steps:
the method for establishing the observation model of the sensor comprises the following steps:
y k ＝Cx k +v k ，
wherein ,yk For the observation value of the dissolved oxygen of the sewage measured by a sensor at the moment k, C is the observation matrix of the system, x k Is the system state at time k, v k Noise term, v, which is an observation of time k k N (0, R), where R is the covariance matrix of the observed noise term.
4. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 3, characterized in that: the estimating a state value of the system using an unbiased finite impulse response filter within a window interval includes:
batch processing of unbiased finite impulse response filters using standard least squares in window intervals [ m, k]The internal estimated value to obtain the estimated value of the system state at k timeThe method comprises the following steps:
wherein ,Hm,k Representing a mapping matrix, Y m,k Representing an extended observation vector, T representing a transpose; m represents the initial time, m=k-n+1, N being the window size;
transforming the system state estimation value to obtain a recursive Kalman system state estimation valueThe method comprises the following steps:
wherein ,Ak Representing the state transition matrix of the system, G k Representing generalized noise power gain matrix, C k Representing an observation matrix of the system;
prediction prior state estimationThe method comprises the following steps:
wherein ,is a priori state estimate, l is an auxiliary variable, iterating from l=m+k to l=k ending to obtain an estimate at the time scale K;
iteratively updating the prior state estimation to obtain the posterior state estimationThe method comprises the following steps:
wherein ,Kl Is the offset correction gain and,G l is the generalized noise power gain,/->z l Is the updated measurement residual,/>
5. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 4, characterized in that: the obtaining the estimated value of the system state of the unbiased finite impulse response filter under different window sizes, selecting the corresponding window size with the highest estimated value precision according to the mahalanobis distance, including:
randomly selecting an estimated value of a system state of the unbiased finite impulse response filter in a preset range of the window size, and primarily evaluating the basic set by using a mahalanobis distance to obtain an initial optimal window size;
expanding the initial optimal window size, re-evaluating the estimated value of the system state corresponding to the expanded initial optimal window size by using the mahalanobis distance, and selecting the window size corresponding to the highest estimated value precision at the moment as the final optimal window size.
6. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 5, characterized in that: the method comprises the steps of randomly selecting an estimation value of a system state of an unbiased finite impulse response filter in a preset range of window sizes, and obtaining an initial optimal window size by using a primary estimation basic set of a mahalanobis distance, wherein the method comprises the following steps of:
setting the maximum value N of the window size N preset range max And a minimum value N min Interval [ N ] min ,N max ]Aliquoting into c+1 parts to obtain a set wherein ,/> The method comprises the steps of equally dividing a set into c+1 intervals, wherein i is { 1..c }, and c is a preset natural number;
at the collectionIn each of the c+1 intervals of (2) takes a value +.>i.e { 1..c+1 }, run +.>The corresponding unbiased finite impulse response filter generates respective state estimation values, and the mahalanobis distance value of each state estimation value is calculated as:
wherein ,yk Is the measurement vector of the data set,is an estimated vector, C k Is a measurement matrix, R k Is the measurement noise covariance;
and selecting the estimated value with the minimum Markov distance value as an initial optimal estimated value, and taking the window size corresponding to the initial optimal estimated value as the initial optimal window size.
7. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 6, characterized in that: expanding the initial optimal window size, re-evaluating the estimated value of the system state corresponding to the expanded initial optimal window size by using the mahalanobis distance, and selecting the window size corresponding to the highest estimated value precision at the moment as the final optimal window size, wherein the method specifically comprises the following steps:
at an initial optimal window sizeThe basic filter bank is obtained by expanding the basic filter bank to the left and the right sides respectively, and the basic filter bank is integrated as follows:
wherein ,respectively correspond to->d is the expansion distance of window size;
operationThe corresponding unbiased finite impulse response filter generates respective state estimation values, and the corresponding mahalanobis distance is calculated to obtain +.>According to->Is selected to the final optimal window size.
8. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 7, characterized in that: said basis isWhen the final optimal window size is selected, three conditions are selected, specifically:
first case: if it isThe final optimal window size +.>
Second case: if it isAn unbiased finite impulse response filter with a smaller window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->Namely, isThe expansion of the optimal window size at the current moment is stopped, and the final optimal window size is +.>
If it isThen the "add window smaller unbiased finite impulse response filter with window size +.>Continuing to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance;
third case: if it isAn unbiased finite impulse response filter with a larger window is added, the window size is +.>Continuing to evaluate the Markov distance to obtain +.>
If it isThen consider->Corresponding->The optimal window size at the current moment is the optimal window size, the expansion is stopped, and the final optimal window size is +.>
If it isThen the "add window larger unbiased finite impulse response filter with window size +.>The mahalanobis distance is continued to be used to evaluate the operation until the current mahalanobis distance is greater than the last mahalanobis distance.
9. The adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation method as claimed in claim 8, characterized in that: when the unbiased finite impulse response filter with smaller adding window and the unbiased finite impulse response filter with larger adding window are used, the following conditions are satisfied at the same time:
at time scale k, the window size of the unbiased finite impulse response filter must be smaller than k;
the mahalanobis distance calculated by adding a state estimation of an unbiased finite impulse response filter with a smaller or larger window needs to be smaller than the mahalanobis distance calculated from the previous state estimation;
the newly added window size must be within the initial preset range.
10. An adaptive unbiased finite impulse response filtered sewage dissolved oxygen concentration estimation system, comprising:
the modeling module is used for establishing a mathematical model of the sewage dissolved oxygen concentration system according to a sewage treatment mechanism, monitoring the sewage dissolved oxygen concentration by using a sensor, and establishing an observation model of the sensor of the sewage dissolved oxygen concentration;
the estimated value acquisition module is used for estimating the state value of the system by using an unbiased finite impulse response filter in a window interval by combining the mathematical model of the system and the observation model of the sensor;
the optimal window size acquisition module is used for acquiring estimated values of the system states of the unbiased finite impulse response filter under different window sizes, and selecting the window size corresponding to the highest estimated value precision according to the mahalanobis distance as the optimal window size;
and the real-time estimation module is used for estimating the state value of the dissolved oxygen concentration of the sewage in real time by using an unbiased finite impulse response filter in a window interval corresponding to the optimal window size.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310772649.4A CN116679026B (en) | 2023-06-27 | 2023-06-27 | Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310772649.4A CN116679026B (en) | 2023-06-27 | 2023-06-27 | Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116679026A true CN116679026A (en) | 2023-09-01 |
| CN116679026B CN116679026B (en) | 2024-07-12 |
Family
ID=87789089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310772649.4A Active CN116679026B (en) | 2023-06-27 | 2023-06-27 | Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116679026B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116947196A (en) * | 2023-09-20 | 2023-10-27 | 深圳市盘古环保科技有限公司 | Advanced oxidation industrial wastewater treatment system based on digitalization |
Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070185693A1 (en) * | 2003-08-11 | 2007-08-09 | Japan Science And Technology Agency | System estimation method and program, recording medium, and system estimation device |
| CN101179707A (en) * | 2007-09-21 | 2008-05-14 | 清华大学 | Wireless network video image multi-visual angle cooperation target tracking measurement method |
| CN102033933A (en) * | 2010-12-17 | 2011-04-27 | 南方医科大学 | Distance metric optimization method for maximizing mean average precision (MAP) |
| CN103136513A (en) * | 2013-02-05 | 2013-06-05 | 山东神思电子技术股份有限公司 | Improved automatic storage management (ASM) facial feature point locating method |
| CN103700114A (en) * | 2012-09-27 | 2014-04-02 | 中国航天科工集团第二研究院二O七所 | Complex background modeling method based on variable Gaussian mixture number |
| CN106228539A (en) * | 2016-07-12 | 2016-12-14 | 北京工业大学 | Multiple geometric primitive automatic identifying method in a kind of three-dimensional point cloud |
| CN106321318A (en) * | 2016-08-18 | 2017-01-11 | 河南职业技术学院 | Starting control system and control method of automobile engine |
| CN107291813A (en) * | 2017-05-19 | 2017-10-24 | 南京邮电大学 | Exemplary search method based on semantic segmentation scene |
| CN107332691A (en) * | 2017-06-08 | 2017-11-07 | 西安电子科技大学 | A kind of method for detecting failure nodes of wireless sensor network |
| US20180269955A1 (en) * | 2017-03-17 | 2018-09-20 | Aireon Llc | Receiver with non-coherent matched filter |
| CN108804784A (en) * | 2018-05-25 | 2018-11-13 | 江南大学 | A kind of instant learning soft-measuring modeling method based on Bayes's gauss hybrid models |
| CN110262431A (en) * | 2019-07-03 | 2019-09-20 | 江南大学 | A kind of sewage disposal process multiobjective optimization control method based on operating mode's switch |
| CN110455517A (en) * | 2019-08-22 | 2019-11-15 | 苏州旋械感知信息科技有限公司 | A kind of tower health monitor method of wind power generating set |
| CN113918873A (en) * | 2021-10-28 | 2022-01-11 | 江南大学 | Estimation method of sewage dissolved oxygen concentration, storage medium, electronic device and system |
| CN114048602A (en) * | 2021-11-11 | 2022-02-15 | 杭城智慧科技研究院(杭州)有限公司 | Aeration quantity estimation method and system based on finite impulse response filtering algorithm |
| CN115048985A (en) * | 2022-05-17 | 2022-09-13 | 国网浙江省电力有限公司嘉兴供电公司 | Electrical equipment fault discrimination method |
| CN115235769A (en) * | 2022-07-11 | 2022-10-25 | 北京华能新锐控制技术有限公司 | Fan bearing fault diagnosis method based on Mahalanobis distance compensation factor |
| CN115290080A (en) * | 2022-07-08 | 2022-11-04 | 杭州天眼信息咨询有限公司 | Pipeline robot positioning method based on unbiased finite impulse response filter |
| CN115290079A (en) * | 2022-07-08 | 2022-11-04 | 杭州天眼信息咨询有限公司 | Robot positioning method and system based on minimum variance unbiased finite impulse response |
| CN115435845A (en) * | 2022-10-21 | 2022-12-06 | 江苏大学 | Intelligent wide-area water quality monitoring and analyzing system and method based on unmanned ship and integrated learning and water quality analyzing method |
| CN115869065A (en) * | 2022-11-23 | 2023-03-31 | 中国科学院深圳先进技术研究院 | Method for estimating biomechanical parameters of in-vivo organ tissues |
| CN116774259A (en) * | 2023-07-10 | 2023-09-19 | 上海海事大学 | Satellite navigation signal processing method based on improved center difference Kalman filter, filter and receiver |
| CN116796110A (en) * | 2023-07-05 | 2023-09-22 | 云南南滇古茶供应链有限公司 | Tea leaf storage environment state estimation method, system and storage medium |
-
2023
- 2023-06-27 CN CN202310772649.4A patent/CN116679026B/en active Active
Patent Citations (23)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070185693A1 (en) * | 2003-08-11 | 2007-08-09 | Japan Science And Technology Agency | System estimation method and program, recording medium, and system estimation device |
| CN101179707A (en) * | 2007-09-21 | 2008-05-14 | 清华大学 | Wireless network video image multi-visual angle cooperation target tracking measurement method |
| CN102033933A (en) * | 2010-12-17 | 2011-04-27 | 南方医科大学 | Distance metric optimization method for maximizing mean average precision (MAP) |
| CN103700114A (en) * | 2012-09-27 | 2014-04-02 | 中国航天科工集团第二研究院二O七所 | Complex background modeling method based on variable Gaussian mixture number |
| CN103136513A (en) * | 2013-02-05 | 2013-06-05 | 山东神思电子技术股份有限公司 | Improved automatic storage management (ASM) facial feature point locating method |
| CN106228539A (en) * | 2016-07-12 | 2016-12-14 | 北京工业大学 | Multiple geometric primitive automatic identifying method in a kind of three-dimensional point cloud |
| CN106321318A (en) * | 2016-08-18 | 2017-01-11 | 河南职业技术学院 | Starting control system and control method of automobile engine |
| US20180269955A1 (en) * | 2017-03-17 | 2018-09-20 | Aireon Llc | Receiver with non-coherent matched filter |
| CN107291813A (en) * | 2017-05-19 | 2017-10-24 | 南京邮电大学 | Exemplary search method based on semantic segmentation scene |
| CN107332691A (en) * | 2017-06-08 | 2017-11-07 | 西安电子科技大学 | A kind of method for detecting failure nodes of wireless sensor network |
| CN108804784A (en) * | 2018-05-25 | 2018-11-13 | 江南大学 | A kind of instant learning soft-measuring modeling method based on Bayes's gauss hybrid models |
| CN110262431A (en) * | 2019-07-03 | 2019-09-20 | 江南大学 | A kind of sewage disposal process multiobjective optimization control method based on operating mode's switch |
| CN110455517A (en) * | 2019-08-22 | 2019-11-15 | 苏州旋械感知信息科技有限公司 | A kind of tower health monitor method of wind power generating set |
| CN113918873A (en) * | 2021-10-28 | 2022-01-11 | 江南大学 | Estimation method of sewage dissolved oxygen concentration, storage medium, electronic device and system |
| CN114048602A (en) * | 2021-11-11 | 2022-02-15 | 杭城智慧科技研究院(杭州)有限公司 | Aeration quantity estimation method and system based on finite impulse response filtering algorithm |
| CN115048985A (en) * | 2022-05-17 | 2022-09-13 | 国网浙江省电力有限公司嘉兴供电公司 | Electrical equipment fault discrimination method |
| CN115290080A (en) * | 2022-07-08 | 2022-11-04 | 杭州天眼信息咨询有限公司 | Pipeline robot positioning method based on unbiased finite impulse response filter |
| CN115290079A (en) * | 2022-07-08 | 2022-11-04 | 杭州天眼信息咨询有限公司 | Robot positioning method and system based on minimum variance unbiased finite impulse response |
| CN115235769A (en) * | 2022-07-11 | 2022-10-25 | 北京华能新锐控制技术有限公司 | Fan bearing fault diagnosis method based on Mahalanobis distance compensation factor |
| CN115435845A (en) * | 2022-10-21 | 2022-12-06 | 江苏大学 | Intelligent wide-area water quality monitoring and analyzing system and method based on unmanned ship and integrated learning and water quality analyzing method |
| CN115869065A (en) * | 2022-11-23 | 2023-03-31 | 中国科学院深圳先进技术研究院 | Method for estimating biomechanical parameters of in-vivo organ tissues |
| CN116796110A (en) * | 2023-07-05 | 2023-09-22 | 云南南滇古茶供应链有限公司 | Tea leaf storage environment state estimation method, system and storage medium |
| CN116774259A (en) * | 2023-07-10 | 2023-09-19 | 上海海事大学 | Satellite navigation signal processing method based on improved center difference Kalman filter, filter and receiver |
Non-Patent Citations (3)
| Title |
|---|
| IGOR MELNYKOV等: "On K-means algorithm with the use of Mahalanobis distances", STATISTICS AND PROBABILITY LETTERS, vol. 84, 31 December 2014 (2014-12-31), pages 88 - 95, XP028781964, DOI: 10.1016/j.spl.2013.09.026 * |
| 吴晓东等: "基于多块信息提取和马氏距离的k近邻故障监测", 《信息与控制》, vol. 50, no. 03, 25 December 2020 (2020-12-25), pages 287 - 296 * |
| 朱占龙;王军芬;: "基于自适应模糊C均值与后处理的图像分割算法", 激光与光电子学进展, no. 01, 17 August 2017 (2017-08-17) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116947196A (en) * | 2023-09-20 | 2023-10-27 | 深圳市盘古环保科技有限公司 | Advanced oxidation industrial wastewater treatment system based on digitalization |
| CN116947196B (en) * | 2023-09-20 | 2023-12-08 | 深圳市盘古环保科技有限公司 | Advanced oxidation industrial wastewater treatment system based on digitalization |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116679026B (en) | 2024-07-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106873379B (en) | Sewage treatment optimal control method based on iterative ADP algorithm | |
| CN108469507B (en) | Effluent BOD soft measurement method based on self-organizing RBF neural network | |
| CN111291937A (en) | Method for predicting quality of treated sewage based on combination of support vector classification and GRU neural network | |
| Jeppsson et al. | Reduced order models for on-line parameter identification of the activated sludge process | |
| Pisa et al. | LSTM-based wastewater treatment plants operation strategies for effluent quality improvement | |
| CN111650834B (en) | Sewage treatment process prediction control method based on extreme learning machine | |
| CN116679026B (en) | Self-adaptive unbiased finite impulse response filtering sewage dissolved oxygen concentration estimation method | |
| CN111762958A (en) | Deep well aeration process optimization method and device for sewage treatment plant based on ASM2D model | |
| CN108536106B (en) | Aeration system dissolved oxygen online regulation and control method based on Kalman filtering-extreme learning machine | |
| CN104914227B (en) | Sewage quality flexible measurement method based on many gaussian kernel self-optimizing Method Using Relevance Vector Machine | |
| CN110542748B (en) | Knowledge-based robust effluent ammonia nitrogen soft measurement method | |
| CN103809436A (en) | Method for intelligent modeling sewage disposal process using activated sludge process | |
| CN115356930A (en) | Multi-objective optimization control system and method in sewage treatment process | |
| CN113918873B (en) | Method for estimating dissolved oxygen concentration in sewage, storage medium, electronic device, and system | |
| CN113608444A (en) | Sewage treatment control method based on self-adaptive prediction control | |
| CN113721468B (en) | Self-healing control method and system for additive faults of sewage treatment system actuator | |
| Civelekoglu et al. | Modeling Carbon and Nitrogen Removal in an Industrial Wastewater Treatment Plant Using an Adaptive Network‐Based Fuzzy Inference System | |
| CN111399558B (en) | Knowledge selection-based multi-objective optimization control method for sewage treatment process | |
| Qiao et al. | Recurrent neural network-based control for wastewater treatment process | |
| Kim et al. | Genetic algorithms for the application of Activated Sludge Model No. 1 | |
| Sun et al. | Enhanced MPC based on unknown state estimation and control compensation | |
| Frącz | Nonlinear modeling of activated sludge process using the Hammerstein-Wiener structure | |
| CN112967761B (en) | Sewage dephosphorization and dosing calculation method and medium based on self-organizing fuzzy neural network | |
| Borowa et al. | Modelling of wastewater treatment plant for monitoring and control purposes by state-space wavelet networks | |
| CN114781166A (en) | Sewage treatment process soft measurement method based on weighted probability slow characteristic model |
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 |