1--灰色系统理论建模学案.pptVIP

即或 （1） 在式（1）中， 称为灰导数，a称为发展系数， 称为白化背景值，b称为灰作用量。 将时刻表 代入（1）式有 引入矩阵向量记号： 于是GM（1，1）模型可表示为 现在问题归结为求a,b在值。用一元线性回归，即最小二乘法求它们的估计值 为 注：实际上回归分析中求估计值是用软件计算的，有标准程序求解，如matlab 等。 GM（1，1）的白化型 对于GM（1，1）的灰微分方程（1），如果将灰导数 的时刻 视为连续变量t，则 视为时间t函数 ，于是 对应于导数量级 ，白化背景值 对应于导数 。于是GM（1，1）的灰微 分方程对应于的白微分方程为 （2） 1.数据的检验与处理 为了保证GM（1，1）建模方法的可行性，需要对已知数据做必要的检验 处理。 设原始数据列为了 ，计算数列的级比 如果所有的级比都落在可容覆盖区间 内，则数据列 可以建立GM（1，1）模型且可以进行灰 色预测。否则，对数据做适当的变换处理，如平移变换： 取C使得数据列 的级比都落在可容覆盖内。 GM（1，1）灰色预测的步骤 2. 建立GM（1，1）模型 不妨设 满足上面的要 求，以它为数据列建立GM（1，1）模型 用回归分析求得a,b的估计值，于是相应的白化模型为 解为 （3） 于是得到预测值 从而相应地得到预测值： 3. 检验预测值 （1）残差检验：计算相对残差 如果对所有的 ，则认为达到较高的要求：否则，若 对所有的 ，则认为达到一般要求。 （2）级比偏差值检验：计算 如果对所有的 ，则认为达到较高的要求；否则 若对所有的 ，则认为达到一般要求。 例 北方某城市1986～1992 年道路交通噪声平均声级数据见表 灰色预测计算实例 1 1986 71.1 2 1987 72.4 3 1988 72.4 4 1989 72.1 5 1990 71.4 6 1991 72.0 7 1992 71.6 第一步: 级比检验 建立交通噪声平均声级数据时间序列如下： =(71.1, 72.4, 72.4, 72.1, 71.4, 72.0, 71.6) 序号 年份 Leq （1）求级比λ (k) =(0.982,1,1.0042,1.0098, 0.9917,1.0056) （2）级比判断 由于所有的λ (k)∈[0.982,1.0098]，k = 2,3,……,7，故可以用x(0)作满意的GM（1，1）建模。 第二步: GM（1，1）建模 （1）对原始数据 作一次累加，即 = (71.1,143.5,215.9,288,359.4,431.4,503） （2）构造数据矩阵B 及数据向量Y （3）计算 ?? 于是得到a =