CN102708240A - Ecological elasticity index model method of town and country ecological planning - Google Patents

Abstract

The invention discloses a method for quantitative partitioning based on an ecological elasticity index model. The method comprises the following steps of: investigating and surveying in a target region, and drawing the result into thematic maps; using a predetermined grading standard to respectively divide single factor grades; using an analytic hierarchy method to respectively calculate the weight value; mutually superposing and combining the thematic maps by a GIS (Geographic Information System); and performing weighted averaging on the weather single factor, the soil single factor, the hydrology single factor, the ground surface physical single factor and the vegetation condition single factor to obtain the analysis result of the ecological elasticity index model, thereby finishing the quantitative partitioning of the target region. According to the embodiment of the invention, by the GIS (Geographic Information System), the information in various aspects of the ecological elasticity recovery capability is comprehensively superposed, the analytic hierarchy method is optimized, a quantitative evaluation system is established, and comprehensive and simple research on the internal conditions of the ecological elasticity recovery capability is expected and practiced in space distribution of city planning.

Description

The ecological elasticity indexes model method of a kind of town and country ecological planning

Technical field

The present invention relates to a kind ofly utilize ecological elasticity indexes model to carry out the quantification partitioned method.

Background technology

Ecological elasticity indexes model method is special a kind of space ecological analysis method to the various ecological resources system resilience abilities in the ecosystem in the ecological planning of town and country.Ecological elasticity indexes is the ability that refers to can the oneself the keeping of the ecosystem, self-control and various pressure of opposing and disturbance, also is ecological elastic force.

Basically determined the ecosystem features situation at ecosystem mesoclimate condition, vegetation state, soil, geology and geomorphology and moisture condition.The various combination of these factors has constituted the different ecological environmental form; The recovery capability or the situation of the different ecosystem-types and the corresponding ecosystem have also just been determined; The antijamming capability that comprises system with disturbed after recovery capability etc.; Weather, the hydrology and soil are the principal elements of recovery capability after decision ecosystem-type and system are interfered; Vegetation state is a principal element of keeping ecosystem stability and antijamming capability, and the combination of these factors can reflect ecosystem features basically.

The defective of the prior art scheme of present ecological elasticity indexes: the situation of wanting the recovery of concentrated expression local area ecological to city planning and construction; Not only to estimate, simultaneously to city planning practitioner's importantly operability with how to implement to the planning management and the enforcement aspect of city planning its ecological elastic recovery capability.

Summary of the invention

The present invention one of is intended to solve the problems of the technologies described above at least to a certain extent or provides a kind of useful commerce to select at least.For this reason, one object of the present invention is to propose a kind of ecological elasticity indexes model of target area that utilizes and carries out the quantification partitioned method.

Carry out the quantification partitioned method according to the ecological elasticity indexes model of the utilization of the embodiment of the invention; May further comprise the steps: step S100; In the target area; Investigate, survey to meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor; And, set up unified coordinate system with GIS vector quantization as a result, be depicted as meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps; Step S200 divides said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial grade of vegetation state respectively with predetermined grade scale; Step S300, the application level analytic approach is calculated meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial weighted value of vegetation state respectively; Step S400 adopts GIS with the merging that superposes mutually of said meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps; Step S500 obtains the analysis result of ecological elasticity indexes model to said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor weighted mean, thereby accomplish the quantification of target area is distinguished.

Carry out the quantification partitioned method according to the ecological elasticity indexes model of the utilization of the embodiment of the invention; Through the informix stack of GIS (geographical information technology) with ecological elastic recovery capability each side; The Optimization Layer fractional analysis; Set up the quantification appraisement system, the inherent situation of ecological elastic recovery capability is studied in expectation comprehensively easily, and implements in the space issue of city planning.

In addition, the ecological elasticity indexes model of usefulness according to the above embodiment of the present invention carries out the quantification partitioned method and can also have following additional technical characterictic:

According to one embodiment of present invention; Step S100 comprises, step S101 is with meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor vector quantization; Polar plot is continuously nonoverlapping unique territory; Unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies, need meet the following conditions

S1∪S2∪...∪Sn＝S

Si∩Sj＝0

S is certain monofactorial scope, and S1, S2, Si, Sj, Sn are meant a plurality of small area that certain single-factor the inside is divided, and n is the zonule sum; Step S102; Set up unified coordinate system system, the meteorological single-factor figure behind the vector quantization, soil single-factor figure, hydrology single-factor figure, face of land physics single-factor figure, vegetation state single-factor figure are unified in define projection among the GIS, conversion is projected as the Gauss projection coordinate system.

According to one embodiment of present invention, said step S300 comprises, step S301, relatively n single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz, thus confirm proportion shared in z, get two factor y at every turn _iAnd y _jUse a _{The ij table}Show y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _Ij, use numeral 1～9 and inverse thereof as scale, obtain weighted value by the monofactorial difference of ecological elasticity shown in the table 1,

Table 1

Judgment value	Comparison	Intensity
			1	Y i＝Y j	Equate
3	Y i＞∞Y j	Y iBe better than Y slightly j
			5	Y i＞＞Y j	Y iObviously be better than Y j
7	Y i＞＞＞Y j	Y iCompare Y jMuch better
			9	Y i＞Y j	Y iExtremely be better than Y j
1/3	Y i＜Y j	Y iSlightly inferior to Y j
			1/5	Y i＜＜Y j	Y iObviously inferior to Y j
1/7	Y i＜＜＜Y j	Y iCompare Y jMuch inferior
			1/9	Y i＜∞Y j	Y iDefinitely inferior to Y j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank

Step S302, n is constituted one by element relatively and compares the judgment matrix A=(a of (in pairs relatively) in twos _Ij) _{N * n}. obviously, judgment matrix has character,

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& ...& a_{1n}\\ a_{21}& a_{22}& ...& a_{2n}\\ .& .& & .\\ .& .& & .\\ .& .& & .\\ a_{n1}& a_{n2}& ...& a_{\mathrm{nn}}\end{array}\right)$

a _ij＞0，

a _ii＝1(i，j＝1，2，...n)

Step S303, weight coefficient calculates.The product Z of each row of judgment matrix _iThen press Formula calculates

Z _iN th root W _iTo vectorial W _i=[W ₁, W ₂... W _n] ^TCarry out normalization and calculate, obtain proper vector, calculate maximum characteristic root λ max: $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Step S304 carries out consistency check to each the ecological elasticity single-factor weighted value that obtains, and the check formula is CI=(λ max-n)/(n 1);

Step S305 obtains the monofactorial weight of each ecological elasticity w as a result _iMean value,

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a certain single-factor weight, and m is for calculating the sum of weight, and di is the monofactorial some weighted values of this ecology load pressure.

Additional aspect of the present invention and advantage part in the following description provide, and part will become obviously from the following description, or recognize through practice of the present invention.

Description of drawings

Above-mentioned and/or additional aspect of the present invention and advantage are from obviously with easily understanding becoming the description of embodiment in conjunction with figs, wherein:

Fig. 1 is the process flow diagram that carries out the quantification partitioned method according to the ecological elasticity indexes model of the utilization of the embodiment of the invention; With

Fig. 2 carries out the synoptic diagram that superposes in the quantification partitioned method according to the ecological elasticity indexes model of the utilization of the embodiment of the invention.

Embodiment

Describe embodiments of the invention below in detail, the example of said embodiment is shown in the drawings, and wherein identical from start to finish or similar label is represented identical or similar elements or the element with identical or similar functions.Be exemplary through the embodiment that is described with reference to the drawings below, be intended to be used to explain the present invention, and can not be interpreted as limitation of the present invention.

In description of the invention; It will be appreciated that; The orientation of indications such as term " " center ", " vertically ", " laterally ", " length ", " width ", " thickness ", " on ", D score, " preceding ", " back ", " left side ", " right side ", " vertically ", " level ", " top ", " end " " interior ", " clockwise ", " counterclockwise " or position relation are for based on orientation shown in the drawings or position relation; only be to describe with simplifying for the ease of describing the present invention; rather than the device or the element of indication or hint indication must have specific orientation, with specific azimuth configuration and operation, therefore can not be interpreted as limitation of the present invention.

In addition, term " first ", " second " only are used to describe purpose, and can not be interpreted as indication or hint relative importance or the implicit quantity that indicates indicated technical characterictic.Thus, one or more a plurality of these characteristics can be shown or impliedly comprised to the characteristic that is limited with " first ", " second " clearly.In description of the invention, except as otherwise noted, the implication of " a plurality of " is two or more, only if clear and definite qualification is arranged in addition.

In the present invention, only if clear and definite regulation and qualification are arranged in addition, broad understanding should be done in terms such as term " installation ", " linking to each other ", " connection ", " fixing ", for example, can be to be fixedly connected, and also can be to removably connect, or connect integratedly; Can be mechanical connection, also can be to be electrically connected; Can be directly to link to each other, also can link to each other indirectly through intermediary, can be the connection of two element internals.For those of ordinary skill in the art, can understand above-mentioned term concrete implication in the present invention as the case may be.

In the present invention; Only if clear and definite regulation and qualification are arranged in addition; First characteristic second characteristic it " on " or D score can comprise that first and second characteristics directly contact, can comprise that also first and second characteristics are not directly contacts but through the contact of the additional features between them.And, first characteristic second characteristic " on ", " top " and " above " comprise first characteristic directly over second characteristic and oblique upper, or only represent that the first characteristic level height is higher than second characteristic.First characteristic second characteristic " under ", " below " and " below " comprise first characteristic directly over second characteristic and oblique upper, or only represent that the first characteristic level height is less than second characteristic.

Describe in detail according to the ecological elasticity indexes model of the utilization of the embodiment of the invention with reference to accompanying drawing below and carry out the quantification partitioned method.

Carry out the quantification partitioned method according to the ecological elasticity indexes model of the utilization of the embodiment of the invention, may further comprise the steps:

Step S100; In the target area; Investigate, survey to meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor; And, set up unified coordinate system with GIS vector quantization as a result, be depicted as meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps.

Particularly, collect all kinds of of target area with ecological for information about flexible.The information content should contain: weather conditions (comprising accumulated temperature, rainfall amount, temperature, meteorological disaster etc.), soil (comprising soil grade, land type, the soil erosion etc.), hydrologic condition (comprising underground water, surface water etc.), face of land physical state (comprising elevation, the gradient etc.), vegetation state (vegetation pattern, afforestation rate etc.).Investigate measurement in the target area, be depicted as thematic maps.

Select the main ecological resources single-factor in weather conditions, soil, the hydrology, surface condition, the vegetation state, the ecological resources single-factor is used the GIS vector quantization, set up unified coordinate system, draw ecological resources elasticity single-factor figure.

With each ecological resources elasticity single-factor vector quantization, it is continuously nonoverlapping unique territory that polar plot requires.Unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies.This be since in city planning every land used only corresponding to a kind of soil utilize character be determine.Need meet the following conditions

S1∪S2∪...∪Sn＝S

Si∩Sj＝0

S is certain monofactorial scope, and S1, S2, Si, Sj, Sn etc. are meant the arbitrarily small regional scope that certain single-factor the inside is divided.N is the zonule sum.

Set up unified coordinate system system, the ecological resources elasticity single-factor figure behind the vector quantization is unified in defines projection among the GIS, conversion is projected as China's Gauss projection coordinate system commonly used.

Step S200 divides said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial grade of vegetation state respectively with predetermined grade scale.Each ecological resources elasticity single-factor of choosing is utilized present situation and load-bearing capacity, divided rank according to what existing code requirement was estimated resource.Grade is low more, and the ecological resources elasticity indexes is low more, and the self-recovery ability is good more; Otherwise grade is high more, and the ecological resources elasticity indexes is strong more, and the self-recovery ability is poor more, needs control ecological disruption behavior.

Step S300, the application level analytic approach is calculated meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial weighted value of vegetation state respectively.

In other words, compare n ecological resources elasticity single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz (ecological resources elasticity indexes), thus confirm their shared proportion in z, get two factor y at every turn _iAnd y _jUse a _IjExpression y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _IjWith numeral 1～9 and reciprocal, obtain weighted value by the difference of monofactorial importance value of the different ecological of quantification shown in the table 1 resource elasticity and influence property as scale.

Table 1

Judgment value	Comparison	Intensity
			1	Y i＝Y j	Equate
3	Y i＞∞Y j	Y iBe better than Y slightly j
			5	Y i＞＞Y j	Y iObviously be better than Y j
7	Y i＞＞＞Y j	Y iCompare Y jMuch better
			9	Y i＞Y j	Y iExtremely be better than Y j
1/3	Y i＜Y j	Y iSlightly inferior to Y j
			1/5	Y i＜＜Y j	Y iObviously inferior to Y j
1/7	Y i＜＜＜Y j	Y iCompare Y jMuch inferior
			1/9	Y i＜∞Y j	Y iDefinitely inferior to Y j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank

N is constituted one by element relatively and compares the judgment matrix A=(a of (in pairs relatively) in twos _Ij) _{N * n}. obviously, judgment matrix has character:

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& ...& a_{1n}\\ a_{21}& a_{22}& ...& a_{2n}\\ .& .& & .\\ .& .& & .\\ .& .& & .\\ a_{n1}& a_{n2}& ...& a_{\mathrm{nn}}\end{array}\right)$

a _Ij＞0, a _Ii=1 (i, j=1,2 ... n) weight coefficient calculates.The product Zi of each row of judgment matrix; Then calculate the n th root Wi of Zi by

formula.Vectorial Wi=[W1, W2...Wn] T is carried out normalization calculate, obtain proper vector.

Calculate maximum characteristic root λ max: $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Ecological resources elasticity single-factor weighted value to obtaining carries out consistency check.

In order to check the consistency problem of judgment matrix; The consistance that needs to calculate coincident indicator CI judgment matrix also has randomness; This consistance at random can represent that the value of RI is relevant with the dimension size of matrix with mean random coincident indicator RI, relative uniformity index CR=CI/RI.It is generally acknowledged that when CR≤0.1, judgment matrix meets the crash consistency condition basically; When CR＞0.1, given judgment matrix does not meet the crash consistency condition, need adjust and revise.

CI＝(λmax-n)/(n-1)

Popular, if 1 pair of single-factor of ecological resources elasticity single-factor, 2 influence value a＞1 (being that single-factor 1 is bigger than single-factor 2), 2 pairs of single-factors of single-factor, 3 influence value b＞1, then in logic, 1 pair of single-factor 3 influence value of single-factor are c=ab＞1 certainly.But it is impossible when being configured to the comparison matrix, requiring to satisfy c=ab.Therefore move back and require paired comparator matrix that certain consistance is arranged.If c＜1, then the result is inconsistent, must recomputate weighted value.

Obtain the monofactorial weight of each ecological resources elasticity mean value of wi as a result.

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a certain single-factor weight, and m is for calculating the sum of weight, and di is the monofactorial some weighted values of this ecological resources elasticity.

Step S400 adopts GIS with the merging that superposes mutually of said meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps.

Utilize the GIS geographical information technology to merge stack each other each ecological resources elasticity single-factor; Stacked system is as shown in Figure 2; The key element of input is represented two different individual ecological resources elasticity single-factors respectively, and the figure layer of output has kept monofactorial all characteristics of two each and every one ecological resources elasticity.

Step S500 obtains the analysis result of ecological elasticity indexes model to said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor weighted mean, thereby accomplish the quantification of target area is distinguished.

In specific words, each ecological resources elasticity single-factor weighted mean is obtained ecological resources flexibility analysis result.

${\mathrm{CSI}}^{\mathrm{eco}}=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_i^{\mathrm{eco}}\×W_i^{\mathrm{eco}}$

In the formula, CSI ^Eco, ecological elasticity indexes; Si ^Eco: each characteristic element of the ecosystem; Wi ^Eco: the corresponding weighted value of each key element of the ecosystem.

According to the analysis result criteria for classification of forcing normal distribution method with the ecological elasticity factor, accomplish ecological resources elasticity indexes analysis chart, thus auxiliary Fu Zaiqu, equilibrium area, the border, overload district of confirming ecological resources.

Forcing normal distribution is according to the normal distribution principle; The regularity of distribution of " broad in the middle, two is little " that promptly is commonly called as; Confirm opinion rating and each grade shared number percent in sum in advance, confirm Fu Zaiqu, equilibrium area, border, overload district according to ecological resources elasticity indexes analysis result is auxiliary then.

Be that example explains that the utilization ecology elasticity indexes model according to the embodiment of the invention carries out the quantification partitioned method below with the Huairou District, Beijing City.

Utilization with the embodiment of the invention is carried out the analysis of ecological resources elasticity indexes to Huairou, Beijing, is weather (S1) in the selected principal element of definite ecosystem elasticity indexes, vegetation (S2), soil (S3), landform (S4) and the hydrology (S5).Choose 11 indexs altogether as the indicator layer evaluation index, specific as follows:

S1={I1, I2, I3}={ >=10 ℃ accumulated temperature, rainfall amount, frostless season }

S2={I4, I5}={ vegetation pattern, forest cover degree }

S3={I6, I7, I8}={ soil grade, soil types, soil erosion situation }

S4={I9}={ height above sea level }

S5={I10, I11}={ underground water, surface water }

Adopt analytical hierarchy process to carry out weight calculation, the result sees table 2.Calculate through carrying out the spatial analysis stack in the GIS-Geographic Information System, ask for evaluation result.Formulate the suitable classification of evaluation set, construction ecological resources elasticity single-factor evaluation district comes that the oneself of ecological resources in the reflecting regional keeps, capacity of self-regulation.

Table 2 ecological resources elasticity single-factor weight table

Ecological resources elasticity single-factor	Weight
		Accumulated temperature	0.035684
Quantity of precipitation	0.028451
		Frostless season	0.022664
Vegetation pattern	0.08552
		Coverage rate	0.12828
The soil grade	0.189775
		Soil types	0.120723
The soil erosion	0.051202
		Height above sea level	0.0775
Underground water	0.15618
		Surface water	0.10412

According to evaluation result, the ecological resources self-recovery ability in most of mountain area is better than the Plain.Ecological resources elasticity indexes such as Huairou, the middle and south, Huairou town town, cities and towns, mausoleum, Yang Song town, north house town are higher, and the self-recovery ability is poor more, need control ecological disruption behavior.Northern ecological resources elasticity indexes is lower, and the ecological resources elasticity indexes is low more, and the self-recovery ability is good more.

Through using the present invention; Overcome lack in the urban ecological planning to all kinds of ecological resources oneselfs keep, the evaluation of self-control and can't correctly estimate the defective of the mining inetesity of various resources; A kind of resource elasticity indexes computing method to the town and country ecological planning of quantification have been proposed, for the city planning developing direction provides reference.

In the description of this instructions, the description of reference term " embodiment ", " some embodiment ", " example ", " concrete example " or " some examples " etc. means the concrete characteristic, structure, material or the characteristics that combine this embodiment or example to describe and is contained at least one embodiment of the present invention or the example.In this manual, the schematic statement to above-mentioned term not necessarily refers to identical embodiment or example.And concrete characteristic, structure, material or the characteristics of description can combine with suitable manner in any one or more embodiment or example.

Although illustrated and described embodiments of the invention above; It is understandable that; The foregoing description is exemplary; Can not be interpreted as limitation of the present invention, those of ordinary skill in the art can change the foregoing description under the situation that does not break away from principle of the present invention and aim within the scope of the invention, modification, replacement and modification.

Claims (3)

1. one kind is utilized ecological elasticity indexes model to carry out the quantification partitioned method, it is characterized in that, may further comprise the steps:

Step S100; In the target area; Investigate, survey to meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor; And, set up unified coordinate system with GIS vector quantization as a result, be depicted as meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps;

Step S200 divides said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial grade of vegetation state respectively with predetermined grade scale;

Step S300, the application level analytic approach is calculated meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, the monofactorial weighted value of vegetation state respectively;

Step S400 adopts GIS with the merging that superposes mutually of said meteorological single-factor thematic maps, soil single-factor thematic maps, hydrology single-factor thematic maps, face of land physics single-factor thematic maps, vegetation state single-factor thematic maps;

Step S500 obtains the analysis result of ecological elasticity indexes model to said meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor weighted mean, thereby accomplish the quantification of target area is distinguished.

2. according to claim 1ly a kind ofly utilize ecological elasticity indexes model to carry out the quantification partitioned method, it is characterized in that, step S100 comprises,

Step S101; With meteorological single-factor, soil single-factor, hydrology single-factor, face of land physics single-factor, vegetation state single-factor vector quantization, polar plot is continuously nonoverlapping unique territory, and unique territory is meant that each key element of vector quantization interconnects but non-overlapping copies; Need meet the following conditions

S1∪S2∪...∪Sn＝S

Si∩Sj＝0

S is certain monofactorial scope, and S1, S2, Si, Sj, Sn are meant a plurality of small area that certain single-factor the inside is divided, and n is the zonule sum;

Step S102; Set up unified coordinate system system, the meteorological single-factor figure behind the vector quantization, soil single-factor figure, hydrology single-factor figure, face of land physics single-factor figure, vegetation state single-factor figure are unified in define projection among the GIS, conversion is projected as the Gauss projection coordinate system.

3. according to claim 1ly a kind ofly utilize ecological elasticity indexes model to carry out the quantification partitioned method, it is characterized in that, said step S300 comprises,

Step S301, relatively n single-factor yz=(yz ₁, yz ₂... yz _n) to the influence of target zz, thus confirm proportion shared in z, get two factor y at every turn _iAnd y _jUse a _IjExpression y _iWith y _jTo the ratio of the influence degree of z, measure a in 1～9 ratio scale _Ij, use numeral 1～9 and inverse thereof as scale, obtain weighted value by the monofactorial difference of ecological elasticity shown in the table 1,

Table 1

Judgment value Comparison Intensity

1 Y _i＝Y _j Equate 3 Y _i＞∞Y _j Y _iBe better than Y slightly _j 5 Y _i＞＞Y _j Y _iObviously be better than Y _j 7 Y _i＞＞＞Y _j Y _iCompare Y _jMuch better 9 Y _i＞Y _j Y _iExtremely be better than Y _j 1/3 Y _i＜Y _j Y _iSlightly inferior to Y _j 1/5 Y _i＜＜Y _j Y _iObviously inferior to Y _j 1/7 Y _i＜＜＜Y _j Y _iCompare Y _jMuch inferior 1/9 Y _i＜∞Y _j Y _iDefinitely inferior to Y _j

Annotate: 2,4,6,8 and 1/2,1/4,1/6,1/8 expression intensity is between corresponding adjacent rank

Step S302, n is constituted one by element relatively and compares the judgment matrix A=(a of (in pairs relatively) in twos _Ij) _{N * n}. obviously, judgment matrix has character,

$A=\left(\begin{array}{cccc}{a}_{11}& a_{12}& ...& a_{1n}\\ a_{21}& a_{22}& ...& a_{2n}\\ .& .& & .\\ .& .& & .\\ .& .& & .\\ a_{n1}& a_{n2}& ...& a_{\mathrm{nn}}\end{array}\right)$

a _ij＞0，

a _ii＝1(i，j＝1，2，...n)

Step S303, weight coefficient calculates.The product Z of each row of judgment matrix _iThen press

Formula calculates Z _iN th root W _iTo vectorial W _i=[W ₁, W ₂... W _n] ^TCarry out normalization and calculate, obtain proper vector, calculate maximum characteristic root λ max: $\mathrm{\λ}\mathrm{Max}={\mathrm{\Σ}}_{i=1}^n\frac{\left(\mathrm{KW}\right)i}{\mathrm{NWi}}$

Step S304 carries out consistency check to each the ecological elasticity single-factor weighted value that obtains, and the check formula is CI=(λ max-n)/(n-1);

Step S305 obtains the monofactorial weight of each ecological elasticity w as a result _iMean value,

$\mathrm{Wi}=\frac{{\mathrm{\Σ}}_{i=1}^m\mathrm{Di}}{m}$ And w1+w2+...+wn=1

Wherein wi is a certain single-factor weight, and m is for calculating the sum of weight, and di is the monofactorial some weighted values of this ecology load pressure.

Images