矢量风场合成的F检验

0708011

做合成场时一般用t检验，但是矢量风场不可以用这种，参照气象家园上这个帖子 使用F检验 ftest ，自己编写了ncl版本的~~~
原帖http://bbs.06climate.com/forum.php?mod=viewthread&tid=19518
感谢原作者麦田_smile好分享。
公式见附件的文献。
应该是对的。这是原帖程序用grads画出来的出的F值，注意最小值是10开始

这是我用ncl画出来的，这里竖线区最小值是4.7，所以区域大一些。

zongqing123

naq=5   ;强年的数目
nar=5   ;弱年的数目
m=2
ucha=0.  ;u强弱年差的平均值
vcha=0.
uvqs=new((/2,naq/),float)  ;用来计算协方差矩阵
uvrs=new((/2,nar/),float)
您的程序是强年与弱年的数目相同，如果不同，后面定义的uvqs与uvrs是不是不合适，
一般来讲强年与弱年的数目是不一定相同的，请问楼主有何良策？

pdikang

谢谢楼主111111111111111

西西

额。很给力啊！自己的原创，不错不错

gxlbxzx

感觉不大对  用你这个画施能文章里面的那个例子 就不一样，希望是我弄错

0708011

gxlbxzx 发表于 2014-12-11 12:32
感觉不大对  用你这个画施能文章里面的那个例子 就不一样，希望是我弄错

你可以贴个图看看。我自己试了是没错的。通过信度的关键区就是我选择的区域上。不过第一次编这种 我也不大确定。

gxlbxzx

检验的区域还是有差别的

麦田_smile

赞赞赞{:5_213:}{:5_213:}

0708011

gxlbxzx 发表于 2014-12-12 12:21
检验的区域还是有差别的

我也照着画了 确实不一样 但大致区域差不多  我又把程序检查了几遍，应该是对的。原因不清，可能是资料不一样。

gxlbxzx

0708011 发表于 2014-12-15 11:21
我也照着画了 确实不一样 但大致区域差不多  我又把程序检查了几遍，应该是对的。原因不清，可能是资料 ...

数据是一样的   都是ncep的   你看我的那个风场和他文章里面是完全一致的

0708011

gxlbxzx 发表于 2014-12-15 19:25
数据是一样的   都是ncep的   你看我的那个风场和他文章里面是完全一致的

嗯 是一样的。可我又仔细检查了  确实没有发现错误。我觉得是不是ncep再分析资料有变动过 他这个文章04年的。一些小改动使得检验的结果出现不一致。我仔细检查这个程序应该是对的。