環(huán)境稅對(duì)經(jīng)濟(jì)發(fā)展和能源強(qiáng)度影響的動(dòng)力學(xué)分析

陳 乾

(揚(yáng)州工業(yè)職業(yè)技術(shù)學(xué)院 基礎(chǔ)科學(xué)部,江蘇 揚(yáng)州 225000)

環(huán)境稅是為了保護(hù)和改善環(huán)境,對(duì)直接向環(huán)境排放污染物的企業(yè)或單位征收的稅種,上述污染物主要指的是大氣污染物,水污染物和固體廢物.征收環(huán)境稅可以很好地從自由經(jīng)濟(jì)市場(chǎng)角度處理環(huán)境問(wèn)題,是調(diào)整經(jīng)濟(jì)結(jié)構(gòu),構(gòu)建低碳經(jīng)濟(jì)的重要突破口,也是實(shí)現(xiàn)綠色發(fā)展的重要舉措.隨著中國(guó)經(jīng)濟(jì)的快速發(fā)展,大量化石能源被使用,化石能源大量使用過(guò)程中產(chǎn)生的二氧化碳、二氧化硫以及各種固體廢棄物給環(huán)境造成了極大的危害[1].為了保護(hù)環(huán)境,治理污染,中國(guó)于2018年正式征收環(huán)境稅,但由于中國(guó)化石能源在總能源中的占比高達(dá)到86%[2],環(huán)境稅的征收無(wú)疑將對(duì)中國(guó)經(jīng)濟(jì)發(fā)展產(chǎn)生影響.因此,研究環(huán)境稅,化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展之間的關(guān)系具有重要意義.本文基于環(huán)境稅、化石能源消耗和經(jīng)濟(jì)發(fā)展之間復(fù)雜的相互作用關(guān)系,建立一個(gè)三維動(dòng)力系統(tǒng)模型,并利用非線性相關(guān)理論,研究環(huán)境稅、化石能源消耗和經(jīng)濟(jì)發(fā)展之間的關(guān)系.

由于現(xiàn)階段環(huán)境稅的實(shí)施還處于初級(jí)階段,因此,環(huán)境稅對(duì)經(jīng)濟(jì)發(fā)展的影響尚不是很明確.Aloi等[3]基于內(nèi)生增長(zhǎng)模型研究發(fā)現(xiàn)環(huán)境稅能通過(guò)加速企業(yè)研發(fā)促進(jìn)經(jīng)濟(jì)增長(zhǎng).牛歡等[4]通過(guò)構(gòu)建環(huán)境稅的世代交替模型,發(fā)現(xiàn)環(huán)境稅的實(shí)施能對(duì)經(jīng)濟(jì)發(fā)展起到一定紅利作用.陳靜茹[5]依據(jù)相關(guān)性檢驗(yàn)分析發(fā)現(xiàn)環(huán)境稅可以從促進(jìn)綠色生產(chǎn),帶動(dòng)環(huán)保業(yè),充實(shí)地方收入等方面帶動(dòng)經(jīng)濟(jì)增長(zhǎng).但是王佳鄧等[6]通過(guò)構(gòu)建包含江蘇省環(huán)境稅和經(jīng)濟(jì)發(fā)展的可計(jì)算均衡模型發(fā)現(xiàn)環(huán)境稅的征收會(huì)對(duì)江蘇省的經(jīng)濟(jì)發(fā)展起到一定的抑制作用.徐欣欣[7]從稅收政策和稅制改革出發(fā),構(gòu)建環(huán)境稅和經(jīng)濟(jì)發(fā)展的回歸模型,研究結(jié)果發(fā)現(xiàn)環(huán)境稅的征收雖然能降低污染排放,但并不能促進(jìn)經(jīng)濟(jì)增長(zhǎng).

大量研究及論證表明,環(huán)境稅的征收在降低化石能源使用和加速綠色能源創(chuàng)新方面具有重要作用.王良斌等[8]認(rèn)為環(huán)境稅的征收短期會(huì)迫使高耗能和高污染的小型企業(yè)放棄部分生產(chǎn)線,長(zhǎng)期將會(huì)使大中型企業(yè)加速能源使用效率等創(chuàng)新.Chien等[9]通過(guò)截面相關(guān)性分析以及AMG和CCEMG穩(wěn)健性檢驗(yàn),發(fā)現(xiàn)環(huán)境稅將對(duì)遏制二氧化碳的排放起到顯著作用.朱新玲等[10]以長(zhǎng)江經(jīng)濟(jì)帶的準(zhǔn)環(huán)境稅數(shù)據(jù)為基礎(chǔ),研究發(fā)現(xiàn)環(huán)境稅的實(shí)施將有效降低單位GDP的能源消耗,從側(cè)面減少化石能源的使用.

前人在上述研究中主要利用統(tǒng)計(jì)學(xué)方法來(lái)模擬預(yù)測(cè)因素之間的線性關(guān)系,而現(xiàn)實(shí)中,不同因素之間都可能存在某種聯(lián)系,因此,為了研究多種因素之間的關(guān)系,本文引入非線性動(dòng)力系統(tǒng).Sun等[11]建立關(guān)于能源供需的動(dòng)力系統(tǒng)模型來(lái)研究中國(guó)中西部之間的能源供需問(wèn)題.Wang等[12]通過(guò)建立關(guān)于能源價(jià)格,能源供給和經(jīng)濟(jì)增長(zhǎng)的三維動(dòng)力系統(tǒng)模型來(lái)研究政府不同控制策略對(duì)能源強(qiáng)度的影響.寧成[13 ]基于資源、環(huán)境和經(jīng)濟(jì)發(fā)展之間的三維動(dòng)力系統(tǒng)關(guān)系,研究發(fā)現(xiàn)要想在發(fā)展經(jīng)濟(jì)的同時(shí)保護(hù)環(huán)境,就一定要控制資源消耗.

綜上所述,為了更準(zhǔn)確地描述環(huán)境稅、化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展之間的關(guān)系,本文將建立三者之間的三維動(dòng)力系統(tǒng)模型,研究環(huán)境稅的征收對(duì)經(jīng)濟(jì)發(fā)展和能源強(qiáng)度的影響.

1 模型的建立

環(huán)境稅的征收無(wú)疑會(huì)對(duì)化石能源消費(fèi)起到一定抑制作用,而且會(huì)大幅減少污染性投資,化石能源消費(fèi)和污染性投資的減少又會(huì)影響經(jīng)濟(jì)的發(fā)展,但征收環(huán)境稅的同時(shí),政府通過(guò)政策調(diào)控促進(jìn)企業(yè)技術(shù)升級(jí),提供更多的綠色能源就業(yè)機(jī)會(huì),又會(huì)帶動(dòng)經(jīng)濟(jì)發(fā)展,此外征收的環(huán)境稅也會(huì)直接充實(shí)國(guó)家收入.因此,環(huán)境稅和化石能源消耗,經(jīng)濟(jì)發(fā)展之間存在著復(fù)雜的關(guān)系.本文基于三者之間這種復(fù)雜的關(guān)系,建立了一個(gè)含有環(huán)境稅、化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展的非線性動(dòng)力系統(tǒng)——環(huán)境稅系統(tǒng),并借助非線性系統(tǒng)相關(guān)理論研究他們之間的關(guān)系.可以用下面的微分方程描述環(huán)境稅系統(tǒng)

(1)

其中，x(t)為隨時(shí)間變化的化石能源消耗,y(t)為隨時(shí)間變化的經(jīng)濟(jì)發(fā)展(GDP),z(t)為隨時(shí)間變化的環(huán)境稅，ai,bi,ci,M為正常數(shù),t∈I,I為一個(gè)經(jīng)濟(jì)時(shí)期，a1為化石能源自身消耗速度,a2為經(jīng)濟(jì)發(fā)展對(duì)化石能源消耗的促進(jìn)系數(shù),a3為經(jīng)濟(jì)發(fā)展和環(huán)境稅共同對(duì)化石能源消耗的影響系數(shù)，b1為化石能源消耗和經(jīng)濟(jì)發(fā)展之間的關(guān)系,b2為經(jīng)濟(jì)發(fā)展對(duì)自身的促進(jìn)作用,b3為適當(dāng)?shù)恼{(diào)控對(duì)環(huán)保產(chǎn)業(yè)的促進(jìn)作用以及直接征收環(huán)境稅對(duì)經(jīng)濟(jì)的帶動(dòng)作用,M為環(huán)境稅對(duì)經(jīng)濟(jì)發(fā)展的閾值，c1為化石能源消耗對(duì)環(huán)境稅的稅率,c2表示市場(chǎng)經(jīng)濟(jì)的通貨膨脹情況,c3表示政府對(duì)環(huán)境稅的控制作用.

能源強(qiáng)度是指單位GDP所需要消耗的能源,一個(gè)國(guó)家或地區(qū)的能源強(qiáng)度可以用下面的公式表示

2 模型的分析

本節(jié)將依據(jù)耗散性、平衡點(diǎn)、相圖、李雅普諾夫指數(shù)圖和分岔圖來(lái)分析環(huán)境稅動(dòng)力系統(tǒng)的行為.

2.1 耗散性分析

環(huán)境稅系統(tǒng)的耗散度為

(2)

根據(jù)勞斯-霍爾維茨理論[14],當(dāng)耗散度為負(fù)時(shí),稱系統(tǒng)是耗散的,此時(shí)環(huán)境稅系統(tǒng)可能會(huì)出現(xiàn)混沌現(xiàn)象.

2.2 平衡點(diǎn)分析

環(huán)境稅動(dòng)力系統(tǒng)描述了環(huán)境稅和化石能源消費(fèi),經(jīng)濟(jì)發(fā)展之間復(fù)雜的動(dòng)力學(xué)關(guān)系,基于MATLAB數(shù)值模擬,發(fā)現(xiàn)當(dāng)環(huán)境稅系統(tǒng)取不同參數(shù)時(shí),系統(tǒng)將出現(xiàn)不同的動(dòng)力學(xué)特征.經(jīng)過(guò)大量的數(shù)值模擬發(fā)現(xiàn)當(dāng)環(huán)境稅系統(tǒng)選定下列參數(shù)時(shí),環(huán)境稅系統(tǒng)將出現(xiàn)混沌行為.

表1 環(huán)境稅系統(tǒng)模擬參數(shù)

此時(shí),環(huán)境稅系統(tǒng)的耗散度為

(3)

所以環(huán)境稅系統(tǒng)是耗散的,并且根據(jù)指數(shù)形式e-(a1-b2+c3)收斂.即,對(duì)于一個(gè)初始體積為V(0)的體積元來(lái)說(shuō),在時(shí)刻t這個(gè)體積元將變成V(0)e-(a1-b2+c3)t,所以當(dāng)t趨于無(wú)窮時(shí),此體積元將以指數(shù)速率-7.4收縮到零.所以在無(wú)窮時(shí)間點(diǎn)時(shí)系統(tǒng)的所有軌跡都將收縮到一個(gè)體積為零的極限點(diǎn)集上,并且它的動(dòng)力學(xué)行為都會(huì)被固定在一個(gè)吸引子上,這就說(shuō)明了吸引子的存在性.

通過(guò)穩(wěn)定性理論的研究方法,令環(huán)境稅系統(tǒng)的左邊微分項(xiàng)為零,并將表1參數(shù)代入系統(tǒng)(1)可得

(4)

求得系統(tǒng)(1)的平衡點(diǎn)為

(5)

在S0(0，0，0)，根據(jù)系統(tǒng)的雅各布矩陣

(6)

得到特征值為λ1=5.968，λ2=-2.024，λ3=-3.044,由于λ1為正實(shí)數(shù),λ2和λ3為負(fù)實(shí)數(shù),因此根據(jù)勞斯-赫爾維茨理論平衡點(diǎn)S0(0，0，0)為鞍點(diǎn),也是不穩(wěn)定的點(diǎn).

在S1(1.162,-0.302,-4.893)，根據(jù)系統(tǒng)的雅各布矩陣

(7)

得到特征值為λ1=-1.740，λ2=1.320+7.254i，λ3=1.320-7.254i，由于λ1為負(fù)實(shí)數(shù),λ2和λ3為有實(shí)部的共軛虛根,因此根據(jù)勞斯-赫爾維茨理論平衡點(diǎn)S1(1.162,-0.302,-4.893)為不穩(wěn)定的鞍點(diǎn).

在S2(1.224,0.281,4.582),根據(jù)系統(tǒng)的雅各布矩陣

(8)

得到特征值為λ1=-1.619，λ2=1.260+7.327i，λ3=1.260-7.327i，由于λ1為負(fù)實(shí)數(shù),λ2和λ3有實(shí)部的共軛虛根,因此根據(jù)勞斯-赫爾維茨理論平衡點(diǎn)S2(1.224,0.281,4.582)為不穩(wěn)定的鞍點(diǎn).

在S3(-4.751,3.105,-2.227)，根據(jù)系統(tǒng)的雅各布矩陣

(9)

得到特征值為λ1=-9.340，λ2=5.120+4.368i，λ3=1.260+7.327i，由于λ1為負(fù)實(shí)數(shù),λ2和λ3有實(shí)部的共軛虛根,因此根據(jù)勞斯-赫爾維茨理論平衡點(diǎn)S3(-4.751,3.105,-2.227)為不穩(wěn)定的鞍點(diǎn).

在S4(-4.778,-2.360,1.875)，根據(jù)系統(tǒng)的雅各布矩陣

(10)

得特征值為λ1=-7.926，λ2=4.413+4.631i，λ3=1.260+7.327i，由于λ1為負(fù)實(shí)數(shù),λ2和λ3有實(shí)部的共軛虛根,因此，根據(jù)勞斯-赫爾維茨理論平衡點(diǎn)S4(-4.778,-2.360,1.875)為不穩(wěn)定的鞍點(diǎn).

根據(jù)以上分析,本系統(tǒng)的五個(gè)平衡點(diǎn)都是不穩(wěn)定的鞍點(diǎn).

由于現(xiàn)實(shí)世界是一個(gè)復(fù)雜多變的世界,因此,混沌的出現(xiàn)是驗(yàn)證系統(tǒng)是否符合現(xiàn)實(shí)世界的重要條件.接下來(lái)將分別從環(huán)境稅系統(tǒng)的相圖,李雅普諾夫指數(shù)圖和分叉圖來(lái)依次判斷系統(tǒng)的混沌特性.由于環(huán)境稅的目的是為了降低能源強(qiáng)度,減少環(huán)境污染,因此,本文將直接選取與能源強(qiáng)度有關(guān)的參數(shù)b2作為研究參數(shù).

2.3 系統(tǒng)的相圖

將系統(tǒng)(1)的初始值選為(0.2,0.4,0.8),應(yīng)用MATLAB模擬分析可以得到系統(tǒng)的三維相圖如圖1所示.

圖1 三維相圖

從上述相圖可以看到環(huán)境稅、化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展隨時(shí)間發(fā)展變化是無(wú)序的,也就是說(shuō)環(huán)境稅動(dòng)力系統(tǒng)可能存在混沌現(xiàn)象.

2.4 李雅普諾夫指數(shù)圖和分叉圖

李雅普諾夫特征指數(shù)指的是混沌系統(tǒng)對(duì)初值敏感判斷的一個(gè)指標(biāo),它表示相鄰軌線間的平均發(fā)散(分離)率,如果非線性系統(tǒng)中的李雅普諾夫指數(shù)大于零，那么此系統(tǒng)可能存在混沌行為,系統(tǒng)的分叉圖代表系統(tǒng)隨時(shí)間演化可能出現(xiàn)的情況.系統(tǒng)(1)的李雅普諾夫指數(shù)圖如圖2,分叉圖如圖3.

根據(jù)此系統(tǒng)李雅普諾夫指數(shù)圖可以發(fā)現(xiàn),該系統(tǒng)在1.44

最后,通過(guò)分岔圖來(lái)觀察系統(tǒng)的混沌現(xiàn)象,經(jīng)觀察系統(tǒng)在b2=1.44時(shí)發(fā)生突變分岔,即b2=1.44是系統(tǒng)穩(wěn)定和不穩(wěn)定的臨界點(diǎn).因此,當(dāng)b2>1.44時(shí)環(huán)境稅系統(tǒng)可能出現(xiàn)混沌現(xiàn)象.綜上所述,當(dāng)1.44

3 數(shù)值模擬

本節(jié)將通過(guò)實(shí)際數(shù)據(jù)展示環(huán)境稅、化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展在此系統(tǒng)下的演化情況.其中數(shù)據(jù)來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》并通過(guò)歸一法對(duì)數(shù)據(jù)進(jìn)行預(yù)處理,然后利用遺傳算法計(jì)算出實(shí)際系統(tǒng)的參數(shù),并研究當(dāng)政策發(fā)生改變從而引起環(huán)境稅系統(tǒng)的參數(shù)改變時(shí),系統(tǒng)可能會(huì)發(fā)生的變化.

3.1 數(shù)據(jù)來(lái)源及預(yù)處理

根據(jù)中華人民共和國(guó)環(huán)境保護(hù)稅法,環(huán)境稅的征收主要依據(jù)應(yīng)稅大氣污染物的排放量,應(yīng)稅水污染物的排放量以及應(yīng)稅固體污染物的排放量.所以,依據(jù)每種污染物的稅率,本節(jié)可以將征收的環(huán)境稅簡(jiǎn)化為1.4×廢水的排放量+1.14×廢氣的排放量+5×固體廢物的排放量.由于我國(guó)自2018年才實(shí)際征收環(huán)境稅,因此2018年之前的環(huán)境稅數(shù)據(jù)采用推測(cè)數(shù)據(jù),化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展數(shù)據(jù)直接來(lái)自中國(guó)統(tǒng)計(jì)年鑒.

經(jīng)歸一化以后得到表2所示數(shù)據(jù).

表2 歸一化數(shù)據(jù)

3.2 遺傳算法和參數(shù)辨識(shí)

遺傳算法(GA)是一種用來(lái)識(shí)別系統(tǒng)參數(shù)的程序算法,在給定系統(tǒng)模型和實(shí)際數(shù)據(jù)的情況下,可以很好地識(shí)別出實(shí)際系統(tǒng)的參數(shù).因此,本文應(yīng)用遺傳算法來(lái)識(shí)別環(huán)境稅動(dòng)力系統(tǒng)的參數(shù).為了使用GA,將遺傳算法重寫成以下向量形式

(11)

其中，X(t)代表系統(tǒng)狀態(tài)向量,α代表系統(tǒng)的參數(shù)向量.將它離散化得到

X(k+1)=X(k)+f(X(k))=F(X(k),α).

(12)

因此,環(huán)境稅系統(tǒng)的參數(shù)識(shí)別就等同于下面的非線性最優(yōu)化問(wèn)題

(13)

其中，α=(α1,α2,…,α10)，i是需要辨識(shí)的參數(shù)個(gè)數(shù).

本系統(tǒng)采用MATLAB對(duì)參數(shù)進(jìn)行識(shí)別,人口總數(shù)設(shè)置為200,交叉率為0.4,突變率為0.1,容錯(cuò)率為10-6.辨識(shí)出的參數(shù)集如表3所示,此參數(shù)不同于原來(lái)設(shè)置的參數(shù),因此,在此參數(shù)下環(huán)境稅系統(tǒng)可能會(huì)發(fā)生不同的行為.

表3 實(shí)際系統(tǒng)參數(shù)

將表3參數(shù)作為環(huán)境稅系統(tǒng)的參數(shù),并設(shè)置2011年數(shù)據(jù)作為系統(tǒng)的初始值,經(jīng)MATLAB模擬得到環(huán)境稅系統(tǒng)的相圖如圖4所示.

圖4 實(shí)際系統(tǒng)的相圖

根據(jù)相圖所示,此系統(tǒng)最終趨于穩(wěn)定,因此,可以認(rèn)定環(huán)境稅動(dòng)力系統(tǒng)以及模擬出的參數(shù)符合實(shí)際情況.

3.3 情景分析

本節(jié)采用靈敏度分析法研究當(dāng)一個(gè)變量發(fā)生改變時(shí)系統(tǒng)的變化.由于本系統(tǒng)主要研究環(huán)境稅變化時(shí)對(duì)經(jīng)濟(jì)發(fā)展和能源強(qiáng)度的影響,因此,本節(jié)將研究與環(huán)境稅和能源強(qiáng)度有關(guān)的參數(shù)b3和c1改變時(shí)環(huán)境稅系統(tǒng)的變化.

3.3.1 政府調(diào)控對(duì)實(shí)際系統(tǒng)的影響

在本環(huán)境稅系統(tǒng)中,參數(shù)b3代表當(dāng)稅收增加,政府通過(guò)適當(dāng)政策調(diào)控,引導(dǎo)綠色產(chǎn)業(yè)發(fā)展,對(duì)我國(guó)經(jīng)濟(jì)發(fā)展的影響.

固定環(huán)境稅系統(tǒng)其他參數(shù)如表3,當(dāng)政府增加征收環(huán)境稅時(shí),部分企業(yè)會(huì)主動(dòng)轉(zhuǎn)型,進(jìn)行技術(shù)產(chǎn)業(yè)升級(jí),政府也會(huì)采取適當(dāng)措施鼓勵(lì)綠色產(chǎn)業(yè)發(fā)展,此時(shí)參數(shù)b3將增加；反之,當(dāng)政府對(duì)環(huán)境監(jiān)管不嚴(yán),減少鼓勵(lì)綠色產(chǎn)業(yè)發(fā)展的措施時(shí),參數(shù)b3將適當(dāng)降低.經(jīng)MATLAB模擬研究發(fā)現(xiàn),當(dāng)0.0263≤b3≤0.3243時(shí),環(huán)境稅系統(tǒng)均會(huì)保持穩(wěn)定,圖5表示b3在不同數(shù)值時(shí)經(jīng)濟(jì)發(fā)展演化趨勢(shì).從圖5可以看出當(dāng)政府采取適當(dāng)措施鼓勵(lì)綠色產(chǎn)業(yè)發(fā)展時(shí),可以促進(jìn)經(jīng)濟(jì)發(fā)展；圖6表示b3在不同數(shù)值時(shí)的能源強(qiáng)度演化趨勢(shì),從圖6可以看出當(dāng)政府采取適當(dāng)措施鼓勵(lì)綠色產(chǎn)業(yè)發(fā)展時(shí),能源強(qiáng)度也會(huì)適當(dāng)?shù)亟档?這說(shuō)明加強(qiáng)政府在綠色產(chǎn)業(yè)方面的投入,鼓勵(lì)環(huán)保產(chǎn)業(yè)發(fā)展不僅可以發(fā)展經(jīng)濟(jì)還可以降低能源強(qiáng)度.

3.3.2 環(huán)境稅率對(duì)實(shí)際系統(tǒng)的影響

本環(huán)境稅系統(tǒng)中,參數(shù)c1代表化石能源消耗對(duì)環(huán)境稅的綜合稅率,當(dāng)化石能源消耗增多時(shí),政府通過(guò)提高環(huán)境稅稅率c1來(lái)控制環(huán)境污染.現(xiàn)固定環(huán)境稅系統(tǒng)其他參數(shù)如表3,并改變參數(shù)c1,經(jīng)MATLAB模擬研究發(fā)現(xiàn),當(dāng)0.0153≤c1≤0.2672時(shí),環(huán)境稅系統(tǒng)均會(huì)保持穩(wěn)定.圖7表示c1在不同參數(shù)時(shí)的經(jīng)濟(jì)發(fā)展演化趨勢(shì),根據(jù)圖7可以發(fā)現(xiàn),當(dāng)政府提高環(huán)境稅稅率即增加c1時(shí),將抑制我國(guó)經(jīng)濟(jì)的發(fā)展,這是因?yàn)榛茉茨壳叭栽谖覈?guó)經(jīng)濟(jì)發(fā)展中占重要地位；圖8表示c1在不同參數(shù)時(shí)能源強(qiáng)度的變化,根據(jù)圖8可以發(fā)現(xiàn),當(dāng)政府提高環(huán)境稅率即增加c1時(shí),將適當(dāng)降低化石能源強(qiáng)度.總而言之,當(dāng)政府提高環(huán)境稅稅率即增加c1時(shí),將抑制我國(guó)經(jīng)濟(jì)的發(fā)展,但能夠降低能源強(qiáng)度.

4 總結(jié)

建立了一個(gè)關(guān)于環(huán)境稅、化石能源消費(fèi)和經(jīng)濟(jì)發(fā)展的三維動(dòng)力系統(tǒng)模型,通過(guò)非線性系統(tǒng)分析,證明了系統(tǒng)符合中國(guó)情況,并利用遺傳算法模擬出符合實(shí)際情況的系統(tǒng)參數(shù).通過(guò)數(shù)值模擬,研究得出當(dāng)政府提高對(duì)綠色產(chǎn)業(yè)的扶持力度,鼓勵(lì)環(huán)保產(chǎn)業(yè)發(fā)展時(shí),不僅可以促進(jìn)經(jīng)濟(jì)發(fā)展,而且可以降低我國(guó)發(fā)展的能源強(qiáng)度；而如果僅僅只是提高環(huán)境稅的征收稅率,雖然可以降低我國(guó)發(fā)展的能源強(qiáng)度,但會(huì)抑制經(jīng)濟(jì)的發(fā)展.