Documentation of corr_LP_RAW

Global Index (short | long) | Local contents | Local Index (short | long)

Cross-Reference Information

Listing of script corr_LP_RAW

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1-7_EOF_yr101-1000.mat; pcsr = pcs; lamr = lam;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; pcsl = pcs; laml = lam;
load LP10_detrend_L1-7_yr101-1000.mat; pcsl = pcs; laml = lam;
[b, a] = butter(6, 2/10);
tempc = pcsr - filtfilt(b, a, pcsr);
a = NaN * ones(2);
for num = 1:2;
  a(num, 1) = corr(pcsr(:,num), pcsl(:,num), 0);
  a(num, 2) = corr(tempc(:,num), pcsl(:,num), 0);
end
LOWPASS:
a
a =
    0.8019    0.9824
    0.7618    0.9160
a.^2
ans =
    0.6430    0.9651
    0.5803    0.8390
HIGHPASS:
a
a =
    0.5798    0.9476
    0.4094    0.6884
a.^2
ans =
    0.3361    0.8979
    0.1676    0.4739