Utilizing Recursive Least Squares Method Computer Science Essay

Published: Last Edited:

This essay has been submitted by a student. This is not an example of the work written by our professional essay writers.

Non linear restricted element analyses of structure that are subject to seismic proceedings have necessary of high quality accelerogram data. Raw accelerogram data needs to be familiar to remove the pressure of the transfer process of the gadget itself. This techinique is known as upgrading. Unfortunately, data about the recording gadget is often unknown or undependable. This is above all often the case for older analogue recordings. This paper utilizes a recursive least squares (RLS) algorithm to recognize the utensil characteristics even when completely unknown. The fallout presented in the paper applies a modern approach to de noising the accelerogram by employing the wavelet transform. This system removes only those workings of the signal whose amplitudes are below a confident doorsill and is not therefore occurrence selective. It takes over from to some extent predictable band pass filter which requires a careful mixture of cut-off frequencies, now needless.


Seismology (is derived from the Greek words seisms, (earthquake) and logia, as a whole "Talking about earthquakes") is the technical study of earthquakes and the broadcast of elastic waves from side to side the Earth. The study of this meadow also contain study of earthquake belongings, like tsunamis as well as diverse seismic sources like volcanic, tectonic, oceanic, atmospheric, and artificial processes (like explosions). Pale seismology is a related field that utilizes geology to infer data regarding past earthquakes. Seismogram is identified a recording of earth movement as a principle of this time.

Seismic waves

Various types of seismic waves are produced by Earthquakes, and other sources, which travel from side to side, astound, and facilitate an effective way to image both sources and structures deep within the Earth. In solids there are three usual types of seismic waves: P-waves, S-waves (both are called body waves) and interface waves. There are other two usual kinds of surface waves (Rayleigh and Love) which travel along a solid-air interface, can be fundamentally explained in terms of interacting P- and/or S-waves.


The earliest brand of body wave is the P-wave, which is also called primary wave. This is the earliest to 'arrive' at a seismic situation, and, as a output, the best brand of seismic wave. The P wave can move from side to side solid astound and fluids, which contains water or the liquid layers of the earth. It shoves and drags the astound it moves from side to side immediately like sound waves push and pull the air. Consist you ever heard a big clap of thunder and heard the windows rattle at the same time? The windows rattle due to the sound waves were pushing and pulling on the window glass much like P waves push and pull on astound. At times animals can perceive sound the P waves of an earthquake. Generally people can only feel the bump and clatter of these waves.

P waves generally known as compressional waves, because utilization of their pushing and pulling behavior.

Figure 1 - A P wave travels from side to side a medium by means of firmness and dilation. Subdivision is standing for by cube in this model.


An S wave is slower compared to a P wave and can only move from side to side solid astound, not from side to side any liquid medium. Based on this property of S waves seismologists conclude that the Earth's outer core is a liquid. S waves move astound particles in the up and down direction, or side-to-side, perpendicular to the direction that the wave is traveling in (the direction of wave propagation).

Figure 2 - An S wave travels from side to side a medium. Particles are represented by cubes in this model.

Surface Waves

Surface waves will have low frequency than body waves, which travels only from side to side the crust, and are easily distinguished on a seismogram as an output. Surface waves are nearly absolutely accountable for the damage and demolition related with earthquakes; however they turn up after body waves. In the deeper earthquakes this damage and the strength of the surface waves are condensed.

Love Waves

The earliest brand of surface wave is called a Love wave, named after a British Mathematician A.E.H. Love, who worked out the mathematical model for this brand of wave in 1911. Restrained to the float up of the crust, Love waves manufacture absolutely horizontal motion.

Figure 3 - A Love wave travels from side to side a medium. Particles are represented by cubes in this model.

Rayleigh Waves

The other brand of surface wave is the Rayleigh wave, named for John William Strutt, Lord Rayleigh, who mathematically predicted the existence of this brand of wave in 1885. Primarily of the trembling felt from an earthquake is due to the Rayleigh wave, which can be much larger than the other waves.

Figure 4 - A Rayleigh wave travels from side to side a medium. Particles are represented by cubes in this model.

For large earthquakes, the normal modes of the Earth can be observed. The earliest observations were made in the 1960s as the advent of higher fidelity instruments coincided with two of the largest earthquakes of the 20th century - the 1960 Great Chilean earthquake and the 1964 Great Alaskan earthquake.

One of the earliest significant discoveries (suggested by Richard Dixon Oldham in 1906 and definitively shown by Harold Jeffreys in 1926) was that the outer core of the Earth is liquid. Pressure waves (P-waves) pass from side to side the core. Transverse or shear waves (S-waves) that vibrate side-to-side requires rigid material so that they do not pass from side to side of the outer core. Therefore, the liquid core cause a "shadow" on the side of the planet opposite of the earthquake where no direct S-waves are observed. When the P-wave velocity of the outer core also is reduced, it causes a considerable delay for P waves penetrating the core from the (seismically faster velocity) mantle.

Seismic waves which are produced by explosions or vibrating controlled sources are one of the primary methods of underground exploration in geophysics (in addition to lot of various electromagnetic methods like induced polarization and magnetotellurics). Controlled source seismology has been utilized to map salt domes, faults, anticlines and other geologic traps in petroleum-bearing rocks, geological faults, astound types, and long-buried giant meteor craters. For instance, the Chicxulub impactor, which is believed to consist killed the dinosaurs, was localized to Central America by analyzing ejecta in the cretaceous boundary, and then physically proven to exist utilizing seismic maps from oil exploration.

Utilizing seismic tomography with earthquake waves, the interior of the Earth has been absolutely mapped to a resolution of several hundred kilometers. This technique has enabled scientists to identify convection cells, mantle plumes and other large-scale features of the inner Earth.

Seismograph is an instrument that sense and record the motion of the Earth. Nowadays, networks of seismographs continuously monitor the seismic environment of the planet, allowing for the monitoring and analysis of global earthquakes and tsunami warnings, as well as recording a variety of seismic signals arising from non-earthquake sources ranging from explosions (nuclear and chemical), to pressure variations on the ocean floor induced by ocean waves (the global microseism), to cryospheric events related with large icebergs and glaciers. Above-ocean meteor strike as large as ten kilotons of TNT, (equivalent to about 4.2 Ã- 1013 J of effective explosive force) consist been recorded by seismographs. The monitoring of nuclear testing is a major motivation for the global instrumentation of the Earth.

One of the earliest attempts at the scientific study of earthquakes followed the 1755 Lisbon earthquake. Other notable earthquakes especially that spurred major developments in the science of seismology contain the 1857 Basilicata earthquake, 1906 San Francisco earthquake, the 1964 Alaska earthquake and the 2004 Sumatra-Andaman earthquake. An extensive list of famous earthquakes can be identified on the List of earthquakes page.

Forecasting a probable timing, location, magnitude and other significant features of a forthcoming seismic event is called earthquake prediction. Primarily seismologists do not believe that a system to facilitate timely warnings for individual earthquakes have yet been developed and lot of believe that such a system would be unlikely to provide significant warning of impending seismic events. More general forecasts, however, are routinely utilized to establish seismic hazard. Such general forecasts estimate the probability of an earthquake of a particular size affecting a particular location within a particular time span and they are routinely utilized in earthquake engineering. Various attempts consist been made by seismologists and others to create effective systems for precise earthquake predictions, including the VAN technique.


Seismic Data Improvement

The main point of seismic processing is to use the acquired data into an image that can be utilized to infer the sub-surface structure. Processing contain the application of a series of computer routine to the obtain information directed by the hand of the processing geophysicist. The interpreter must be involved at all stages to check that processing decisions do not radically alter the interpretability of the results in a detrimental manner.

In general, processing routines fall into one of the following categories:

* enhancing signal at the expense of noise

* facilitating velocity data

* breaking up diffractions and placing dipping events in their true subsurface locations

* increasing resolution.

Correction of Seismic information records

Correction of Seismic information records involves standard techniques of de convolving instrument and structural responses from seismic accelerograms. Instrument and structural responses are a result of convolution of ground motion with the transfer process of the recording instrument and structure on which the instrument is mounted. To recover an approximation of the ground motion of the seismic event the instrument response can be de-convolved in the time or frequency domain, which is usually followed by a band-pass filtering of the information records. A band-pass filter removes any ground motion outside the band and hence the accelerograms may not effectively represent an approximation of true ground motion. In order to overcome this it is proposed therefore to utilize wavelet de noising as an substitute to band-pass filtering. The key concepts of de-noising are it removes minimum and high frequency corrupting signals but holds relevant information records providing a better approximation of true ground motion. Utilizing the stationary wavelet transform (SWT) some data collected as sample seismic signals are threshold de noised and match up to with the more standard band-pass filtering techniques.

De noising utilizing a band-pass filtering techinique

Digital band-pass filter are designed to separate the minimum and high frequency noise. Ormsby filter is the widely utilized digital band-pass filter. In order to reach the essential band-pass design specifications the Ormsby filter essential large values of N and at the time utilizes a considerable computational effort. In order to reach the essential stringent band-pass specification and reduce computation time Infinite impulse response (IIR) digital filters are now normally employed. The phase response of these filters is non linear whereas a zero phase response is essential. While it is still superior to the computational effort essential for the FIR filter there is a loss of computational optimality. While, by executing the information records in both the forward and reverse direction a zero phase response is also secured. In this techinique the information records is first filtered in the forward direction, then the filtered information records is reversed and run back through the filter. This techinique is computationally more efficient than former ones.

De noising utilizing wavelets

In this technique the discrete wavelet transform (DWT) of a signal will be taken and passed the transform through a threshold which removes the coefficients below a certain value and then taking the inverse DWT in order to rebuild a de noised time signal. After minimum pass filtering with the appropriate filter weights depending on the choice of a wavelet basis the DWT is able to concentrate most of the energy of the signal into a small number of wavelet coefficients. Match up to to those of the noise coefficients obtained after high pass filtering the dimensions of the wavelet coefficients will be large. Hence thresholding or shrinking the wavelet transform will eliminate the minimum amplitude noise in the wavelet domain and the inverse DWT will recover the preferred signal with minute loss of detail.

The stationary wavelet transform (SWT) (translation invariant DWT)

The key problem with the DWT approach is that it is not translation invariant. The signal is no longer orthogonal to most of the basis functions due to the coefficients of the DWT do not shift with a signal. To describe the signal more coefficients would be necessary and the co-efficient dimensions would also be much smaller, which reduces the effectiveness of any de noising scheme.

While from the outcome obtained it can be clearly stated that the differences among utilizing a DWT and SWT for de noising are insignificant; output plots like power spectral densities and the acceleration response spectra are virtually the same. While discontinuities in seismic records do occur therefore the SWT seems a more effective choice than DWT and digital band-pass filters for de noise purposes. The SWT only removes "noise" that has minimum energy content, which is independent of frequency. SWT de noising is computationally efficient and obviates the requirement to adjust filter cut off's to fit particular seismic events. For different groups of researchers around the world it is evident that choice of filter cut off frequencies varies. Rather unsurprisingly, at the minimum and high frequency range of the spectrum, the differences among band-pass filtering and SWT techniques exists. In the design of large dams or tall building structures the minimum frequency or long period end is of importance.

Different techniques of seismic information records correction:

Recursive least squares (RLS) Algorithm

This is an Adaptive filter algorithm utilized to find the filter coefficients that relate to recursively producing the least squares of the error signal. RLS algorithm is contrast to other algorithms, which aims to reduce the mean square error. The key difference among RLS filters and MSE filters is that RLS filters are dependent on the signals themselves, whereas MSE filters are dependent on their statistics.

Recursive Least Squares and Wavelet De Noising based Seismic correction demonstrates inverse filtering utilizing the RLS algorithm, which yields acceptable outcome when match up to the standard 2nd order type de convolution. In order to de-convolve the instrument response the follow-on inverse filter is then implemented to the information records. Power spectral plots and the total acceleration response spectra are equivalent of two sets of seismic events in this procedure. The key advantage of RLS over LMS is that RLS algorithm is dependent on the incoming information records samples rather than the statistics of the ensemble average as in the case of the LMS algorithm, this is the reason why the RLS algorithm was chosen in predilection to the least mean squares (LMS) adaptive algorithm.

Total Least Squares

When instrument parameters are not represented the TLS algorithm is a useful tool for correcting seismic information records. The basic advantage of a least squares based techinique of de convolution of an instrument response from seismic information records is that it does not necessitate any data regarding the instrument; it only necessitates the recorded seismic information records. The only requirement is the original recording from the seismograph and the algorithm can then produce the inverse filter with which to de convolve the instrument response. Utilizing information records from 4 instrument types the algorithm was gone through testing and was found to be in good agreement in all cases except with the QR RLS and the 2nd order SDOF responses. From this we can conclude that the TLS may not be as robust as the QR RLS in securing the instrument response, while the TLS performance in some cases is better than that of the QR RLS and the 2nd order SDOF with a standard filter, due to it reflects the anti alias filter whose details were in this case available in the record.. A particular problem in seismic correction techniques is that some times the transfer process of the recording instruments is not represented. In order to de couple the instrument response where instrument parameters are afforded, a 2nd order SDOF transfer process is implemented in either the time or frequency domain.

On the contrary when instrument parameters are not available, then if we don't make some assumptions regarding instrument parameter, it is not possible to de convolve (de couple) the instrument response from the seismic information records utilizing a 2nd order SDOF expression. An estimation of the instruments inverse frequency response was afforded with which to de convolve the information records to get an approximation of the ground motion.

Partial Total Least Squares approach

The problem in the Total Least Squares (TLS) AX ≈ B is solved by The Partial Total Least Squares (PTLS) subroutine by utilizing a Partial Singular Value Decomposition (PSVD), hereby improving the computational efficiency of the classical TLS algorithm (a factor 2 approximately). The TLS solution will be computed from any orthogonal basis of the right singular subspace corresponding to the smallest singular values of the augmented matrix [A;B] can be found based on a clear study, the computational complexity can be improved. For the smallest singular values the dimension of this subspace may be specified by the rank of the TLS approximation of [A;B] or depend on a specified upper bound (e.g. a perturbation level) .More than one vector in B is supported in the routine PTLS. By computing the least amount norm solution it also solves specified and undetermined sets of equations.

It has been demonstrated that the PTLS algorithm produces good outcome when match up to to the standard TLS. While the inverse responses are not representative of the instrument response when more rows and/or columns unmodified. Therefore applying the PTLS as opposed to the TLS isn't a net gain.

SDOF techniques

Seismic correction techniques generally apply a 2nd order, single degree of freedom (SDOF) instrument process with which to inverse filter or de convolute the accelerometer response. Unfortunately the significant instrument parameters are either not specified or are not representing in many of seismic time histories. Therefore in order to de convolute instrument response researchers utilize an approximation for the instrument parameters, or they do not instrument de convolute at all. Recently more modern inverse system identification techniques contain been implemented to de-convolute the instrument response. The main important techniques include the QR. Recursive Least Squares (QR RLS), the Total least Squares with wavelet de-noising and wavelet decomposition.

QR Recursive Least Squares (QR RLS) method:

QR decomposition based RLS algorithm is developed based on the square root Kalman filter counterpart. The algorithm is derived depending on the utilization of an orthogonal triangulation procedure represent as QR decomposition. Based on a Givens rotation or a Householder reflection the unitary matrix utilized in the algorithm is derived, which zero's out the elements of the input information records vector and changed the square root of the inverse correlation matrix. In order to get an approximation of the true ground motion these inverse filter weights are an approximation of the inverse transfer process of the instrument and these are then convoluted with the original seismic information records.

In this paper we developed a seismic data analysis technique based on the RLS algorithm, which can effectively correct the data of differing spectra; provided they fit the model assumed by the RLS method. We examined the RLS prediction error as an event analyzed signal and in addition developed the coefficient change signal for seismic data. Our experiments indicate that this technique is very effective about the variations in the event ordering, event positions or starting position of the algorithm consist only minor effects on the filter coefficients. Therefore, this technique must be useful in those situations where the events are close to all-values and the noise level is not high.