testmakeh

 $Id: testmakeh.m,v 1.4 2004/05/11 15:08:06 stuver Exp $

0001 function [timeSeries, hij] = testmakeh(simId, varargin)
0002 
0003 % $Id: testmakeh.m,v 1.4 2004/05/11 15:08:06 stuver Exp $
0004 
0005 addpath ../
0006 
0007 hij = makeh(simId, 1e-21, 1, 16384);
0008 % Assumes a sampling rate of 16384 samples per second
0009 
0010 timeSteps = size(hij, 3);
0011 
0012 if length(varargin) == 0
0013     
0014     for iTimeStep = 1:timeSteps
0015         timeSeries(1,iTimeStep) = hij(1,1,iTimeStep);
0016         timeSeries(2,iTimeStep) = hij(1,2,iTimeStep);
0017         timeSeries(3,iTimeStep) = hij(1,3,iTimeStep);
0018         timeSeries(4,iTimeStep) = hij(2,1,iTimeStep);
0019         timeSeries(5,iTimeStep) = hij(2,2,iTimeStep);
0020         timeSeries(6,iTimeStep) = hij(2,3,iTimeStep);
0021         timeSeries(7,iTimeStep) = hij(3,1,iTimeStep);
0022         timeSeries(8,iTimeStep) = hij(3,2,iTimeStep);
0023         timeSeries(9,iTimeStep) = hij(3,3,iTimeStep);
0024     end
0025     
0026     figure(1)
0027     clf
0028     subplot(331)
0029     plot((1:timeSteps)./16384,timeSeries(1,:))
0030     xlabel('Seconds (assumes 16384 samples per second)')
0031     ylabel('Amplitude')
0032     grid on
0033     subplot(332)
0034     plot((1:timeSteps)./16384,timeSeries(2,:))
0035     title(['h_{ij} Timeseries  (' upper(simId) ')'])
0036     xlabel('Seconds (assumes 16384 samples per second)')
0037     ylabel('Amplitude')
0038     grid on
0039     subplot(333)
0040     plot((1:timeSteps)./16384,timeSeries(3,:))
0041     xlabel('Seconds (assumes 16384 samples per second)')
0042     ylabel('Amplitude')
0043     grid on
0044     subplot(334)
0045     plot((1:timeSteps)./16384,timeSeries(4,:))
0046     xlabel('Seconds (assumes 16384 samples per second)')
0047     ylabel('Amplitude')
0048     grid on
0049     subplot(335)
0050     plot((1:timeSteps)./16384,timeSeries(5,:))
0051     xlabel('Seconds (assumes 16384 samples per second)')
0052     ylabel('Amplitude')
0053     grid on
0054     subplot(336)
0055     plot((1:timeSteps)./16384,timeSeries(6,:))
0056     xlabel('Seconds (assumes 16384 samples per second)')
0057     ylabel('Amplitude')
0058     grid on
0059     subplot(337)
0060     plot((1:timeSteps)./16384,timeSeries(7,:))
0061     xlabel('Seconds (assumes 16384 samples per second)')
0062     ylabel('Amplitude')
0063     grid on
0064     subplot(338)
0065     plot((1:timeSteps)./16384,timeSeries(8,:))
0066     xlabel('Seconds (assumes 16384 samples per second)')
0067     ylabel('Amplitude')
0068     grid on
0069     subplot(339)
0070     plot((1:timeSteps)./16384,timeSeries(9,:))
0071     xlabel('Seconds (assumes 16384 samples per second)')
0072     ylabel('Amplitude')
0073     grid on
0074     shg
0075 elseif length(varargin) == 2
0076     [row, col] = deal(varargin{1}, varargin{2});
0077     
0078     for iTimeStep = 1:timeSteps
0079         timeSeries(iTimeStep) = hij(row, col, iTimeStep);
0080     end
0081     
0082     figure(1)
0083     clf
0084     plot((1:timeSteps)./16384,timeSeries)
0085     title(['h_{' num2str(row)  num2str(col) '} Timeseries  (' upper(simId) ')'])
0086     xlabel('Seconds (assumes 16384 samples per second)')
0087     ylabel('Amplitude')
0088     grid on
0089     shg
0090 else
0091     error('testmakeh: length(varargin) must be 0 or 2')
0092 end    
0093 
0094 return