Models for the description of continuous systems

Ч exp(- p Чt ) or:

p Ч T_Б Ч(1+ T Ч p)

Ч exp(- p Чt )

The linear differential equations (DE) it is the most general form of the description of dynamics, however their practical use and engineering practice often can be difficult.

The most general view of record DE:

a		d ⁿ x	+ ...+ a		dx	+ a x = b	d ^mu	+ ...+ b	du	+ b u (2.3).
ⁿ⁺¹ _dt n	² dt
			^m⁺¹ _dt m	² dt

Transfer Functions (TF). In engineering practice MMD in ви are often used by PF. For objects with self-alignment of PF in a general view:

_W ₍ _p₎ ₌ _C _Ч ^b1 ⁺ ^b2 ^Ч ^p ⁺ ⁺^bm+1 ^Ч ^p^m _Ч _exp₍_- _p _Чt ₎ _{. (n≥m) (2.4)}

a₁+ a₂Ч p + + a_n₊₁ Ч pⁿ

When using "a method of the areas" (b1=a1=1) []:

W ( p)= C Ч

b₁+ b₂Ч p

Ч exp(- p Чt ) or (2.5)

+ a

Ч p + a

Ч p ²

W ( p)= C Ч

b₁

Ч exp(- p Чt ) (2.6)

+ a

Ч p + a

Ч p ²

+ a

Ч p³

Other types, for example:

W ( p)=

C Чexp(- p Чt )

etc. (2.7)

Х(1 + T_i Ч p)

i =1

Very often apply the simplified expression of PF in the form of an aperiodic link with delay:

W ( p)=	C	Ч exp(- p Чt )	(2.8)
1 + T Ч p

For objects without self-alignment transfer function in the general case is used used in a look:

й_D	b + b Ч p + + b	_Ч_p^m щ
W ( p)= к	- C Ч		m+1		ъ	Чexp(- p Чt) (2.9)

к ^p	a + a	Ч p + + a	n+1	Ч pⁿ ъ
л			ы

(we pay attention that there is an integrated component - size 1/p) or the simplified expressions for objects without self-alignment:

The concrete type of W(p) gets out of a condition of ensuring adequacy and convenience of calculations.

Adequacy of mathematical modelcan be estimated, for example, on a formula (2.12). Ifthe value d found on this formula doesn't exceed 3-7%, then the model is considered adequate:

			N			-	Щ
	100%	X	i	X i
d =			^Че
	N			X	i
		i =1

For the equation of a look (2.8) an exit X i example in the figure 8.7):

Щ		й	t -t	щ
X (t)= C	Ч DU Ч_к1	- exp(-	)_ъ	(for

Щ		л	T	ы
		<t ) (2.13A)
X (t)=0(for	t

(2.12)

is determined analytically by a formula (see an

t іt ) (2.13)

and for W(p) of a look (2.7) at n=2: W ( p)

Щ	ж		T		t -t
X (t)= C Ч DU Ч з1	-		Ч exp(-

	з	T₁- T₂		T₁
	и

=		C Чexp(- p Чt )

(1 + T₁ Ч p) Ч (1 + T₂ Ч p)

		T		t -t	ц
) +		Ч exp(-		) ч (2.14)

	T₁- T₂		T₂	ч
		ш

In more general case as criterion of compliance at the solution of this task take criterion of a look:

n _ж	Щ ц²
з	ч	(2.15)
min(Fai), F (ai ) = е_з X i - X i _ч
j =1и	ш

where X _i – experimental value;

X i –calculated value.

For finding of coefficients of ai work out the equations:

¶F	= 0	(2.16)

¶a_i

Thus, turns out the system of the equations, solving which it is possible to define ai. Adequacy of the received mathematical models of a statics can be checked was checked by

Fischer [1] criterion (see also next lecture). In addition so-called criteria of suitability of approach [5] the R-square (determination coefficient) used for Estimate of accuracy of nonlinear models can be for this purpose used more convenient. Criterion a R-square he can accept values only from zero to unit and the closer to unit, the better the parametrical model brings closer basic data.

For his definition the criterion of SSE (Sum of squares due to error) - the sum of squares of errors on a formula is calculated in the beginning:

SSE =_еw_k Ч( y_k - y€_k )²

k =1

where w_k - weight (at us they aren't set, and are considered equal to unit); y_k - experimental (initial) values of data for each experience;

y€_k - the calculated (predicted) values of data for each experience, are received on a

formula (1);

n - amount of experimental values (for example, n=20).

The criterion the R-square (designated below as R) is defined as the relation of the sum of squares concerning regression of SSR to the full sum of squares (SST), i.e.:

	n	n		SSR	SSE
SSR =еw_k Ч( y€_k - y_k )²;	SST =еw_k Ч( y_k - y_k )²;	R =	=1 -	,
SST	SST
	k =1	k =1

where ^y^k - average value experimental (initial) data of value.
The proximity of the received values of criterion tells a R-square to unit about high precision
of the description of an experiment, for example expression of a look (2.2). Usually accepted for
practice consider values of criterion a R-square higher than 0,9.
This indicator is a statistical measure of a consent. The coefficient of determination changes
in the range from 0 to 1. If he is equal 0, it means that communication between variables of
regression model is absent. On the contrary, if the coefficient of determination is equal 1, it
corresponds to ideal model when all points of observations lie precisely on the line of regression,
i.e. the sum of squares of their deviations is equal to 0.
Sometimes indicators of narrowness of communication can give quality standard (Cheddok's
scale):
Quantitative measure of narrowness of communication Qualitative characteristic of force of
communication
0,1	- 0,3	Weak
0,3	- 0,5	Moderate
0,5	- 0,7	Noticeable
0,7	- 0,9	High
0,9	- 0,99	Very high

Functional communication arises at value equal 1, and lack of communication — 0. At values of indicators of narrowness of communication less than 0,7 size of coefficient of determination will always be lower than 50%.

Frequency Characteristics (FC). When giving on an entrance of linear system of a signal:

U (t)= A_u sin(wt)= A_u Чexp( jwt)

(2.17)

at the exit there will be a signal:

X (t)= A_x sin(wt + j)= A_x Чexp( jwt + j)	(2.18)
and APFC has an appearance:
	b + b Ч jw + + b		Ч jw^m
W ( jw)= C Ч			m+1		Чexp(- jw Чt )	=
			Ч jwⁿ
	a + a		Ч jw + + a	n+1

= A(w)Чexp( jf(w))=ReW (w)+ImW (w)
(2.19),
Where:

A(w)=	[ReW (w)]² + [ImW (w)]²		(2.20)
		[ImW (w)]
j(w) = arctan _[_ReW ₍_w₎_]	(2.21)

Also weight functions are used:

t	t
x(t)=_т g(t -t ) y(t )dt =_т g(t ) y(t -t )dt
		(2.22)

And systems of matrix linear equations of space of a status:

Model in space of parameters of a status

x = Ax + bU

y = C^T x	(2.23)

where

U - input vector; x - a vector of state variables; y - system output vector;

A - a matrix of dynamics of system; B - a control matrix; C^T - a matrix of measurement (sensors)

x = AЧ x(t)+ B Чu(t);
y(t)= C Ч x(t)+ D Чu(t)	(2.24)

1 234

Дата добавления: 2017-05-18; просмотров: 447;