Documentation of plot_pcs

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

Help text

         tim(timind), pcs(timind, 2), '--k');

Cross-Reference Information

Listing of script plot_pcs

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs3 = pcs;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
load RAW_detrend_L1-7_EOF_yr101-1000.mat; pcs3 = pcs;
tim = 101:1000;
figure(1);
for i = 1:3;
  timind = 300*(i-1) + [1:300];
  subplot(4,1,i);
    plot(tim(timind), pcs(timind, 2), '-k');%, ...
    axis([min(tim(timind))-1 max(tim(timind))+1 -3 3]);
    grid on
    set(gca, 'YTick', -3:3, ...
        'XTick', [(min(tim(timind))-1):25:(max(tim(timind))+1)]);
    ylabel('STD')
end
subplot(4,1,1);
  title(['Normalized Unfiltered PC' num2str(num) ' of 0-270m Heat Content']);


%  Get filter response

[b, a] = butter(6, 2/10);
[h, w] = freqz(b, a, 128);
r = abs(h).^2;
f2 = w/(2*pi);

subplot(4,2,7);
  h1 = plot(f2, r, '-k', f2, (1-r), '--k');
  axis([0 .25 -.1 1.1]);
  grid on
  set(gca, 'YTick', 0:.25:1, 'XTick', 0:.05:.5);
  title('Butterworth Filter Response');
  legend(h1, 'LP10', 'HP10');

%  Get spectrum from compare_cts.m routine

subplot(4,2,8);
       semilogy(f2, p(:,1), 'b-', f2, rn, 'b-', ...
	      f2, rner5, 'b--', f2, rner1, 'b-.')
     set(gca, 'YTick', [.01 .1 1 10 100]);
     axis([0 0.25 .1 15])
     grid
     title(['Power Spectrum for UNFILT PC' num2str(num)]);
     xlabel(['Frequency:  yr^-^1'])
     set(gca, 'XTick', 0:.05:.5);

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap1/Plots