mft.m

来自「气候分析中检验突变的经典程序」· M 代码 · 共 16 行

function f=mft(a,b)

%MFT moving-f test
%  load ltzhshshl
%  x=z;
%  for i=1:55
%  a=x(1:i);
%  b=x(i:55);
%  c(i)=mft(a,b);
%  end

n1=length(a);
n2=length(b);
s1=sum((a-mean(a)).^2)/(n1-1);
s2=sum((b-mean(b)).^2)/(n2-1);;
f=s1/s2;