function m4(tt)
t=size(tt);
n=t(3);
for
i=1:n
for b=1:9
r=tt(:,b,i);
t=tt(:,b+9,i);
deglich(r,t);
load R
load T
f(i,b)=rtf1(R,T);
R_T=R-f(i)*T;
RnT=R+f(i)*T;
for h=1:6
for k=1:4
vd(k)=R_T(h,k,13)/var(R_T(h,k,1:R_T(h,k,13)));
vs(k)=RnT(h,k,13)/var(RnT(h,k,1:RnT(h,k,13)));
ivd(k)=1/vd(k);
ivs(k)=1/vs(k);
if k==2 | k==3
md(k)=
-sum(R_T(h,k,1:R_T(h,k,13)))/R_T(h,k,13);
ms(k)=
-sum(RnT(h,k,1:RnT(h,k,13)))/RnT(h,k,13);
else
md(k)=
sum(R_T(h,k,1:R_T(h,k,13)))/R_T(h,k,13);
ms(k)= sum(RnT(h,k,1:RnT(h,k,13)))/RnT(h,k,13);
end
end
for k=1:4
wd(k)=1/4;
ws(k)=1/4;
%
wd(k)=vd(k)/sum(vd);
%
ws(k)=vs(k)/sum(vs);
end
difff(h,b,i)=wd*md';
sigdif(h,b,i)=sqrt(sum(ivd)/16);
summm(h,i,b)=ws*ms';
isigsummm2(h,i,b)=16/sum(ivs);
end
end
end
save f
save difff
save sigdif
%%%%%%%%%
summ %%%%%%%%%
for b=1:9
wi=sum(sum(isigsummm2(:,:,b)));
summ1=isigsummm2(:,:,b)*summm(:,:,b)'/wi;
summ(b)=0;
for
k=1:6
summ(b)=summ(b)+summ1(k,k);
end
subplot(3,3,b)
j=0;
for i=1:9
for h=1:6
j=j+1;
ssum(j)=summm(h,i,b);
hold on;
set(gca,'XLim',[0 n*6]);
quiver(j,summm(h,i,b),0,sqrt(1/isigsummm2(h,i,b)),'r*');
quiver(j,summm(h,i,b),0,-sqrt(1/isigsummm2(h,i,b)),'r*');
hline(summ(b),'b-');
hold off;
end
end
sigsum2(b)=var(ssum)/54;
end
summ
save
summ
save sigsum2
%^^^^^^^^^^^^^^^^^^^^^^^^^^^^