Applied Mathematics

p-ISSN: 2163-1409    e-ISSN: 2163-1425

2012;  2(5): 146-156

doi: 10.5923/j.am.20120205.01

Suppression of Random Noise by the Separation of Frequencies

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:

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
(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
(2)
where - the measure generated by the distribution of the random corner frequency on a subset of transcendental numbers .
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
(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]
(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
(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 - the number of input - output number. Instead of accurate baseline data and can be observed or synchronous writes only approximate data, and 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]
(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]
(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
(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
(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
(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 , where- number of implementation, - set number, shows the weight function , where on the set of processes , 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 exact spectra of the first implementation of the input at the inlet and the exact components of the projection of the second implementation , 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,. 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 . Such a system of sets of input actions denoted as
Where
- 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
(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, fordimension 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 and the output signal). 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
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.