计量地理学6趋势面分析方法.pptx

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第六章 趋势面分析方法,趋势面分析普通原理,趋势面模型适度检验,趋势面分析应用实例,1/39,先看一个例子,流域降水量及观察点地理位置数据,序号,降水量,Z/,mm,横坐标,x,/,10,4,m,纵坐标,y,/,10,4,m,1,2,3,4,5,6,7,8,9,10,11,12,27.6,38.4,24,24.7,32,55.5,40.4,37.5,31,31.7,53,44.9,0,1.1,1.8,2.95,3.4,1.8,0.7,0.2,0.85,1.65,2.65,3.65,1,0.6,0,0,0.2,1.7,1.3,2,3.35,3.15,3.1,2.55,2/39,a,Y,Z,3/39,X,Y,Z,a,4/39,某流域降水量二次多项式趋势面,5/39,某流域降水量三次多项式趋势面,6/39,7/39,第九章 趋势面分析方法,趋势面分析普通原理,趋势面模型适度检验,趋势面分析应用实例,8/39,一、趋势面分析普通原理,趋势面分析，是利用数学曲面模拟地理系统要素在空间上分布及改变趋势一个数学方法。,它实质上是经过回归分析原理，利用最小二乘法拟合一个二维非线性函数，模拟地理要素在空间上分布规律，展示地理要素在地域空间上改变趋势。,9/39,趋势面分析方法经常被用来模拟资源、环境、人口及经济要素在空间上分布规律，它在空间分析方面含有主要应用价值。,趋势面是一个抽象数学曲面，它抽象并过滤掉了一些局域随机原因影响，使地理要素空间分布规律显著化。,10/39,通常把实际地理曲面分解为趋势面和剩下面两部分，前者反应地理要素宏观分布规律，属于确定性原因作用结果；而后者则对应于微观局域，是随机原因影响结果。,趋势面分析一个基本要求，就是所选择趋势面模型应该是剩下值最小，而趋势值最大，这么拟合度精度才能到达足够准确性。空间趋势面分析，正是从地理要素分布实际数据中分解出趋势值和剩下值，从而揭示地理要素空间分布趋势与规律。,11/39,（一）,建立趋势面模型,设某地理要素实际观察数据为,，趋势面拟合值为,，则有,式中：,i,即为剩下值（残差值）。,显然，当（,x,i,，y,i,）在空间上变动时，（3.6.1）式就刻画了地理要素实际分 布曲面、趋势面和剩下面之间互动关系。,（3.6.1）,12/39,趋势面分析关键,从实际观察值出发推算趋势面，普通采取回归分析方法，使得残差平方和趋于最小，即,这就是在最小二乘法意义下趋势面拟合。,13/39,用来计算趋势面数学方程式有多项式函数和傅立叶级数，其中最为惯用是多项式函数形式。因为任何一个函数都能够在一个适当范围内用多项式来迫近，而且调整多项式次数，可使所求回归方程适合实际问题需要。,14/39,多项式趋势面形式,一次趋势面模型,二次趋势面模型三次趋势面模型,（3.6.2）,（3.6.3）,（3.6.4）,15/39,（二）预计趋势面模型参数,实质,依据观察值,z,i,，,x,i,，,y,i,（,i,=1,2,n,）确定多项式系数,a,0,，,a,1,，,a,p,，使残差平方和最小。,过程,将多项式回归（非线性模型）模型转化为多元线性回归模型。,16/39,令,则,其残差平方和为,（3.6.5）,17/39,求,Q,对,a,0,，,a,1,，,a,p,偏导数，并令其等于0，得正规方程组(式中 为,p,+1个未知量）,（3.6.6）,18/39,用矩阵形式表示,则,（3.6.6）,式变为,（3.6.7）,19/39,对于二元二次多项式有,其正规方程组为,20/39,21/39,由式（3.6.7）求解，可得,（3.6.8）,22/39,二、趋势面模型适度检验,趋势面拟适当度,R,2,检验,趋势面拟适当度显著性,F,检验,趋势面适度逐次检验,趋势面分析拟合程度与回归模型效果直接相关，所以，对趋势面分析进行适度性检验是一个关系到趋势面能否在实际研究中加以应用关键问题，也是趋势面分析中不可缺乏主要步骤。,这能够经过以下检验来完成：,23/39,（一）,趋势面拟适当度,R,2,检验,趋势面与实际面拟合度系数,R,2,是测定回归模型拟合优度主要指标。,普通用变量,z,总离差平方和中回归平方和所占比重表示回归模型拟合优度。,总离差平方和等于回归平方和与剩下平方和之和。即,24/39,越大（或,越小）就表示因变量与自变量关系越亲密，回归规律性越强、效果越好。,记,越大，趋势面拟合度就越高。,为剩下平方和，,它表示随机原因对离差影响，,为回归平方和，,它表示自变量对因变量离差总影响。,（3.6.9）,25/39,（二）趋势面拟适当度显著性,F,检验,趋势面适度,F,检验，是对趋势面回归模型整体显著性检验。,方法：是利用变量,z,总离差平方和中剩下平方和与回归平方和比值，确定变量,z,与自变量,x,、,y,之间回归关系是否显著。即,结果分析：在显著性水平,下，查,F,分布表得,F,，若计算,F,值大于临界值,F,，则认为趋势面方程显著；反之则不显著。,（3.6.10）,26/39,（三）趋势面适度逐次检验,方法,(1)求出较高次多项式方程回归平方和与较低次多项式方程回归平方和之差;,(2)将此差除以回归平方和自由度之差，得出因为多项式次数增高所产生回归均方差;,(3)将此均方差除以较高次多项式剩下均方差，得出相继两个阶次趋势面模型适度性比较检验值,F,。,若所得,F,值是显著，则较高次多项式对回归作出了新贡献，若,F,值不显著，则较高次多项式对于回归并无新贡献。对应方差分析表见表3.6.1。,27/39,表3.6.1 多项式趋势面由,K,次增高至（,K,+1）次 回归显著性检验,离差起源,平方和,自由度,均方差,F,检验,（,K,+1,）次回归,p,（K+1）次剩下,n,p,1,K,次回归,q,K次剩下,n,q,1,由K次增高至,（K+1）次回归,p,q,总离差,SS,T,28/39,需要注意是，在实际应用中，往往用次数低趋势面迫近改变比较小地理要素数据，用次数高趋势面迫近起伏改变比较复杂地理要素数据。次数低趋势面使用起来比较方便，但详细到某点拟合较差；次数较高趋势面只在观察点附近效果很好，而在外推和内插时则效果较差。,29/39,三、趋势面分析应用实例,某流域1月份降水量与各观察点坐标位置数据如表3.6.2所表示。下面，我们以降水量为因变量,z,，地理位置横坐标和纵坐标分别为自变量,x,、,y,，进行趋势面分析，并对趋势面方程进行适度,F,检验。,30/39,表3.6.2 流域降水量及观察点地理位置数据,序号,降水量,Z/,mm,横坐标,x,/,10,4,m,纵坐标,y,/,10,4,m,1,2,3,4,5,6,7,8,9,10,11,12,27.6,38.4,24,24.7,32,55.5,40.4,37.5,31,31.7,53,44.9,0,1.1,1.8,2.95,3.4,1.8,0.7,0.2,0.85,1.65,2.65,3.65,1,0.6,0,0,0.2,1.7,1.3,2,3.35,3.15,3.1,2.55,31/39,建立趋势面模型,(1),首先采取二次多项式进行趋势面拟合，用最小二乘法求得拟合方程为,解题步骤,32/39,图3.6.1 某流域降水量二次多项式趋势面,33/39,(2),再采取三次趋势面进行拟合，用最小二乘法求得拟合方程为,34/39,图3.6.2 某流域降水量三次多项式趋势面,35/39,(1)趋势面拟适当度,R,2,检验:依据,R,2,检验方法计算，结果表明，二次趋势面判定系数为,R,2,2,=0.839，三次趋势面判定系数为,R,3,2,=0.965，可见二次趋势面回归模型和三次趋势面回归模型显著性都较高，而且三次趋势面较二次趋势面含有更高拟合程度。,模型检验,36/39,(2),趋势面适度显著性,F,检验:依据,F,检验方法计算，结果表明，二次趋势面和三次趋势面,F,值分别为,F,2,=6.236和,F,3,=6.054。在置信水平,=0.05下，查,F,分布表得,F,3,=。,显然 ，,故二次趋势面回归方程显著而三次趋势面不显著。所以，,F,检验结果表明，用二次趋势面进行拟合比较合理。,。,37/39,趋势面比较：在二次和三次趋势面检验中，对两个阶次趋势面模型适度进行比较，对应方差分析计算结果见表3.6.3。,(3)趋势面适度逐次检验：,离差起源,平方和,自由度,均方差,F,检验,三次回归,三次剩下,二次回归,二次剩下,由二次增高至,三次回归,1 129.789,41.474,982.244,189.018,147.545,9,12,-,9,-,1,5,12,-,5,-,1,4,125.532,20.737,196.449,31.503,36.886,6.054,6.236,1.779,表3.6.3 二次和三次趋势面回归模型逐次检验方差分析表,38/39,分析：从二次趋势面增加到三次趋势面，,F,32,=1.779。在置信度水平,=0.05下，查,F,分布表得,F,0.05,(4,2)=6.94，因为,F,32,F,0.05,(4,2)=6.94，故将趋势面拟合次数由二次增高至三次，对回归方程并无新贡献，因而选取二次趋势面比较适当。这也深入验证了趋势面拟适当度显著性,F,检验结论。,39/39,