fortran 求趋势系数

believefuture

river 发表于 2014-8-31 21:55
首先还有一个赋初值的问题，d(i)=0.0为什么要放在别的循环里啊，直接放在
do i=1,662
d(i)=p(i)/sqrt(q ...

就是看的这个公式，我就是按这个编的，我的原始数据就是662个站点50年的数据

believefuture

言深深 发表于 2014-8-31 20:44
p(i)=p(i)+(a(i,j)-m(i))*(j-25.5)
q(i)=q(i)+(a(i,j)-m(i))*(a(i,j)-m(i))*(j-25.5)*(j-25.5)

就是用的这个趋势系数的公式

cx6689942

赞一个，希望继续为大家提供自己的体会

言深深

believefuture 发表于 2014-8-31 23:31
就是用的这个趋势系数的公式

你再看看和你的一样么？无语了···
下面的是x序列和t序列的方差，即两条序列分别的距平平方累计乘积的开放，为何到你这边就变成了距平平方乘积的累积了呢？这差别很大的好么？

believefuture

言深深 发表于 2014-9-1 08:26
你再看看和你的一样么？无语了···
下面的是x序列和t序列的方差，即两条序列分别的距平平方累计乘积的 ...

谢谢你，自己太笨了

言深深

believefuture 发表于 2014-9-1 08:40
谢谢你，自己太笨了

这是我用最小二乘法计算的趋势，其实就是斜率。我不知道和你的趋势系数是不是一个东西，不过话说回来，你给的公式其实是相关系数不知道楼主到底需要的是什么呢？
贴上我自己编写的一段计算序列x的斜率的子程序。
subroutine trend(x,n,a)
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!        【程序说明】                                                                !
!        此程序段用最小二乘法求解线性趋势                        !
!        使用公式：y=a*x+b                                                        !
!                a=(nΣxy-ΣxΣy)/(nΣx^2-(Σx)^2)                !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
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)
end

river

believefuture 发表于 2014-8-31 23:30
就是看的这个公式，我就是按这个编的，我的原始数据就是662个站点50年的数据

你在哪看到的这个公式？

believefuture

言深深 发表于 2014-9-1 08:46
这是我用最小二乘法计算的趋势，其实就是斜率。我不知道和你的趋势系数是不是一个东西，不过话说回来，你给 ...

万分感谢你，我其实就是求变量随时间的趋势系数，
我求的趋势系数：样本长度为n的要素的时间序列与自然数1,2,3，。。。，n之间的相关系数

言深深

believefuture 发表于 2014-9-1 11:16
万分感谢你，我其实就是求变量随时间的趋势系数，
我求的趋势系数：样本长度为n的要素的时间序列与自然 ...

好