灰色预测方法.pptVIP

第五步：求关联度计算结果表明，运输业和工业的关联程度大于农业、商业和工业的关联程度。为参考序列时，计算类似，这里略去。第31页，共72页，星期日，2025年，2月5日aita2=1.00000.50380.37050.3333aita3=1.00000.83900.52540.5035aita4=1.00000.63500.49730.3542loadgldr=0.5aita2=(u+r*v)./(y2+r*v)aita3=(u+r*v)./(y3+r*v)aita4=(u+r*v)./(y4+r*v)r12=(1/4)*sum(aita2)r13=(1/4)*sum(aita3)r14=(1/4)*sum(aita4)r12=0.5519r13=0.7170r14=0.6216第32页，共72页，星期日，2025年，2月5日8.3灰色预测模型一、GM（1，1）模型的建立设时间序列有n个观察值：通过累加生成新序列令为的紧邻均值生成序列第33页，共72页，星期日，2025年，2月5日则GM（1，1）的灰微分方程模型为：其中：α称为发展灰数；μ称为内生控制灰数。设为待估参数向量，可利用最小二乘法求解。解得：其中第34页，共72页，星期日，2025年，2月5日第35页，共72页，星期日，2025年，2月5日称微分方程：为灰色微分方程的白化方程，也称影子方程。（1）白化方程的解也称时间响应函数，其为（2）GM（1，1）灰微分方程的时间响应序列为第36页，共72页，星期日，2025年，2月5日（3）取（4）还原值注：（1）原始数据不一定全部用来建模，数据不同，模型不同。（2）数据的取舍应保证建模序列等时矩、相连，不得有跳跃出现。（3）一般建模数据应当有最新的数据及其相邻数据构成。第37页，共72页，星期日，2025年，2月5日二、模型检验（1）残差检验按预测模型计算并将累减生成然后计算原始序列与的绝对误差序列及相对误差序列。灰色预测检验一般有残差检验、关联度检验和后验差检验。给定a,当成立,称模型为参差合格模型.a取0.01、0.05、0.1所对应的模型分别为：优、合格、勉强合格第38页，共72页，星期日，2025年，2月5日（2）关联度检验根据前面所述关联度的计算方法算出与原始序列的关联系数，然后计算出关联度，根据经验，当ρ=0.5时，关联度大于0.6便满意了。第39页，共72页，星期日，2025年，2月5日（3）后验差检验a.计算原始序列标准差:b.计算绝对误差序列的标准差：c.计算方差比：第40页，共72页，星期日，2025年，2月5日d.计算小误差概率：令:，则P0.950.800.70≤0.70C0.350.500.65≥0.65好合格勉强合格不合格第41页，共72页，星期日，2025年，2月5日例某矿某年3-7月份的轻伤事故情况如表所示:原始数据列为：累加生成数列为：月份34567轻伤人次2629313334表1轻伤事故人次第42页，共72页，星期日，2025年，2月5日第43页，共72页，星期日，2025年，2月5日所以其中α称为发展灰数；μ称为内生控制灰数。第44页，共72页，星期日，2025年，2月5日所以，即事故预测公式为：第45页，共72页，星期日，2025年，2月5日生成数列的预测值、原始数列的还原值分别如表所示。为了得到原始数列的预测值，需要将生成数列的预测值作累减还原为原始值，即根据下式求得：表2生成数列的预测值与误差检验-0.03153.0315340.30118.701193-0.1486.14862-0.2755.27551026260k第46页，共72页，星期日，2025年，2月5日表3原始数