全球专业中文经管百科,由121,994位网友共同编写而成,共计436,047个条目

邏輯斯蒂方程

用手机看条目

出自 MBA智库百科(https://wiki.mbalib.com/)

邏輯斯蒂方程(Logistic Equation)

目錄

邏輯斯蒂方程的推導

  當一種新產品剛面世時,廠家和商家總是採取各種措施促進銷售。他們都希望對這種產品的推銷速度做到心中有數,這樣廠家便於組織生產,商家便於安排進貨。怎樣建立數學模型描述新產品推銷速度呢?

  首先要考慮社會的需求量.社會對產品的需求狀況一般依如下兩個特性確定:

  1. 對產品的需求有一個飽和水平.當產品需求量達到一定數量時,對這種產品的需求也飽和了,設飽和水平為a;

  2. 假設在時刻t,社會對產品的需求量為x=x(t),需求的增長速度dx/dt正比於需求量x(t)與需求接近飽和水平的程度a-x(t)之乘積,記比例繫數為k ;

  根據上述實際背景的兩個特征,可建立如下微分方程:

  \boldsymbol{\frac{dx}{dt}=kx(a-x)}.......................(1)

  分離變數,得:

  \boldsymbol{\frac{x}{x(a-x)}=kdt}

  兩邊積分,得:

  \boldsymbol{\frac{dx}{a-x}=e^{ac}e^{akt}=Ae^{bt}}

  其中:\boldsymbol{A=e^{ac}}

  從而,通解為:

  \boldsymbol{x(t)=\frac{aAe^{bt}}{Ae^{bt}+1}=\frac{a}{1+Be^{-bt}}}......(2)

  其中\boldsymbol{B=\frac{1}{A}=e^{-x}},B和b為正常數,可由初始條件確定。式(1)稱為邏輯斯蒂方程(1ogistic equation),式(2)稱為邏輯斯蒂曲線。

邏輯斯蒂方程的基本性質

  1.當t=O時,x(t)的值為:X_0=\frac{a}{1+B}

  2.x(t)的增長率\frac{dx(t)}{dt}=\frac{abBe^{-bt}}{(1+B^{-bt})^2}>0,因此,x(t)是增函數;

  3.當B值較大而t較小時,\boldsymbol{Be^{-bt}}將很大,1+Be^{-bt} \approx Be^{-bt},於是

  x(t) \approx \frac{a}{Be^{-bt}}=\frac{a}{b}e^{bt}

  x(t)近似於依指數函數增大,銷售速度不斷增大;

  4.當t增大以後,\boldsymbol{Be^{-bt}}越來越接近於零,分母越來越接近於1,銷售速度開始下降,x(t)的值接近於a(飽和值)。

邏輯斯蒂方程的應用

  1.人口限制增長問題

  人口的增長不是呈指數型增長的,這是由於環境的限制、有限的資源和人為的影響,最終人口的增長將減慢下來。實際上,人口增長規律滿足邏輯斯蒂方程。

  2. 信息傳播問題

  所謂信息傳播可以是一則新聞,一條謠言或市場上某種新商品有關的知識,在初期,知道這一信息的人很少,但是隨時間的推移,知道的人越來越多,到一定時間,社會上大部分人都知道了這一信息.這裡的數量關係可以用邏輯斯蒂方程來描述。若以t表示從信息產生算起的時間,P表示已知信息的人口比例,則邏輯斯蒂方程變為:

  P=\frac{1}{1+Be^{-bt}}...................(3)

  例如,當某種商品調價的通知下達時,有10%的市民聽到這一通知,2小時以後,25%的市民知道了這一信息,由邏輯斯蒂方程可算出有75%的市民瞭解這一情況所需要的時間。

  在方程(3)中,由t=0時,P=10%可得 B=9;再由t=2時,P=25%可得,b=\frac{1}{2}ln3

  當P=75%時,有:

  75%=\frac{1}{1+9e^{-\frac{1}{2}ln3}}

  解得t=6,即6小時後,全市有75%的人瞭解這一通知。

  3.商品銷售預測問題

  例如,某種商品的銷售,開始時,知道的人很少,銷售量也很小。當這種商品信息傳播出去後,銷售量大量增加,到接近飽和時銷售量增加極為緩慢。比如,這種商品飽和量估計a=500(百萬件),大約5年可達飽和,常數b經測定為b=lnl0,B=100。下麵我們來預測一下第3年末的銷售量是多少。

  由x=\frac{a}{1+Be^{-bt}},有:

  x=\frac{500}{1+100e^{-ln10*3}}=\frac{500}{1+100(\frac{1}{10})^3}=454.5(百萬件)

  所以第三年末的市場銷售量大約為454.5百萬件,這樣可以做到有計劃地生產。

  邏輯斯蒂方程的應用比較廣泛。如果問題的基本數量特征是:在時間t很小時,呈指數型增長,而當t增大時,增長速度就下降,且越來越接近於一個確定的值,這類問題可以用邏輯斯蒂方程加以解決。

邏輯斯蒂方程案例分析

案例一:基於邏輯斯蒂方程的垃圾郵件過濾特征方法的研究[1]

  1.邏輯斯蒂(Logistic)方程

  Y=\frac{K}{1+ae^{bX}}(1)

  Logistic曲線其特點是開始增長緩慢,而在以後的某一範圍內迅速增長,達到某限度後,增長又緩慢下來∞3.曲線略呈拉長的“S”型.Logistic增長的速度函數圖像為單峰曲線,這表明Logistic曲線增長過程分為三個階段:慢一快一慢。

  在實際中,特征項也根據重要程度的不同分為三部分,並且與Logistic曲線的三個階段相對應:第一部分是在垃圾郵件分類中貢獻甚微的特征項,特征權值不需要明顯的變化,全部是不予選擇的特征項,對應Logistic曲線變化幅度較小的第一階段;第二部分是要進行甄別的特征項,需要區別特征項重要程度強弱,特征權值變化幅度越大越有利於特征項的選擇,與曲線的變化幅度最快的第二階段對應;第三個部分是必須選取的高特征權值的特征項,對應曲線第三階段,這部分特征項同樣不需要顯著的差別,全部是要選擇的特征項.2.參數設定。

  2.參數設定

  以影響特征選擇方法的三個要素為因數,設邏輯斯蒂(Logistic)方程的參數X為:

  X=\alpha+\beta+\sqrt{\gamma}(2)

  根據(8)式測定垃圾郵件中特征項集的X值分佈情況:P(X<100)>90%,即特征權值小於100的特征項占待選擇特征項集的95%以上,可判定需要甄別的特征項集中在X值小於lOO的特征項集合內,由此設定一組樣本值(X_1,Y_1),(X_2,Y_2),\cdots(\cdotsX_10,Y_10),,其數值如下表所示.以此分析此結果是否能以Logistic方程描述。

  表:X,Y參數值表

X10203040506070809095
y9121630506577849093

  在(1)式中,當X→∞,\lim_{x \to \infty}=K,所以特征項重要程度的無限延長的終極量是K.假定這組值的K值為100(%),即特征項的最重要程度為100%.將(1)式移項得:\frac{K-Y}{Y}=ae^{-bX}(3)

  令Y^\prime=\ln(\frac{K-Y}{Y}),則由最小平方法估計的Y^\prime

  Y^\prime=\ln a-bX(4)

  因此,Y和X對於Logistic方程的符合度,可由Y^\prime和X的線性相關係數\gamma Y^\prime X給出:

  \gamma Y^\prime X=\frac{\sum XY^\prime-\frac{1}{n}(\sum X)(\sum Y^\prime)}{\sqrt{\left[\sum X^2-\frac{1}{n}(\sum X^2)\right]\left[\sum(Y^\prime)^2-\frac{1}{n}(\sum Y^\prime)^2\right]}}=\frac{\sum xy^\prime}{\sqrt{\sum X^2\cdot\sum(y^\prime)^2}}(5)

  (5)中,

  y^\prime=\ln(\frac{K-Y}{Y})-\overline{y^\prime},\overline{y^\prime}=\sum\ln(\frac{K-Y}{Y})/n(6)

  x=X-\overline{x},\overline{x}=\sum X/n(7)

  以為樣本容量。回歸參數a和b的樣本估計值為:

  -b=\frac{\sum XY^\prime-\frac{1}{n}(\sum X)(\sum Y^\prime)}{\sum X^2-\frac{1}{n}(\sum X)^2}=\frac{\sum xy^\prime}{\sum X^2}(8)

  \ln a=\overline{y^\prime}+b\overline{x}(9)

  根據表2以及式(6)(7)求得:

  \sum x^2=7822.5,\sum(y^\prime)^2=28.6894,\sum xy^\prime=-471.7887

  故根據(11)得:\gamma Y^\prime X=-0.9958為極顯著.故知X和Y的關係以Logistic方程配合是合適的.進一步根據(6)和(7)求得b值為0.0603和a值為23.1316.那麼所求Logistic方程為:

  Y=\frac{100}{1+23.1316e^{-0.0603\cdot X}}

  式(10)的曲線圖如圖所示

  Image:Logistic方程曲线图.jpg

  3.實驗和結果

  採用的是CCERT的中文郵件樣本集,該樣本集由中國教育科研網收集並維護.該數據集在國內和教育網內是常用的一個衡量反中文垃圾郵件技術的數據集.採用數據集的8000封垃圾郵件以及8000封正常郵件.選擇貝葉斯方法進行垃圾郵件分類,分類結果利用查全率(Recall)、準確率(Precision)和F值評測過濾系統性能。

  查全率(R)=正確過濾掉的郵件數/應該過濾掉的垃圾郵件數;

  準確率(P)=正確過濾掉的郵件數/實際過濾掉的郵件數;

  \frac{2RP}{R+P};值實際是查全率和準確率的調和平均,綜合成一個指標.實驗結果如下表。

  表:統評測結果

信息熵互信息幾率比期望交叉熵CHI統計量Logistic方程
Recall0.84460.86150.89540.87020.87370.9366
Precision0.96980.95740.94210.93630.93160.9493
F0.90280.90370.91810.90200.90170.9429

  實驗結果表明,利用Logistic方程進行分類時的F值有顯著的提高,說明過濾系統整體性能有所提高,其中查全率有明顯的提高,說明這種方法檢出垃圾郵件的能力較高,使漏網的垃圾郵件減少,但是準確率沒有提升,說明合法郵件誤判為垃圾郵件的可能性沒有明顯降低。

參考文獻

  1. 妥曉娜,吳承勇.基於邏輯斯蒂方程的垃圾郵件過濾特征方法的研究[J].內蒙古大學學報(自然科學版),2010,41(4):450-455
本條目對我有幫助106
MBA智库APP

扫一扫,下载MBA智库APP

分享到:
  如果您認為本條目還有待完善,需要補充新內容或修改錯誤內容,請編輯條目投訴舉報

評論(共8條)

提示:評論內容為網友針對條目"邏輯斯蒂方程"展開的討論,與本站觀點立場無關。
111.78.92.* 在 2010年4月21日 08:32 發表

謝謝你的講解,講解很詳細!

回複評論
218.19.34.* 在 2010年4月27日 12:28 發表

謝謝

回複評論
123.150.197.* 在 2010年12月28日 10:33 發表

微分方程求解過程有點問題吧

回複評論
125.73.52.* 在 2011年3月25日 12:01 發表

非常感謝,,多重改良邏輯斯蒂,都改些什麼呢 ?

回複評論
58.254.93.* 在 2011年7月3日 23:35 發表

用傳統方法推啊 我看不懂你的方法啊

回複評論
112.80.183.* 在 2012年3月15日 14:34 發表

有圖嗎

回複評論
乔巴 (討論 | 貢獻) 在 2012年3月15日 18:08 發表

112.80.183.* 在 2012年3月15日 14:34 發表

有圖嗎

附圖並添加了新的案例,希望對您有幫助!

回複評論
121.232.238.* 在 2012年3月29日 19:48 發表

微分方程順序錯了

回複評論

發表評論請文明上網,理性發言並遵守有關規定。

打开APP

以上内容根据网友推荐自动排序生成

官方社群
下载APP

闽公网安备 35020302032707号