如何计算线性趋势

向晚_clouds

把西北太平洋分为2°*2°的格点后统计了每个格点中35年来的台风的頻数，现在想看台风頻数在哪些区域有增加趋势，哪些区域有减少趋势，要怎么计算每个格点的线性趋势？

言深深

趋势计算比较简单，我写了一个去趋势的子程序，里面有趋势项的求解，可以使用：


subroutine dtrend(x,n,y)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!        【程序说明】                                                                !
!        此程序段用最小二乘法去掉线性趋势                        !
!                                                                                                !
!        对输入数组x(i)，默认横坐标是自然数；                !
!                输出数组y(i)是去掉线性趋势之后的数组。        !
!        使用公式：y=a*x+b                                                        !
!                a=(nΣxy-ΣxΣy)/(nΣx^2-(Σx)^2)                !
!                b=y(平均)-a*x（平均）                                        !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
real x(n),y(n)
        s1=0;s2=0;s3=0;s4=0
        do i=1,n
        s1=s1+i*x(i)
        s2=s2+i
        s3=s3+x(i)
        s4=s4+i**2
        enddo
                a=(s1*n-s2*s3)/(s4*n-s2**2)
                b=(s3-a*s2)/n
        do i=1,n
        y(i)=x(i)-i*a-b
        enddo
end

其中a就是趋势！

river

言深深 发表于 2015-3-30 15:39
趋势计算比较简单，我写了一个去趋势的子程序，里面有趋势项的求解，可以使用：

你又出现了，感觉你好久没活动了

言深深

river 发表于 2015-3-30 16:35
你又出现了，感觉你好久没活动了

额，一直都在，并未远去。

river

言深深 发表于 2015-3-30 16:37
额，一直都在，并未远去。

看来现在很多人都是这个状态···

向晚_clouds

言深深 发表于 2015-3-30 15:39
趋势计算比较简单，我写了一个去趋势的子程序，里面有趋势项的求解，可以使用：

谢谢~学渣我研究一下。。。

小赵4403

想问下能在什么软件中运行啊，MATLAB中能使用吗

lqouc

小赵4403 发表于 2015-3-30 22:16
想问下能在什么软件中运行啊，MATLAB中能使用吗

这是fortran语言，matlab有现成函数。

小赵4403

lqouc 发表于 2015-3-31 09:54
这是fortran语言，matlab有现成函数。

是什么函数怎么用啊

小赵4403

lqouc 发表于 2015-3-31 09:54
这是fortran语言，matlab有现成函数。

是什么函数怎么用啊