預(yù)測(cè)模型

1、灰色預(yù)測(cè)
適用范圍：它是一種基于小樣本數(shù)據(jù)進(jìn)行預(yù)測(cè)的模型。灰色系統(tǒng)理論它主要針對(duì)時(shí)間序列少、統(tǒng)計(jì)數(shù)據(jù)少、信息不完全的分析和建模。

灰色系統(tǒng)：我們將信息完全未確定的系統(tǒng)為黑色系統(tǒng)，信息完全確定的為白色系統(tǒng)，既有已知信息又有未知信息的為灰色系統(tǒng)

灰色系統(tǒng)是對(duì)離散序列建立的微分方程，它通過(guò)對(duì)原始數(shù)據(jù)的整理來(lái)尋求其變化規(guī)律，也就是灰色序列的生產(chǎn)。

灰色生成：累加生成、累減生成、均值生成、級(jí)比生成

1數(shù)據(jù)檢驗(yàn)
計(jì)算數(shù)列的級(jí)比，所有級(jí)比都落在可覆蓋區(qū)間（e^-2/n+1,e2/n+1)內(nèi)，則可以建立GM（1，1）模型

2構(gòu)建模型
生成x1的鄰值生成序列z1，定義模型為：

d(k)+az1(k)=b

（a為發(fā)展系數(shù)，z1(k)為白化背景值，b為灰作用量）

1.1預(yù)測(cè)檢驗(yàn)
殘差檢驗(yàn)：

s(k)=x0(k)-x^0(k)/x0(k)

所有的|s(k)|<0.1，認(rèn)為達(dá)到較高的要求，

|s(k)|<0.2認(rèn)為達(dá)到一般要求

級(jí)比偏差值檢驗(yàn)：

p(k)=1-(1-0.5a)/(1+0.5a)*s(k)/s^(k)

|p(k)|<0.1認(rèn)為達(dá)到較高要求

|p(k)|<0.2認(rèn)為達(dá)到一般要求

2、馬爾可夫鏈

狀態(tài)空間從一個(gè)狀態(tài)到另一個(gè)狀態(tài)轉(zhuǎn)換的隨機(jī)過(guò)程。它必須滿足兩個(gè)條件：

1 任意時(shí)刻可以用有限個(gè)可能狀態(tài)之一描述

2 無(wú)后效性

狀態(tài)向量：x=（x1，x2，x3,…xk)

各個(gè)元素分別表示第幾次觀測(cè)時(shí)第i個(gè)狀態(tài)的概率

系統(tǒng)的可能狀態(tài)為k個(gè)

p 轉(zhuǎn)移概率矩陣

x(n+1)=x(0)P(n)

3、線性回歸

線性回歸在因變量和一個(gè)或多個(gè)自變量之間建立一種關(guān)系 a為截距，b表示斜率，然后衡量自變量與因變量的相關(guān)程度 較多適用于數(shù)據(jù)多、計(jì)算簡(jiǎn)單。

其表達(dá)形式為：y=w’x+e.其中e為誤差服從均值為0的正態(tài)分布

最小二乘法：以殘差平方和最小確定直線的位置

調(diào)用python的sklearn可以實(shí)現(xiàn)線性回歸算法，訓(xùn)練模型，從而求出截距和斜率

//灰色預(yù)測(cè) import numpy as np import pandas as pd def ratio(x0):n=len(x0)ratio=[x0[i]/x0[i+1] for i in range(n-1)]print(f"級(jí)比:{ratio}")min_r=min(ratio)max_r=max(ratio)r=[np.exp(-2/(n+2)),np.exp(2/(n+2))]print(r)if min_r >= r[0] and max_r <= r[1]:print("級(jí)比滿足要求")else:print("級(jí)比超過(guò)灰色模型預(yù)測(cè)范圍")return ratio,r def predict(x0):n=len(x0)x1=np.cumsum(x0) #累加序列z=(x1[:len(x1)-1]+x1[1:])/2 z1=z.reshape(len(z),1)B=np.append(-z1,np.ones_like(z1),axis=1)Y=x0[1:].reshape(n-1,1)u=np.dot(np.dot(np.linalg.inv(np.dot(B.T,B)),B.T),Y)x12=np.zeros(n)x02=np.zeros(n)x12[0]=x02[0]=x0[0]for i in range(1,n):x12[i]=(x0[0]-u[1]/u[0])*np.exp(-u[0]*i)+u[1]/u[0]for i in range(1,n):x02[i]=x12[i]-x12[i-1]print(f"原始值：{x0}")print(f"預(yù)測(cè)值：{x02}")return x02,x12,u def accuracy(x0,x02,ratio,u):n=len(x0)deltas=[]for i in range(n):eps=x0[i]-x02[i]delta=abs(eps/x0[i])deltas.append(delta)print(f"殘差:{deltas}")ratio2=[x02[i]/x02[i+1] for i in range(n-1)]rhos=[]for i in range(n-1):rho=1-(1-0.5*u[0]/u[1])/(1+0.5*u[0]/u[1])*(ratio[i]/ratio2[i])rhos.append(rho)print(f"級(jí)比偏差：{rhos}") if __name__=='__main__':x0=np.array([71.1,72.4,72.4,72.1,71.4,72.0,71.6])ratio,r=ratio(x0)x02,x12,u=predict(x0)accuracy(x0,x02,ratio,u) //馬爾可夫鏈 import numpy as np import random as rm def markov():start=np.array([1,0])trans=np.array([[0.286,0.714],[0.455,0.545]])restmp=startfor i in range(15):res=np.dot(restmp,trans)print(res)restmp=res markov() //線性回歸 import pandas as pd import numpy as np from sklearn.linear_model import LinearRegression Model=LinearRegression() data={"學(xué)習(xí)時(shí)間":[0.50,0.75,1.00,1.25,1.50,1.75,2.00,2.25,2.50],"成績(jī)":[10,20,13,43,20,22,33,50,62]} df=pd.DataFrame(data) x=df['學(xué)習(xí)時(shí)間'] y=df['成績(jī)'] r=df.corr() print(f"相關(guān)系數(shù):{r}") x=x.values.reshape(-1,1) y=y.values.reshape(-1,1) Model.fit(x,y) a=Model.intercept_[0] print(f"截距：{a}") b=Model.coef_[0][0] print(f"斜率:{b}")

總結(jié)

以上是生活随笔為你收集整理的预测模型算法的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。