CN103868198B - Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise - Google Patents

Abstract

The invention discloses Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise, the room air conditioner of different Energy Efficiency Ratio grades Energy Efficiency Ratio test result under cooling condition can be modified: gather 14 groups of room air conditioner Energy Efficiency Ratio data under operating condition of test, by each constant coefficient value in Bessel Formula certain Energy Efficiency Ratio regression equation, so that it is determined that Energy Efficiency Ratio drift correction formula; Gather model machine Energy Efficiency Ratio under any wet and dry bulb temperature in national standard, determine Energy Efficiency Ratio drift correction value by Energy Efficiency Ratio drift correction formula, thus obtaining the actual measurement operating mode recurrence correction value to nominal condition. The present invention provides a kind of Energy Efficiency Ratio Acquisition Error modification method being applied to performance of room air conditioners test and appraisal, the method can improve test and appraisal accuracy, reflect performance of room air conditioners really, it is to avoid manual control room air conditioner efficiency grade occurs predictable mistake to divide; Can instructing the setting of the air-conditioning equipment operational factors such as convertible frequency air-conditioner, the adjustment of space air state is more accurate simultaneously.

Description

Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise

Technical field

The present invention relates to a kind of Energy Efficiency Ratio Acquisition Error modification method, a kind of method being applied to room air conditioner Energy Efficiency Ratio Acquisition Error correction specifically.

Background technology

Room air conditioner Energy Efficiency Standard is the Energy Efficiency Standard of the electrical equipment that China works out the earliest, and front and back experienced by three revisions, and the related specifications of the U.S., European Union and Japan has been followed in the formulation of Energy Efficiency Standard and revision. The revision of 2000 and twice Energy Efficiency Standard in 2004 is mainly reflected in the rise of energy efficiency market, and existing room air conditioner Energy Efficiency Standard, its efficiency grade is changed to three grades of efficiencies of existing GB GB12021.3-2010 by the Pyatyi efficiency of former GB GB12021.3-2004: existing Energy Efficiency Standard directly delete three in former Energy Efficiency Standard, four and Pyatyi efficiency grade, and former Energy Efficiency Standard I and II efficiency grade is newly defined as the two of existing Energy Efficiency Standard, three grades of efficiency grades. Based on this when unique step, the one-level first order as working standard that two grades of efficiencies of existing Energy Efficiency Standard are boosted.

In order to be appropriately determined the efficiency grade of room air conditioner, introducing thermodynamics consummating degree room air conditioner is carried out efficiency evaluation in " the thermodynamics consummating degree analysis of room air conditioner " literary composition, this is a kind of room air conditioner energy efficiency analysis method for air merited attention. Meanwhile, room air conditioner Energy Efficiency Standard updates revision, and accordingly, its performance standard is also required to revision, and matches, and namely avoids because measuring the former mistake division thus resulting in room air conditioner efficiency grade on to greatest extent.

Specified refrigerating capacity and Energy Efficiency Ratio need to reflect room air conditioner performance under nominal condition point, and this is also the foundation of room air conditioner efficiency ranking. When Cooling Capacity For The Room Air Conditioner and Energy Efficiency Ratio are surveyed, the impact of tested person condition and anthropic factor, the actual condition point of performance of room air conditioners test often offsets in tolerance. The franchise respectively ± 1 DEG C of correlation standard wet and dry bulb temperature and ± 0.5 DEG C.This test result often leading to refrigerating capacity and Energy Efficiency Ratio is not at nominal condition point, but records under a certain operating point in wet and dry bulb tolerance. Based on this test result, room air conditioner is carried out efficiency ranking, it is possible to the mistake division to room air conditioner efficiency grade can be caused. On the one hand, can reach under nominal condition (not up to) room air conditioner of three grades of energy efficiency market, it is likely to be due to actual measurement condition and is positioned at wet and dry bulb temperature tolerance lower limit (upper limit), and be divided into by mistake and can not reach (reaching) three grades of energy efficiency market, and three grades of energy efficiency market are the threshold values of the room air conditioner market access. On the other hand, along with instrument and meter precision improves constantly, the actual measurement condition point of manual control is positioned at wet and dry bulb temperature high tolerance or lower limit, it is also possible to cause that room air conditioner can occur predictable mistake to divide.

Summary of the invention

Technical problem: for the problems referred to above, the present invention provides one to be applied to room air conditioner, can obtain wet and dry bulb temperature in national standard tolerance under any operating point Energy Efficiency Ratio to the recurrence correction value of nominal condition, improve Energy Efficiency Ratio modification method in the wet and dry bulb temperature franchise of test and appraisal accuracy.

Technical scheme: Energy Efficiency Ratio modification method in the wet and dry bulb temperature franchise of the present invention, comprises the following steps:

1) the room air conditioner Energy Efficiency Ratio data under operating condition of test are gathered, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition is Δ E_i(i=1,2,14), i represents different measurement condition:

Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1～5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6～10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;

3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;

2) by described step 1) in gather each floor data substitute in following regression equation respectively, obtain the over-determined systems of 14 hexa-atomic linear function compositions, then solve described equation group with Bessel Formula for criterion, obtain every constant coefficient E of regression equation₀, a, b, c, d, e value:

$\mathrm{\Δ}E(t_1,t_2)=E_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t₁Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode₂Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;

E₀Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature₁Coefficient, b represent temperature difference item t of indoor wet bulb temperature₂Coefficient, c be indoor dry-bulb temperature secondary temperature difference item t² ₁Coefficient, d represent indoor wet bulb temperature secondary temperature difference item t² ₂Coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t₁t₂Coefficient;

3) the Energy Efficiency Ratio drift correction amount in wet and dry bulb temperature tolerance, under any practical running operating point, it is possible to calculate according to following formula:

ΔE(t₁,t₂)=E₀+ζΛQt+(Qt)^TΛ Qt;

Wherein, t is temperature difference vector, t=(t₁,t₂)^T, T represents the transposition of vector;

Λ is weight matrix, Λ=diag (θ₁, θ₂), diag represents diagonal matrix, θ₁Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ₂Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\α}}_{11}& {\mathrm{\α}}_{12}\\ {\mathrm{\α}}_{21}& {\mathrm{\α}}_{22}\end{array}\right),$ α₁₁Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₁Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α₁₂Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₂Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;

ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ₁, ζ₂); ζ₁Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ₂Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;

4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains₁,t₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:

E_t=E_m-ΔE(t₁,t₂);

Wherein, E_mFor the Energy Efficiency Ratio measured value under any operating mode of room air conditioner tested in wet and dry bulb tolerance; Δ E (t₁,t₂) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, E_tCorrection value is returned for Energy Efficiency Ratio under this operating mode.

Energy Efficiency Ratio measured value under operating mode any in wet and dry bulb tolerance is modified by the present invention, obtain the recurrence correction value to nominal condition of the Energy Efficiency Ratio under this operating mode, can be used for the actual operation parameters adjustment of convertible frequency air-conditioner: when Energy Efficiency Ratio returns correction value more than Energy Efficiency Ratio measured value, suitably turn down compressor operating frequency; When Energy Efficiency Ratio returns correction value less than Energy Efficiency Ratio measured value, suitably heighten compressor operating frequency; When Energy Efficiency Ratio returns correction value equal to Energy Efficiency Ratio measured value, it is not necessary to adjust.

In the inventive method, described step 2) in Bessel Formula be:

ε=k σ, wherein ε represents monitoring residual error, and k represents monitoring coefficient, and σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation, $\mathrm{\σ}=\sqrt{\frac{\underset{i=1}{\overset{14}{\Σ}}\[{\mathrm{\ΔE}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\Σ}}{\mathrm{\ΔE}}_i\]}{13}},(i=1,2,\mathrm{...},14);$

Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has at least the Energy Efficiency Ratio actual measurement deviation of 50% to meet relational expression ε_i≤ ε; As k=2, the Energy Efficiency Ratio actual measurement deviation of 95% is had at least to meet relational expression ε_i≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had at least to meet relational expression ε_i≤ ε;

ε_iRepresent the residual error of Energy Efficiency Ratio measurement error and equation deviation, ε under different operating point_i=| Δ E_i-ΔE_i(t₁,t₂) |; Δ E_iRepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ E_i(t₁,t₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different measurement condition, i=1,2, and, 14.

In the inventive method, described step 3) in,

${\mathrm{\α}}_{11}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\α}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\α}}_{21}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\α}}_{22}=(c-d+\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\ζ}}_1=\frac{({\mathrm{a\α}}_{22}-{\mathrm{b\α}}_{12})}{{\mathrm{\θ}}_1({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})},{\mathrm{\ζ}}_2=\frac{({\mathrm{b\α}}_{11}-{\mathrm{a\α}}_{21})}{{\mathrm{\θ}}_2({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})};$

Wherein, λ₁、λ₂For the eigenvalue of secondary temperature difference term coefficient matrix B, p₁、p₂For the characteristic vector of secondary temperature difference term coefficient matrix B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),{\mathrm{\&lambda;}}_1=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d+\sqrt{{(c-d)}^2+e^2}}{2},$ $\left|{\mathrm{\&lambda;}}_1\right|<\left|{\mathrm{\&lambda;}}_2\right|;\left|\right|p_1\left|\right|={(e^2+{\&lsqb;c-d-\sqrt{{(c-d)}^2+e^2}\&rsqb;}^2)}^{\frac{1}{2}},\left|\right|p_2\left|\right|={(e^2+{\&lsqb;c-d+\sqrt{{(c-d)}^2+e^2}\&rsqb;}^2)}^{\frac{1}{2}}).$

Beneficial effect: the present invention compared with prior art, has the advantage that

Energy Efficiency Ratio measured value in the wet and dry bulb temperature tolerance of the inventive method application finite discrete, describes the Energy Efficiency Ratio continuous print variation characteristic about wet and dry bulb temperature.

The standard deviation of deviation function deviation is surveyed for object with Energy Efficiency Ratio, with the numerical relation of test value and the residual sum standard deviation of equation value for target, each measurement condition is clicked on Mobile state monitoring, to determine the specific features and accuracy thereof that Energy Efficiency Ratio deviation changes with wet-bulb depression seriality.

The inventive method can obtain wet and dry bulb temperature, and in national standard tolerance, under any operating point, Energy Efficiency Ratio is to the recurrence correction value of nominal condition, thus instructing the setting of the air-conditioning equipment carrying out practically parameters such as convertible frequency air-conditioner, equipment runs more accurate;

The inventive method can obtain wet and dry bulb temperature in national standard tolerance under any operating point Energy Efficiency Ratio to the correction value of nominal condition, improve test and appraisal accuracy, reflect performance of room air conditioners more really, and provide reference to the revision of national standard;

Energy Efficiency Ratio drift correction formula is once determined, the room air conditioner Energy Efficiency Ratio measured value of follow-up identical Energy Efficiency Ratio grade is determined either directly through this correction formula, it is to avoid the huge consumption of the human and material resources that repeatedly repeated trials brings;

Along with instrument and meter precision improves constantly, it is to avoid the actual measurement condition point of manual control is positioned at wet and dry bulb temperature high tolerance or lower limit, thus causing that room air conditioner efficiency grade occurs predictable mistake to divide.

Accompanying drawing explanation

Fig. 1 is the process step figure of the inventive method.

Detailed description of the invention

Below in conjunction with drawings and Examples, the present invention will be further described:

In the wet and dry bulb temperature franchise of the present invention, the room air conditioner of different Energy Efficiency Ratio grades Energy Efficiency Ratio test result under cooling condition can be modified by Energy Efficiency Ratio modification method. It is modified to embodiment with the Energy Efficiency Ratio measured result under three class room air conditioner cooling conditions below to illustrate, specifically comprises the following steps that

Step 1) gather the room air conditioner Energy Efficiency Ratio data under operating condition of test, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition is

ΔE₁=-0.03267kW`kW<sup>-1</sup>、ΔE<sub>2</sub>=-0.01971kW`kW^-1、ΔE₃=0kW`kW^-1、

ΔE₄=0.0329kW`kW<sup>-1</sup>、ΔE<sub>5</sub>=0.05315kW`kW^-1、ΔE₆=-0.02574kW`kW^-1、

ΔE₇=-0.0094kW`kW<sup>-1</sup>、ΔE<sub>8</sub>=0kW`kW^-1、ΔE₉=0.03773kW`kW^-1、

ΔE₁₀=0.05315kW`kW<sup>-1</sup>、ΔE<sub>11</sub>=-0.06462kW`kW^-1、ΔE₁₂=0.01385kW`kW^-1、

ΔE₁₃=-0.05007kW`kW<sup>-1</sup>、ΔE<sub>14</sub>=0.10285kW`kW^-1:

Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1～5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6～10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;

3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;

Step 2) by described step 1) in gather each floor data substitute into following regression equation respectively:

$\mathrm{\&Delta;}E(t_1,t_2)=E_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t₁Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode₂Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;

E₀Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature₁Coefficient, b represent temperature difference item t of indoor wet bulb temperature₂Coefficient, c be indoor dry-bulb temperature secondary temperature difference itemCoefficient, d represent indoor wet bulb temperature secondary temperature difference itemCoefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t₁t₂Coefficient;

The 3rd, under 6 operating modes, t₁、t₂、t₁t₂、ΔE_i(t₁,t₂) value be 0, therefore obtain the over-determined systems of 12 hexa-atomic linear functions composition:

E₀-a+c=-0.02574;

E₀-0.5a+0.25c=-0.01278;

E₀+ 0.5a+0.25c=0.0329;

E₀+ a+c=0.05315;

E₀-0.5b+0.25d=-0.02574;

E₀-0.2b+0.04=-0.0094;

E₀+ 0.3b+0.09d=0.03773;

E₀+ 0.5b+0.25d=0.05315;

E₀-a-0.5b+c+0.25d+0.5e=-0.06462;

E₀-a+0.5b+c+0.25d-0.5e=0.01385;

E₀+ a-0.5b+c+0.25d-0.5e=-0.05007;

E₀+ a+0.5b+c+0.25d+0.5e=0.10285

Solve above-mentioned equation group, obtain every constant coefficient E of regression equation₀, a, b, c, d, e value respectively 0.002,0.044,0.088,0.0049,0.019,0.0069;

Meanwhile, the convergence of Bessel Formula checking Energy Efficiency Ratio regression equation is adopted:

The residual error of test value and equation value is ε_i=| Δ E_i-ΔE_i(t₁,t₂) |, Δ E_iRepresent the difference of Energy Efficiency Ratio measured value and nominal value, Δ E_i(t₁,t₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; It is respectively as follows: ε according to the residual error that above-mentioned formula tries to achieve each operating point corresponding₁=0.00245, ε₂=0.00121, ε₃=0.000059, ε₄=0.00106, ε₅=0.001304, ε₆=0.00357, ε₇=0.00372, ε₈=0.00023, ε₉=0.00094, ε₁₀=0.00128, ε₁₁=0.00275, ε₁₂=0.0024, ε₁₃=0.00118 ε₁₄=0.00132

σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation, $\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\&Sigma;}}{\&lsqb;{\mathrm{\&Delta;E}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\&Sigma;}}{\mathrm{\&Delta;E}}_i\&rsqb;}^2}{13}}=0.002,(i=1,2,\mathrm{...},14);$

Bessel Formula is ε=k σ, respectively when monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has the Energy Efficiency Ratio actual measurement deviation of 57% to meet relational expression ε_i≤ ε, (i represents different measurement condition, i=1,2, and, 14); As k=2, the Energy Efficiency Ratio actual measurement deviation of 100% is had to meet relational expression ε_i≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had to meet relational expression ε_i≤ ε. Therefore the convergence of Energy Efficiency Ratio regression equation, W₀, a, b, c, d, e value respectively 0.0019,0.19,0.37,0.0027,0.011,0.0098;

$\mathrm{\&Delta;}W(t_1,t_2)=W_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2$

$=0.002+0.044t_1+0.088t_2+0.0049{t_1}^2+0.019{t_2}^2+0.0069t_1{t}_2$

Step 3) in wet and dry bulb temperature tolerance, a room air conditioner model machine (specified refrigerating capacity and Energy Efficiency Ratio respectively 12kW, 3.32kW kW of the present invention^-1) actual measurement Energy Efficiency Ratio E_m=3.114kW`kW^-1(operating condition of test point: indoor dry-bulb temperature is 26.7 DEG C, wet bulb temperature is 18.9 DEG C), then the Energy Efficiency Ratio drift correction amount under this operating mode, it is possible to calculate according to following formula:

ΔE(t₁,t₂)=E₀+ζΛQt+(Qt)^TΛ Qt;

Wherein, t is temperature difference vector, t=(t₁,t₂)^T, T represents the transposition of vector;

Λ is weight matrix, Λ=diag (θ₁, θ₂), diag represents diagonal matrix, θ₁Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ₂Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α₁₁Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₁Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α₁₂Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₂Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;

ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ₁, ζ₂); ζ₁Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ₂Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;

θ₁、θ₂、α₁₁、α₁₂、α₂₁、α₂₂、ζ₁、ζ₂Determine by the following method:

Coefficient c, d, e of secondary temperature term in Energy Efficiency Ratio regression equation are carried out quadratic form and changes to obtain Second Order with Constant Coefficients matrix $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),$ Determine its eigenvalue λ₁、λ₂And characteristic of correspondence vector p₁、p₂:

${\mathrm{\&lambda;}}_1=\frac{c+d+\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2},$

$p_1={(e,c-d-\sqrt{{(c-d)}^2+e^2})}^T,p_2={(e,c-d+\sqrt{{(c-d)}^2+e^2})}^T;$

Then, ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},$ Weight matrix Λ=diag (θ₁, θ₂); ${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,$ ${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d+\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,$ Temperature transition matrix $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right);{\mathrm{\&zeta;}}_1=\frac{({\mathrm{a\&alpha;}}_{22}-{\mathrm{b\&alpha;}}_{12})}{{\mathrm{\&theta;}}_1({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},{\mathrm{\&zeta;}}_2=\frac{({\mathrm{b\&alpha;}}_{11}-{\mathrm{a\&alpha;}}_{21})}{{\mathrm{\&theta;}}_2({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},$ Temperature difference item is relative to high order temperature difference item relative effect coefficient vector ζ=(ζ to Energy Efficiency Ratio₁, ζ₂);

Namely Energy Efficiency Ratio drift correction formula is:

$\begin{array}{c}\mathrm{\&Delta;}E(t_1,t_2)=0.002+(0.84,0.77)\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\\ {\left(\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\right)}^T\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.16& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\end{array}$

The Energy Efficiency Ratio drift correction amount under this operating mode is calculated according to above-mentioned solving result:

$\begin{array}{c}\mathrm{\&Delta;}E(26.7,18.9)=0.002+(0.84,0.77)\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\\ +{\left(\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)\right)}^T\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)=-0.01916lW`{\mathrm{kW}}^{-1};\end{array}$

Step 4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains₁,t₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:

E_t=E_m-ΔE(t₁,t₂)=3.0953kW`kW<sup>-1</sup>+0.01916kW`kW^-1=3.1145kW`kW^-1

Namely dry-bulb temperature be 26.7 DEG C, wet bulb temperature be that at 18.9 DEG C, Energy Efficiency Ratio returns correction value be 3.1145kW`kW^-1。

Wherein, E_mFor the Energy Efficiency Ratio measured value under any operating mode of tested room air conditioner, for 3.0953kW`kW<sup>-1</sup>; Δ E (t<sub>1</sub>,t<sub>2</sub>) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, for-0.01916kW`kW^-1; E_tCorrection value is returned for Energy Efficiency Ratio under this operating mode.

Claims (3)

1. Energy Efficiency Ratio modification method in a wet and dry bulb temperature franchise, it is characterised in that the method comprises the following steps:

1) the room air conditioner Energy Efficiency Ratio data under operating condition of test are gathered, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition point is Δ E_i(i=1,2,14), i represents different measurement condition:

Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1～5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6～10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;

3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;

2) by described step 1) in gather each floor data substitute in following regression equation respectively, obtain the over-determined systems of 14 hexa-atomic linear function compositions, then solve described equation group with Bessel Formula for criterion, obtain every constant coefficient E of regression equation₀, a, b, c, d, e value:

$\mathrm{\&Delta;}E(t_1,t_2)=E_0+{\mathrm{at}}_1+{\mathrm{bt}}_2+{\mathrm{ct}}_1^2+{\mathrm{dt}}_2^2+{\mathrm{et}}_1{t}_2,$

Wherein, t₁Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode₂Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;

E₀Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature₁Coefficient, b represent temperature difference item t of indoor wet bulb temperature₂Coefficient, c be indoor dry-bulb temperature secondary temperature difference item t² ₁Coefficient, d represent indoor wet bulb temperature secondary temperature difference item t² ₂Coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t₁t₂Coefficient;

3) the Energy Efficiency Ratio drift correction amount in wet and dry bulb temperature tolerance, under any practical running operating point, it is possible to calculate according to following formula:

ΔE(t₁,t₂)=E₀+ζΛQt+(Qt)^TΛ Qt;

Wherein, t is temperature difference vector, t=(t₁,t₂)^T, T represents the transposition of vector;

Λ is weight matrix, Λ=diag (θ₁, θ₂), diag represents diagonal matrix, θ₁Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ₂Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α₁₁Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₁Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α₁₂Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α₂₂Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;

ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ₁, ζ₂); ζ₁Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ₂Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;

4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains₁,t₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:

E_t=E_m-ΔE(t₁,t₂);

Wherein, E_mFor the Energy Efficiency Ratio measured value under any operating mode of room air conditioner tested in wet and dry bulb tolerance; Δ E (t₁,t₂) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, E_tCorrection value is returned for Energy Efficiency Ratio under this operating mode.

2. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, it is characterised in that described step 2) in Bessel Formula be: ε=k σ, wherein ε represents monitoring residual error, k represents monitoring coefficient, and σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation $\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\&Sigma;}}{\&lsqb;{\mathrm{\&Delta;E}}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\&Sigma;}}{\mathrm{\&Delta;E}}_i\&rsqb;}^2}{13}},(i=1,2,...,14);$

Described Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has at least the Energy Efficiency Ratio actual measurement deviation of 50% to meet relational expression ε_i≤ ε; As k=2, the Energy Efficiency Ratio actual measurement deviation of 95% is had at least to meet relational expression ε_i≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had at least to meet relational expression ε_i≤ ε;

ε_iRepresent the measurement error of Energy Efficiency Ratio under different operating mode and the residual error of equation deviation, ε_i=| Δ E_i-ΔE_i(t₁,t₂) |; Δ E_iRepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ E_i(t₁,t₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different measurement condition, i=1,2,14.

3. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, it is characterised in that described step 3) in, ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|};$

${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d-\sqrt{{(c-d)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\&zeta;}}_1=\frac{\left({\mathrm{a\&alpha;}}_{22}-{\mathrm{b\&alpha;}}_{12}\right)}{{\mathrm{\&theta;}}_1\left({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21}\right)},{\mathrm{\&zeta;}}_2=\frac{\left({\mathrm{b\&alpha;}}_{11}-{\mathrm{a\&alpha;}}_{21}\right)}{{\mathrm{\&theta;}}_2\left({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21}\right)};$

Wherein, λ₁、λ₂For the eigenvalue of secondary temperature difference term coefficient matrix B, p₁、p₂For the characteristic vector of secondary temperature difference term coefficient matrix B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),{\mathrm{\&lambda;}}_1=\frac{c+d-\sqrt{{(c-d)}^2+e^2}}{2},{\mathrm{\&lambda;}}_2=\frac{c+d+\sqrt{{(c-d)}^2+e^2}}{2},$ |λ₁| < | λ₂|; $\left|\right|p_1\left|\right|={(e^2+{\&lsqb;c-d-\sqrt{{(c-d)}^2+e^2}\&rsqb;}^2)}^{\frac{1}{2}},\left|\right|p_2\left|\right|={(e^2+{\&lsqb;c-d+\sqrt{{(c-d)}^2+e^2}\&rsqb;}^2)}^{\frac{1}{2}}).$