shape - E. A. Robinson’s subroutine to shape one signal to another
call shape(lb, b, ld, d, la, a, lc, c, ase, space, wht)
Basic Seismic
Utilities (BSU) subroutine originally published by Robinson
(1967) and modified by Michaels to include a
"whitening" stabalizaton factor. This factor adds
a fraction of the zero lag autocorrelation back to the zero
lag autocorrelation in the usual way to avoid unstable
results from inverting spectral notches. The routine calls
eureka which solves the normal equations by the Levinson
recursion algorithm.
Fortran Version.
Arguments
lb (int)
Number of
samples in signal b
b (float)
Signal to be shaped
ld (int)
Number of
samples in signal d
d (float)
Desired output signal from shaping b
la (int)
Number of
samples in signal a
a (float)
The filter which does the shaping, output from subroutine
lc (int)
Number of
samples in signal c
c (float)
The actual output from applying the filter to b, should approximate d.
ase (float)
Averge square error. The smaller the better the shaping.
space (float)
Work space vector (size 3*la)
wht (float)
The amount of whitening, fraction of zero lag autocorrelation to add back in.
NOTE:
Calls cross, eureka, dot, and fold
subroutines.
cross(3), fold(3), dot(3), eureka(3)
Robinson, E.A.,
1967, "Multichannel Time Series Analysis
with Digital Computer Programs", Holden-Day, 298p.
Levinson, N.,
1946, "The Wiener RMS (root mean square) error
criterion
in filter design and prediction", J. Math. Phys.
25, 261-278.
No known bugs.
Copyright © 1967 by Holden-Day Inc.
Holden-Day Inc., Liquidation Trust, 1259 S.W. 14th Street Boca Raton, FL 33486 Phone: 561.750.9229 Fax: 561.394.6809
This subroutine is included for distribution with Basic Seismic Utilities (BSU) by joint permission of the author, Enders Robinson, and Holden-Day Inc. Liquidation Trust, Frederick H. Murphy, Trustee. Being part of the BSU package, it is governed by the terms of the GPL license, and is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
E. A. Robinson