Global Index (short | long) | Local contents | Local Index (short | long)
c =
0.9975 0.9972 0.9981
0.9954 0.9929 0.9965
0.9663 0.9680 0.8842
0.9510 0.9692 0.8838
0.9715 0.9796 0.8202
0.9019 0.9562 0.8204
0.8747 -0.9577 0.9182
0.8049 -0.8268 -0.7883
0.6526 0.8392 -0.5245
0.8167 -0.7615 0.4552
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs/data
load heat_pcs.mat
rpcs1 = rpcs; rlds1 = rlds; rper1 = rper; rlam1 = rlam;
hpcs1 = hpcs; hlds1 = hlds; hper1 = hper; hlam1 = hlam;
lpcs1 = lpcs; llds1 = llds; lper1 = lper; llam1 = llam;
cd /home/disk/tao/dvimont/Papers/hbcsiro/matlab/data
load heat_pcs.mat
for num = 1:10;
c(num,:) = [corr(rpcs1(:,num), rpcs(:,num)) ...
corr(hpcs1(:,num), hpcs(:,num)) ...
corr(lpcs1(:,num), lpcs(:,num))];
end
for num = 1:10;
c(num,:) = [corr(rlds1(:,num), rlds(:,num)) ...
corr(hlds1(:,num), hlds(:,num)) ...
corr(llds1(:,num), llds(:,num))];
end
%c =
%
% 0.9950 0.9955 0.9949
% 0.9932 0.9902 0.9948
% 0.9628 0.9644 0.8402
% 0.9466 0.9655 0.9122
% 0.9627 0.9771 0.7989
% 0.9065 0.9553 0.8292
% 0.8575 -0.9562 0.9135
% 0.8040 -0.8199 -0.7917
% 0.6539 0.8392 -0.5424
% 0.8107 -0.7578 0.4445
c = [rper1(1:10); rper(1:10); hper1(1:10); ...
hper(1:10); lper1(1:10); lper(1:10)];
%c =
%
% Columns 1 through 7
%
% 23.0987 14.5977 6.5529 5.4953 4.3051 3.3947 3.0556
% 25.5371 13.9148 6.4043 5.9417 4.0821 3.3908 2.7992
% 24.3872 10.3137 8.0885 5.9591 4.3519 3.1088 2.7630
% 27.3167 10.0151 7.5675 6.0246 4.4364 3.3145 2.5804
% 27.0578 16.4205 6.4654 5.6393 4.4124 3.7083 3.4323
% 29.3011 15.7046 6.4292 5.6882 4.3123 3.5863 3.2916
%
% Columns 8 through 10
%
% 2.7843 2.3224 2.1713
% 2.5232 2.3092 2.1930
% 2.5973 2.4224 2.0055
% 2.4645 2.1879 1.8516
% 2.6298 2.4116 2.3510
% 2.5819 2.3258 2.1829
c = [rlam1(1:10); rlam(1:10); hlam1(1:10); ...
hlam(1:10); llam1(1:10); llam(1:10)];
%c =
%
% 1.0e+19 *
%
% Columns 1 through 7
%
% 1.4200 0.8974 0.4029 0.3378 0.2647 0.2087 0.1878
% 1.2880 0.7018 0.3230 0.2997 0.2059 0.1710 0.1412
% 0.5806 0.2455 0.1926 0.1419 0.1036 0.0740 0.0658
% 0.5434 0.1992 0.1505 0.1198 0.0882 0.0659 0.0513
% 1.0200 0.6190 0.2437 0.2126 0.1663 0.1398 0.1294
% 0.8957 0.4800 0.1965 0.1739 0.1318 0.1096 0.1006
%
% Columns 8 through 10
%
% 0.1712 0.1428 0.1335
% 0.1273 0.1165 0.1106
% 0.0618 0.0577 0.0477
% 0.0490 0.0435 0.0368
% 0.0991 0.0909 0.0886
% 0.0789 0.0711 0.0667
ct = getct;
ct = detrend(ct);
[b1, a1] = butter(9, 2/10);
[b2, a2] = butter(9, 2/9);
hpct = ct - filtfilt(b1, a1, ct);
lpct = filtfilt(b2, a2, ct);
c = [corr(rpcs1(:,1), ct) corr(rpcs(:,1), ct); ...
corr(hpcs1(:,1), hpct) corr(hpcs(:,1), hpct); ...
corr(lpcs1(:,1), lpct) corr(lpcs(:,1), lpct)];
%c =
%
% -0.8858 -0.9025
% 0.8702 0.8867
% -0.8814 -0.8997