趋势面分析.ppt

1 趋势面分析方法 第1节趋势面分析的一般原理第2节趋势面模型的适度检验第3节趋势面分析应用实例 2 趋势面分析是利用数学曲面模拟地理系统要素在空间上的分布及变化趋势的一种数学方法 第1节趋势面分析的一般原理 趋势面分析的实质 通过回归分析原理 运用最小二乘法拟合一个二维非线性函数 模拟地理要素在空间上的分布规律 展示地理要素在地域空间上的变化趋势 3 趋势面分析方法常被用来模拟资源 环境 人口及经济要素在空间上的分布规律 它在空间分析方面具有重要的应用价值 举例 粮食产量与气温和降雨量等自然因素的关系是非线性关系 可以利用趋势分析来拟合回归模型从而预测粮食产量 4 趋势面是一种抽象的数学曲面 它抽象并过滤掉一些局域随机因素的影响 使地理要素的空间分布规律明显化 所谓的空间趋势面并不是地理要素的实际分布面 而是一个模拟地理要素空间分布的近似曲面 5 变量的实测数据Mi xi yi zi 分布在趋势面上或趋势面上下 如图所示 实测点在趋势面上的投影 实测点 趋势面 趋势面示意图 6 实际曲面 趋势面 剩余曲面 趋势面反映地理要素的宏观分布规律 属于确定性因素作用的结果 剩余面反映局部性变化特点 它受局部因素和随机因素的影响 局部异常 随机干扰和模型本身的误差 7 观测值 确定性函数值 随机性函数值 趋势值 剩余值我们通常用回归分析的方法来求趋势值和剩余值 趋势面分析的基本要求 使剩余值最小 趋势值最大 这样拟合度精度才能达到足够的准确性 空间趋势面分析 正是从地理要素分布的实际数据中分解出趋势值和剩余值 从而揭示地理要素空间分布的趋势与规律 8 1 1建立趋势面模型 设某地理要素的实际观测数据为趋势面拟合值为则有 3 6 1 式中 i即为剩余值 残差值 显然 当 xi yi 在空间上变动时 3 6 1 式就刻画了地理要素的实际分布曲面 趋势面和剩余面之间的互动关系 9 趋势面分析的核心 从实际观测值出发推算趋势面 采用回归分析方法 使得残差平方和趋于最小 即 以此来估计趋势面参数 这就是在最小二乘法意义下的趋势曲面拟合 10 用来计算趋势面的数学方程式有多项式函数和傅立叶级数 其中最为常用的是多项式函数形式 因为任何一个函数都可以在一个适当的范围内用多项式来逼近 而且调整多项式的次数 可使所求的回归方程适合实际问题的需要 11 多项式趋势面的形式 一次趋势面模型 二次趋势面模型 三次趋势面模型 3 6 2 3 6 3 3 6 4 12 在实际的空间趋势面模拟中 按照对事物认识由易到难的规律 应首先考虑用公式 3 6 2 表达的倾斜平面去拟合 然后再用公式 3 6 3 描述的二次抛物趋势面去模拟 如果还不能满足研究需求 则需选用三次趋势面 四次趋势面甚至更高次趋势面进行拟合 13 实质 根据观测值zi xi yi i 1 2 n 确定多项式的系数a0 a1 ap 使残差平方和最小 过程 将多项式回归 非线性模型 模型转化为多元线性回归模型 令则 1 2估计趋势面模型的参数 14 求Q对a0 a1 ap的偏导数 并令其等于0 得正规方程组 式中a0 a1 ap为p 1个未知量 其残差平方和为 3 6 5 3 6 6 15 用矩阵形式表示 则 3 6 6 式变为 3 6 7 对于二元二次多项式有其正规方程组为 16 由式 3 6 7 求解 可得 3 6 8 17 2 1趋势面拟合适度的R2检验2 2趋势面拟合适度的显著性F检验2 3趋势面适度的逐次检验 第2节趋势面模型的适度检验 18 趋势面与实际面的拟合度系数R2是测定回归模型拟合优度的重要指标 一般用变量z的总离差平方和中回归平方和所占的比重表示回归模型的拟合优度 总离差平方和等于回归平方和与剩余平方和之和 即 2 1趋势面拟合适度的R2检验 19 为剩余平方和 它表示随机因素对Z的离差的影响 为回归平方和 它表示p个自变量对因变量z的离差的总影响 越大 或越小 就表示因变量与自变量的关系越密切 回归的规律性越强 效果越好 记越大 趋势面的拟合度就越高 3 6 9 20 趋势面适度的F检验是对趋势面回归模型整体的显著性检验 方法 利用变量z的总离差平方和中剩余平方和与回归平方和的比值 确定变量z与自变量x y之间的回归关系是否显著 即结果分析 在显著性水平 下 查F分布表得F 若计算的F值大于临界值F 则认为趋势面方程显著 反之则不显著 3 6 10 2 2趋势面拟合适度的显著性F检验趋势 21 步骤 1 求出较高次多项式方程的回归平方和与较低次多项式方程的回归平方和之差 2 将此差除以回归平方和的自由度之差 得出由于多项式次数增高所产生的回归均方差 3 将此均方差除以较高次多项式的剩余均方差 得出相继两个阶次趋势面模型的适度性比较检验值F 若所得的F值是显著的 则较高次多项式对回归作出了新贡献 若F值不显著 则较高次多项式对于回归并无新贡献 相应的方差分析表见表9 1 2 3趋势面适度的逐次检验 22 表9 1多项式趋势面由K次增高至 K 1 次的回归显著性检验 23 注意 在实际应用中 往往用次数低的趋势面逼近变化比较小的地理要素数据 用次数高的趋势面逼近起伏变化比较复杂的地理要素数据 次数低的趋势面使用起来比较方便 但具体到某点拟合较差 次数较高的趋势面只在观测点附近效果较好 而在外推和内插时效果较差 24 实例1 某流域1月份降水量与各观测点的坐标位置数据如表3 6 2所示 下面 我们以降水量为因变量z 地理位置的横坐标和纵坐标分别为自变量x y进行趋势面分析 并对趋势面方程进行适度F检验 第3节趋势面分析应用实例 25 表9 2流域降水量及观测点的地理位置数据 26 1 建立趋势面模型 解题步骤 二次趋势面模型 采用二次多项式进行趋势面拟合 用最小二乘法求得拟合方程为 27 求取拟合方程的具体过程 将相关数据输入Excel软件中 令x1 x x2 y x3 x2 x4 xy x5 y2 并根据原始数据进行运算 结果如下 28 将前述准备好的数据复制到SPSS软件中 如下图 29 按图示进行操作 利用SPSS软件求取方程系数 30 31 运算结果 32 33 图3 6 1某流域降水量的二次多项式趋势面 34 三次趋势面模型 按照同样方法采用三次趋势面进行拟合 用最小二乘法求得拟合方程为 35 基于SPSS的运算结果 36 37 图3 6 2某流域降水量的三次多项式趋势面 38 趋势面拟合适度的R2检验根据R2检验方法计算 结果表明 二次趋势面的判定系数R22 0 839 三次趋势面的判定系数R32 0 965 可见二次趋势面回归模型和三次趋势面回归模型的显著性都较高 而且三次趋势面较二次趋势面具有更高的拟合程度 2 模型检验 39 趋势面适度的显著性F检验根据F检验方法计算 结果表明 二次趋势面和三次趋势面的F值分别为F2 6 236和F3 6 054 在置信水平 0 05下 查F分布表得F2 F0 05 5 6 4 39 F3 F0 05 9 2 19 38 显然 F2 F2 而F3 F3 故二次趋势面的回归方程显著而三次趋势面不显著 因此 F检验的结果表明 用二次趋势面进行拟合比较合理 40 趋势面比较 在二次和三次趋势面检验中 对两个阶次趋势面模型的适度进行比较 相应的方差分析计算结果见表3 6 3 趋势面适度的逐次检验 表3 6 3二次和三次趋势面回归模型的逐次检验方差分析表 41 42 结果分析 从二次趋势面增加到三次趋势面 F3 2 1 779 在置信度水平 0 05下 查F分布表得F0 05 4 2 19 25 由于F3 2 F0 05 4 2 19 25 故将趋势面拟合次数由二次增高至三次 对回归方程并无新贡献 因而选取二次趋势面比较合适 这也进一步验证了趋势面拟合适度的显著性F检验的结论 43 利用Matlab软件实现制图 绘制上述二次趋势面图形 可以直接调用如下函数命令 x y meshgrid 0 0 25 4 z 5 988 17 438 x 29 787 y 3 558 x 2 0 357 x y 8 070 y 2 surface x y z 具体过程如后图所示 44 45 46 47 48 绘制上述三次趋