- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Mathematical Problems in Engineering

Volume 2012 (2012), Article ID 927821, 13 pages

http://dx.doi.org/10.1155/2012/927821

## Advanced Signal Processing and Command Synthesis for Memory-Limited Complex Systems

Faculty of Applied Sciences, Politechnica University, Hagi-Ghita 81, 060032 Bucharest, Romania

Received 30 June 2011; Revised 3 September 2011; Accepted 3 September 2011

Academic Editor: Carlo Cattani

Copyright © 2012 Cristian Toma. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper presents advanced signal processing methods and command synthesis for memory-limited complex systems. For accurate measurements performed on limited time interval, some specific methods should be added. For signal processing, a robust filtering and sampling procedure performed on a specific working interval is required, so as the influence of low-amplitude and high-frequency fluctuations to be diminished. This study shows that such a signal processing method for the case of memory-limited complex systems requires the use of certain differentiation/integration procedures performed by oscillating systems, so as robust results suitable for efficient command synthesis to be available. A brief comparison with uncertainty aspects in modern physics (where quantum aspects can be considered as features of complex systems) is also presented.

#### 1. Introduction

As it is known, an important aspect in observing and modeling dynamic environmental phenomena consists in measuring with higher accuracy some physical quantities corresponding to changes in the environment. Yet for accurate measurements performed on limited time interval, for memory-limited complex systems, some specific methods should be used. Sudden (sharp) changes in the environment require a pair of consecutive values for the measured quantity so as any difference to be detected as soon as possible. Moreover, any value taken into consideration by the complex system should be established using a robust filtering and sampling procedure performed on a specific working interval, so as the influence of low-amplitude and high-frequency fluctuations to be decreased in a significant manner. Being quite possible for sharp (sudden) changes in the environment to appear during such a working interval (on which filtering and sampling procedures are performed), it results that specific signal processing methods based on the values achieved on a set of succesive time intervals are necessary.

Filtering and sampling devices usually consist of asymptotically stable systems, sometimes an integration of the output over a certain time interval being added. Yet such structures are very sensitive at random variations of the integration period, being recommended for the signal which is integrated to be approximately equal to zero at the end of the integration period. For this reason, oscillating systems for filtering the received signal should be used, so as the filtered signal and its slope to be approximately zero at the end of a certain time interval (at the end of an oscillation). For avoiding instability of such oscillating systems on extended time intervals, certain electronic devices (gates) controlled by computer commands should be added, so as to restore the initial null conditions for the oscillating system before a new working cycle to start [1].

The filtering performances of asymptotically stable systems are determined by their transfer function. a Filtering and sampling devices consisting of low-pass filters of first or second order having the transfer function (for a first-order system) and (for a second-order system) attenuate an alternating signal of angular frequency about times (for a first-order system) or about times (for a second-order system). The response time of such systems at a continuous useful signal is about ( for the first-order system and for the second-order system). If the signal given by the first- or second-order system is integrated over such a period, a supplementary attenuation for the alternating signal of about can be obtained.

But such structures are very sensitive at the random variations of the integration period (for unity-step input, the signal which is integrated is equal to unity at the sampling moment of time), and the use of oscillators with a very high accuracy cannot solve the problem due to switching phenomena appearing at the end of the integration period (when an electric current charging a capacitor is interrupted).

These random variations cannot be avoided if we use asymptotically stable filters. For robustness, the signal processing structure based on an integration procedure should provide a null value for the integrating signal at the end of a certain working interval. This property is similar to wavelets aspects presented in [2, 3].

Mathematically, an ideal solution could consist in using an extended Dirac function for multiplying the received signal before the integration (see [1]) but is very hard to generate such extended Dirac functions (a kind of acausal pulses) using nonlinear differential equations for (i) symmetrical pulses (see [4]) or (ii) asymmetrical pulses (see [5] for more details).

A heuristic algorithm for generating practical test functions using MATLAB procedures was presented in [4]. First, it has been shown that ideal test functions cannot be generated by differential equations, being emphasized the fact that differential equations can only generate functions similar to test functions (defined as practical test functions). Then a step-by-step algorithm for designing the most simple differential equation able to generate a practical test function was presented, based on the invariance properties of the differential equation and on standard MATLAB procedures. The result of the algorithm consists in an oscillating second-order system working at the stability limit from initial null conditions, on a limited working interval corresponding to the period of the generated oscillations.

It was shown that the simplest structure possessing such properties is represented by an oscillating second-order system having the transfer function receiving a step input and working on the time interval . For initial conditions equal to zero, the response of the oscillating system at a step input with amplitude will have the form By integrating this result on the time interval , we obtain the result , and we can also notice that the quantity which is integrated and its slope are equal to zero at the end of the integration period. Thus, the influence of the random variations of the integration period (generated by the switching phenomena) is practically rejected.

This oscillating system attenuates about times such an input, and the influence of the integrator consists in a supplementary attenuation of about times. The oscillations having the form generated by the input alternating component have a lower amplitude and give a null result after an integration over the time interval .

These results have shown that such a structure provides practically the same performances as a structure consisting of an asymptotically stable second-order system and an integrator (response time of about , an attenuation of about times for an alternating component having frequency ) moreover being less sensitive at the random variations of the integration period. For restoring the initial null conditions after the sampling procedure (at the end of the working period), some electronic devices must be added. Yet the previous analysis is valid for step inputs which are active on the whole working interval (the integration period).

In [6] has been performed the analysis of this structure by considering that the input is represented by a unity short-step pulse (instead of a unity step pulse) which differs to zero on the time interval . It was shown that certain free oscillations of the second-order oscillating system are generated for (when the action of the external short-step command has ceased). These free oscillations have the angular velocity , the amplitude and the initial phase

Thus, the output corresponding to the free oscillations of the system for (when the action of the external short-step command has ceased) can be written as

However, we must notice that, usually, such a filtering and sampling structure receives an electronic signal presenting possible step changes from an already measured value to a final unknown value. Since the previously measured value can be substracted from the received signal during subsequent working intervals, the analysis of sudden (sharp) changes in the environment could start by considering that the input of the second-order oscillating system is represented by a null signal for (the first part of the working interval) and by a signal with amplitude for (the second part of the working interval).

#### 2. The Oscillating Signal Processing System for the Case of Short-Step Inputs

We will continue the analysis of this structure by considering that the input is represented by a short-step pulse which differs to zero on the time interval . This means that the input can be represented under the form or using the Heaviside function where corresponds to the function if we apply the Laplace transformation.

The transfer function of the second-order oscillating system is On the time interval , the output of the second-order oscillating system equals zero On the time interval , the output of the second-order oscillating system is represented (using the Laplace transformation) as which corresponds to the output which can be written as where .

By denoting with the integral of (considering as initial moment the zero moment of time), it results at the time moment the set of values:

It can be easily noticed that which can be written also as

This result shows that the sampled values for (the integral of ) and for at the time moment (the end of the working interval) can be used in a simple manner for obtaining the quantity

For this purpose, we can divide the sampled value for by and add this result to the sampled value for (the integral of ). All these operations can be performed electronically (using analog devices) in an accurate manner. It can be easily noticed that the quantity represents the active time (on which the step input acts upon the second-order oscillating system).

On the subsequent working interval, we can consider that the input of the second-order oscillating system equals on the whole time interval . As a consequence, the integral of the generated output equals which allows a robust estimation of the amplitude of the step change for the input as where is the sampled quantity for the integral over a period for the oscillating system output (starting to work from initial null conditions on the next working interval), and is a constant value. This operation can be also performed electronically in an easy manner.

However, this result is far of being useful for practical applications. Since the differential equation of the second-order oscillating system is it results that where represents the amplitude of the received signal. So the algorithm previously presented performs the integral of the received signal on a period, without any filtering procedure. The influence of low-amplitude high-frequency alternating components of the received signal is not diminished in a significant manner, the advantages of a filtering procedure based on second-order systems being lost.

The previously presented algorithm could be accepted if the integral can be performed on an extended time interval. In this case, we can simply estimate the quantity on a time period several times greater that (the period of the oscillating system) so as requirements regarding filtering performances (for rejecting the influence of low-frequency high-amplitude components) to be fulfilled. Thus, the influence of an alternating component with time constant and angular frequency is decreased about times. The sampling procedure is not robust any more (the signal which is integrated differs to zero at the end of the integration interval, in case of step changes) but for extended intervals, the relative error generated by switching phenomena can be neglected (the switching interval is very narrow as related to the integration interval, so the integral of stochastic switching phenomena represents a low value as related to the integral performed on the whole working interval).

#### 3. Synthesis of Quick Command for Compensating Dynamic Environmental Changes

In case of electric drives, an extended integration interval for quantity (as presented in the second section of this paper) could be allowed if it corresponds not to the controlled quantity (the angular frequency of a shaft, for example) but to its derivative (supposing that the resistive torque can vary). As a consequence, can be written as ( being a constant). It results that so it corresponds to the variation of the angular frequency on the time interval which should be adjusted by a supplementary active torque represented by where corresponds to the active torque and corresponds to the resistive torque. The active and resistive torque should be equal in the stationary regime, so the supplementary torque should act on a limited time interval, so as to compensate the difference previously detected.

For this purpose, we could notice that quantity is proportional to this difference, so a supplementary active torque can be transmitted in a limited time interval as ( being a constant). When , both function and its derivative have null values, so the action of the supplementary torque can be stopped avoiding the influence of any switching phenomena (which can generate errors).

Since it results that the action of this supplementary torque on this time interval can be represented as

While is proportional to the difference (as has been shown), it finally results that

By adjusting the relation between and , we can set quantity to unity, so as the action of the step change in the environment upon the angular velocity to be compensated in a subsequent finite time interval.

This intuitive model is also valid for any complex (biological) system which should maintain its position or the velocity of certain components at a specific value.

#### 4. Efficient Signal Processing Methods Based on Two State Variables

The signal processing method presented in previous paragraph is based on sampled values for three successive working intervals of the second-order oscillating system. Considering that a step change for the input is detected on a certain working interval, the previous value for the input is determined on an initial working interval (on which the input equals a certain value ), and the final value for the input is determined on a final working interval (when the input equals a final value ). During the middle working interval, the quantity is substracted from the received signal, and the quantity is available at the end of the interval (using amplifying and sampling procedures). Finally, the step change and the time moment are determined (using substracting and dividing procedures).

However, the linearity of this second-order oscillating system allows a more efficient and robust procedure to be used. The algorithm presented in the previous paragraph requires two identical oscillating second-order systems working at the same time: one for processing the input (so as to determine the estimated values for an input considered to be constant on that interval) and another for processing the difference between the received signal and the previously sampled value (so as to detect a possible step change during this interval by determining the quantity , where stands for the step change and stands for the active time). Moreover, filtering aspects could require extended time intervals for processing the input signal (as was shown), and the absence of any control action during such an interval could allow significant changes for the output of the complex system from the desired value.

In previous paragraph, it was shown that filtering properties and robustness require an extended time interval for processing the input signal (received from transducers). For this purpose, we can use either a second-order oscillating system with a period equal to the working interval (which means ), or a second-order oscillating system with a period corresponding to a submultiple of the working interval (which means ). The last choice is far more convenient, as it will be shown in this paragraph. Let us suppose that on the whole working interval the input is represented just by the constant value . In this case, the output is represented by . Supposing that quantity was determined on previous working interval, we can determine in a very simple manner the expected values which would be sampled for when phase equals (this means , the effect of high-frequency low-amplitude fluctuations being neglected). At these time moments quantity would be zero, so the sampling procedure would be robust as related to the constant input .

This suggests the possibility of sampling at these time moments (this means when equals ) and subtracting the expected values previously presented. Due to the linearity of the second-order oscillating system, if a step change with amplitude is detected during the working interval, the sampled values at the moments of time will be represented by a sum of values determined by step input and by a step change starting to act at moment . Thus, the result is

Thus, the successive values of the subtracting procedure will correspond just to the influence of the step change .

A quick analysis for first pair of values for and after the step change starts to act (when significant differences from expected values are detected) can be performed by comparing quantities where

The phase difference between and equals (as it was shown). It results that

It can be easily noticed that the ratio can be written as

A quick comparison of and is useful for an approximation of . For high values for and low values for , can be approximated as , and can be approximated as . For low values for and high values for , can be approximated as 0 and can be approximated as . For similar values for and , can be approximated as , and can be approximated as . Using , can be determined as where stands for the first time moment in the set determined by , for which significant differences from expected values are detected for sampled values . This algorithm is efficient because it allows a preliminary command to be transmitted to the controlled system before the working interval (on which the received signal is integrated and sampled) to come to an end. It is more effective if the working interval includes several periods of the second-order oscillating system (thus, more moments for estimative sampling are established inside the working interval, and the preliminary command is transmitted faster).

#### 5. Aspects Connected with the Uncertainty Principle

In the third section of this paper has been presented an algorithm for a preliminary estimation of a step change in the environment based on two state variables (two successive values for a function sampled at two successive time moments) with the phase difference corresponding to these time moments determined as

It can be noticed that measurements performed for sinusoidal functions at time moments when the phase difference equals imply some aspects regarding opposite requirements for sampling moments of time. Let us suppose that we are sampling just a sine function for phase

It results that sampled values correspond to the pair

This pair can be considered also as a value for a sine function and a value for its derivative for the same phase . By taking into account the fact that any sampling procedure requires a nonzero time interval, it results that it is desirable for these functions to be almost constant on a very short time interval necessary for this sampling procedure. Mathematically, this would imply the necessity of both functions to have null values at sampling moments of time. However, it is well known that a maximum/minimum value for sine function (when the differential of sine function equals zero) corresponds to a null value for its derivative, the cosine function (at these time moments, the differential of the cosine function has the greatest value for its slope, considered as modulus). Quite similar, a maximum/minimum value for cosine function (when the differential of cosine function equals zero) corresponds to a null value for its derivative, the function (at these time moments, the differential of the sine function has the greatest value for its slope, considered as modulus). So it is impossible to select a suitable phase so as the measurement accuracy for both sine and cosine functions to be the best possible if nonzero sampling intervals are taken into account.

This aspect is similar to the uncertainty principle in physics, when measurements corresponding to a certain physical variable and to its conjugated variable (usually corresponding to a derivative of a function in respect to the previous variable) are involved. Moreover, it should be noticed that for a cosine function the action of the operator would correspond to It can be noticed that another cosine function has been generated.

Quite similar, the action of the operator would correspond to

It can be noticed that the cosine function has been annihilated. This aspect is similar to creation/annihilation of particles in advanced quantum mechanics, where such operators are derived using the decomposition of certain fields in plane waves. However, aspects connected with momentum and position operators are hard to be noticed at this stage of research.

As a conclusion for memory-limited complex systems, we can notice that the use of second-order oscillating systems allows ot just robust sampling procedures on extended time intervals for certain quantities corresponding to step changes in the environment, but also the use of just two state variables (corresponding to sampled values for sine and cosine function at certain moments of time) for a preliminary estimation of such step changes during the working interval. For an extended working interval which includes several alternances of the oscillations generated by the received signal, this implies the possibility of transmitting quick preliminary commands towards the actuators, a final adjustment being determined at the end of the whole working interval, based on the difference between required action and the already-performed action. Similar to aspects presented in previous section, these preliminary commands should be better transmitted as a pulse defined by for (a period of the oscillating signal considered from the second moment of time when the preliminary sampling procedure for was performed, and the pair of values sampled at and has been analyzed). Thus, subsequent expected values for and at time moments determined by could be computed in an easy manner by the signal processing device, and any differences could be analyzed (for a preliminary conclusion) in an easy manner on next periods of the oscillating system until the working interval comes to an end.

For this reason, the algorithm previously presented is suitable for memory-limited complex systems since it performs both a preliminary analysis of signal received from environment for detecting step changes (with preliminary commands transmitted towards actuators) and a final accurate estimation for the required action on next extended time intervals (computed as a difference between the whole action required by the step change and the action already performed by preliminary commands).

As in case of biological systems, this algorithm is based on values sampled at some successive moments of time. It generates a sequence of certain commands towards the environment as a sequence of pulses, analyzes the difference between expected values and real values for the signal received from the environment, and adjusts the command with higher accuracy after an extended time interval (a kind of multilevel control and command). Another important similarity between this algorithm and behaviour of biological systems should be noticed; the sampling moments of time (when the processed signal is recorded) differ to the time moments when the filtered (processed) signal has a great slope (considered as modulus) so as to allow a robust estimation using just two sampled (recorded) values. The fact that less memory is involved is essential for complex systems which have to survey a great number of parameters (motion parameters, for instance) in the environment and to check the effect of commands transmitted towards a great number of actuators, see the case of vision processing studied in [7, 8].

#### 6. Conclusions

This study has presented advanced signal processing methods and command synthesis for memory-limited complex systems. It was shown that for observing, modeling, and controlling dynamic environmental phenomena in case of memory-limited complex systems, some specific methods based on accurate measurements performed on limited time intervals are required. Starting from the necessity of a set of consecutive measurements performed in a robust manner for detecting step changes in the environment, it was shown that an extended time interval for processing the input signal is necessary. For this reason, the use of second-order oscillating systems was improved by adding a supplementary algorithm so as preliminary values for step changes in the environment to be available for control and command during the signal processing interval. A method for generating a robust command towards the control equipment on a limited time interval in order to compensate the step changes detected on previous working interval was also presented. Finally, similarities between measurements of a certain quantity and of its derivative for a sine function, by one side, and the uncertainty principle in physics (by the other side) were briefly mentioned.

#### References

- C. Toma, “An extension of the notion of observability at filtering and sampling devices,” in
*Proceedings of the International Symposium on Signals, Circuits and Systems Iasi (SCS '01)*, p. 233, Iasi, Romania, July 2001. - C. Cattani and J. J. Rushchitsky,
*Wavelet and Wave Analysis as applied to Materials with Micro or Nanostructure*, vol. 74 of*Series on Advances in Mathematics for Applied Sciences*, World Scientific, Toh Tuck Link , Singapore, 2007. - C. Cattani, “Connection coefficients of Shannon wavelets,”
*Mathematical Modelling and Analysis*, vol. 11, no. 2, pp. 117–132, 2006. View at Google Scholar · View at Zentralblatt MATH - G. Toma, “Practical test-functions generated by computer algorithms,”
*Lecture Notes Computer Science*, vol. 3482, pp. 576–584, 2005. View at Publisher · View at Google Scholar - S. Pusca, “Invariance properties of practical test-functions used for generating asymmetrical pulses,”
*Lecture Notes Computer Science*, vol. 3980, pp. 763–770, 2006. View at Publisher · View at Google Scholar - A. Toma and C. Morarescu, “Detection of short-step pulses using practical test-functions and resonance aspects,”
*Mathematical Problems in Engineering*, vol. 2008, Article ID 543457, 15 pages, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. Y. Chen and Y. F. Li, “Vision sensor planning for 3-D model acquisition,”
*IEEE Transactions on Systems, Man and Cybernetics, Part B*, vol. 35, no. 5, pp. 894–904, 2005. View at Publisher · View at Google Scholar - S. Y. Chen, Y. F. Li, and J. Zhang, “Vision processing for realtime 3-D data acquisition based on coded structured light,”
*IEEE Transactions on Image Processing*, vol. 17, no. 2, pp. 167–176, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at MathSciNet