CN105913199A - Energy, economy and environment coordination degree calculation method and device based on ellipsoidal model - Google Patents

Abstract

The present invention provides a method and device for calculating energy, economy, and environmental coordination based on an ellipsoid model, and relates to the field of regional information evaluation. The method provided by the present invention establishes an ellipsoid analysis model in advance and uses the ellipsoid The model adjusts for the degree of coordination between the environment and the energy-economy in the case of coordination between the energy and the economy, the degree of coordination between the economy and the energy-economy in the case of coordination between the energy and the environment system, and the degree of coordination between the economy and the environment system Harmony of energy and economy-environment. And finally make the calculated energy, economy and environment coordination degree more accurate.

Description

Calculation method and device for energy, economy and environment coordination degree based on ellipsoid model

technical field

The invention relates to the field of regional information evaluation, in particular to a calculation method and device for energy, economy and environment coordination degree based on an ellipsoid model.

Background technique

With the rapid development of social economy and the advancement of industrialization, the problems of resource shortage and environmental pollution have become increasingly prominent, and people are paying more and more attention to the coordinated development of energy-economy-environment (Energy-Economy-Environment, 3E). Energy is the basis and guarantee for national economic and social development, and it is also the focus of international politics, economy, military and diplomacy; the economic base determines the superstructure and is the center of all development issues; the environment is an important basis for key to sustainable development.

Among them, how to determine the coordination degree of energy, economy and environment more reasonably and accurately is the basis of other theoretical researches. Therefore, how to accurately determine the coordination degree of energy, economy and environment has become the focus.

Since the 3E system is a large and complex system, it is difficult to describe the real development level of each subsystem with a single index linear structure. Therefore, some specific mathematical models are used in related technologies to assist in the analysis and calculation of the degree of coordination. Specifically, for example, in related technologies, a cube model is used to analyze and calculate the degree of coordination, but this analysis and calculation method is still not accurate enough.

That is, in the related art, the method of calculating the coordination degree of energy, economy and environment is not accurate enough.

Contents of the invention

The object of the present invention is to provide a method and device for calculating the degree of coordination of energy, economy and environment based on an ellipsoid model, so as to improve the accuracy of calculating the degree of coordination of energy, economy and environment.

In the first aspect, an embodiment of the present invention provides a calculation method for energy, economy and environment coordination degree based on an ellipsoid model, including:

According to the statistical table of historical data of the target area, analyze and obtain the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area;

Calculate the energy, economy and environment coordination degree of the target region according to the following formula,

$\mathrm{<span\; class="notranslate">\; E.</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&alpha;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&beta;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&gamma;</span>}}{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span>}<span\; class="notranslate">\; ;</span>$

in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region.

In combination with the first aspect, the embodiment of the present invention provides the first possible implementation mode of the first aspect, wherein, according to the statistical tables of historical data of the target region, the energy development speed V _a and the economic development speed V a of the target region are analyzed and obtained. _b and the environmental development velocity V _c include:

Calculate the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area according to the following formula:

${\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; A</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; C</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ;</span>$

Among them, t is the statistical year of the historical data statistics table, and A, B, and C are the curve fitting equations of the comprehensive development level values of energy, economy, and environment in the target area, respectively.

In combination with the first aspect, the embodiment of the present invention provides a second possible implementation of the first aspect, wherein, according to the historical data statistics table of the target region, the energy development speed V _a and the economic development speed V a of the target region are analyzed and obtained. _b and the speed of environmental development V _c also include: using the principal component analysis method to determine the economic standard index data in the economic system of the target region and the weight value of the economic standard index data, and to determine the energy standard index data in the energy system of the target region and the above-mentioned The weight value of the energy standard index data, and determine the weight value of the environmental standard index data and the environmental standard index data in the environmental system of the target region;

Using the method of weighted calculation, the value of the comprehensive economic development level of the target area is calculated according to the value of the determined economic standard index data and the weight value of the economic standard index data, and the value of the comprehensive economic development level of the target area is calculated by a fitting method The curve fitting equation of ;

Using the method of weighted calculation, according to the value of the determined energy standard index data and the weight value of the energy standard index data, the value of the comprehensive energy development level of the target area is calculated, and the value of the comprehensive development level of energy is calculated by a fitting method The curve fitting equation of ;

Using the method of weighted calculation, according to the value of the determined environmental standard index data and the weight value of the environmental standard index data, the value of the comprehensive environmental development level of the target area is calculated, and the value of the comprehensive development level of the environment is calculated by a fitting method The curve fitting equation for .

In combination with the first aspect, the embodiment of the present invention provides a third possible implementation manner of the first aspect, wherein, the energy development speed V _a and the economic development rate of the target region are analyzed and obtained according to the historical data statistics table of the target region. Velocity V _b and ambient development velocity V _c also include:

According to the following formulas, the values of the comprehensive development levels of energy, economy and environment are calculated respectively,

F=Z×T;

Among them, F is an n×1 matrix, each row represents the value of the comprehensive development level of energy, economy and environment in the corresponding year, n is the number of years obtained through the historical data of the target area; The energy, economic or environmental standard index data calculated from the statistical table of regional data over the years, T is a p×1 matrix, each row represents the weight value under the corresponding index of energy, economy or environment, and p is obtained through the historical data of the target region The number of indicators.

In combination with the first aspect, the embodiment of the present invention provides a fourth possible implementation manner of the first aspect, wherein the energy development speed V _a , the economic development Velocity V _b and ambient development velocity V _c also include:

Calculate the weight value of energy, economic or environmental standard index data according to the following formula,

$\mathrm{<span\; class="notranslate">\; T</span>}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; K</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}<span\; class="notranslate">\; \&times;</span>\frac{{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}}{\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; K</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}<span\; class="notranslate">\; ,</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}\mathrm{<span\; class="notranslate">\; K</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; ;</span>$

${\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; j</span>}}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; p</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; j</span>}}}{\sqrt{{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}}}<span\; class="notranslate">\; ,</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; p</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}\mathrm{<span\; class="notranslate">\; K</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; ;</span>$

Among them, T _j =(T _1j ,T _2j ,...T _ij ...T _pj )(j=1,2...K) means that among the components extracted from the historical data statistics table of the target region, the first K principal components The weight of j principal components, K the number of principal components, T is the weight value of energy, economic or environmental standard index data, λ _k is the characteristic root of the principal components, p=(p _1j ,p _2j ,…p _ij , …p _pj )(j=1,2,…K) is the factor loading matrix derived from the statistical tables of historical data in the target region.

In combination with the first aspect, the embodiment of the present invention provides a fifth possible implementation of the first aspect, wherein, according to the historical data statistics table of the target region, the energy development speed V _a and the economic development speed V a of the target region are analyzed and obtained. _b and the environmental development velocity V _c also include:

Calculate the economic standard indicator data through the following steps;

Extract the matrix containing the data in the original indicators from the historical data statistics table of the target region;

Standardize the data in the original indicators according to the following formula, Among them, X=(X _i1 ,X _i2 ,…X _ip )(i=1,2…,n) is a matrix containing original index data of energy, economy or environment, Z=(Z _i1 ,Z _i2 ,…Z _ip ) (i=1,2,...n) is a standardized matrix of energy, economic or environmental data;

Statistically standardize each data in the standardized matrix to obtain energy, economic or environmental standard index data.

In combination with the first aspect, the embodiment of the present invention provides a sixth possible implementation manner of the first aspect, wherein, according to the statistical tables of historical data of the target region, the energy development speed V _a and the economic development speed V a of the target region are analyzed and obtained. _b and the environmental development velocity V _c also include:

According to the standardized matrix, calculate its correlation matrix;

Calculate the characteristic root of the correlation matrix;

Among the characteristic roots of the correlation matrix, the characteristic roots whose values are not negative are selected as the values of the principal components.

With reference to the first aspect, the embodiment of the present invention provides a seventh possible implementation manner of the first aspect, which further includes:

Establish a three-dimensional coordinate system with an ellipsoid, and bring the values of energy development speed, economic development speed and environmental development speed into the three-dimensional coordinate system. The X, Y and Z axes of the three-dimensional coordinate system represent energy Values for speed of development, speed of economic development and speed of environmental development;

Adjust the calculation methods of α, β, and γ according to the parameters of the three-dimensional coordinate system brought in, where α represents the degree of coordination between the environment and energy-economy under the condition that energy and economy are coordinated, and β represents the degree of coordination between energy and the environment system In the case of coordination, the coordination degree of economy and energy-economy, γ indicates the coordination degree of energy and economy-environment in the case of coordination of economy and environment system.

In the second aspect, the embodiment of the present invention also provides an energy, economy and environment coordination degree calculation device based on an ellipsoid model, including:

The acquisition module is used to analyze and obtain the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region according to the historical data statistics table of the target region;

The calculation module calculates the energy, economy and environmental coordination degree of the target region according to the following formula,

$\mathrm{<span\; class="notranslate">\; E.</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&alpha;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&beta;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&gamma;</span>}}{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span>}<span\; class="notranslate">\; ;</span>$

in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region.

With reference to the second aspect, the embodiment of the present invention provides a first possible implementation manner of the second aspect, wherein the acquisition module includes:

The calculation unit is used to calculate the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region according to the following formula:

${\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; A</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; C</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ;</span>$

Among them, t is the statistical year of the historical data statistics table, and A, B, and C are the curve fitting equations of the comprehensive development level values of energy, economy, and environment in the target area, respectively.

The method for calculating the energy, economy and environment coordination degree based on the ellipsoid model provided by the embodiment of the present invention adopts the analysis method of the ellipsoid model, and uses the cube model in the prior art to analyze and calculate the coordination degree of energy, economy and environment , making the calculated degree of coordination inaccurate, compared with the established ellipsoid analysis model in advance, and using the ellipsoid analysis model to adjust the coordination degree of environment and energy-economy in the case of energy and economy coordination, in The degree of coordination of economy and energy-economy in the case of coordinated energy and environmental systems, and the degree of coordination of energy and economy-environment in the case of coordinated economic and environmental systems. And finally make the calculated energy, economy and environment coordination degree more accurate.

In order to make the above-mentioned objects, features and advantages of the present invention more comprehensible, preferred embodiments will be described in detail below together with the accompanying drawings.

Description of drawings

In order to illustrate the technical solutions of the embodiments of the present invention more clearly, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of the present invention, and thus It should be regarded as a limitation on the scope, and those skilled in the art can also obtain other related drawings based on these drawings without creative work.

Figure 1 shows a modular schematic diagram of the coupling relationship between energy, economy and environment in related technologies;

Fig. 2 shows a schematic diagram of the trajectory of the energy-economic evolution relationship provided by the embodiment of the present invention;

Fig. 3 shows a schematic diagram of the energy-economy-environment evolution track provided by the embodiment of the present invention.

detailed description

The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. The components of the embodiments of the invention generally described and illustrated in the figures herein may be arranged and designed in a variety of different configurations. Accordingly, the following detailed description of the embodiments of the invention provided in the accompanying drawings is not intended to limit the scope of the claimed invention, but merely represents selected embodiments of the invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without making creative efforts belong to the protection scope of the present invention.

In related technologies, when analyzing and calculating the coordination degree of energy, economy and environment, a cube model is adopted. However, the applicant found that the degree of coordination calculated by using this model for analysis and calculation is not accurate enough. Further, the applicant provides a method for analyzing and calculating the degree of coordination using the ellipsoid model.

In the following, the formula used to calculate the value of the comprehensive development level will be explained first.

Step 1. In general, the basis for calculating the degree of coordination is the statistical table of historical data in the target region. The matrix containing the data in the original index can be directly extracted from the statistical table, and then the matrix should be standardized. Among them, the matrix containing the data in the original indicator is shown in the following matrix: Each element in the matrix is data from the original indicator.

Specifically, standardization can be performed according to the following formula:

Among them, n is the number of samples (plans), p is the evaluation index X ₁ , X ₂ ,...,X _p , the matrix X represents the data in the original index, and the matrix Z represents the standardized data. The resulting standardized matrix is as follows:

Step 2 is to calculate the correlation matrix R based on the standardized matrix. It is calculated as follows:

R=(r _jk ) _p×p , j,k=1,2,...,p Formula 3; where, r _jk is the correlation coefficient between the jth index and the kth index.

further, And there are r _jj = 1, r _jk = r _kj formula 4;

In step 3, the characteristic root of the correlation matrix R needs to be solved. Among them, the characteristic equation of R is

|R-λI|＝0 Formula 5;

λ _i (i=1,2,...p) represents the non-negative characteristic root of R. And confirm the principal components by the following formula:

λ _i (i=1,2,…p)>1 Formula 6;

The number of characteristic roots greater than 1 is taken as K, and the first K characteristic roots are taken as the eigenvalues of the main components.

Step 4, calculate the factor loading matrix p _ij , and find the eigenvector matrix:

_Wherein , the corresponding feature vector T _j =(T _1j , T _2j , . . . , T _ij , .

Afterwards, the composite weight of the standard indicator data is calculated:

Where T represents the comprehensive weight, K is the number of main components,

Finally, after calculating the standard indicator data in step 1 and calculating the comprehensive weight in step 8, the value of the comprehensive development level can be calculated as follows:

F=Z×T Formula 9;

Among them, F is an n×1 matrix, each row represents the value of the comprehensive development level of energy, economy and environment in the corresponding year, n is the number of years obtained through the historical data of the target area; The energy, economic or environmental standard index data calculated from the statistical table of regional data over the years, T is a p×1 matrix, each row represents the weight value under the corresponding index of energy, economy or environment, and p is obtained through the historical data of the target region The number of indicators.

Corresponding to the above steps, an example of combining actual data is also provided below.

Step 1, use formula (1) and formula (2) to standardize the original data. (The original data X is shown in Table 1). In the example, 12 indicators in the economy were selected, namely n=12, p=12, and the standardized matrix data (matrix Z) was calculated.

Table 1 Raw data (matrix X)

Table 2 Standardized Data (Matrix Z)

Step 2, calculate the sample correlation matrix R by formula (3) and formula (4). This step can be performed with SPSS software.

Table 3 Correlation coefficient matrix R

Step 3. Use SPSS software to find the characteristic root of correlation matrix R according to formula (5). It can be seen from the formula (6) that two principal components are taken, that is, K=2.

Table 4 Characteristic root λ _i of correlation matrix R

Step 4. Use SPSS software to obtain the factor loading matrix p _ij as shown in Table 5, use the formula (7) to find the eigenvector matrix, which is the weight of the first principal component and the second principal component, and use the formula (8) to obtain the comprehensive weight , as shown in Table 6.

Table 5 initial factor loading matrix p _ij

Table 6 Economic system comprehensive principal component weight

Step 5. According to the weights in Table 4, use formula (9) to find the scores of the first principal component and the second principal component, and use formula (10) to calculate the comprehensive development level value, as shown in Table 7.

Table 7 Comprehensive principal component values of the economic system

Similarly, the weights and development levels of energy and environment can be obtained as shown in the table below.

Table 8 Energy system comprehensive principal component weights

Table 9 Environmental System Weights

Table 10 Comprehensive principal component values of the energy system

years First Principal Component F1 The second principal component F2 overall ratings 2003 -4.304 -0.943 -3.830 2004 -3.767 -0.067 -3.245 2005 -2.691 0.941 -2.178 2006 -1.816 1.527 -1.344 2007 -0.408 0.791 -0.239 2008 0.867 -2.003 0.461 2009 1.168 -1.403 0.806 2010 1.451 -0.105 1.231 2011 2.653 -0.607 2.193 2012 3.151 0.853 2.827 2013 3.696 1.018 3.318 2014 -4.304 -0.943 -3.830

Table 11 Comprehensive principal component value of environmental system

After calculating the value of the comprehensive development level, the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area can be calculated by means of differentiation:

Among them, t is the statistical year of the statistical table of historical data, and A, B and C are the curve fitting equations of the comprehensive development level values of energy, economy and environment in the target area respectively.

Finally, put the calculated energy development speed V _a , economic development speed V _b and environmental development speed V _c of the target region into the formula obtained by using the ellipsoid model to obtain the target region’s energy, Economic and environmental coordination.

The specific calculation formula is as follows:

in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region.

Next, the process of obtaining formula 11 is described:

According to the coupled system theory, the coupling relationship of the 3E system is an objective representation of the mutual influence and interdependence among the energy, economic and environmental systems and among the internal elements. situation. For example, in the 3E system structure, the carrying capacity and supply capacity of the ecological environment system, as well as the excavation capacity and driving capacity of the energy system are very limited to some extent, which determines that the excessive development of the economic system will inevitably lead to energy Transition loss and damage to the ecological environment system. Conversely, the energy system and the environmental system will also restrict and slow down the development of the economic system through resource shortages, environmental pollution, and government intervention. The coupling relationship between energy, economy and environment is shown in Figure 1.

On the basis of analyzing the relationship between the constituent elements of the 3E system, referring to the dynamic development idea in general systematics, the nonlinear system evolution equation of energy, economy and environment is defined (formula 12, namely Among them, in Indicates the change speed of the comprehensive development index of energy, economy and environment, x ₁ , x ₂ ,...,x _n (i=1,2,...n) represent the index of energy, economy and environment subsystem), (Formula 13, namely f(x ₁ ,x ₂ ,…,x _n )＝f(0)+a ₁ x ₁ +a ₂ x ₂ +…+a _n x _n +θ(x ₁ ,x ₂ ,…,x _n )) is the derivation of (Equation 12), combining Equations 12 and 13 to obtain an approximate linear system (Equation 14, namely ).

According to the idea of describing nonlinear systems by differential equations, construct general functions of energy, economy, and environment nonlinear systems (formulas 15, 16, 17, namely Among them, f(Ener) represents the energy system, f(Eco) represents the economic system, and f(Envi) represents the environmental system. Among them, x _i , y _j , z _k are indicators in energy, economy, and environment systems respectively, and a _i , b _j , c _k are the influence weights of each indicator respectively).

As a unified whole, the 3E system has only three elements f(Ener), f(Eco), and f(Envi), so the evolution equation of the whole 3E system can be established in the form of (formula 18, 19, 20, which is

Among them, A, B, and C represent the comprehensive development level curves of energy, economy, and environmental systems affected by their own system and external systems, respectively.

The energy, economy, and environmental evolution development systems are closely related to the passage of time. Therefore, taking time t as a variable to define the evolution speed of each subsystem, the formula (10) is obtained.

V _a is the evolution and development speed of the energy system, V _b is the evolution and development speed of the economic system, and V _c is the evolution and development speed of the environmental system. Then V _3e can be regarded as a function of V _a , V _b , and V _c , that is, V _3e =g(V _a , V _b , V _c ). It is assumed that changes in the economic system are cyclical, so the energy system and environmental system related to the economic system will also be cyclical. Due to the existence of the development relationship of V _3e =g(V _a , V _b , V _c ), the relationship between the three changes can be analyzed in a three-dimensional Cartesian coordinate system from the perspective of system development speed. Because the development of the economic system has a certain restrictive effect on the energy system and the environmental system, and the development of the economic system has a greater range of change than the energy system and the environmental system, and the energy system has a greater range of change than the environmental system, so the relationship between the three subsystems Two pairs of elliptical change tracks form an ellipsoidal action relationship track in three-dimensional space.

The projection of the ellipsoid on each plane forms the trajectory of the evolution and development relationship of the binary system, that is, the ellipse. For example, the projection of the ellipsoid on the V _a -V _b plane is an ellipse, that is, the trajectory diagram of the relationship between energy and economic evolution, as shown in Figure 2.

Figure 2 reflects the trajectory diagram of the evolution and development relationship between energy and economy. Assuming that point A is the point on the plane where the development speed of energy economy is on the plane, the evolution and development status of energy system and economic system and the coupling degree of coordinated development can be expressed by angle a, from It can be seen in the figure (i.e. Equation 11), then (ie Equation 11). Angle a can also be regarded as a function of V _a and V _b , denoted as m(V _a , V _b ). Then the angle between V _a and V _c is denoted by b, the angle between V _b and V _c is denoted by c, and the coupling degree functions of energy and environment, economy and environment can be recorded as m(V _a , V _c ), m (V _b , V _c ).

The three-dimensional coordinate system is shown in the figure below.

Figure 3 reflects the evolution trajectory of energy-economy-environment. Assuming that point A is a point on the three-dimensional coordinates of energy, economy, and environment development speed, it can be seen from the figure that the angle between OA and plane V _a -V _b is α, and α represents the coordination between energy and economy. Under the circumstances, the degree of coordination between environment and energy-economy. Therefore, the angle between OA and plane V _a -V _c and the angle between plane V _b -V _c are recorded as β and γ, and β represents the coordination degree of economy and energy-economy under the condition of energy and environmental system coordination , γ indicates the coordination degree of energy and economy-environment under the condition that economy and environment system are coordinated. It can be seen from the figure the same way. (Formula 11).

Therefore, it is concluded that in the case of energy and economic system coordination, the degree of coordination between the environment and energy-economy: In the case of energy and environmental system coordination, the degree of coordination of economy and energy-economy: In the case of coordinated economic and environmental systems, the coordination degree of energy and economy-environment: (Formula 11)

Since the overall coordination degree of 3E is inseparable from the three sub-elements, not only the coordination relationship between energy, economy and environment, but also the interaction among them should be considered. Therefore, this paper takes the energy-economy, energy-environment, and economy-environment coordination degree as the weight, and finds the weighted average of α, β, and γ as the 3E system coordination degree. The constructed coordination formula is

Next, a specific example is used to illustrate the rationality of calculation using the method provided in this application.

Step 1, basic parameter acquisition,

First of all, according to the scientific, systematic, comprehensive, hierarchical and regional principles of the establishment of the index system, from the two directions of depth and breadth, drawing on the existing research results of the coordinated development index system, establish the Beijing energy, economy, environment Coordinate the development of the indicator system. According to the availability of indicator data, the data of Beijing from 2003 to 2014 are selected for calculation. The data are all from the "National Statistical Yearbook", "Beijing Statistical Yearbook" and "National Bureau of Statistics". Tables 12, 13, and 14 are the indicators and data of energy, economy, and environment systems respectively.

Table 12 Indicators and data of energy system in Beijing from 2003 to 2014

Table 13 Indicators and data of Beijing's economic system from 2003 to 2014

Table 14 Indicators and data of Beijing's environmental system from 2003 to 2014

Step 2, calculation, uses the principal component analysis method to determine the weight of the indicators, and calculates the comprehensive development level of energy, economy, and environment according to formulas 4, 5, and 6. The results are shown in Table 15.

Table 15 Comprehensive development level of energy, economy and environment

The three curves of comprehensive development levels of energy, economy and environment from 2003 to 2014 were respectively polynomially fitted to obtain the evolution fitting expressions of energy, economy and environment:

f(Ener)=0.9056t+0.0019t ² -0.0016t ³ -4.8909, fitting degree R ² =0.99;

f(Eco)=-0.2208t+0.1915t ² -0.0103t ³ -3.4066, the fitting degree R ² =0.99;

f(Envi)=-0.1235t-0.1477t ² +0.0094t ³ -3.7978, fitting degree R ² =0.98;

Using Equation 10, calculate:

V _α =0.9056+0.0038t-0.0048t ² ;

V _β =-0.2208+0.3830t-0.0309t ² ;

V _γ =-0.1235-0.2954t+0.0282t ² ;

The range of t is 1 to 12, corresponding to 2003 to 2014.

Then, formula 11 can be used to obtain the value of the energy-economy-environment ternary coupling degree. The following is a list of energy-economy-environmental coupling, and only the results from 2007 to 2014 are listed in tabular form.

Table 16 List of Energy-Economy-Environmental Coupling Degrees

Corresponding to the method provided above, the embodiment of the present application also provides an ellipsoid model-based energy, economy and environmental coordination degree calculation device, including:

The acquisition module is used to analyze and obtain the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region according to the historical data statistics table of the target region;

The calculation module calculates the energy, economy and environmental coordination degree of the target region according to the following formula,

$\mathrm{<span\; class="notranslate">\; E.</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&alpha;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&beta;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&gamma;</span>}}{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span>}<span\; class="notranslate">\; ;</span>$

in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region.

Those skilled in the art can clearly understand that for the convenience and brevity of the description, the specific working process of the above-described system, device and unit can refer to the corresponding process in the foregoing method embodiment, which will not be repeated here.

The above is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Anyone skilled in the art can easily think of changes or substitutions within the technical scope disclosed in the present invention. Should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims (10)

1. The calculation method of energy, economy and environment coordination degree based on ellipsoidal model, is characterized in that, comprises: According to the statistical table of historical data of the target area, analyze and obtain the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area; Calculate the energy, economy and environment coordination degree of the target region according to the following formula, $\mathrm{<span\; class="notranslate">\; E.</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&alpha;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&beta;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&gamma;</span>}}{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span>}<span\; class="notranslate">\; ;</span>$ in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region. 2. The method according to claim 1, characterized in that, according to the historical data statistics table of the target area, the energy development speed V _a , the economic development speed V _b and the environment development speed V _c of the target area are analyzed and include : Calculate the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area according to the following formula: ${\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; A</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; C</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ;</span>$ Among them, t is the statistical year of the statistical table of historical data, and A, B and C are the curve fitting equations of the comprehensive development level values of energy, economy and environment in the target region respectively. 3. The method according to claim 2, characterized in that, according to the statistical table of historical data of the target area, the energy development speed V _a , the economic development speed V _b and the environment development speed V _c of the target area are analyzed and obtained. Including: using the principal component analysis method to determine the economic standard index data in the economic system of the target region and the weight value of the economic standard index data, and determine the energy standard index data in the energy system of the target region and the weight of the above energy standard index data value, and determine the weight value of the environmental standard index data and the environmental standard index data in the environmental system of the target region; Using the method of weighted calculation, the value of the comprehensive economic development level of the target region is calculated according to the value of the determined economic standard index data and the weight value of the economic standard index data, and the value of the comprehensive economic development level of the target area is calculated by a fitting method The curve fitting equation of ; Using the method of weighted calculation, according to the value of the determined energy standard index data and the weight value of the energy standard index data, the value of the comprehensive energy development level of the target area is calculated, and the value of the comprehensive development level of energy is calculated by a fitting method The curve fitting equation of ; Using the method of weighted calculation, according to the value of the determined environmental standard index data and the weight value of the environmental standard index data, the value of the comprehensive environmental development level of the target area is calculated, and the value of the comprehensive development level of the environment is calculated by a fitting method The curve fitting equation for . 4. The method according to claim 3, characterized in that, according to the historical data statistics table of the target area, the energy development speed V _a , the economic development speed V _b and the environment development speed V _c of the target area are analyzed and obtained. include: According to the following formulas, the values of comprehensive development levels of energy, economy and environment are calculated respectively, F=Z×T; Among them, F is an n×1 matrix, each row represents the value of the comprehensive development level of energy, economy and environment in the corresponding year, n is the number of years obtained through the historical data of the target area; The energy, economic or environmental standard index data calculated from the statistical table of regional data over the years, T is a p×1 matrix, each row represents the weight value under the corresponding index of energy, economy or environment, and p is obtained through the historical data of the target region The number of indicators. 5. The method according to claim 4, characterized in that, according to the statistical table of historical data of the target area, the energy development speed V _a , the economic development speed V _b and the environment development speed V _c of the target area are analyzed or obtained. include: Calculate the weight value of energy, economic or environmental standard index data according to the following formula, $\mathrm{<span\; class="notranslate">\; T</span>}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; K</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}<span\; class="notranslate">\; \&times;</span>\frac{{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}}{\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; K</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}<span\; class="notranslate">\; ,</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}\mathrm{<span\; class="notranslate">\; K</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; ;</span>$ ${\mathrm{<span\; class="notranslate">\; T</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; j</span>}}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; p</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; j</span>}}}{\sqrt{{\mathrm{<span\; class="notranslate">\; \&lambda;</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}}}<span\; class="notranslate">\; ,</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; p</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; ...</span>}\mathrm{<span\; class="notranslate">\; K</span>}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; ;</span>$ Among them, T _j ＝(T _1j ,T _2j ,...T _ij ...T _pj )(j=1,2...K) means that among the components extracted from the statistical tables of historical data in the target region, the first K principal components The weight of j principal components, K the number of principal components, T is the weight value of energy, economic or environmental standard index data, λ _k is the characteristic root of principal components, p=(p _1j ,p _2j ,…p _ij , …p _pj )(j=1,2,…K) is the factor loading matrix derived from the statistical tables of historical data in the target area. 6. The method according to claim 5, characterized in that, the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target area are analyzed according to the historical data statistical tables of the target area. include: Calculate the economic standard indicator data through the following steps; Extract the matrix containing the data in the original indicators from the historical data statistics table of the target region; Standardize the data in the original indicators according to the following formula, Among them, X=(X _i1 ,X _i2 ,…X _ip )(i=1,2…,n) is a matrix containing original index data of energy, economy or environment, Z=(Z _i1 ,Z _i2 ,…Z _ip ) (i=1,2,...n) is a standardized matrix of energy, economic or environmental data; Statistically standardize each data in the standardized matrix to obtain energy, economic or environmental standard index data. 7. The method according to claim 6, characterized in that, the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region are analyzed according to the historical data statistical table of the target region. include: According to the standardized matrix, calculate its correlation matrix; Calculate the characteristic root of the correlation matrix; Among the characteristic roots of the correlation matrix, the characteristic roots whose values are not negative are selected as the values of the principal components. 8. The method according to claim 7, further comprising: Establish a three-dimensional coordinate system with an ellipsoid, and bring the numerical values of energy development speed, economic development speed and environmental development speed into the three-dimensional coordinate system, and the X, Y and Z of the three-dimensional coordinate system The axes represent the values of the energy development speed, the economic development speed and the environmental development speed respectively; Adjust the calculation method of α, β, and γ according to the parameter content of the three-dimensional coordinate system brought in, wherein, α represents the coordination degree of the environment and energy-economy under the condition that energy and economy are coordinated, and β represents the degree of coordination between the energy and the economy. In the case of energy and environmental system coordination, the coordination degree of economy and energy-economy, γ indicates the coordination degree of energy and economy-environment in the case of coordination of economic and environmental systems. 9. The energy, economy and environment coordination degree calculation device based on the ellipsoid model, characterized in that it includes: The acquisition module is used to analyze and obtain the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region according to the historical data statistics table of the target region; The calculation module calculates the energy, economy and environmental coordination degree of the target region according to the following formula, $\mathrm{<span\; class="notranslate">\; E.</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&alpha;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&beta;</span>}<span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; *</span>\mathrm{<span\; class="notranslate">\; \&gamma;</span>}}{\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; +</span>\mathrm{<span\; class="notranslate">\; m</span>}<span\; class="notranslate">\; (</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; )</span>}<span\; class="notranslate">\; ;</span>$ in, V _a is the energy development speed of the target region, V _b is the economic development speed of the target region, and V _c is the environmental development speed of the target region. 10. The device according to claim 9, wherein the acquiring module comprises: The calculation unit is used to calculate the energy development speed V _a , the economic development speed V _b and the environmental development speed V _c of the target region according to the following formula: ${\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; a</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; A</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; b</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; B</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; V</span>}}_{\mathrm{<span\; class="notranslate">\; c</span>}}<span\; class="notranslate">\; =</span>\frac{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; C</span>}}{\mathrm{<span\; class="notranslate">\; d</span>}\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; ;</span>$ Among them, t is the statistical year of the statistical table of historical data, and A, B and C are the curve fitting equations of the comprehensive development level values of energy, economy and environment in the target region respectively.