CC BY
55
MSC 60H25, 60H40
ON THE MEASUREMENT OF THE «WHITE NOISE»
A.L. Shestakov, South Ural State University (Chelyabinsk, Russian Federation), G.A. Sviridyuk, South Ural State University (Chelyabinsk, Russian Federation)
In the framework of the Leontieff type equations theory we consider the mathematical model of the measuring transducer, demonstrating the mechanical lag effect. In studying of the model with deterministic external signal the methods and results of the Sobolev type equations theory and degenerate groups of operators are very useful, because they helped to create an efficient computational algorithm. Now, the model assumes a presence of white noise along with the deterministic signal. Since the model is represented by a degenerate system of ordinary differential equations, it is difficult to apply existing nowadays approaches such as Ito - Stratonovich - Skorohod and Melnikova - Filinkov - Alshansky in which the white noise is understood as a generalized derivative of the Wiener process. Instead of it, we propose a new concept of the «white noise», which is equal to the symmetric mean derivative (in the paper - the derivative of the Nelson - Gliklikh) of the Wiener process, and in the framework of the Einstein - Smoluchowsky coincides with the «ordinary» derivative of Brownian motion. The first part of the paper contains the basic facts of the Nelson - Gliklikh derivative theory adapted to this situation. The second part deals with the weakened Showalter - Sidorov problem and gives exact formulas for its solution. As an example, we present a concrete model of a measuring transducer.
Keywords: Leontieff type equations, weakened Showalter - Sidorov problem, symmetric mean derivative, Wiener process.
Introduction
One of the authors [1] proposed a mathematical model of the measuring transducer (MT), and the co-author [2] proposed to consider it as a part of Leontieff type equations
Lx = Mx + y, (1)
titled for obvious similarity with the balance Leontieff equations. Model (1) proved to be adequate to a wide range of measured phenomena [3], [4]; obtained from it data agreed well with the results of field experiments. Later, on the basis of the model (1) the authors developed a method of signal restoration, distorted by mechanical lag [5], and the resonance in the circuits of MT [6]. On the basis of this method a numerical algorithm for reconstructing a signal distorted by lag effect of MT [7] has been constructed.
In (1) L, M are square matrixes of order n, and det L = 0. If we use the well-known theory of Kronecker - Weierstrass (see for example [8], Chapter 12), then in the case of a regular pencil L + AM system (1) can be reduced to the equivalent system
Lx = Mi + y, (2)
where matrices L = diag{Nv1, NV2 ,...,NVk, Ij}, M = diag{Im ,S}, NVj are Jordan boxes of
order Vj, j = 1,k, with zeros on main diagonals; Ij and Im are identity matrices, l = n — m,
k
m = Vj; S is a square matrix of order l. In (2) m components of the vector function y = y(t)
j=i
correspond to the output signal and the remaining components characterize the condition the
MT: the vector-function y = y(t) models the input signal. The same must be said about their pre-images from (1). Algorithm of reducing of (1) to (2) in the theory is very simple, but very unstable in the numerical implementation. Therefore, to investigate (1) we used the results of the degenerate operator groups theory [9], Chapter 4. Such approach let to create not only a numerical algorithm for solving of the problems with inertia and resonances of MT, but also to construct complexes of programs [10, 11].
In all our previous works it was tacitly assumed that the input signal y = y(t) is determined, that is, it does not contain random disturbances, such as additive white noise. So far as it is difficult to imagine the MT without the white noise (WN), we intend to study this phenomenon in the model (1), where the input signal is not only an useful signal, but it’s also the WN. However, before we start this project, we want to make few remarks about the WN as a phenomenon.
It is believed that the history of the study of white-noise is based on the theory of Brownian motion of A. Einstein and M. Smoluchowsky. From this theory it follows that the displacement of a particle in Brownian motion is proportional to the \/t, where t is time. Therefore, the particle velocity is proportional to the (2\/t)-1 and therefore it is not defined at the instant time t = 0. The next step in this direction was made by N. Wiener, who suggested that the displacement of a particle is determined by a random process, which later got his name. Thus, the Wiener process is the random process w(t), and has the following properties:
(w1) w(0) = 0 almost surely (a.s.), and sample paths w(t) are a.s. continuous;
(w2) mathematical expectation E(w(t)) = 0, and autocorrelation function E((w(t) — w(s))2) = |t — s|;
(w3) sample paths w(t) are a.s. nondifferentiable for all t € [0, +ro) and on any arbitrarily small interval have unbounded variation.
Usually, under the white noise we understand the generalized derivative of the Wiener process (since «:regular> derivative does not exist due to (w3)). In this sense the white noise is presented, for example, in linear stochastic differential equation
dx = (Sx + y)dt + A5w. (3)
Here, in the right side, the symbol 5w means the generalized differential of the Wiener process
w(t), that is the WN. The equations of the form (3) were studied for the first time by K. Ito [12],
and then R.L. Stratonovich [13] and A.V. Skorokhod [14] joined the research. Their approaches
differ mainly in the interpretation of the integral / A5w(t), which appears in the right-hand side
J 0
of (3) after the integration. At the present time this approach is extended to partial differential equations [15]. In addition, recently in I.V. Melnikova school [16] a new direction of stochastic differential equations theory was formed and is actively developing. Here, the equation (3) is understood in the form
X = Sx + y + w, (4)
where all derivatives are considered in the Schwartz space. Thus, the generalized derivative of the Wiener process in the right-hand side of (4) is an additive WN.
However, the approaches of Ito - Stratonovich - Skorohod and Melnikova - Filinkov -Alshansky can hardly be applied to the equations of the form (1) research. In the first approach as it follows from the Kronecker - Weierstrass theory equations (1) (or in the equivalent form (2)) are divided into two parts; one part can be solved by the integration, such as (3), but the other part can be solved only by the repeated differentiation. And if the first derivative of the Wiener process has caused so much debate, what can we say about the second, the third and etc.?! The second approach is inapplicable due to the fact that the authors developed a theory of optimal dynamic measurements based on the theory of the optimal control of equations (1) solutions, but
it exists now only in the framework of Hilbert spaces, and nothing is heard about its distribution on a locally convex space (as in (4)).
The way out is seen by authors in the usage of the mean derivative instead of the generalized derivative of the Wiener process. E. Nelson [17], laid the fundamentals of the derivatives theory, and then Yu.E. Gliklikh [18] developed the theory till to its current state. One of the most important objects of this theory is the symmetric mean derivative of the random process, which is also called as the current rate of this process. Further for brevity, exactly this derivative will be called as the derivative of Nelson - Gliklikh, but we will keep the author’s notations. Exactly, the Nelson - Gliklikh derivative of the Wiener process w(t), for example, we will designate by symbol Dsw(t).
Let’s enumerate the advantages of such substitution. Firstly, in the case of deterministic (i.e. non-random) smooth process, the Nelson - Gliklikh derivative coincides with the «ordinary> derivative just as the generalized derivative coincides with the «ordinary> derivative of the smooth function. Secondly, the Nelson - Gliklikh derivative of the Wiener process w(t) has calculated and has been the following form: Dsw(t) = (2t)-1 w(t). Exactly this random process we call the «white noise> (<WN>), paying attention to the quotes. As well as the generalized derivative of the Wiener process, our «WN> due to (w2) has zero expectation. Finally, thirdly, if the Wiener process w(t) simulates the displacement of a particle in Brownian motion, then, according to Einstein - Smoluchowsky theory, its sample paths are a. s. equivalent to \<rt. From this Dsw(t) is a. s. equivalent to (2\/i)-1, which coincides with the «ordinary> derivative of Brownian motion.
This paper is organized as follows. The first part contains the basic facts of Nelson - Gliklikh derivative theory adapted to our situation. Here we follow mainly [18]. In the second part we investigate the Showalter - Sidorov problem for the equations of the form (1) with the additive «WN> in the right side. Moreover, this problem is taken as a «weakened> in the sense of S.G. Krein. We show the existence of this problem solutions and give the exact formula. As an example, the specific model of the WN is given. As we repeatedly noted the benefits of the Showalter - Sidorov problem compared to the Cauchy problem for Leontieff type equations (and even more general Sobolev type), refer the interested reader to [7, 19].
Finally, the authors consider their pleasant duty to express their gratitude to Yu.E. Gliklih for his strict but constructive criticism.
1. The Nelson — Gliklikh derivative
Let Q = (Q, A, P) be complete probability space, Rn - dimensional vector space endowed with the Borel a-algebra. Let’s call the mapping £ : Q ^ Rn as the random value, and denote the set of random variables by the symbol V(Q; Rn). In this set we select a Lebesgue space Lq(Q; Rn), q € [0, +ro), and note that the inclusion Lq(Q;Rn) ^ Lr(Q;Rn) is dense and continuous, if q > r, and the set Q is bounded. Suppose further that Ao is the certain a-algebra on Q, and Ao C A; by Qo = (Q, Ao, P) denote corresponding complete probability space. It is easily seen that, L2(Qo; Rn) is closed space in the L2(Q; Rn). By n : L2(Q; Rn) ^ L2(Qo; Rn) denote the orthoprojector. In the case of boundeaness of Q (further this requirement will be assume as the default setting) we uniquely continue the orthoprojector n to the projector on Lq(Q; Rn), which we denote by the same symbol n : Lq(Q; Rn) ^ Lq(Qo; Rn), q € [1, 2].
Definition 1. Let £ € Li(Q; Rn). The random value n£ € Li (Qo; Rn) is called the conditional mathematic expectation £ relatively to Ao and is denoted by E(£|Ao).
Note that if E(£|A) = £ and E(£|Ao) is mathematic expectation, if Ao = (0, Q}. Let’s also note that E(£|Ao) is unique (up to values on the set of zero probability) random variable from
Li(Q; Rn) such that / £dP = / E(£|Ao)dP for any A € Fo- This consequence of the Radon -JA JA
Nikodym theorem can be considered the equivalent definition E(£|Ao).
Let la C R be some interval, —to < a < b < +to. Consider the following mappings: f : Iba ^ V(Q; Rn), for each t € Iba putting in correspondence £ € V(Q; Rn) and g : V(Q; Rn) x Q ^ Rn, for each pair (£, u) putting in correspondence £(u) € Rn- Under the random process we call the mapping n : Iba x Q ^ Rn which has the form n = n(t,u) = g(f (t), u). Thus, at each fixed t € Ií the random process n = n(t, ■) is the random variable, i.e. n(t, ■) € V(Q; Rn), and at each fixed u € Q the random process n = n(',u) is called as the (sample) path- Set of random processes we denote by symbol P (Iba x Q; Rn).
With each £ € V(Q; Rn) we connect a-algebra A^ C A, i.e. minimal a-subalgebra A, relatively to which £ is measurable. A^ is called a-algebra, generated by £. The equivalent definition of A^ is that it’s minimal a-algebra, containing prototypes of all Borel sets in Rn under mapping £ : Q ^ Rn. With each n € P(1Í x Q; Rn) we connect three families of a-subalgebras of the a-algebra A:
- the past , generated by random variables n(S, ■) at each S € (a, t);
- the future Fn, generated by random variables n(S, ■) at each S € (t, b);
- the present , generated by the random variable n(t, ■).
All a-algebras we consider complete, i.e. containing sets of zero probability.
Let’s rename Lq(Q) = Lq(Q; Rn) and by Lq(Iba x Q) we denote the set {n € P(Iba x Q; Rn) : n(t, ■) € Lq(Q) at all t € 3^}, Q € [1, +to). Remind that n € P(Iba x Q; Rn) is called random process with a.s. continuous paths, if for P-a.s. all u € Q paths n(',u) are continuous. The set of
random processes from the Lq(ia x Q), whose paths are a.s. continuous, donote by L^(3ba x Q).
For brevity, we also rename En = E(n|Nn).
Definition 2. Let n € L0(ia x Q), q € [1, 2], the mean right Dn(t, ■) (left D*n(t, ■)) derivative of the random process n at the point t € Iba is called as the random variable
Dn(t, ■ )= lim En (n(t + Ai’ > - n(t' ■ >
/v ’ 7 Ai^0+ * \ At
(D.n(t,) = ALm- E?( > - A - >))
if the limit exists in the sense of the uniform metric on R. The random process n is called mean right (left) differentiable on the Iba, if in each point t € Iba the mean right (left) derivative exists.
So, let n € L0(ia x Q), q € [1, 2], be mean right (left) differentiable on la. Its mean right (left) derivative is also the random process, which we denote by the symbol Dn (D*n). If n € L0(ia x Q), q € [1, 2] is mean right and left differentiable on Iba, we can determine the symmetrical (antisymmetrical) mean derivative Dsn = 1 (D + D*)n (D^n = 1 (D - D*)n). Further for the sake of brevity symmetrical mean derivative Ds will be called the Nelson - Gliklikh derivative. Note that if all random variables n(t, ■ ), t € Iba, of the process n € L0(ia x Q) take values with the probability 1, and all paths n( ■ ,u) are a.s. differentiable, the Nelson - Gliklikh derivative coincides with the partial derivative with respect to the first argument.
Consider the Wiener processes and their Nelson - Gliklikh derivative. Assume Iba = R+ (= [0, +to)) and consider the family B = {Bt C A : t € R+} of a-subalgebras of a-algebra A on Q. Assume the family B be never-decreasing, i.e. D for all s > t. The process n € P(R+ x Q; Rn) is call the martingale with respect to family B, if for any t > s we have the equality E(n(t, ■ )|BS) = n(s, ■ ).
Definition 3. The random process w € L2(R+ x Q) is called the Wiener process, subaltern to B, if the Wiener process is the martingale with respect to B, w(0, ■ ) = 0 and E((w(t, ■ )— w(s, ■ ))2|Bt) = t — s for t > s.
Let’s formulate some interesting and useful consequence of Definition 3 as properties of the Wiener process.
(Cw1) (B. Levi) the Wiener process w has stationary independent Gaussian increments, moreover E(w(t, ■ ) — w(s, ■ )) = 0 and E((w(t, ■ ) — w(s, ■ ))2) = |t — s|.
In other words, the increment w(t, ■ ) — w(s, ■ ) of the Wiener process does not depend on family B and has the same (Gaussian) distribution, as the process w(t — s, ■ ) for t > s. The same can be said about the Wiener process.
(Cw2) Paths w( ■ , u) of the Wiener process are a.s. not differentiable at any point t € R+ and have unbounded variation on any interval Iba C R+.
As already mentioned, this property is the main obstacle in the mathematical research of the WN. Now list the properties of the Wiener process, which will be useful in the investigation of the «WN>.
Theorem 1. (J.E. Glikhlih) (i) DSw(t, ■ ) = (2t)-1w(t, ■ ), D|w(t, ■ ) = —(2t2)-1w(t, ■ );
t
(ii) J ||DSw(s)||ds < +to.
0
Here || ■ || - is the norm in L1(Q). Define L^((0,t) x Q) = {n € Lk-1((0,T) x Q) : there exist DSn, and all its paths are a.s. continuous on (0,t) for any t > 0}, k € N, q € [1, +to). From the theorem 1 (i) follows that w € L2((0,t) x Q), q € [1, 2]. There is stronger result:
Theorem 2. (Yu.E. Glikhlih) w € L^((0,t) x Q), k € N, q € [1, 2].
2. The Weakened Showalter — Sidorov problem for Leontieff type equations with the additive «white noise»
Let’s start with the simplification of the definitions. First of all, we rename P(Ia x Q) = P(iaxQ; Rn) the set of random processes with the values in Rn. The random process n € P(ia xQ) will be denoted by symbols n = n(t), considering its dependence on the second variable u € Q taking place by default. Nelson - Gliklikh derivative Ds (if it exists) of the random process n
o o o
will be denoted by symbol n, i.e. Dsn = Dsn(t) =n=n (t). The subset (the Lebesgue space) of
random processes P(ia x Q), with a.s. continuously differentiable up to order k inclusive (in the
sense of Nelson - Glikhlih) paths, denote by L^(3^ x Q), k € {0} U N.
Let further L h M be square matrixes of order n, det L = 0, and the pencil ^L — M be p-regular [2]. Consider Leontieff type equations
oo
L n= Mn+ w, (5)
where tv = (2t)-1w(t) - «WN>. The random process n € L1((0,t) x Q) is called a solution of the equation (5), if a.s. all its paths satisfy (5) at any t € (0, +to). A weakened (in the sense of
S.G. Krein) Showalter -Sidorov problem we call the following problem
tlim[fiL(M )]p+1(n(i) — io) = 0, (6)
t—^0+
where RL(M) = (aL — M)-1 L is the right L-resolvent of the matrix M, a € pL(M) - L-resolvent set of the matrix M. We note immediately that since ker[RL(M)]p+1 = ker P, im[RL(M )]p+1 =
imP, where
P = 2^ JY Rl(M)dp
is projector on Rn, so the problem (6) is equivalent to the problem
Um p(n(t) — Co) = 0- (7)
t—0+
So, the solution n = n(t) of the equations (5) will be called the solution of the problem (5), (6) (or, that is also equivalently, - (5), (7)), if it a.s. satisfies (6).
Along with the projector P we introduce for consideration the projector
1
Q = 2^ / LL(M W'
where LL(M) = L(jL — M) 1 is the left L-resolvent of the matrix M, and y € C, as above, is the closed contour, bounding the L-spectrum aL(M) of the matrix M.
Lemma 1. Let the pencil jL — M be p-regular, then dimker P = dimker Q, LP = QL, MP = QM.
In view of p-regularity of the pencil jL — M without loss of generality, we can assume det M = 0. Indeed, by making in (5) replacement n(t) = &(t)eat, where a € pL(M), we will reach at a system of the form (5), and in the right side will be the matrix M' = M — aL. Obviously, det M' = 0. Let’s construct the matrix (In — P)M-1(In — Q)L(In — P) = H.
Lemma 2. Let the pencil jL — M be p-regular, det M = 0, then the matrix H is nilpotent with degree not above p.
Lemma 3. Let the pencil jL — M be p-regular, then there exists the square matrix A of order n and AQL = LP A = diag{Om, Ii} is fulfilled up to permutations of rows, where m = dimker P,
l = n — m.
Assume S = AQM and construct
1
etS = 2- J R^(M)e^dp.
It is easy to see that the family {etS : t € R} forms a group, and its unit etS|t=o = P.
Theorem 3. Let the pencil pL — M be p-regular, det M = 0. Then for any £0 € V(Q; Rn) and
t € (0, +to) there exists the solution n € L^(0,t) x Q) of the problem (5), (6), and all solutions
a.s. have the form
p r t
n(t) = — V HkM-1(I„ — Q)Dks w (t) + etS^0 + / e(t-s)SAQ W (s)ds. (8)
fc=0 j0
Here H0 = \n — P by construction. From theorem 1 (ii) and theorem 2 it follows that n € L1((0,t) x Q) for all t € (0, +to). By direct substitution we see that (8) is a solution of (5). Similarly [9], Chapter 4, we see that there are a.s. no solutions of (5) different from (8). Since by action of the projector P one the first term (in fact, subtrahend) in (8), we obtain identically zero at all t € (0,t), it is clear that (7), and thus (6) are also satisfied.
Example 1. One of the simplest models of the MT, considered in [6], [7], has the following form (applied to the situation considered here)
a = Aa+ u, /3 = Ca. (9)
Here the random process a = a(t) simulates condition of the MT, matrixes A h C - its devices, 3 = 3(t) is the observed process, in this case, «WN> u =u (t). Assuming n = (a1,..., am,31, ...,3i), L = diag{Ii, Om},
M = ( A O M V C —Im
we will reach the system (5) w = (u 1,... ,u i, 0,..., 0). As it is easy to see, firstly, det M = 0,
l
since det M = (—1)m det A and det A = 0 by construction. And secondly, the pencil pL — M is 0-regular. Projectors Q = L,
P = ( I O
\C Om
operator A = L, the group of solution operators
etM = ( etA O
e V CetA Om
where etA = 2—: J(pli — A)-1e^tdp, and the contour 7 limits the region containing the spectrum
Y
a (A) of the matrix A (Incidentally, here a (A) = aL(M)). In the model (9) we are interested just
in observation 3, and in addition, for technical reasons, the initial random variable £0 (in (6),
(7)) can be set equal to zero a.s. Therefore, by Theorem 3 is valid
Corollary 1. All observations in the model (9) a.s. are given by formula
3(t) = C I e(t-s)A u (s)ds, t € R+.
0
This raises two important conclusions: firstly, the random process of observation 3(0) =0 is
a.s., and secondly, paths 3(t) are a.s. continuous.
References
1. Shestakov A.L. Dynamic Accuracy of the Transmitter with a Corrective Device in the Form of Sensor Model. Measurement Techniques, 1987, no. 2, pp. 26-34. (in Russian)
2. Sviridyuk G.A., Brychev S. V. Numerical Solution of Systems of Equations of Leontief Type. Russian Mathematics, 2003, vol. 47, no. 8, pp. 44-50.
3. Shestakov A.L. Dynamic Error Correction Transducer Linear Filter-based Sensor Model. Izvestiya VUZ. Priborostroenie, 1991, vol. 34, no. 4. pp. 8-13. (in Russian)
4. Shestakov A.L. Modal Synthesis of the Transmitter. J. of Computer and Systems Sciences International, 1995, no. 4, pp. 67-75. (in Russian)
5. Shestakov A.L., Sviridyuk G.A. A New Approach to Measuring Dynamically Distorted
Signals. Vestnik Yuzhno-Ural’skogo gosudarstvennogo universiteta. Seriya «Matematicheskoe modelirovanie i programmirovanie> - Bulletin of South Ural State University. Seria
«Mathematical Modelling, Programming & Computer Software», 2010, no. 16 (192), issue 5, pp. 116-120. (in Russian)
6. Shestakov A.L., Sviridyuk G.A. Optimal Measurement of Dynamically Distorted
Signals. Vestnik Yuzhno-Ural’skogo gosudarstvennogo universiteta. Seriya «Matematicheskoe modelirovanie i programmirovanie» - Bulletin of South Ural State University. Seria
«Mathematical Modelling, Programming & Computer Software», 2011, no. 17 (234), issue. 8.
- pp. 70-75.
7. Shestakov A.L., Keller A.V., Nazarova E.I. Numerical Solution of the Optimal Measurement Problem. Automation and Remote Control, 2012, vol. 73, no. 1, pp. 97-104.
8. Gantmacher F.R. The Theory of Matrices. AMS Chelsea Publishing: Reprinted by American Mathematical Society, 2000. 660 p.
9. Sviridyuk G.A., Fedorov V.E. Linear Sobolev Type Equations and Degenerate Semigroups of Operators. Utrecht, Boston, Koln, Tokyo, VSP, 2003.
10. Keller A.V., Nazarova E.I. The Problem of Optimal Measurement: Numerical Solution, Algorithm for Programs. News of Irkutsk State University. Series «Mathematics», 2011, vol. 4, no. 3, pp. 74-82.
11. Showolter - Sidorov Problem (Shosid Problem): Certificate 2010616865 / Keller A.V.(RU); possessor of the right(s) - South Ural State University (Chelyabinsk, Russian Federation).
- 210615137; application form 16.08.2010; registered 14.10.2010, Register of computer programs.
12. Ito K. Essentials of Stochastic Processes (Translations of Mathematical Monographs, V. 231). American Mathematical Society, 2006.
13. Stratonovich R.L. Conditional Markov Processes and Their Applications to the Theory of Optimal Control, N.-Y., Elsevier, 1968.
14. Skorokhod A.V. Markov Processes and Probabilistic Applications in Analysis. Probability theory - 1, Itogi Nauki i Tekhniki. Ser. Sovrem. Probl. Mat. Fund. Napr., 43, VINITI, Moscow, 1989, pp. 147-188.
15. Kovdcs M., Larsson S. Introduction to Stochastic Partial Differential Equations. Proceedings of -<New Directions in the Mathematical and Computer Sciences», National Universities Commission, Abuja, Nigeria, October 8-12, 2007. Publications of the ICMCS, 2008, vol. 4, pp. 159-232.
16. Melnikova I.V., Filinkov A.I., Alshansky M.A. Abstract Stochastic Equations II. Solutions In Spaces Of Abstract Stochastic Distributions. J. of Mathematical Sciences, 2003, vol. 116, no. 5, pp. 3620-3656.
17. Nelson E. Dynamical Theories of Brownian Motion. Princeton: Princeton University Press, 1967.
18. Gliklikh Yu.E. Global and Stochastic Analysis with Applications to Mathematical Physics. London, Dordrecht, Heidelberg, N.-Y., Springer, 2011.
19. Sviridyuk G.A., Zagrebina S.A. The Showalter - Sidorov Problem as Phenomena of the Sobolev-type Equations. News of Irkutsk State University. Series «Mathematics», 2010, vol.
3, no. 1, pp. 51-72.
УДК 517.9 + 681.2.08
ОБ ИЗМЕРЕНИИ «БЕЛОГО ШУМА»
А.Л. Шестаков, Г.А. Свиридюк
В рамках теории уравнений леонтьевского типа рассмотрена математическая модель измерительного устройства, демонстрирующая эффект механической инерционности. При изучении модели с детерминированным внешним сигналом очень полезными оказались методы и результаты теории уравнений соболевского типа и вырожденных групп операторов, поскольку они позволили создать эффективный вычислительный алгоритм. Теперь в модели предполагается наряду с детерминированным сигналом наличие белого шума. Поскольку модель представлена вырожденной системой обыкновенных дифференциальных уравнений, то к ней трудно применимы существующие ныне подходы Ито - Стратоновича - Скорохода и Мельниковой - Филинкова -Альшанского, в которых белый шум понимается как обобщенная производная винеров-ского процесса. Вместо этого предлагается новая концепция «белого шума», равного симметрической производной в среднем (в статье - производной Нельсона - Гликлиха) винеровского процесса, причем подмечено, что в рамках теории Эйнштейна - Смолу-ховского данная производная совпадает с «обычной» производной броуновского движения. В первой части статьи собраны основные факты теории производной Нельсона
- Гликлиха, адаптированные к рассматриваемой ситуации. Во второй - рассмотрена ослабленная задача Шоуолтера - Сидорова и даны точные формулы ее решения. В качестве примера приведена конкретная модель измерительного устройства.
Ключевые слова: уравнения леонтьевского типа, ослабленная задача Шоуолтера
- Сидорова, симметрическая производная в среднем, виннеровский процесс.
Литература
1. Шестаков, А.Л. Динамическая точность измерительного преобразователя с корректирующим устройством в виде модели датчика / А.Л. Шестаков // Метрология. - 1987. -№ 2. - С. 26-34.
2. Свиридюк, Г.А. Численное решение систем уравнений леонтьевского типа / Г.А. Свиридюк, С.В. Брычев // Изв. высш. учеб. заведений. Математика. - 2003. - № 8. - С. 46-52.
3. Шестаков, А.Л. Коррекция динамической погрешности измерительного преобразователя линейным фильтром на основе модели датчика / А.Л. Шестаков // Изв. высш. учеб. заведений. Приборостроение. - 1991. - Т. 34, № 4. - С. 8-13.
4. Шестаков, А.Л. Модальный синтез измерительного преобразователя / А.Л. Шестаков // Известия РАН. Теория и системы управления. - 1995. - № 4. - С. 67-75.
5. Шестаков, А.Л. Новый подход к измерению динамически искаженных сигналов / А.Л. Шестаков, Г.А. Свиридюк // Вестн. Юж.-Урал. гос. ун-та. Сер.: Мат. моделирование и программирование. Челябинск. - 2010. - №16 (192), вып. 5. - C. 116-120.
6. Shestakov, A.L. Optimal measurement of dynamically distorted signals / A.L. Shestakov, G.A. Sviridyuk // Вестн. Юж.-Урал. гос. ун-та. Сер.: Мат. моделирование и программирование. Челябинск. - 2011. - № 17 (234), вып. 8. - С. 70-75.
7. Шестаков, А.Л. Численное решение задачи оптимального измерения / А.Л. Шестаков, А.В. Келлер, Е.И. Назарова // Автоматика и телемеханика. - 2012. - № 1. - C. 107-115.
8. Гантмахер, Ф.Р. Теория матриц / Ф.Р. Гантмахер. - М: Физматлит, 2004.
9. Sviridyuk, G.A. Linear Sobolev Type Equations and Degenerate Semigroups of Operators / G.A. Sviridyuk, V.E. Fedorov. - Utrecht; Boston; Köln; Tokyo: VSP, 2003.
10. Келлер, А.В. Задача оптимального измерения: численное решение, алгоритм программ / А.В. Келлер, Е.И. Назарова // Известия Иркут. гос. ун-та. Сер. Математика. -Иркутск, 2011. - Т.4, №3. - С. 74-82.
11. Showolter - Sidorov problem (shosid problem): свидетельство 2010616865 / Келлер
A.B.(RU); правообладатель ГОУ ВПО «Южно-Уральский государственный университета - 210615137; заявл. 16.08.2010; зарегестр. 14.10.2010, Реестр программ для ЭВМ.
12. Ito, K. Essentials of Stochastic Processes (Translations of Mathematical Monographs, V. 231)/ K. Ito. - American Mathematical Society, 2006.
13. Stratonovich, R.L. Conditional Markov Processes and Their Applications to the Theory of Optimal Control / R.L. Stratonovich. - N.-Y.: Elsevier, 1968.
14. Скороход, А.В. Марковские процессы и вероятностные приложения в анализе / А.В. Скороход // Итоги науки и техники. Сер. Совр. проблемы математики. Фундамент. направления. - Т.43. - М.:ВИНИТИ, 1989. - С. 147-188.
15. Kovacs, M. Introduction to stochastic partial differential equations / M. Kovacs, S. Larsson // Proceedings of «New Directions in the Mathematical and Computer Sciences >, National Universities Commission, Abuja, Nigeria, October 8-12, 2007. Publications of the ICMCS. -V. 4. - 2008. - P. 159-232.
16. Melnikova, I.V. Abstract Stochastic Equations II. Solutions In Spaces Of Abstract Stochastic Distributions / I.V. Melnikova, A.I. Filinkov, M.A. Alshansky // J. of Mathematical Sciences.
- 2003. - V. 116, №5. - P. 3620-3656.
17. Nelson, E. Dynamical Theories of Brownian Motion / E. Nelson. - Princeton: Princeton University Press, 1967.
18. Gliklikh, Yu.E. Global and Stochastic Analysis with Applications to Mathematical Physics / Yu.E. Gliklikh. - London; Dordrecht; Heidelberg; N.-Y.: Springer, 2011.
19. Свиридюк, Г.А. Задача Шоуолтера - Сидорова как феномен уравнений соболевского типа / Г.А. Свиридюк, С.А. Загребина // Известия Иркут. гос. ун-та. Сер. Математика.
- Иркутск, 2010. - Т.3, №1. - С. 51-72.
Александр Леонидович Шестаков, доктор технических наук, профессор, кафедра «Информационно-измерительная техника>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), admin@susu.ac.ru.
Георгий Анатольевич Свиридюк, доктор физико-математических наук, профессор, кафедра «Уравнения математической физики>, Южно-Уральский государственный университет (г. Челябинск, Российская Федерация), Georgy_Sviridyuk@mail.ru.
Поступила в редакцию 15 июня 2012 г.
