请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 22519|回复: 27

[求助] 如何计算线性趋势

[复制链接]

新浪微博达人勋

发表于 2015-3-30 12:45:20 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
把西北太平洋分为2°*2°的格点后统计了每个格点中35年来的台风的頻数,现在想看台风頻数在哪些区域有增加趋势,哪些区域有减少趋势,要怎么计算每个格点的线性趋势?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2015-3-30 15:39: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就是趋势!
密码修改失败请联系微信:mofangbao
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-3-30 16:35:48 | 显示全部楼层
言深深 发表于 2015-3-30 15:39
趋势计算比较简单,我写了一个去趋势的子程序,里面有趋势项的求解,可以使用:

你又出现了,感觉你好久没活动了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 0
发表于 2015-3-30 16:37:18 | 显示全部楼层
river 发表于 2015-3-30 16:35
你又出现了,感觉你好久没活动了

额,一直都在,并未远去。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-30 17:08:18 | 显示全部楼层
言深深 发表于 2015-3-30 16:37
额,一直都在,并未远去。

看来现在很多人都是这个状态···
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-3-30 18:52:12 | 显示全部楼层
言深深 发表于 2015-3-30 15:39
趋势计算比较简单,我写了一个去趋势的子程序,里面有趋势项的求解,可以使用:

谢谢~学渣我研究一下。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-30 22:16:52 | 显示全部楼层
想问下能在什么软件中运行啊,MATLAB中能使用吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-31 09:54:39 | 显示全部楼层
小赵4403 发表于 2015-3-30 22:16
想问下能在什么软件中运行啊,MATLAB中能使用吗

这是fortran语言,matlab有现成函数。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-31 12:01:51 | 显示全部楼层
lqouc 发表于 2015-3-31 09:54
这是fortran语言,matlab有现成函数。

是什么函数怎么用啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-31 12:02:04 | 显示全部楼层
lqouc 发表于 2015-3-31 09:54
这是fortran语言,matlab有现成函数。

是什么函数怎么用啊
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表