Documentation of compare_cts3

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

Cross-Reference Information

Listing of script compare_cts3

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load NP_RAW_T1_EOF_yr101-550.mat;  pcs1 = pcs;
load NP_RAW_T1_EOF_yr551-1000.mat; pcs2 = pcs;
load NP_RAW_T1_EOF_yr101-1000.mat; pcs3 = pcs;
num = 1;
for i = 1:3;
  if i == 1;
    ctann = pcs1(:,num);
    tit = 'Years 101 - 550';
  elseif i == 2;
    ctann = pcs2(:,num);
    tit = 'Years 551 - 1000';
  elseif i == 3;
    ctann = pcs3(:,num);
    tit = 'Years 101 - 1000';
  end;
  nx = length(ctann);
  nfft = 128; 
  noverlap = nfft/2;
  [p, f] = spectrum(ctann, nfft, noverlap);
  n2 = nfft/2;
  f2 = 0.5 * (0:n2)/n2;
  rho = (corr(ctann, ctann, 1) + sqrt(corr(ctann, ctann, 2)))/2;% + ...
  rn = (1-rho^2) ./ (1-2*rho*cos((pi/n2)*(0:n2))+rho^2);
  rn = (rn / mean(rn))*mean(p(:,1));
  if i == 3;
    rner1 = rn * 2.05;
    rner5 = rn * 1.68;
  else
    rner1 = rn * 2.64;
    rner5 = rn * 2.01;
  end
  dofx = nx / n2; 
  figure(num); figure_orient;
  subplot(3,1,i)
     semilogy(f2, p(:,1), 'b-', f2, rn, 'b-', ...
	      f2, rner5, 'b--', f2, rner1, 'b-.')
     axis([0 0.5 .05 15])
     grid
     set(gca, 'XTick', [0:.05:.5]);
     title(['Power Spectrum for PC' num2str(num) ', 20N - 65N, ' tit]);
     xlabel(['Frequency:  yr^-^1;  Degrees of Freedom:  ' ...
              num2str(round(dofx)) ';' ...
             '  Bandwidth:  7.8 x 10^-^3 yr^-^1'])
     l=legend(['PC1 Spectrum'],...
              'Red Noise', '95% Confidence', '99% Confidence');
end
cd /home/disk/tao/dvimont/matlab/CSIRO/Temperature/Plots_Spectra