15.08.2013 Views

Journal of Accident Investigation

Journal of Accident Investigation

Journal of Accident Investigation

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

JOSEPH GREGOR<br />

When the correlated data for this case were plotted as in<br />

figure 9, the resulting curve fit turned out to be linear, but with<br />

b = 0.9831. This indicated that the time-base in one unit<br />

was operating at a constant but different rate when compared<br />

with the other. If the CVR had been a tape-based unit, the<br />

assumption would have been that the tape speed was <strong>of</strong>f, and<br />

the resulting transform applied to adjust the CVR data to match<br />

the FDR data. However, since both units were solid-state in<br />

this case, a different assumption was required.<br />

Researching these correlation problems with the manufacturer<br />

turned up an anomaly in the FDR recording system. In a typical<br />

installation, data are fed to the FDR in groups identified by SRN.<br />

Each group represents a repeating set <strong>of</strong> data points for the<br />

various parameters being measured. These groups are sent to<br />

the FDR at a rate <strong>of</strong> 1 per second – and each SRN is assumed to<br />

represent 1 second <strong>of</strong> flight data. In this particular installation,<br />

however, the data stream was being fed to the FDR at a rate<br />

<strong>of</strong> 61 groups per minute. As a result, each SRN represented<br />

roughly 60/61 = 0.9836 seconds <strong>of</strong> data, consistent with<br />

transform results found for b. Once this transform was applied,<br />

the corrected FDR on-key/<strong>of</strong>f-key data were found to overlay<br />

clearly with the CVR data at the <strong>of</strong>fset originally determined by<br />

the cross-correlation function.<br />

CONSIDERATIONS FOR MORE COMPLEX<br />

TRANSFORMS<br />

Valid correlations may always be obtained by assuming<br />

a linear transform as given in Eq. (2), provided that the<br />

time-bases in both units are operating at a fixed, constant rate.<br />

This will not be the case, however, if the time-base in one or<br />

both units changes rate during the recording. This may occur in<br />

the case <strong>of</strong> a solid-state recorder if power is momentarily lost to<br />

one <strong>of</strong> the boxes. In this instance, a linear fit as in Eq. (2) would<br />

still be valid on each side <strong>of</strong> the discontinuity. This may also<br />

occur for a tape-based CVR due to power supply fluctuations<br />

or mechanical changes in the tape drive mechanism. Either<br />

<strong>of</strong> these factors could change the tape speed, and hence the<br />

amount <strong>of</strong> information corresponding to 1 second <strong>of</strong> digitized<br />

data. When played back at a constant speed, the effect will be<br />

an apparent change in the time-base for that audio recording.<br />

Note that, since a tape-based CVR is an analog device, there<br />

is no reason to assume that the tape speed must change from<br />

one discrete constant value to another discrete constant value.<br />

Instead, the tape speed – and hence the time-base for the<br />

playback – may in principle take on any analog value. A linear<br />

transform as in Eq. (2) will not suffice in this case.<br />

In the past, such correlations have been performed by<br />

employing a piecewise linear fit to subsets <strong>of</strong> the recorded<br />

data. These sections are then reassembled to obtain a single<br />

contiguous data set. This is essentially equivalent to employing<br />

Eq. (2) over segments <strong>of</strong> the data, and performing a curve fit to<br />

obtain the values for b and C corresponding to a unique solution<br />

within each segment. Matching solutions at each end, one<br />

could obtain a transform for the entire data set <strong>of</strong> the form,<br />

tcvr= b1 * tfdr + C1 t0 < tfdr ≤ t1<br />

tcvr= b2 * tfdr + C2 t1 < tfdr ≤ t2<br />

tcvr= b3 * tfdr + C3 t2 < tfdr ≤ t3<br />

tcvr= b4 * tfdr + C4 t3 < tfdr ≤ t4<br />

Eq. (5)<br />

This method is not without its disadvantages. Picking the<br />

limits required to obtain a correlation using a minimum number<br />

<strong>of</strong> segments for the desired accuracy involves a significant<br />

amount <strong>of</strong> labor. If the flow <strong>of</strong> time in one <strong>of</strong> the units varies<br />

continuously, as may be expected for an analog device, the<br />

number <strong>of</strong> segments required to obtain a good fit could become<br />

extremely large. In addition, while the magnitude <strong>of</strong> the data<br />

may match at the boundaries, the derivative <strong>of</strong> the resulting<br />

curve will necessarily be discontinuous at the ends <strong>of</strong> each<br />

segment – an unphysical result.<br />

We may, however, extend this method to obtain a more<br />

accurate result with far less labor by recalling from basic<br />

calculus that any continuous function can be represented by<br />

a number <strong>of</strong> straight line segments. The larger the number <strong>of</strong><br />

segments, the closer the resulting representation will come to<br />

the true curve. The original curve can be reproduced exactly<br />

in the limit that an infinite number <strong>of</strong> segments are employed.<br />

This is nothing more than the description <strong>of</strong> a generalized curve<br />

fit. An accurate correlation to any data set not yielding to a<br />

linear fit may then be obtained by cross-correlating the data<br />

to obtain a first-order time alignment, and then calculating a<br />

cubic spline fit or similar generalized curve fit. This procedure<br />

may be easily implemented on any desktop computer using a<br />

wide variety <strong>of</strong> s<strong>of</strong>tware tools. The results may be plotted as in<br />

figure 9 to check the validity <strong>of</strong> the final solution.<br />

CONCLUSIONS<br />

A new method has been described for performing the<br />

correlation <strong>of</strong> time as represented in the CVR with time<br />

as represented in the FDR. This method employs the use<br />

2 NTSB JOURNAL OF ACCIDENT INVESTIGATION, SPRING 2006; VOLUME 2, ISSUE 1

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!