Moore指数计算时间序列的Mt以及theta角度值Matlab源代码:
%计算Wit2值
for i = 0:34;
i=i+1;
for j = 1:3;
j=j+1;
Wit2(i,j) = Wit(i,j)*Wit(i,j);
Wit2(i,1) = Wit(i,1);
end
end
%计算Mt值
%先算一下分子
a=0;
b=0;
i=0;
for i = 0:33;
i=i+1;
for k = 1:3;
k = k+1;
a = Wit(i,k)*Wit(i+1,k);
b = b+a;
end
LuPengyaosuandefenzi(i,1) = b;
a = 0;
b = 0;
end
%再算一下分母
LuPengyaoxiedefenmu = [];
i=0;
k=0;
c1=0;
d1=0;
e1=[];
for i = 0:33;
i = i+1;
for k = 1:3;
k = k+1;
c1 = Wit2(i,k)*Wit2(i,k);
d1 = d1+c1;
end
e1(i,1) = sqrt(d1)
c1=0;
d1=0;
end
i=0;
k=0;
c2=0;
d2=0;
e2=[];
for i = 0:33;
i = i+1;
for k = 1:3;
k = k+1;
c2 = Wit2(i+1,k)*Wit2(i+1,k);
d2 = d2+c2;
end
e2(i,1) = sqrt(d2)
c2=0;
d2=0;
end
LuPengyaoxiedefenmu = e1.*e2;
%计算Mt的值
i = 0;
Mt = [];
for i = 0:33;
i = i+1;
Mt(i,1) = LuPengyaosuandefenzi(i,1)/LuPengyaoxiedefenmu(i,1);
end
%计算theta的值
theta = [];
i=0;
for i = 0:33;
i = i+1;
theta(i,1) = acos(Mt(i,1));
end
|