`Mathematical Problems in EngineeringVolume 2010 (2010), Article ID 324818, 11 pageshttp://dx.doi.org/10.1155/2010/324818`
Research Article

## Specific Differential Equations for Generating Pulse Sequences

Department of Physics, Politechnica University, 313 Spl. Independentei, 060042 Bucharest, Romania

Received 17 September 2009; Accepted 6 October 2009

Copyright © 2010 Ghiocel 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 study presents nonlinear differential equations capable to generate continuous functions similar to pulse sequences. First are studied some basic properties of second-order differential equations with time-dependent coefficients generating bounded oscillating functions similar to test-functions (the function and its derivative being equal to zero at the same time moments). The necessary intercorrelations between the phase of generated oscillations and the time-dependent coefficients is presented, being shown also that the external command function should be set to a constant value at these time moments so as to determine the amplitude and the sign of generated oscillations. Then some possibilities of using previous differential equation for generating positive-definite functions with null values for the function and its derivative at the same time moments and with constant slope for its amplitude are presented, being shown that the corresponding external command function should present also alternating components. Finally all previous results are used for determining a set of second-order differential equations with time dependent coefficients and a set of external command and corrective functions for generating a pulse sequence useful for modelling time series.

#### 1. Introduction

As was presented in [1], practical test-functions are important in signal analysis due to the fact that such a function and a finite number of its derivatives are equal to zero at the beginning and the end of a certain working interval and sampling procedures can be performed in a robust manner. In the same paper were studied invariance properties of differential equations able to generate such practical test-functions on a limited time interval. However, aspects presented in that initial study should be extended to functions defined on extended time intervals, so as to investigate possibilities of generating sequences of such practical test functions. Stress should be laid upon sequence of sharp mathematical functions similar to pulses, while in the limit case pulse sequences can be considered as time series.

#### 2. Differential Equations for Generating Bounded Oscillating Functions

For a limited time interval it was shown [2] that a second-order system working at the limit of stability and an integrator can generate a function represented by

for an external constant input (the integration time interval being represented by a period of the oscillating second-order system, while ) in a robust manner. At the end of the integration, the value can be sampled in a robust manner, while both and are equal to zero. Moreover, filtering properties of such a structure are similar to those of a structure consisting of an asymptotically stable second-order system and an integrator.

Let us check the possibility of obtaining a more general positive-definite oscillating function as a result of a differential equation. We consider the function written as

(its average value is equal to unity on an unlimited time interval). This implies

We are looking for a second-order differential equation under the form

(similar to linear differential equations with time-invariant coefficients generating damped oscillations, it is considered that and are equal to zero when equals its average value, this means when ).

By substituting in the previous equation , respectively, , with , respectively, it results

For a robust solution we should equate functions multiplying , respectively on both sides. It results

By considering it results

Thus the differential equation generating the output (and its derivative ) can be written as

or

where equals . Since is a positive-definite function, at time moments when both and vanishes (implying , ) and the previous equation becomes

so for at this time moments the second derivative and thus function can be approximated as zero on a large time interval.

Next step consists in studying the possibilities of generating a function represented by

It results that

We consider a slightly changed differential equation—the term in previous differential equations being substituted by , where is an external command function (it can be noticed that the amplitude of the external command should be proportional to the amplitude of function ). Thus, the differential equation becomes

Substituting the mathematical expressions of and , previously obtained in this modified differential equation, it results

By equating the terms which do not depend on , and the functions multiplying on both sides of this equation, it results

(i)(ii)

which implies

(by substituting from previous equation), and

then

(similar to the analysis of a possible function presented at the beginning of this paragraph), and

(iii)

Since , it results and by substituting with this mathematical expression in last equation it results

(the right-hand side of the equation was reduced by dividing both the numerator and the denominator with ) and

with possible solutions: ( is a constant function) and/or ( is a constant function). However, we should notice that cannot be a constant function, while function cannot be equal to zero (this would correspond to a differential equation which does not depend on the external command ).

Substituting in equation previously obtained (see (i)) it results

(the right-hand side of the equation was reduced by dividing both the numerator and the denominator with ) and

which can be written also as

Substituting (according to the previous conclusion; see (iii)) it results

with possible solutions and/or .

By studying possible solutions for the equations presented at (i), (ii) and (iii), it results that (), and ( and/or ) and ( and/or ). This implies that (as in the previous case), and , so should be a constant function.

A more general second-order differential equation generating a bounded positive-definite oscillating function is represented by

(it can be noticed that the external command denoted by could be different from the function corresponding to the amplitude of ). The analysis is similar to the one performed in the previous case (when the external command and the amplitude of the generated function were considered to be equal). In the same manner, and are computed and substituted in the second-order differential equation. By equating the functions multiplying and on both sides of this mathematical expression, it results the same equations as in the previous case (see (ii), (iii)), with the same consequences: (by equating functions multiplying ), and (by equating functions multiplying ).

By equating the terms which do not depend on on both sides of this mathematical expression (for a robust solution) it results

(which differs to equation obtained in the previous case; see (i)).

Since it results . Substituting , in previous equation it results

which implies

(it was shown that cannot be equal to zero). So for a robust solution it is necessary for the external command to be equal to the amplitude of , both being constant functions (, so ).

It is also possible to equate just the functions multiplying on both sides of the mathematical expression (this implies and, consequently, , as has been shown). Performing a reduction of similar terms on both sides of equation, it results

Substituting ,   it results

The first fraction on right side has been obtained by dividing both numerator and denominator by (a nonzero expression, since , as was shown) and the function multiplying on right side has been obtained by dividing both numerator and denominator by (a nonzero value, since as was shown). It results

so the function (the external command) can be written as

This result can be used for determining the external command function if the function (the amplitude of the function to be generated) is known. However, it can be noticed that contains an oscillating part , and so the differential equation is far of being robust.

It can be noticed that the function has an important property: the function and its derivative are equal to zero at the same time moments , when and (supposing to be a constant function, as has been shown). Inspecting the second-order differential equation (where and is a constant function—as has been shown) it results that at this time moments the differential equation can be written as

So at this time moments (when the state variables and of the second-order differential equations are equal to zero) the second derivative is determined by the external command . This suggests that at this time moments a new function:

can be generated by a differential equation:

where the time moment has become the origin of time for the new working interval, with the external command function represented by the constant function . The amplitude and the sign of the external command at the time moments when and are equal to zero determine the amplitude and sign of the function on next working time interval (until both become zero again). The phase of the function and the function which determines the differential equations on each working interval are intercorrelated by , as has been shown.

It can be noticed that the previous differential equations can be easily implemented (using either analogue circuits if a high working frequency is necessary or digital circuits if more accurate results are required). In case of analogue circuits, mathematical operations as multiplying and/or dividing functions are available due to high performances of operational amplifiers. The errors between the required function and the function generated by analogue circuits cannot be avoided; however, the cumulative errors can be easily set to zero at the end of each working interval , by simply setting to zero the functions at these time moments (when they should present a null value, according to theory). By performing this operation before applying a new external command function , a robust implementation of the mathematical model can be acquired.

The same aspect is valid also if digital circuits are used. A higher accuracy can be obtained, but numerical errors cannot be avoided and thus numerical values of and should be set to zero at the end of each working interval (as in case of analogue circuits).

#### 3. Functions Similar to Pulse Sequences Generated by Second-Order Differential Equations

The mathematical expression for the second-order equation presented in previous paragraph should be used for generating functions similar to pulse sequences or time series. For this purpose, at first step it will be presented a differential equation able to generate a function with a constant slope for its amplitude, this means a function which can be written as

(the slope of the amplitude has been considered to be equal to unity). It results

By substituting the previous expression for into the left side of the differential equation (where , as was shown), by substituting the previous expressions for and on the right side of the differential equation and by considering that the external command function can be represented by a sum:

it results

By equating the terms which do not depend on , and the functions multiplying on both sides of this equation, it results

(i)(ii)(iii)

Since , (as was shown) from previous three equations it results that

So the external command function necessary for generating a function by a differential second-order equation is represented by

However, pulse sequences are best represented by functions as

which is similar to a constant time series for (the functions differs to zero just when ).

For , it results

By substituting in the second term of right-hand side, it results

The first term on right-hand side can be written as

(where ).

The second term on right-hand side can be written as

where , , .

The third term on right side can be written under the form

where and as was shown.

So the second derivative can be written under the form

It can be noticed that the function can be generated by a second-order differential equation similar to the one presented in previous paragraph. The time-dependent coefficient (multiplying ) equals the time-dependent coefficient of the differential equation presented in the previous paragraph for any ; the time-dependent coefficient (multiplying the difference between the external command and the function to be generated) is represented by time-dependent coefficient of the differential equation presented in the previous paragraph divided by .

Two important aspects should be emphasized.

(i)The first term of the external command function is represented by (since ).(ii)The second term of the external command function is represented by a corrective function which can be written as: ; so, for the function and thus the external command function can be approximated by .

This suggests that a set of differential equations

can be used for generating the complete set of functions starting from . First and are used for generating , then and are used for generating and so on.

Such a set of differential equations can be also implemented using analogue or digital circuits for multiplying and/or dividing functions. The errors between the required set of functions and the functions generated by analogue or digital circuits become larger as increases, because the number of functions which should be previously generated also increases (it is necessary to generate functions corresponding to , , before generating a function corresponding to a certain and the errors are cumulative). For this reason, the operation of setting to zero all generated functions should be done more often and the working interval should be shortened (it can include just a few alternating components of generated signal).

#### 4. Conclusions

This study has presented nonlinear differential equations capable to generate continuous mathematical functions similar to pulse sequences. First were studied some basic properties of second-order differential equations with time-dependent coefficients , generating bounded oscillating functions similar to test-functions (the function and its derivative being equal to zero at the same time moments) and similar to wavelets also [3]. The necessary intercorrelations between the phase of generated oscillations and the time-dependent coefficients have been obtained as , , being shown also that the external command function should be set to a constant value at these time moments so as to determine the amplitude and sign of generated oscillations. Then some possibilities of using previous differential equation for generating positive-definite functions with null values for the function and its derivative at the same time moments and with constant slope for its amplitude are presented (similar to solitary waves, see also [4]), being shown that the corresponding external command function should present also alternating components as . Finally all previous results are used for determining a set of second-order differential equations with time dependent coefficients and a set of external command and corrective functions for generating a pulse sequence useful for modelling time series (similar to time series application to surface analysis [5] or to fracture phenomena [6]) or for computing autocorrelation function [7] with extension to network traffic [8]). The study is also suitable for extending wavelets dynamical aspects (see [9]) to correlation-based computational methods (see [10]).

#### References

1. G. Toma, “Practical test-functions generated by computer algorithms,” in Proceedings of the International Conference on Computational Science and Its Applications (ICCSA '05), vol. 3482 of Lecture Notes in Computer Science, pp. 576–585, Singapore, May 2005.
2. C. I. Toma, “Filtering possibilities based on oscillating systems for optoelectronic signals,” in ROMOPTO 2000: Sixth Conference on Optics, vol. 4430 of Proceedings of SPIE, pp. 842–845, Bucharest, Romania, September 2000.
3. C. Cattani, “Harmonic wavelets towards the solution of nonlinear PDE,” Computers & Mathematics with Applications, vol. 50, no. 8-9, pp. 1191–1210, 2005.
4. J. J. Rushchitsky, C. Cattani, and E. V. Terletskaya, “Wavelet analysis of the evolution of a solitary wave in a composite material,” International Applied Mechanics, vol. 40, no. 3, pp. 311–318, 2004.
5. V. P. Paun, “Fractal surface analysis of Zircaloy-4 SEM micrographs using the time-series method,” Central European Journal of Physics, vol. 7, no. 2, pp. 264–269, 2009.
6. M. Olteanu, V.-P. Pǎun, and M. Tǎnase, “Fractal analysis of Zircaloy-4 fracture surface,” Revista de Chimie, vol. 56, no. 1, pp. 96–100, 2005.
7. M. Li, “Modeling autocorrelation functions of long-range dependent teletraffic series based on optimal approximation in Hilbert space—a further study,” Applied Mathematical Modelling, vol. 31, no. 3, pp. 625–631, 2007.
8. M. Li and S. C. Lim, “Modeling network traffic using generalized Cauchy process,” Physica A, vol. 387, no. 3, pp. 2584–2594, 2008.
9. C. Cattani, “Wavelet analysis of dynamical systems,” Electronics and Communications, vol. 17, pp. 115–124, 2002.
10. M. Li and C.-H. Chi, “A correlation-based computational model for synthesizing long-range dependent data,” Journal of the Franklin Institute, vol. 340, no. 6-7, pp. 503–514, 2003.