• Nie Znaleziono Wyników

O pewnym rachunku operatorowym dla nieliniowych równań dynamiki

N/A
N/A
Protected

Academic year: 2022

Share "O pewnym rachunku operatorowym dla nieliniowych równań dynamiki"

Copied!
11
0
0

Pełen tekst

(1)

ZESZYTY NAUKOWE POLITECHNIKI ŚLĄSKIEJ Seria: ELEKTRYKA, z. 37

_______1972 Nr kol. 356

Magdalena Umińska-Bortliczek Instytut Podstawowych Problemów Elektrotechniki i Energoelektroniki

0 PEWNYM RACHUNKU OPERATOROWYM DLA NIELINIOWYCH RÓWNAŃ DYNAMIKI

Streszczenie. W pracy przedstawiono pewien nowy komplet operacji w dziedzinie nieliniowych równań dynamiki nazwany transformacją Cauchy-Taylor - Cauchy’ego.

A. Stosowany powszechnie w elektrotechnice rachunek operatorowy oparty o transformację Laplace’a dotyczy w zasadzie równań różniczkowych linio­

wych. Zastosowaniem rachunku operatorowego do pewnych szczególnych nie­

liniowych równań dynamiki zajmował się P. Nowacki [i] . Problematyka ta podjęta została przez autorkę w pracy [2]. Korzystając z wyników uzyska­

nych w [2] można mianowicie twierdzić, że równanie różniczkowe nieliniowe Idalej rrn) jest transformowalne wg Laplace’a wraz z rozwiązaniami i wszy­

stkimi jego kolejnymi przybliżeniami, jeżeli posiada postać następującą:

f-j^ (x,X, . . *x) + fjj iX]X, •• *x) 3 f(t), ^1/

gdzie

t - zmienna niezależna, czas

flt) - zastępcza funkcja wymuszająca typu wykładniczego, rzędu pierwszego

UI-1)

fHU,i,..-x) - nieliniowa część rrn rzędu co najwyżej UI-1/spełniająca zmodyfikowany warunek Lipschitza [5] oraz

U J -1)

fN lx,x,...x) x=x dla każdego t > 0 o

przy M, q rzeczywistych In)

fL^x,x,...x) - liniowa część rrn dowolnego N rzędu pochodnej x It) - szukane rozwiązanie.

(2)

W dalszym ciągu przedstawiona zostanie pewna nowa transformacja całko­

wa, nazwana przez autorkę transformacją Cauchy-Taylor-Cauchy’ego [5]. Za­

stosowanie transformacji Cauchy-Taylor-Cauchy’ego do rrn typu U ) z nie- liniowościami wielomianowymi pozwala zestawić w postaci tablic słownik transformacji, umożliwiający poszukiwanie rozwiązania drogą rekurencyjną.

Stanowi to istotną zaletę omawianej metody.

*B. Przyjmując ograniczenia i własności rrn U ) takie jak dla metody splotowej [3] (4] oraz zakładając podział naturalny rrn przy najwyższej pochodnej liniowej, równanie to możemy przepisać w postaci następującej:

i l t ) = A U ) . <N) . < N -1)

''' ' " \x,x,...x)

gdzie

A U ) = flt) - fj,U,x,...x; - fL (x,x,...x) ^2)

x U ) - zmienna zależna zmiennej rzeczywistej t

f U ) - jednowartościowa, ciągła ograniczona funkcja wymusza­

jąca klasy C 00, typu wykładniczego, w szczególności liniowa kombinacja funkcji typu wykładniczego}

vN- 1 ) Uf)

fU x , x , ...x,x) - liniowa część rrn 11) lub (2), o współczynnikach sta­

łych H-tego rzędu pochodnej.posiadająca znane własnoś­

ci w sensie odpowiedzi ipulsowej:

t 0

f

st at

y U ) = 2T3

J eS ds = °^e

a rzeczywiste», U - 1 )

fj}(x,x,...x) - funkcjonał zaburzenia; w szczególności są to nielinio wości typu wielomianowego [4] ,

A U ) - zmodyfikowana wymuszająca funkcja stanu będąca funk­

cją analityczną zmiennej rzeczywistej t w związku z czym współczynników jej można poszukiwać drogą reku­

rencyjną.

Zgodnie z powyższym rozwiązaniem rrn (2) jest:

t t t

xU) = J ... J J AU) dt ( 3 )

o 0 0

N-krotnie

X ) ,

W oparciu o zasadę Riemanna-Schwartza można teraz twierdzić, że ist­

nieje operacja Cauchy’ego, która jest:

Zasada Riemanna-Schwartza mówi, że: jeżeli w części wspólnej dwóch ob­

szarów w naszym przypadku jest to półoś rzeczywista funkcja x(t) oraz Wlfc) są identyczne, wówczas funkcje te określają w obszarze

Ż

Uy\.jeó- ną funkcję analityczną.

(3)

O pewnym rachunku operatorowym.. 291

1° operacją przedłużenia analitycznego rozwiązania rrn typu 12) z dzie­

dziny rzeczywistej X do dziedziny zespolonej j\. poprzez półoś rzeczywi­

stą t,

2° transformacją addytywną i jednorodną wszystkich własności rrn typu (2) z dziedziny X doj\_.

Napiszemy zatem:

UJ) . I N -1)

flt) = fL (x,X,...x) + fjjlx,x,... x)

U l) . U ł - 1 )

=S>F(X) = PT|(w,...w) + Fw lw,s,...w)

U )

gdzie

x U ) 6 X

w U“) e w(X)

PU»)

IN)

- zmienna rzeczywista X

- zmienna zespolona

A

- funkcja analityczna zmiennej zespolonej a także ostateczna postać poszukiwanego rozwiązania rrn

(2) w dziedzinie A

- obraz funkcji wymuszającej w dziedzinie A

F^lwUO.w (X)... w(X)) - obraz części liniowej rrn (1) w dziedzinieA PN iw^,) ,wU»),. .w) - obraz części nieliniowej rrn (1) w dziedzinieA Wykorzystując własności przedstawionej operacji zajmiemy się teraz rrn zmiennej zespolonej:

UJ) UJ—1)

FlX) = FL iw(X) ,w(X) , ..w(X) + Fjj(w IX), wix),...w (X)) 15)

przekształconym do postaci:

In)

w<X) = AtX)

A(X) - FtX) - PN Cw

, lN-1) U - 1 )

,w,...w) - Ft U»,,...w). 16) gdzie

A(X) - zespolona zmodyfikowana funkcja stanu, której własności pozwala ją założyć, żes

(4)

Pokażemy z kolei, że zach.0o~.1S

A(X) =

^2

wn,N

L n-0

•lm+1)

. t ’

n! IX) n!

£n=0

wn.N X

TUTS / “ Í T <*■ ' STTJ J "' m+i ' dX

n,m

n! w 1 . nN gdzie

i - delta Kronecker’a.n,m Równocześnie:

A ^ (X) n!

2 # j

h

dX

X=0

Porównując wyrażenia 18) i 19) otrzymujemy:

A iH) (X}| - wnjH= ^ - J dX

IX=o C A

Ponieważ zachodzi także:

U) lim A(X) = A(o) » Wio) X— o

lim wn,N * woN

to

IN)

lim A (X) = lim W (X) = lim w X — o X — o n— - o n,J51

A ^ = 1 ™ W n » K X D

na O

iesi szeregiem faylork jednostajnie zbieżnym wewnątrz

K

■jo,

|xj

funkcji będącej B-tą pochodną rozwiązania.

(8)

(9)

(1 0)

(11)

(

1 2

)

(

1 3

)

1 ■ dla

(5)

O pewnym rachunku operatorowym» 293

Następujące wyrażenia stanowią zatem parę transformacyjną transforma­

cji Taylora-Cauchy ego:

lTC U O

W(X) vn,N *

TFZ

» U)

< j > M L

J . n-n+T (14)

% i i * » , » } - » w t i 5 >

Łącząc operację przedłużenia analitycznego w sensie Cauchy’ego 14) z parą transformacyjną Taylora-Cauchy’ego (14) (15) uzyskujemy pewien komplet operacji, który w dalszym ciągu nazywać będziemy transformacją Cauchy-Tay lor-Cauchy’ego.

C. Własności rrn w sensie transformacji Cauchy-Taylor-Cauchy’ego zba­

dano w pracy [5]. W pracy tej obliczono transformaty typowych funkcji wy­

muszających i typowych składników rrn oraz zestawiono tablice transformat przy generalnym założeniu niezerowych warunków początkowych. Tablice te zostały przytoczone (tablice transformat CTC). Omówioną metodę ilustruje poniższy przykład.

D. Przykład.

Dany jest szeregowy dwójnik R,L,C z nieliniową pojemnością Cn (q). Na­

leży określić odpowiedź dwójnika na wymuszenie typu wykładniczego. W roz­

patrywanym układzie przy założeniu zerowych warunków początkowych obowią­

zuje równanie:

R q(t) + L q(t) + uc (q) = e(t) (16) Kładąc arbitralnie:

uc (q) - ( jj- + j} - q) q

TR

1 TTT » b

o

27E - C

%

e(t) = f(t) oraz podstawiając

q(t) = x(t)

(6)

otrzymujemy dla t > 0 następujące rrn rzędu drugiego*

x + ax + bx + cx2 =« f(t) (17)

Równanie (17) je3t transformowalne według Cauchy’ego - por. (4) [3]- z dziedziny rzeczywistej X do dziedziny zespolonej_Ą_ :

2

W(X) + a W(X) + b WU) + c W(X) = P(X) (18) W(X) - F U ) - c W2 U ) - bW(A) - a W(X) (19)

Korzystając z tablic transformat (tablice CTC) obliczamy kolejno:

W(X) = wQ)2

P U ) - f „

a* U ) - a ^2=1*2 + a 5 ^ bw

b W ( X ) + b ¿ n , 1

cW^U) - cV 4- "z %-4-z

/ , U + i j CzV?)' * ( n -'4 -z+ ? ) ( ri-4-'zV2"J'

Z

b

Ó

\ ± S 4 S ) f ł ! . , 2- ,

|_q»1 z»o

n-2-z Cn-g-zTlTTn-?-" ^

_z=o

Współczynniki w Q 2 rozwinięcia Taylora dla drugiej pochodnej poszukiwa­

nego rozwiązania obliczamy według równania rekurencyjnego ułożonego na podstawie definicji (14)*

-..2 - <-’> " £ - 2 <-'> °

£

-2 ^ \ , 1 *

■ ó K - " ) ' ■ J n , l " ’ M S o , o *

n-4 w

S

Tź+TTTż+5) * (n -4 -z+ 1 ) U -z -4 + 2 7wg n - 4 - z z=o

- O n S f e r [ * ¥ * " “ ■>] - « ł f 1 W V ,

(7)

O pewnym rachunku operatorowym.«. 295

Dalsze obliczenia przeprowadzimy dla następujących danych*

a = 2 b - 1 c = 1

f(t) = c~4t - 2 e~2t (21)

x(o) » + 1 x(o) » 0 Podstawiając (21) do (20) otrzymujemy*

w0,2 - - 3

*1,2 + + 7 w 2,2 " - 15/2!

W 3,2 “ + 31/3!

'"

4*2a

o m

- 155/4!

W5,2 " + 311/5!

Dokonując z kolei na równaniu (20) transformacji odwrotnej według defini­

cji (15) otrzymujemy:

W ( X) = - 3 - 7 x - ^ x 2 + ^ j . x 3 - + ^ 5 7 X 5 + . . .

+ W U ? * 1 (2 n+2- l ) i a + . . . (22) t-ni

stąd*

(X) = y i z l f f l ( 2 ^ - 1 ) ^ “ (23) n«o

X, X.

w (X) =

j J

W(X) dX (2n-1) xn (24)

Powracając do dziedziny rzeczywistej X wg Cauchy’ego oraz uwzględniając że x * q otrzymujemy następujące rozwiązanie*

00 \n + 1

q(t) = V " ■ (2 -1) tn - - e“2t-e_t (25)

_ n!

n« o

(8)

Wobec czego prądowa odpowiedź układu szeregowego R,L,C z zadaną nieli­

niowością jest następująca?

ilt) = qlt) = 2

e~2i

+ e-1' (2 6)

LITERATURA

1. Nowacki P.: Die Behandlung von Nichtlinearen Problemen in in Regulung- stechnik. Warszawa 1957.

2. Umińska-Bortliczek M.: 0 własnościach transformacyjnych pewnego typu równań różniczkowych nieliniowych. Zeszyty Naukowe Pol.SI. Elektryka z. 22, 1967.

3. Umińska-Bortliczek M.: 0 pewnej operacji w dziedzinie równań różniczko­

wych nieliniowych. Zeszyty Naukowe Pol.SI. Elektryka z. 36, 1972.

4. Umińska-Bortliczek M.s 0 metodzie splotowej poszukiwania odpowiedzi za­

stępczych dwójników nieliniowych. Zeszyty Naukowe Pol.SI. Elektryka z.

36, 1972 r.

5. Umińska-Bortliczek M.: Transformacja Cauchy-Taylor-Cauchy ego i jej za­

stosowanie do badania stabilności pewnych nieliniowych układów elek­

trycznych. Praca doktorska, Gliwice 1971.

O HEKOTOPOM OÍIEPATOPHOM KUMŁCJIEHMi flJŁi HEJIKHEWHŁE yPA3HHil.il JU.HAi4.ia.

P e 3 u M e

B p a ń o T e n p e » O T a B J i e H a H e i c o T o p a a H O B a s C H C T e a a n p e o Ö p a 3 0 B a H n ä b o b a a c T u H e j i H H e ü H H X y p a B H e j u t k A H H a u H K H , H a s B a H j t a a n p e o 6 p a 3 0 B a H n e u K b n u - T e i i f l e p - K o B H

ON CERTAIN OPERATOR CALCULUS POR NONLINEAR DYNAMIC EQUATIONS

S u m m a r y

In the paper a certain new complete of operations for nonlinear dyna­

mic equations is given. The complete is called Cauchy-Taylor-Cauchy tran­

sformation.

(9)

O pewnym rachunku operatorowym.«» 297

T d B L i C A T R A N S F O R M A T „ C T C -

f(t) C{f(t))‘

-F(A) TC ( T M ) - f n

1 £ « r

m*o

f * H ) " - n r

l A * j j f w r , .

A o * H f (a ' ^

3 sin ort

¿ i ' sin u A u f w2n*ł

( 2 H ) I

4 e ^ s i n o ł ¿“ sin a)A 5 ^ 1 & i ) ! % * ) !

5 m d cos u A

H f - g J -

6 eiosfiJł

m*0 5*0 ¿“ awuA

7 {9 orf

^ r - 22,) ¡

2

>d

łg o)A i ^ ł ł J n *

l i < « ' ' r

12

° ) i

1 Bn- Ium temoulieoo 8 orcłgM)

, nm OJ2"*1 jftH 2 _ . 1 *' M *

m*o a rc fg fa l)

n f -,

2

n

-1

i (2mH)

9 sh d V " !&»♦<

Z _ (&#♦»! *

m*0 s h cjA

0 5 N

M !

<0 ch d

m*0 c łiw A

63 *"

(fn)j

\

!-8 ł Z a T

m*o

<

f- aA

Q n v ł o A |

<

f

12 i ( ( « ) ) £ r * i w

m«0 4 (n)

(10)

tablicatbansfdbhat.C

(11)

O pewnym rachunku operatorowym..» 299

Cytaty

Powiązane dokumenty

W artykule pokazano, jak mo»na rozwi¡za¢ pewne równania ró»niczkowe za po- moc¡ przeksztaªcenia (transformacji) Laplace'a.. Metoda ta zostaªa omówiona na kilku

This is a substantial generalization of the so called isosceles orthogonal ad- ditivity equation in normed spaces, corresponding to the case where tp = || • ||..

Zatem, liczba warunków interpolacyjnych, które nakładamy, jest równa wymiarowi przestrzeni funkcji sklejanych rozpiętej przez nasze funkcje B-sklejane, dzięki czemu warunki brzegowe

[r]

Widzimy, że metoda Newtona sprowadza się do odpowiedniego układu równań (przechodzenie przez dany punkt i równość pochodnych)..

Katedra Ekonomii Politycznej zosta³a w³¹- czona do Wydzia³u Technologii Maszyn, a po zmianie nazwy Wydzia³u na Wydzia³ Mechaniczny Technologiczny – wesz³a w jego sk³ad. kie-

W pracy sprowadzono problem planowania posiadający reprezentację STRIPS do dwóch postaci: programu liniowego oraz układu równań i nierówności.. Następnie porównano

Typowym przykładem równań stanu o okresowo zmiennych współczynnikach jest opis matematyczny maszyny synchronicznej przy obciążeniu