V. N. Tibabishev
Private laboratory identification of dynamic systems, office 151/36, Lugansk, 91004, Ukraine
Correspondence to: V. N. Tibabishev , Private laboratory identification of dynamic systems, office 151/36, Lugansk, 91004, Ukraine.
Email: | ![](email/10.5923.j.am.20120205.01.gif) |
Copyright © 2012 Scientific & Academic Publishing. All Rights Reserved.
Abstract
It is shown that a well-known theory of random stationary processes contain contradictions. Integral representations of correlation functions and random stationary processes are investigated further. The new method of struggle with handicaps is received on the basis of the carried out researches. Method of dealing with noise leads to a new method of identification of dynamic characteristics of control objects in a class of multidimensional linear stationary models. As an example, describe an algorithm for obtaining the differential equation-wire feed control aircraft pitch, taking into account the elastic deformation of the structure of class IL -96 aircraft.
Keywords:
Stationary Random Processes, Almost Periodic Functions, Hilbert Spaces
Cite this paper:
V. N. Tibabishev , "Suppression of Random Noise by the Separation of Frequencies", Applied Mathematics, Vol. 2 No. 5, 2012, pp. 146-156. doi: 10.5923/j.am.20120205.01.
1. Introduction
The initial (primary) data are given in the presence of noise in the solution of applied problems ever. Certain properties of the signals used to control noise. The modern theory of random wide-sense stationary processes with finite variance (the random process) is described in many books, for example, in[1- 4]. In the theory of random processes, the empirical[5] and theoretical methods, for example, published in[2],[6 -7]. The random stationary process and its correlation function are represented by a trigonometric Fourier series of arbitrary frequencies, and therefore belong to the Hilbert space of almost periodic functions of Bohr. This follows from the empirical studies of E. Slutsky, which are summarized in[4]. A stochastic process does not have the ergodic property with a discrete spectrum[2].The first theoretical concept was published by A. Khinchin, according to which the correlation function of random process has two components[6]. One component belongs to the Hilbert space of almost periodic function
. Average value of the time of the square the other component is always zero. This implies that another component belongs to the Hilbert space of square-integrable Lebesgue functions
.Another theoretical concept was published by Wiener, according to which the correlation function of a stationary process belongs only to the Hilbert space of square integrable functions, Lebesgue integral is the Fourier-Plancherel terms of the spectral density[7]. Such a model is a random process has the ergodic property of the theorem, published in[2].In the works of Norbert Wiener and Khintchine representation of random processes they were not considered. This problem was solved empirically first E. Slutsky, then in theory by Kolmogorov[8], further Cramer[2] and other researchers. In all views, except of E. Slutsky, a stationary random process is an integral of the Fourier-Stieltjes.All three concepts of a random stationary process contradict each other. In Wiener[9] shows that if the process belongs to the Hilbert space
, then the correlation function, obtained by averaging over time is identically zero on the entire line. This assertion N. Wiener is contradiction ergodic theorem given in[2].D. Middleton[3] and other researchers[10] pointed out that the integral of the square of a random stationary process for all direct costs, even when the average value of zero. Therefore, a stationary random process can not belong to a Hilbert space
.Known[11], that the space
is a set of functions that decay sufficiently fast at infinity. Such functions are not suitable for description of random stationary processes on the whole line. In many monographs as examples of stationary processes are functions that do not belong to the space
, a elements of almost periodic functions
[12],[2 - 4]. It is still unclear why the Lebesgue-Fourier-Stieltjes in N. Wiener[7] has led to conflicting results.A hypothesis, that the contradictions are due to incorrect assumptions adopted in theoretical studies. First, all theoretical studies suggest that the angular frequency ω is a continuous variable in the integration[8] and differentiation[7]. Second, the probability distribution function contains a step and continuous components[7],[2],[6]. In this connection there arose the first task of assessing the correctness of these assumptions.
2. Investigation of the Known Assumptions
First, we show that a certain assumption about the continuity of the independent variable angular frequency is incorrect assumption. Angular frequency ω is related to the cyclic frequency f numerical function
. Cyclic frequency belongs to a continuous set of real numbers
. The answer to the question, to what set of numbers belongs to the set of corner frequencies, the following theorem.It is theorem 1. Let the real function
defined on the set of real numbers
for all
. If the ratio
is an element of a subset of irrational numbers, transcendental
, for example,
the numerical function that displays the specified set of real numbers
, only a subset of irrational numbers, transcendental
, where
- a subset of rational numbers
- a subset of algebraic irrationals, for example,
where the union
- is an infinite countable set of algebraic numbers.It is proof. Known[13], that there exists
an inverse
for each element such that
. Multiply the left and right sides of a numerical function to return the item to an arbitrarily chosen
. Obtain
. It follows that for
an arbitrary product
of an irrational transcendental number. Of all the possible random numbers
choose an arbitrary rational number
, where m and n are arbitrary nonzero integers. Of all the random numbers m , choose a number
. In this case, we find that
. Theorem 1 is proof.The question arises, under what factors
numerical function
defined on the real line
for all
will display the specified set of real numbers into itself? The answer to this question is the following theorem.It is theorem 2. Let the real function
, defined on the set of real numbers
for all
. If the coefficient
belongs to a subset of rational numbers
, then the numerical function displays the specified set of real numbers
into itself
.It is proof. All rational coefficients
are presented as a fraction
, where
and
are arbitrary nonzero integers. Let
all of the random numbers. In this case
It is consequence of theorem 2. Angle ω and cyclic frequency f are related by the formula
. It is according with theorem 2 the set of numbers
. Since
it belongs to the set of irrational numbers, transcendental, then by theorem 1 the set of corner frequencies
belong to the subset of irrational transcendental number, which is not continuous set of numbers. In an infinite uncountable interval only irrational transcendental numbers contains an infinite countable number of seats are not engaged in irrational, transcendental numbers. On an infinite line of real numbers
, these seats are a countable infinite set of algebraic numbers[13]. It follows that an infinite range of irrational transcendental frequency contains an infinite countable set of discontinuity points. On an infinite line of real numbers, these seats are a countable infinite set of algebraic numbers[13]. It follows that an infinite range of irrational transcendental frequency contains an infinite countable set of discontinuity points. On a subset of irrational numbers, transcendental notions
of mathematical analysis as an infinitely small quantity
, the derivative (spectral density) and the differential of the independent variable angular frequency
does not exist. With respect to the integral sums of Riemann, Riemann-Stieltjes and Fourier and Fourier-Plancherel not exist.Second, the well-known assumption that the distribution function of angular frequency contains a continuous component is an invalid assumption. It is theorem 3. In general, the probability distribution as a function of real random variables consists of three components - a step function, absolutely continuous differentiable function and a singular function. If the probability distribution is a function of random angular frequency, defined on the set of irrational numbers, transcendental, then the probability distribution function of random angular frequency contains only a step function with an infinite countable number of discontinuity points and contains no nonzero continuous components.It is proof. It is known[2], that in general the distribution function of random real variables
consists of three components
where
- a step function of discrete values
,
, continuous variable
, and
can also take an infinite value,
- absolutely continuous differentiable function on a continuous variable
,
- a singular function - a continuous function, which has a derivative almost everywhere equal to zero, for example,
where
- arbitrary constant.From the properties of the probability distribution function defined, for example, random variables
, it follows that the entire function
and its two components
and
. At the same time
. It follows that
. The distribution function is represented as
. From the other properties of the probability distribution of the variable
that
. Hence![](image/10.5923.j.am.20120205.01_079.gif) | (1) |
It is known[14], that the step function can have a finite number of steps, for example,
or an infinite number of steps, but necessarily countable. Therefore, the component
belongs to a countable set of rational numbers and can take values in the interval[0,1]. Let
where
and
- integers with
.Into an infinite number of steps tends to its maximum possible value of the countable
. In this case, a discrete component tends to its maximum value
and the continuous component tends to its minimum value
. Since when the number
of breaks in the step function equal to an infinite countable set of values, for all
the distribution function continuous component
. Theorem 3 is proved.Thus, the probability distribution function of angular frequency has only a discrete component, containing an infinite but countable set of stairs. In this case, it contains no nonzero continuous components.
3. Representation of Correlation Functions and Processes
In this section we repeat the well-known findings of Wiener[7] to represent the correlation functions and the known output representation for a random process and Cramer[2] with the difference that we use is not anticipated by the authors properties for the angular frequency and its distribution function, as proved above properties.Corollary 1 is of theorem 3. If the distribution function of the random variable contains a step function with an infinite countable number of steps, then the Lebesgue measure, generated by this distribution function contains only discrete measure. Therefore, the components of the Lebesgue-Stieltjes measure on a continuous non-existent.Since the angular frequency ω is not a continuous variable, the function
of the variable
is not continuous. Fourier-Stieltjes integral is defined only for continuous functions[14]. Since this condition is not fulfilled, instead of the Riemann-Stieltjes integral must be taken Lebesgue-Stieltjes integral[11], as did Norbert Wiener![](image/10.5923.j.am.20120205.01_094.gif) | (2) |
where
- the measure generated by the distribution of the random corner frequency
on a subset of transcendental numbers ![](image/10.5923.j.am.20120205.01_097.gif)
.Known[11], that in cases with a discrete measure
of the Lebesgue-Stieltjes integral (6) leads to no integral concepts defined in the space of almost periodic functions![](image/10.5923.j.am.20120205.01_100.gif) | (3) |
where
- at a frequency jump
function
.Submission of a random process differs from the representation of the correlation function so that the integral representation (2) instead of the distribution function
uses a random complex continuous with bounded variation with zero mean and orthogonal increments spectral process
[2]![](image/10.5923.j.am.20120205.01_106.gif) | (4) |
It is corollary 2 of Theorem 3. It is known[14], which the function of bounded variation in a certain way is the sum of functions of the jumps and the amount of continuous distribution functions. Therefore Theorem 3 is valid not only for the non-decreasing function of probability apply, but also for functions with limited modifications.It is corollary 3 of Theorem 3. Since any function of bounded variation is the difference between two step and the difference of two continuous distribution functions[14], the function of bounded variation of the angular frequency
, which does not contain continuous components, can generate only a discrete measure
.Function
of the variable ω contains an infinite countable set of discontinuity points. Therefore, for integration into the mapping (4) must be taken Lebesgue-Stieltjes integral. By corollary 3, Embedding theorems 3 Spectral process
, which contains an infinite countable set of discontinuity points, can produce only a discrete measure
. Known[11] that at least a discrete integral representation (4) becomes complex Fourier series, defined in the Hilbert space of almost periodic functions of Bohr ![](image/10.5923.j.am.20120205.01_112.gif)
![](image/10.5923.j.am.20120205.01_113.gif) | (5) |
where
- the point of tears.Our representation for correlation functions (3) and for stationary random processes (5) are differ from the known representations received by Wiener and G. Kramer. First, a random process and its correlation function do not belong to Hilbert space
and the Hilbert space of almost periodic functions
. In this space integral of the square of the random process does not diverge. Second, a stationary random process and its correlation function has only discrete spectrum. In this theoretical study are consistent with empirical studies E. Slutsky.
4. Functional Space for Random Stationary Processes
We shall define a functional space for random stationary processes with discrete spectra on the example of the object control, containing the
inputs and
outputs. It is believed that a multi-dimensional control object allows
an approximate description of the dynamic performance of control channels in a class of linear stationary models. Original raw data presented in the form of single realizations of random nonergodic processes only with discrete spectra, which simultaneously observed (recorded) on the inputs
and output
of control object, where
- number of implementation ![](image/10.5923.j.am.20120205.01_123.gif)
- the number of input
- output number. Instead of accurate baseline data
and
can be observed or synchronous writes only approximate data, ![](image/10.5923.j.am.20120205.01_128.gif)
and ![](image/10.5923.j.am.20120205.01_130.gif)
in which exact components distorted by random noise or additive
and
.Each implementation of approximate initial data belongs to the Hilbert space of almost periodic functions. There are various generalizations of almost periodic functions[15]. We will consider such control objects, whose dynamics obeys the laws of classical mechanics, in particular, second order differential equation of Newton, for example, in the form of
,
- the second derivative, for example, the output signal ,
- force. According to the known existence theorem[16] function
must be continuous. In many cases[16] the function
is the Fourier series. Since the functions
are continuous functions, the Fourier series for these functions are convergent series in the known space[17]. In the Hilbert space of almost periodic functions in the sense of Besicovitch each element appears to be convergent Fourier series[15]. Therefore, we assume that each component of the approximate initial data is elements of the Hilbert space of almost periodic functions in the sense of Besicovitch
.We consider the ergodic properties of random processes in this space. Let the stationary random process
with bounded non-zero expectation
elongs to the Hilbert space of almost periodic functions in the sense of Besicovitch
. It is known[12], that in this space, the notion of average
where i number of implementation.Known ergodic theorem Birkhoff - Khinchin proved that for non-Hilbert space of integrable functions on Riemann[2], that is, for the space
. This theorem is generalized to the Hilbert space of almost periodic functions in the sense of Besicovitch.It is theorem 4. Suppose that, in general, a stationary random process belongs to the set of random processes
with limited non-zero expectation
. If a random stationary process is an element of the Hilbert space of almost periodic functions in the sense of Besicovitch, then with probability 1, the average value of any implementation of this set
equal to the mean stationary random process.It is proof. Every realization of a random process from a variety of
can be represented as
, where
- centered stationary random process. Each cantered stationary random process that belongs to the Hilbert space of almost periodic functions in the sense of Besicovitch is convergent trigonometric series[4]![](image/10.5923.j.am.20120205.01_152.gif) | (6) |
where
. It is
since with probability 1, the average value of a centered random process is zero. It follows that a stationary random process with discrete spectrum has the ergodicity of the first order[1]. Theorem 4 is proved.We show that the second and mixed moments, e.g., centred random stationary processes belonging to the Hilbert space of almost periodic functions in the sense of Besicovitch, Ergodic properties of second order do not possess. Let
, where
- the set of almost periodic functions by generating a correlation function of a certain species[4]![](image/10.5923.j.am.20120205.01_157.gif) | (7) |
where
. The operators display the space
into itself[12], because it is completely continuous operator of convolution type[12]. Let the operator equation of convolution type
kernel is generated by a function of
. In this case, we find that![](image/10.5923.j.am.20120205.01_162.gif) | (8) |
This shows that the stationary processes that belong to the Hilbert space of almost periodic functions in the sense of Besicovitch, do not possess the second-order ergodicity
.The exact component of the input signal generates a forced movement control channels. Noise, distorting the exact component of the input signal, no effect on the forced movement of the control channels. The output signals induced motion is distorted by other noise not associated with precise input to any obstacle that distorts the exact input. Interference and accurate random processes are of different nature. If the components are the exact image and inverse image of the observed signals of linear time-invariant operators, then the interference is linearly independent and therefore uncorrelated processes. In the theory of random processes produce correlated and uncorrelated random processes. Uncorrelated random processes in the sense of the conditions
,
,
,
,
, for all
and
where
- the symbol of averaging over infinite-dimensional set pieces implementations. If, for example, the control channel between the input and output q p is the control channel, the exact components have a nonzero cross-correlation
. In the theory of random processes, these differences in the correlations used to combat interference. To check that these conditions are necessary baseline data in the form of sets containing an infinite number of realizations. It is practically impossible.We consider the conditions under which, firstly, the primary source data given in the form only isolated implementations and, secondly, the single realization of random processes do not possess the ergodic property. Receipt of secondary source data in the form of deterministic correlation functions in such circumstances is impossible. Therefore, we use the condition of linear independence of the individual realizations, under which condition is uncorrelated to be checked on a set of implementations, is always satisfied.Consider the equation for the submission of a centered random process (6). In this expression, the amplitude of random
and
deterministic harmonic oscillations with frequencies
generates lots of random processes. From this expression it follows that all of a given set of random processes is conceived as a whole because it is the linear hull spanned by a deterministic basis of harmonic functions whose frequencies coincide with frequencies of the harmonics of the correlation function
. The set of frequencies of harmonic components of the observed signals is a deterministic secondary source data. The main difference between the proposed methods for dealing with noise from the known method is that as secondary information is not used correlation functions, and set the frequency of the harmonic components. It turned out that the properties of systems of sets of frequency harmonic components of the observed signals depend on the correlation (linear dependence) constitute the observed signals.
5. Properties of Independent Random Processes
Distinguish multidimensional control objects with linearly dependent (correlated) and linearly independent (uncorrelated) input actions. First, consider the control objects with linearly independent input actions. You can specify a number of linearly independent systems of sets of random processes. In particular, a system of linearly independent sets of exact input actions
where
- exact sets of random input processes (impacts). Elements
are linearly independent if the equality
follows from the equality
for all
.System of linearly independent sets of exact input signals and additive noise, distorting the accurate baseline data
, where
,
. Each set of a lot of thought as a whole due to the fact that
averaging over the set of products centered random processes, such as
Generates a deterministic correlation function In systems of sets, for example,
introduce a scalar product for complex centered random functions
and
, for example,
where
and
- complex conjugate functions.It is known[12] that uncountable harmonic basis has a non-separable complex Hilbert space of almost periodic functions
. In this space, each vector is only a countable sum of nonzero orthogonal projections. Therefore, we assume that, for example,
is a separable subspace of complex functions. In this regard, given the above mathematical expression is translated into the complex plane.By hypothesis, any two simultaneously recorded sale
and
must be linearly independent for all
. It is known[11], if the vectors
and
are orthogonal, then they are linearly independent. This implies another condition of linear independence of random realizations of random processes
, where
- complex conjugate functions.From formula (6) shows that for each i th realization of a random process, belonging to a particular set of random processes
, exists regardless of the number of realization of certain deterministic set of frequencies of harmonic components
in the representation of random processes. Relationship between the sets of realizations
and the system of multiple frequencies
, generating a basis of subspace establishes
the following lemma 1.It is lemma 1. Given a finite or countable system of sets, for example,
random stationary in the broad sense of nonergodic processes with discrete spectra
where
number of implementation,
set number. Each set
system of sets
generates a correlation function of the general form
where
,
for everyone q if the system sets
is the union of linearly independent sets of random processes, the system sets the angular frequency
generating orthonormal trigonometric basis for each set of random processes
is a semiring.It is proof. Rewrite the orthogonality condition for the realization of random processes in complex form![](image/10.5923.j.am.20120205.01_220.gif) | (9) |
where
- randomly chosen realization of a random process in a complex form,
- other arbitrarily chosen the complex conjugate of realization of a random process
,
.Condition (9) holds for all values of products of coefficients
, if for all
and all
. This condition is always satisfied if the system of sets of angular frequency
is a semiring, which suppress the pair wise for all
is empty is empty
[11]. Lemma 1 is proved.
6. Properties of Linearly Dependent Random Processes
In the multidimensional control objects on the arbitrarily chosen output synchronously with the input processes
there is an output random process![](image/10.5923.j.am.20120205.01_235.gif) | (10) |
where the operators
are generated by precise linearly independent random components of the input signals
.From equation (10) we find randomly selected dedicated channel controls for which the image and preimage associated operator equation
. Property for a set of orthogonal frequency of the harmonic basis for the dedicated control channel is defined by lemma 2.Lemma 2. If a linear stationary operator of convolution type
, defined in the Hilbert space of almost periodic functions in the sense of Besicovitch
, generated by a stationary nonergodic process with a discrete spectrum ![](image/10.5923.j.am.20120205.01_241.gif)
, where
- number of implementation,
- set number, shows the weight function
, where
on the set of processes ![](image/10.5923.j.am.20120205.01_247.gif)
, the n-dimensional countable set of frequencies
,
generated by a countable-dimensional harmonic basis in a separable subspace
and n-dimensional countable set of frequencies
,
, generated by a countable-dimensional harmonic basis in a separable subspace
are equal
.It is proof. The proof is an obvious consequence of the properties of completely continuous operators, the normal form. It is known[12] that in a Hilbert space
, the operator of convolution type
, which completely continuous normal type operator is mapping
into it.By definition, each basis vector is transferred by the operator of a normal form with a coefficient equal to the eigenvalue[17]. It follows that the orthonormal system of basis functions and transform the image of the same.Known[17] that each eigenspace corresponding nonzero eigenvalue of a completely continuous symmetric operator
is finite. The angular frequency
has n -dimensional size. Lemma 2 is proved.
7. System with Independent Input Actions
It is theorem 5. Given mathematical model of a multidimensional linear stationary control object that contains multiple inputs and multiple outputs. Each output is additively associated with each input through the control channel. The control channel is the operator of convolution type, defined in the Hilbert space of almost periodic functions
accurate input random process and the weighting function given by the number input and output number. Each input accurate and precise output processes are distorted by various additive noise and form approximate the original data. The exact components of the process and distort their noise is uncorrelated between a stationary nonergodic random processes with discrete spectra, belonging to the Hilbert space of almost periodic functions in the sense of Besicovitch
. If the system sets all of the approximate input signals is linearly independent random processes, then all sets of projections on harmonic bases approximate the input and output processes can be approached to allocate the exact components of the projection of input actions and projection of the exact components of the output processes for each control channel.It is proof. Arbitrarily select the input q on which there is an exact realization of a random process, distorted by an additive uncorrelated noise
According to Lemma 1 the set of frequencies of harmonic components (projections) in the representation of this random process
is a semiring
.Each output is connected to the input signals according to the formula (10) in high-dimensional objects of control.The observed random process
, where
,
is an element of a system of linearly independent sets. By Lemmas 1 and 2, the system sets the frequency of the harmonic components (projections) of a random process observed at the output,
is a semiring
for all
,
for all
.We find the source data as the intersection of the sets of frequencies
approximate the original data, the observed input
and
output,
where
. The resulting intersection of the frequencies coincides with the set frequency, which determines the harmonic basis for accurate component of the input to the
-th input. In the limit of infinite observation time, you can find the projection ![](image/10.5923.j.am.20120205.01_281.gif)
exact spectra of the first implementation of the input
at the inlet
and the exact components of the projection of the second implementation ![](image/10.5923.j.am.20120205.01_285.gif)
, observed at the output
, the harmonic basis, determined by a set of frequencies
accurate projections on the background of additive noise
and
. Theorem 5 is proved.Dedicated accurate projection of input and output components of the processes on the set of input data can be used to address a number of applications in the presence of interference. If we solve the problem of identification of dynamic performances, given the relationship
=
, You can find the frequency transfer function of the control channel between the
-th entrance and exit to
-th
on the set of frequencies
. If we solve the problem of indirect measurements, such as the acceleration input
against the background of additive noise
, then filtered from the noise
,
- i implementation of acceleration is given by
.Apparently, the processes with discrete spectra are divided into periodic, non-periodic and almost periodic processes. If the process is periodic, there exist finite values of time n
where we have, for example, process
Periodic process is predictable. It suffices to find the minimum value of the final period T . The frequency of the first harmonic and other frequency harmonics in the Fourier series are comparable (the frequency of the first harmonic is a measure for the multiple frequencies of the other harmonics). The main difference between non-periodic processes is that the frequencies (periods), the harmonics are incommensurable. For each harmonic is necessary to determine the angular frequency, which is expressed in non-periodic infinite decimal. Instrumental and methodological errors limit the accuracy of the frequencies. Find the exact mathematical model of the observed non-periodic process is impossible. Therefore, it is unpredictable throughout the time axis and on this basis is a model of a random process.Almost periodic processes occupy an intermediate value. Frequency harmonic components can be set arbitrarily. These may include harmonic components with multiple and disparate frequencies, i.e. contain periodic and a periodic components. Due to the non-periodic component of almost periodic processes are unpredictable, as well as non-periodic processes.There are various methods for determining the approximate frequency of the harmonic orthogonal basis of the observed random processes, such as input signals
. Among them are well known method of determining the amplitude of the current spectrum[18]. Seeking an array of frequencies is determined by the frequency at which the current amplitude spectrum has a local maximum.
8. System with Correlated Input Actions
Described the frequency control method of noise is applied only in cases where all input actions are linearly independent. Moreover, for all inputs must be executed Conditions
at
,![](image/10.5923.j.am.20120205.01_307.gif)
. If at least one pair of different inputs of this condition is not satisfied, then it indicates that the test object is an object of control with correlated input actions.Degree of connection between the exact components of the input signals in real form observed on different inputs, for example,
and
can be estimated by the cosine of the angle between them
[11]. In this case, there are three cases. If
, that the signals observed at the entrances
and
, coincide up to sign. If after determining the proper dimension of the object controls the number of entries it appears that all remaining inputs, the condition
, this indicates that the condition of theorem 5 holds.In real terms, as a rule, for all inputs or a subset of inputs, the condition
for
, which indicates the correlation of the input signals. In this case, the above method of dealing with noise is not applicable. The described method of interference mitigation can be applied if, instead of the original system with correlated effects using the subsystem with uncorrelated input actions. Based on theorems 6 a subsystem can be obtained.It is theorem 6. If the system is real accurate component of the input signal is a system that satisfies
when
for all inputs or a subset of input signals, and for another subset of the input signals
, then the original system input signals that are distorted correlated and uncorrelated noise, we can distinguish linearly independent subsystem of input signals, distorted only by additive uncorrelated noise.It is proof. Change the previously used model inputs so that would condition
holds for all![](image/10.5923.j.am.20120205.01_322.gif)
. Such a system of sets of input actions denoted as
Where![](image/10.5923.j.am.20120205.01_325.gif)
![](image/10.5923.j.am.20120205.01_326.gif)
- the exact components of the input signal and the random additive noise
at the input ofq,
- i- th realization of the random correlated noise in the form of a stationary random function of communication, which is also observed at different entrances q and p.All the components of each realization of random processes are linearly independent (uncorrelated) with each other. Therefore, only the stationary random processes that occur at different inputs and generate non-trivial cross-correlation function
.The set of frequencies of harmonic components of random processes at each input
,
, for all
. The set of frequencies of harmonic components of the function of communication
observed between the inputs q and p are by definition the intersection of sets
.Obtain the system sets the corner frequency of the harmonic components of all the functions of communication,
for all
. For an arbitrarily selected input q will find the difference[11] systems of sets of angular frequency![](image/10.5923.j.am.20120205.01_340.gif) | (11) |
Performing the operation (11) for all
obtain a system of sets of frequencies
which is a semicircle, as
for all
. Thus, applying the subtraction (11) to the original system input correlated effects
You can select from it a different system
with uncorrelated effects, which is determined by the system sets the frequency of all harmonics
exact components of linearly independent input signals
and all harmonic frequencies of additive noise
. Theorem 6 is proved.
9. An example of Solving the Problem of Identification
The proposed method of dealing with noise is used in solving the problems of identification of dynamic characteristics of the multidimensional control objects in a class of linear stationary models. In particular, the initial data, observed in the regime of one automatic landing aircraft (class IL -96), the mathematical model of a control channel in a class of linear stationary models between input
- a given pitch, at the input-wire control system for pitch, and output
- the actual angle pitch plane. On the pitch is affected by other inputs: the position of thrust lever
(t), the angular position of the flaps
and angular position of the stabilizer
Synchronous records listed the output and input signals were provided in the form of 274 discrete samples in 0, 5 sec. (ZAK_51_IL_V30). Inputs
and
signals are relay type, derivatives of which are described by finite functions on the whole line in the form of rectangular pulses. As discussed above (Section 4), the proposed method is applied to random processes, which have a continuous second derivative. For signal relay type, this condition is not satisfied, making it difficult for the solution of problem identification through
-
and
-
. Known method for solving the problem of identification of the dynamics of multidimensional control systems reduces to the solution of integral equations of the first kind[19]. System of equations can be obtained only for the control objects for which the number of inputs equals the number of outputs. Baseline data are presented in the form of autocorrelation functions of all input signals, cross correlation functions between all the different inputs and cross correlation functions between all input and all output signals. Observed signals are stationary random fragments in the broad sense of processes that do not possess the ergodic property. This fact excludes the receipt of the correlation functions by averaging over time of individual implementations. Obtaining the correlation functions over the set of realizations is hampered for the reason that it is impossible to perform many aircraft landings in the same conditions, such as weather. In addition, all aircraft like objects have a dimension of control over all inputs greater than that of all the exits. For these reasons, the task set cannot be solved by a known method.The task of identifying the dynamic characteristics of the channel management has been solved by a new method proposed in the following sequence. First by local maxima of the current amplitude spectra of the observed signals were determined frequency harmonic components simultaneously observed signals, which are secondary source data? Received the original system sets the frequency of the harmonic components of input and output signals
where
- set the frequency of the harmonic components of input signals on the i-th entry,
- set the frequency of the harmonic components of output signal
. Theorem 5 implies that if the system sets the approximate data observed at all entry points is linearly independent random processes, then all sets of projections on harmonic bases approximate the input and output processes close, you can select the exact components of the projection of input actions and projection of the exact components of output processes. Typically, the system sets the input signal does not satisfy Theorem 5. Therefore, by subtraction of sets (11) of the resulting system sets the frequency of the harmonic components of input signals is a subsystem of the sets of frequencies only independent of the harmonic components of input signals, which are tertiary source data. Software, this operation is performed as follows. Chosen criterion under which two frequencies are coincident. If
, where
T-duration of implementation, for![](image/10.5923.j.am.20120205.01_368.gif)
dimension array of frequencies
,
- the dimension of the array of frequencies
,
then the frequencies are coincident and are removed from each source array frequency input signals. The array of frequencies
it turns out, usually with smaller dimensions in the result. For example, the original dimension of the array of frequencies at the first entrance was
and the dimension of the array of frequencies independent of harmonic components on the same first entry was
.By lemma 2, for a given control channel
is defined by a set of frequencies of harmonic components that generate and describe the forced motion in the channel controls on the criterion that the frequencies of harmonic components of input and output signals of the selected channel management
where
- the dimension of the resulting array of frequencies harmonic output
, where d- the dimension of the matching frequency harmonic components, generators, and describing the forced motion in the channel
In this example, it was found that
. A subset of frequencies is
quaternary source data.For all received frequencies
determined by the Fourier exponents of the input signal ![](image/10.5923.j.am.20120205.01_385.gif)
and the output signal![](image/10.5923.j.am.20120205.01_387.gif)
). At the lowest frequency
sets of frequencies
astatizma determines the order in the system for the control action. For what is the value of the frequency transfer function at the lowest frequency
. Depending on the position of the point
in the complex plane is determined by the order astatizma on the control action. Usually, the specified point can be located in the first, second or third quadrant of the complex plane, which corresponds to the zero, first or second order astatism
[20].Obtained by Fourier exponents and order astatism (in this example
) can make a number of systems of algebraic equations on the 2 nd to the d- th order, which are the Fourier transforms of systems of linear ordinary differential equations with unknown constant coefficients of the 2 to d - order. For example, a system of algebraic equations, given the current q-th order, where
, is represented as![](image/10.5923.j.am.20120205.01_396.gif)
where
From this system of equations should be two systems of equations for real and imaginary components. Solving these equations successively from 2 to d -th order with respect to coefficients
, we find that the maximum order of the system of equations
, at which the system of equations is consistent. To solve this problem it turned out that the dynamic characteristics of the selected channel management, taking into account the elastic strain aircraft structure can be described as a first approximation of ordinary differential equations 9-th order with constant coefficients:
= 0;
=- 2,8531;
=- 14,5651;
=- 5,4709;
=- 3,9165;
=- 1,0721;
=- 0,2584;
=- 0,06054;
=- 0.004845;
=- 0,001023. In aerodynamics there is a rule for a given character and the actual pitch of the aircraft[21]. Therefore, the coefficients have negative signs for the differential equation.In a similar way was obtained a differential equation for the other control channel, which determines the effect of engine power control on the pitch. More information about a new method for identification of dynamic characteristics is published on http://asvt51.narod.ru.The obtained differential equations in two-channel pitch control allows the system to find the control of the aircraft taking into account the elastic deformation the airframe in the space of states.
10. Conclusions
1. In the modern theory of random wide-sense stationary processes are the fundamental results of empirical research on E. Slutsky and the results of theoretical studies of A. Khintchine and Wiener respect to the representation of correlation functions of random stationary processes. Representation of themselves random stationary processes in a broad sense have been derived from empirical studies E. Slutsky and as a result of theoretical research in the papers of Kolmogorov, Cramer and other scientists. In theoretical studies suggest that, firstly, the independent variable is the angular frequency is a continuous variable. Secondly, the probability distribution function of angular frequency includes a step and continuous components. It turned out that the results of empirical studies contradict the theoretical studies. It was conjectured that the contradictions are the result of incorrect assumptions. It turned out that theoretical studies were performed at incorrect assumptions. Following the theoretical studies with reasonable assumptions, it turns out that random stationary in the broad sense has a unique representation as a Fourier series, defined in the Hilbert space of almost periodic functions that have only a discrete spectrum, has the ergodic property of the first order and has no ergodic property of the second order. In this empirical and theoretical studies do not contradict each other. 2. To combat the uncorrelated noise in Wiener instead of primary source data - observable input and output random processes using secondary source data in the form of deterministic correlation functions of random processes. Due to the non-ergodicity of random processes, second order, this method of dealing with noise is limited in the application of those rare cases where there is a possibility of correlation functions by averaging over multiple realizations, and the set of processes generated by the forced motion in the channel of the control object is not distorted by even small interference. In this regard, has been tasked to develop a method for interference mitigation in circumstances where the original data given in the form of single realizations of random nonergodic processes. 3. A new method for interference mitigation, based on the separation of double frequency harmonic components of the sets of bases observed signals of multidimensional control systems. As a secondary deterministic data are encouraged to use the deterministic set of harmonic frequencies of bases observed at the inputs and outputs. In many cases the set of input signals is correlated with each other. This leads to a correlated noise output signals with input signals. To overcome this limitation is proposed tertiary treatment of the original input data. Harmonic frequencies of the sets of input signals are removed coincident frequencies. The task of identification of dynamic characteristics for dependent input signals is reduced to a simpler problem of identifying the dynamic characteristics of the independent input signals. To identify the selected control channel used tertiary data input signals and the secondary data of output signals. The quaternary initial data there is a subset of conterminous frequencies between secondary and tertiary initial source data. Quaternary source data contains only a subset of the frequencies of harmonic components of accurate baseline data. The identification problem and the problem of filtering the input and output signals from noise is solved for the selected control channel only on a subset of the frequencies of Quaternary data. 4. The described method of dealing with noise is used for solving the problems of identification of dynamic characteristics of the control channels of multidimensional control objects in a class of linear stationary models. Is an example of solving the problem of identification of the control channel - given pitch angle, at the input-wire control system pitch - the actual pitch of the aircraft in the class of linear stationary models based on initial data obtained during one automatic landing (class of IL-96). Known method of solving the problem proved to be ineffective for two reasons. Firstly, in the test facility management is not satisfied the restriction on the number of inputs equal to the number of outputs. At one output simultaneously affect signals observed in the four entrances. This eliminates the need to obtain a system of integral equations of the first kind. Secondly, it is impossible to obtain the correlation functions by averaging the random signals to multiple implementations. The proposed method of dealing with noise in the conditions described above allowed us to obtain dynamic response of the control channel in a first approximation, taking into account the elastic strain aircraft structure in the form of ordinary differential equations 9-th order with constant coefficients. The control channel is first order astatism. After that, it was decided to identify the dynamic characteristics of the other control channel in the coordinates - the position of control levers engines - the pitch plane. The dynamics of this control channel is described as a first approximation, taking into account the effect of elastic deformation the airframe ordinary differential equation of fifth order with constant coefficients. The information obtained on the differential equations for the two channels control the pitch control system allows you to build a pitch in the space of states.
References
[1] | S. M. Rytov Introduction to Statistical Radiophysics. Part I. Stochastic processes. Moscow: Nauka, 1974. |
[2] | G. Kramer, M. Litbetter Stationary random processes. Moscow: Mir, 1969. |
[3] | D. Middleton, Introduction to statistical communication theory. Moscow: Soviet Radio, 1961. |
[4] | B. V. Gnedenko, Course in probability theory. Moscow: Nauka, 1969. |
[5] | E. E. Slutsky, Economic and Statistical Works: Music-M.: Penguin Books, 2010. |
[6] | A. Y. Khinchin The theory of correlation of stationary stochastic processes, "Uspekhi Mathematical Sciences, 1938. № 5. 42-51 c. |
[7] | N. Wiener Cybernetics, or Control and Communication in the animal and machine. Moscow: Nauka. |
[8] | A. N. Kolmogorov Statistical Theory of oscillations with a continuous spectrum. Anniversary Collection of the USSR, Part 1, 1947. |
[9] | N. Wiener, The Fourier integral and certain of its applications. Moscow: Fizmatgiz, 1963. |
[10] | V. I. Tikhonov Statistical radio engineering. Moscow: Soviet Radio, 1966. |
[11] | A. N. Kolmogorov, S. V. Fomin Elements of Function Theory and Functional Analysis. Moscow: Nauka, 1972. |
[12] | N. I. Akhiezer, I. M. Glazman The theory of linear operators in Hilbert space. Volume I, Kharkov: Highest School, 1977. |
[13] | G. E. Shilov Mathematical analysis (functions of one variable). Parts 1-2. Moscow: Nauka, 1969. |
[14] | V. I. Smirnov A Course of Higher Mathematics. Volume Five. Moscow: Fizmatgiz, 1959. |
[15] | R. S. Guter, L. D. Kudryavtsev, B. M. Levitan, "Elements of the theory of functions. Moscow: Fizmatgiz, 1963. |
[16] | V. V. Stepanov, Course of differential equations. Moscow: Fizmatgiz, 1958. |
[17] | G. E. Shilov, Mathematical analysis. Special rate. Moscow: Fizmatgiz. 1961. |
[18] | A. A. Kharkevich Spectra and analysis. Moscow: Fizmatgiz. 1962. |
[19] | N. S. Raybman, V. M. Chadeev. Construction of models of production processes. Moscow: Energiya, 1975. |
[20] | Technical Cybernetics. The theory of automatic control. Book 1. Ed. Doctor of Technical Sciences, Professor V. V. Solodovnikova. Moscow: Mashinostroenie, 1967. |
[21] | A. M. Tarasenkov, V. G. Braga, V. T. Taranenko Flight Dynamics. Part II. Moscow: Publishing VVIA them. NE Zhukovsky. 1974. |