求最优气候值方法的实现代码

cuijxin

大家好，刚接触气象这方面的知识，想做一个最优气候值方法的程序，用来求干旱或者降雨的预测值，请大家给点建议，思路

topmad

最优气候值方法? 怎么讲

cuijxin

“从1992年起，美国国家气象中心／气候分析中心建立了一项新的业务。即每年四次发布试验性的长超前预测公报（Experimental long-lead Forecast Bulletin)，汇集若干不同来源的各种长期预报与方法，以便相互交流，促进长期预报工作的发展。公报中的方法之一，就是用“最优”气候值做地面气温的季度预报12。这种方法的中心思想是认为气候变化具有某些内在规律，如阶段性和周期性，因此，简单地用多年(如30年)周期的气候平均值作为来年该要素的估计值的方法有明显的不足，持续性预报也是如此。因此，应力求从历史资料中去获取最优气候值，为统计预报方法的建立提供依据。本文借鉴此方法进行我国短期气候预测试验。”

cuijxin

“二．方法
    我们利用全国160站月(或季)平均温度(或降水)资料，分别依次取前1、2、3、……、20年的平均做为来年的预报值，如以1951～1970年的实况资料为基础，制作71年的预报，分别取70年、69+70、68+69+70、……等20个平均值而得到20个预报值，以次类推；以预测值与实况最接近为标准，从统计上得出每年“最优”的平均年数，然后，取“最优”平均年数频率最高的年数作为计算单站“最优”气候平均值应取的年数，并依次制作下一年的预报。从1971～1996年的26个个例中，挑选出最优气候值，可制作出1997年的预报。

1.对i年可有20个预测值:

前一年的平均  y1(i)=s(i-1)
前二年的平均  y2(i)=(s(i-1)+s(i-2))/2
前三年的平均  y3(i)=(s(i-1)+s(i-2)+s(i-3))/3
前四年的平均  y4(i)=(s(i-1)+s(i-2)+s(i-3)+s(i-4))/4
…
…
前二十年的平均  y20(i)=(s(i-1)+……+s(i-20))/20

2.对这20个预测值与实况比较,取最小值对应的年数n:

Z(n)=min(y1(i)-s(i),y2(i)-s(i),……,y20(i)-s(i))

3.从1971～1996年26个n中，挑选出频次最高的N

4.用前N年的平均来制作1997年的预报:

Y(i)=(s(i-1)+……+s(i-N))/N”

言深深

cuijxin 发表于 2012-3-15 16:25

                               
登录/注册后可看大图

“二．方法
    我们利用全国160站月(或季)平均温度(或降水)资料，分别依次取前1、2、3、……、20年的平均 ...

楼主真有耐心啊！
呵呵，刚刚看完了您写了几点内容，稍微对您提点儿建议，不合适的话还请见谅。
1，您所说的“最优气候值”仅仅是您阅读的一篇文献中的个例，还不能作为整个行业的一个标准，因此没有专门的程序。（但是想诸如EOF，小波分析，或者简单一些的方差，标准差，滑动平均等方式是有固定的方法，也是有相关程序的）；

2，如果您打算复制您目前所阅读的文献（也就是实现文献中的做法，从理解文章内容，到完完整整的计算出来，这也是很多新手熟悉一个新的行业，或者领域的入门必修），我觉得您不妨自己花点儿功夫和事件，研究一下方法的具体实施和操作过程，阅读您所拥有文献的参考资料，弄清楚每一个疑点的来源；

3，学习一些简单的程序知识，在第二点的基础上尝试自己编写程序。因为别人写的一个稍微大点儿的程序，理解起来是非常困难的，尤其对于基础不是很扎实的来说，更加难以理解；

4，最后，我想说的是，在您实现以上三个过程中有一些疑问，可以量化为一些小的问题，具体拿出来问。不然，你这个问题贴的题目太大，要是有谁能够回答完全的话，应该能够写一篇核心期刊的论文了。

祝您好运，说的不妥的地方，望请海涵

mofangbao

楼主是给气象局做外包开发的？

mofangbao

深深说的没错，没有气象基础知识这些东西不是这么快就能掌握的，想做这种程序很吃力。所以楼主如果是接的别人的活儿那可以找一个懂气象的和你一起合作，这样会省力不少。
如果楼主是自己要开始学习气象知识的话请原谅我的猜测。
光看上面的文献，简单的理解就是从后往前依次做平均值，找相差最小的那个作为最优，算法都给你写出来了

cuijxin

言深深 发表于 2012-3-15 17:10

                               
登录/注册后可看大图

楼主真有耐心啊！
呵呵，刚刚看完了您写了几点内容，稍微对您提点儿建议，不合适的话还请见谅。
1，您所 ...

谢谢您的建议 说的很好 我很受用

看了几天文档了 也开始动笔写了些程序，从观测文件里读取了数据 计算了平均值 现在之所以做不下去了 可能

是有些需求还不明确 我这么问 确实有点不妥了 不过还是谢谢您的建议

cuijxin

言深深 发表于 2012-3-15 17:10

                               
登录/注册后可看大图

楼主真有耐心啊！
呵呵，刚刚看完了您写了几点内容，稍微对您提点儿建议，不合适的话还请见谅。
1，您所 ...

是的 我们接了个项目 需要了解的东西还很多啊

言深深

cuijxin 发表于 2012-3-15 18:32

                               
登录/注册后可看大图

是的 我们接了个项目 需要了解的东西还很多啊

按算法走吧，还是可以实现的