(已解决)时间序列画趋势线的问题

LeoBao

最近做论文，画完时间序列图之后，觉得看起来不是很明显，想画趋势线，手册中查到用trege(expr1，expr2,tdim1,tdim2)这个函数可以实现，于是把年份1959-2008放到一个grd里面，并编了ctl,作为变量1，然后原来50年平均降水50个值作为变量2并编写gs如下
'reinit'
'open d:\essay\qs\ttt.ctl'
'open d:\essay\qs\tt.ctl'
'set x 1'
'set y 1'
'set z 1'
'set grid on'
'set grads off'
'enable print d:\essay\qs\rain111.gmf'

'set t 1 50'
'd ot.2'

'set t 1'
'd tregr(st,ot.2,t=1,t=50)'
'define coeff=tregr(st,ot.2,t=1,t=50)'
'define stave=ave(st,t=1,t=50)'
'define otave=ave(ot.2,t=1,t=50)'
'd coeff*(st-stave)+otave'
'print'
'c'
'disable'
;
奇怪的是最后只出了折线图，最后几行编的没用，没有趋势图，为什么啊？

lqouc

shrily521 发表于 2014-6-12 09:33
这个我试了下，只能画直线，不能画曲线，有时候画出来的趋势线并不合适。要是能画二次拟合曲线就好了，不知 ...

单纯二次拟合也不一定能画出很好的效果，可以考虑多项式拟合，当然这个在grads里面不好实现求解线性方程组。可以用fortran计算后在画图。

LeoBao

问题解决了，在'd coeff*(st-stave)+otave'上少了‘set t 1 50’

和清谈

我调了一下，感觉后面几句没什么用，不知道哪里出错了，求指教

易小凯

我每次都会用excel的添加趋势线来完成
没想到grads也有这种功能

LeoBao

和清谈 发表于 2012-5-10 20:16

                               
登录/注册后可看大图

我调了一下，感觉后面几句没什么用，不知道哪里出错了，求指教

后面几句是求线性趋势的，就是一元回归方程

LeoBao

和清谈 发表于 2012-5-10 20:16

                               
登录/注册后可看大图

我调了一下，感觉后面几句没什么用，不知道哪里出错了，求指教

后面几句是求线性趋势的，就是一元回归方程

和清谈

这个我知道，不过还是纠结了好久，最后不知道改了什么就出来了，谢谢哦

nuister

搞不定哎，st指的是年份吗？

nuister

ttt.ctl 和tt.ctl分别描述的grd文件内容是什么呀？谢谢

LeoBao

nuister 发表于 2012-5-12 19:57

                               
登录/注册后可看大图

ttt.ctl 和tt.ctl分别描述的grd文件内容是什么呀？谢谢

tt是要画的参数，ttt是对应年份（时间），gs中还有个小错，是d ot.2