`scatstat1`文档

的scatstat1函数返回的统计值给定的1 d半径内的所有点的值。这类似于将一个移动的意思,但点不等距的,x值也不需要单调递增。

另请参阅scatstat2。

内容

ybar = scatstat1 (x, y,半径)ybar = scatstat1 (x, y,半径,乐趣)ybar = scatstat1(…,选项)

ybar = scatstat1 (x, y,半径)返回所有的意思y值在指定半径每一点x。

ybar = scatstat1 (x, y,半径,有趣)适用于任何函数有趣的来y值,默认有趣的是@mean,但可以@median,@nanstd等。

选择ybar = scatstat1 (…)允许函数接受的任何选项。例如,“omitnan”。

得到当地的中值点在10单位每个x点。首先创建一些随机non-sorted数据与N = 10000分:

N = 10000;x =兰迪(300 N, 1) + 20 + 3 * randn (N, 1);y = 3 *信德(x) + randn(大小(x)) + 3;情节(x, y,“k”。)轴紧盒子从

现在15 x单位内所有点的移动值的每个值:

yb = scatstat1 (@median x, y, 15日,“omitnan”);持有在情节(x, yb,“波”)

注意,这个函数不插入或执行任何间距相等。最后的值对应于每个x点,这需要不等距的或排序。这里是一个时间性,所以你可以看到x值不等距的:

轴((167 173 6))

从1978年到1988年,NSIDC海冰时间序列包含关于一个测量大约每两天。1988年之后,这是一个日常的时间序列。如果你想要一个2年海冰的移动平均时间序列,在这个数据集没有常数时间分辨率?

的现状是:

负载seaice_extent.mat%加载示例数据图绘制(t, extent_N)轴紧盒子从ylabel北半球的海冰范围1 \ times10 ^ 6公里)'

在1988年之前,可用的数据是每隔一天,然后有一个差距的数据,然后它的日常数据。情节之间的时间测量(天)的梯度datenum价值观:

图绘制(t,梯度(datenum (t))) ntitle“天之间测量”盒子从

你看,有一个缺口。

两年移动平均值(半径为1年),将datetime格式t来datenum的格式,并使用scatstat1半径为365.25天。

extent_N_2yr = scatstat1 (datenum (t) extent_N, 365.25);图(1)%的海冰时间序列图持有在情节(t, extent_N_2yr“k”,“线宽”,2)

当然,你可能不想相信任何1内半径(1年2年移动平均值)的目的,这就是为什么奇怪的摆动出现在开始和结束的时间序列。

这个函数是由乍得a。格林德克萨斯大学奥斯汀分校的地球物理研究所(UTIG), 2016年6月。