Matematici Speciale Umc
Matematici Speciale Umc
Matematici Speciale Umc
Capitole de matematici:
Speciale, probabilităţi
şi statistică
E d itu r a
N A U T ICA
2
IOAN-MIRCEA POPOVICI DORINA POPOVICI
MARIANA DUMITRU ALINA COSTEA
Capitole de matematici:
Speciale, probabilităţi
şi statistică
E d itu r a
N A U TICA
3
Editura NAUTICA, 2007
Editură recunoscută de CNCSIS
Str. Mircea cel Bătrân nr.104
900663 Constanţa, România
tel.: +40-241-66.47.40
fax: +40-241-61.72.60
4
Cuvânt introductiv
Autorii
Constanta, 12 martie 2007
5
6
Pe lângă textul foarte direct aplicat pe probleme,
această carte oferă avantajul unor linkuri care fac conexiunea
cu site-ul unei bibliografii uşor accesibile din zona
http://www.mathworld.wolfram.com.
Astfel: Analiza complexă şi celelalte capitole pot fi consultate
accesand următoarele linkuri:
http://mathworld.wolfram.com/ComplexNumber.html
http://mathworld.wolfram.com/EulerFormula.html
şi chiar o trimitere mai acoperitoare în
http://mathworld.wolfram.com/topics/ComplexNumbers.html
http://functions.wolfram.com/ElementaryFunctions
http://mathworld.wolfram.com/Cauchy-RiemannEquations.html
http://mathworld.wolfram.com/topics/ContourIntegration.html
http://mathworld.wolfram.com/CauchyIntegralFormula.html
http://mathworld.wolfram.com/CauchyIntegralTheorem.html
http://mathworld.wolfram.com/FourierTransform.html
http://mathworld.wolfram.com/LaplaceTransform.html
http://mathworld.wolfram.com/ProbabilityAxioms.html
http://mathworld.wolfram.com/topics/ProbabilityandStatistics.ht
mlhttp://mathworld.wolfram.com/topics/Probability.html
http://mathworld.wolfram.com/Probability.html
http://mathworld.wolfram.com/ProbabilityFunction.html
http://mathworld.wolfram.com/Statistics.html
http://mathworld.wolfram.com/topics/GameTheory.html
7
http://www.maths.uq.edu.au/%7Einfinity/Infinity
%2014/Win_Strats.html (un joc)
http://www.maths.uq.edu.au/~infinity/index_news.html (un link
Club Infinity, cu arborescenta in 17 directii )
8
I. ECUAŢII CU DERIVATE PARŢIALE
Motto: Unde formele canonice
sunt ţinta si cheia
10
1.2. Caracteristici
Ecuaţia diferenţială
a ( x , y ) y 2 2 b( x , y ) y + c( x , y) 0 (3)
a ( x , y) 2b( x , y) x + c( x , y) x 2 0 (4)
11
Considerăm ecuaţia (2) pe un domeniu D 0 D şi o
schimbare de variabilă arbitrară r : D 0 Δ 0 r D 0 .
1 ( x , y), 2 ( x , y)
(8)
( x , y), ( x , y)
cu şi funcţii de clasă C 2 pe D 0 .
Fie inversa sa:
x x , y y , (9)
Cu această schimbare de variabilă, vom obţine din (2) o nouă
ecuaţie cu derivate parţiale pentru funcţia U:
U , u x (, ), y(, ) , pe Δ 0 .
Ţinând seama că:
u ( x , y) U ( x , y), ( x, y) ,
avem derivatele funcţiilor compuse
u U U
+
x x x
u U U
+
y y y
2 2
�2 u �2U �
�� �2U � � �2U ���
� �+ 2 + 2 � �+
�x2 � 2 �x
� � �� �
x �x � �x �
�
U � �
� 2
U �
2
+ +
� �
x2 � �x2
12
� �2U �� � � � �
�2 u �2U �
2 + � + �+
��
x y � �x �y ��
��x �y �y �x �
�2U � � �U �2 �U �2
+ + +
� 2 �x �y � ��
x y � ��
x y
2 2
�2 u �2U �� � �2U � � �2U ���
2 � �
+ 2 + 2 � �+
�y 2
��
� y � �� �y �
y � ��
y�
(*)
U � �
� 2
U �
2
+ +
� � y 2
�
� y2
13
U
Ecuaţia (11) se scrie 0
U
Deci f (nu depinde de ) unde f este o
funcţie arbitrară care admite primitive cel puţin local.
Fie una din primitive. Din ultima egalitate rezultă că:
U , +
14
pentru care nu mai putem scrie mulţimea soluţiilor folosind
funcţii reale arbitrare.
1.4. Aplicaţii
Soluţie: b 2 ac x 2 2x 2 + y 3 x 2 + y 3.
După cum 0, 0, 0 ecuaţia este de tip
hiperbolic, parabolic sau eliptic. Se va lua un punct deasupra
parabolei y x 2 + 3, unul pe ea respectiv un punct sub
parabolă.
15
d 2 y dy
2 1 0
dx 2 dx
�1 1
1� 1+ 8 �
2 2 1 0 � 1,2 �� 1
4 �2
� 2
dy dy 1
adică 1, . Soluţiile generale ale acestor
dx dx 2
ecuaţii diferenţiale sunt: y c1 ,2 y + x c 2 .
x
Făcând schimbarea de variabile
x y, x + 2 y cu formulele (a) de derivare obţinem:
2u u
3 0.
u
cu u x ,0 y3 , y.
x0 x x 0
16
2u u
0 0, implică
u
f , u f ()d + h () g + h ()
Deci u ( x , y ) g 3x + y + h ( 2 y + x ).
u
3
Cu condiţiile iniţiale: u (0, y) y , 0, y y
x
obţinem:
�u (0, y) g ( y ) + h(2 y) y 3 (1)
�
�� u
� (0, y ) 3 g � ( y ) + h�
(2 y) y (2)
�� x
Integrând (2) în raport cu y avem:
� 1 y2
�3 g ( y ) + h(2 y ) +c
� 2 2
�g ( y ) + h(2 y ) y 3
�
şi de aici
� 1 1
�g ( y) y 2 y3 + k
� 5 5
� şi
�h(2 y ) 6 y 3 1 y 2 k
� 5 5
17
1 1
3x + y 3x + y + k �
2 3
g (3x + y ) �
5 5 �
� (3)
3 1
h( x + 2 y ) (x + 2 y) (x + 2 y) k �
3 2
20 20 �
18
II. ANALIZA FOURIER
Motto: Nu se putea ca armonia să nu aibă
loc in idealul matematic.
2. SERII FOURIER
2.1. Introducere
19
2.2. Serii trigonometrice Fourier
20
a0
f (x) + (a k cos kx + b k sin kx ) ,
2 k 1
2 +T
unde ak
T f ( x ) cos kxdx ,
2 +T
bk
T f ( x ) sin kxdx .
Prin urmare, în acest caz folosim sistemul trigonometric
1
(T ) : , cos x , sin x ,..., cos nk , sin nx ,...
2
care este ortogonal pe orice interval de forma [, + T ],
unde pulsaţia 2 / T.
In mod analog , folosind sistemul exponenţial
E : {e ikx
}k Z putem ataşa seria sa trigonometrică
+
f (x) c k e ikx
k
, unde
1 +T
ck
T f ( x )e ikx dx .
21
2.3. Aplicaţii
2 2 �x 2 2 �
�x 2 cos kxdx � � sin kx �x sin kxdx � �
0 �k 0 0
�
4 4 �x 1 �
�x sin kxdx � cos kx + �cos kxdx �
k 0 k � k 0 k 0
�
4 4
2 cos k (1) k 2 , b k 0.
k k
2 2 2
a 0 (f )
0 x 2 dx
3
.
Deci
2 cos kx
f (x)
3
+ 4 1 k
k 1
k2
.
22
2. Fie f : [0, ] R , f(x) x. Să se prelungească f
până la o funcţie periodică, pară, de perioadă T 2 şi
apoi să se dezvolte în serie Fourier trigonometrică.
2
2 / 2 , ak � x cos kxdx ,
2
ak
2
0
x cos kxdx
2
k 2
1 k 1 şi b k 0.
4
Deci , a 2 n 0, a 2n -1 , n 1
2n 1 2
şi
4 cos(2n 1) x
f (x)
2
n 1 ( 2n 1) 2
.
23
b) Să se determine seria Fourier numai de cosinusuri asociată
funcţiei pe 0, .
Soluţie. Funcţia f este continuă pe R , este integrabilă pe
orice interval compact , deci problema determinării seriei
Fourier asociate ei pe un anumit interval are sens.
a) Lungimea intervalului este 2l .
În acest caz formulele generale care ne dau coeficienţii
sunt:
2 4( 4n 2 + 1) 4
an
0
x cos x cos 2nx dx 2
( 4n 1) 2
, a0
,
2 8n
bn
0
x cos x sin 2nx dx
4n 2 - 1
, nN .
24
2
�
an x cos x cos nx dx
0
1
�
x[cos(n + 1) x + cos(n 1) x ]dx
0
� 0 dacă n 2m + 1
�
� 4 n +12
� dacă n 2m, a1
� (n 1)
2 2
2
deci , seria cerută este:
2 4 4n 2 + 1
+ cos x
2 4n
n 1
2
2
1
cos 2nx .
Aplicaţii on-line:
a. Identificaţi pe google documente pdf din acest
capitol.
b. Deschideti documentul corespunzător din
mathworld.wolfram.com
c. Trimiteţi tutorelui lucrarea cu atingerea solicitărilor din
a. si b..
25
3. TRANSFORMATA FOURIER
c
in
( x ) n e 1 (2)
1 +T
e
in
cn d, cu 2/T, T perioada
T
26
Să observăm că (1) construieşte o funcţie periodică
de perioadă 2 , ca suprapunere de oscilaţii armonice pure.
Încercând o trecere la limită după 1 cu 1 în (4),
găsim o reprezentare a unei funcţii definită pe toată axa
x ca suprapunere de oscilaţii armonice.
Ţinând cont de argumentul discret
n
n
l
se transformă în argumentul continuu prin trecere
la limită după 1, din (4) obţinem:
� �1 � �
( x) �� � ei x d � d
� 2 �
�
1 � �
� d � e i �
ei x d
2 � �
(5)
1 � �
� ei x d � e i d
2 � �
2
1 �
� ( )ei x d
2 �
2
e i
d (6)
formula căutată pentru dezvoltarea lui în oscilaţii armonice
simple este:
27
1 �
� e
i x
( x) d (7)
2 �
sau ( 6 )
1 �
� x e i x
dx
2 �
respectiv
�
( x) � ( )ei x d
�
sau ( 7 )
1 � i x
( x) � ( )e d
2 �
pentru transformata Fourier ( 6 ) şi inversa sa ( 7 ).
După factorul din faţă este evident cu care din formule se va
lucra.
28
3.3. Transformata Fourier şi diferite operaţii
n n
b k a k să avem f (b k ) f (a k ) .
k 1
k 1
Să observăm că absolut continuitatea este mai tare ca
uniform continuitatea.
Să presupunem că ( x ) este absolut integrabilă şi absolut
continuă în vecinătatea oricărui punct şi este integrabilă
pe R. Datorită integralităţii lui avem:
29
ix
ix
F((x)) (x)e dx (x)e + i (x)e ixdx Cum
() ( ) 0 conform celor relatate anterior,
obţinem :
F() iF() (8)
Cu alte cuvinte, derivării funcţiei ( x ) îi corespunde
înmulţirea funcţiei F cu i . Dacă are
derivate integrabile până la ordinul m, atunci repetând (8)
obţinem:
F x i k F , k 0, m
k (9)
Aplicaţii
Soluţie:
Folosind
30
1
e
i x
x d d , (*)
2
�
2 cos x sin a
�
�
d
31
Soluţie: conform egalităţii (*) şi a funcţiei date vom scrie
� 3
1
x
2
d
�
�
�sin cos ( x ) + i sin ( x ) d
3
� 3
1
2
d
�
�
�sin cos x d +
3
� 3
i
+ d �
� sin sin x d
2
144
�3
4 44 2 4 4 4 4 43
0
� 3
sin sin x sin d
2
� d
0 �
� 3
2 .
� sin xd � sin sin d
0
0 0
32
3 3
1
0
sin sin d
2 cos( 1) cos( + 1)
0
1 sin( 1) sin( + 1) 3 sin3
2 1 +1 0 2 -1
Deci
2 sin x sin 3
x
0 2 1
d
Aplicaţii on-line:
b. Identificaţi pe google documente pdf din acest
capitol.
b. Deschideti documentul corespunzător din
mathworld.wolfram.com
c. Trimiteţi tutorelui lucrarea cu atingerea solicitărilor din
a. si b..
33
III. ANALIZA COMPLEXA
34
4.1. Operaţii cu numere complexe
z1 z 2 x1 x 2 + i y1 y 2 ,
z1 z 2 x 1x 2 y1 y 2 + i x 1 y 2 + x 2 y1 ,
z1 x 1x 2 + y1 y 2 + i x 2 y1 x 1 y 2
,
z2 x 22 + y 22
i 1 + 2
z1 z 2 1 2 e ,
z1 1 i 1 2
e ,
z2 2
+ 2 k
in1 1 i
z1n 1n e ; n z1 e n , k 0, n - 1.
n
36
Dacă z1 z 2 implică f ( z1 ) f ( z 2 ) şi reciproc, pentru
orice z1 şi z 2 D, atunci f(z) este univalentă pe D.
Funcţia f(z) este uniformă pe D dacă îşi conservă
valoarea f ( z 0 ) din punctul z 0 şi la revenirea variabilei z
în z 0 după ce în prealabil a descris un contur din D
pentru orice z 0 D. Dacă nu este uniformă, atunci f(z) este
multiformă . Vezi funcţiile radical şi logaritm.
Funcţia f(z) derivabilă în z 0 se numeşte monogenă în z 0 .
Funcţia f(z) monogenă în orice punct din D se numeşte
olomorfă pe D.
Funcţia radical
znw (3)
este inversa funcţiei putere w z n . (4)
Dacă w e i , atunci soluţii distincte sunt:
+ 2 k
i
zk n e n ,k 0, n 1 . (5)
38
Argumentele lui z k se scriu (pentru )
0
n
2 4 2( n 1)
0 , 0 + , 0 + ,..., 0 + .
n n n
(6)
Atunci, planul (z) va fi împărţit în sectoare prin semidreptele
de ecuaţie
2k
arg z 0 + , k 0, n 1 (7)
n
Toate semidreptele acestea au ca imagine în planul (w)
semidreapta arg w n 0 .
Funcţia este univalentă în sectoarele
2 k 2( k + 1)
I k 0 + , 0 + , k 0, n 1.
n n
w 0 şi w sunt puncte critice algebrice .
39
Funcţii construite cu funcţia exponenţială
Din relaţiile lui Euler
e iu + e iu e iu e iu
cos u ; sinu
2 2i
se obţin extinderile în complex
e iz + e iz e iz e iz
cos z ; sinz (8)
2 2i
Funcţiile hiperbolice
e z + e z e z e z
cosh z ; sinhz . (9)
2 2
Funcţiile acestea sunt olomorfe în orice domeniu care conţine
punctul de la infinit. Se folosesc aceleaşi reguli de derivare ca
în cazul real.
40
Pentru k întreg rezultă că funcţia logaritm este
multiformă cu o infinitate de ramuri şi are ca puncte critice z =
0 şi z =
, numite puncte critice logaritmice.
2
4.5. Aplicaţii cu zz z .
41
4.6. Aplicaţii la olomorfie
Enunţuri
1. a) Stabiliţi domeniul de olomorfie al funcţiilor:
1 z
f ( z) ln z; f(z) ; f (z) .
z 1+ z
b) Reprezentaţi în planul complex domeniul respectiv.
3. Fie funcţiile
a) u ( x , y) (e x e x ) cos y ;
b) u ( x , y) e x ( x cos y y sin y) ;
c) v( x , y) e x + e x sin y .
Demonstraţi existenţa unei funcţii f(z) monogene pe un
domeniu Δ f (care se va determina) şi stabiliţi expresia
funcţiei respective.
Indicaţii şi soluţii
42
dacă z e i .
y
Deci , u ( x, y ) ln x 2 + y 2 ; v(x,y) arctg + 2k .
x
Condiţiile Cauchy-Riemann
��u � v � x 1
�� �x 2 + y 2
�x � y � x2 + y 2
� vor deveni �
��u
� v � y y
�
��y � x �x + y
2 2
x x + y2
2
�
cu soluţiile y = 0, x > 0.
Domeniul de monogenitate este semiaxa pozitivă a axei reale.
1 x + iy
Dacă f(z)= atunci f ( z ) 2
z x + y2
de unde
� x
u
� ( x , y )
� x + y2
2
� .
� y
v ( x, y ) 2
�
� x + y2
Condiţiile Cauchy-Riemann ne vor conduce la sistemul:
43
y 2 x 2 x 2 y 2
2xy 2xy
cu soluţia unică
x 0
y 0
care reprezintă originea sistemului de axe, dar care nu
aparţine domeniului de definiţie al funcţiei.
Deci , funcţia dată nu e monogenă în nici un punct din
planul complex.
Domeniul de monogenitate este mulţimea vidă. Pentru
z x + x 2 y 2 + iy(2 x + 1)
f (z) f ( z) .
1+ z (1 + x ) 2 + y 2
x + x 2 y2
Deci , u ( x , y) şi
(1 + x ) 2 + y 2
y(2 x + 1)
v( x , y) .
(1 + x ) 2 + y 2
Condiţiile Cauchy-Riemann ne conduc la sistemul:
44
(1 + x)(x + x 2 y 2 ) y 2 (2x + 1)
1 + 3x + 2x 2 (1 + x)(2xy + y)
(1 + x)(x + x 2 y 2 ) y 2 (1 + 2x )
sau
.
45
Se verifică uşor.
b) Determinarea funcţiei f(z) se face folosind condiţiile
Cauchy-Riemann, din care se obţine funcţia u(x,y).
x v y v
u ( x , y)
x 0 y
t,y
dt +
y 0 x
x ,t
dt
0
x0
+ e x cos y e x cos y 0 e x cos y e cos y 0 .
Ultimul termen reprezintă o constantă , deci
u ( x , y) e 2 cos y + C.
f (z ) e x cos y + C + ie x sin y .
Din condiţia dată f (0) 1 e 0 + C 1 C 0,
deci f ( z) e x cos y + ie x sin y.
c) f ( z) e x cos y + i sin y e x e iy e x + iy e z .
46
5. SERII COMPLEXE
47
1
seriile de puteri reale că este : R , unde
l
c n +1
, c n 0 sau lim sup c n .
l lim n
n cn n
Teoremă. Oricare ar fi z cu z a r ,
za (z a ) n ( n )
f ( z ) f (a ) + f (a ) + ... + f (a ) + R n ( z )
1! n!
(numită formula lui Taylor a funcţiei f(z) în punctul z = a).
lui f în z 0 luând
f (z 0 ) lim f ( z). Din acest motiv o
z z 0
49
c) Când partea principală are o infinitate de termeni,
punctul z 0 se numeşte punct singular esenţial.
Seria geometrică
1
1 + z + z 2 + ... + z n + ... , cu R = 1.
1 z
eiz eiz z z 3 z 5
sin z + + ...
2i 1! 3! 5!
�
z 2n +1
�(1)n 2n + 1 !
cu R �
n 0
50
5.6. Aplicaţii:
an +1 1
Soluţie : Dacă l lim 1, atunci R 1.
n �� a l
n
cos k 1 sin k 1
2
2
, 2
şi seria armonică este
k k k k2
convergentă. Deci, seria converge pe z 1.
51
n! 3 n 1
2.2.
n 0 3n !
z . R:
e
Aplicaţii on-line:
c. Identificaţi pe google documente pdf din acest
capitol.
b. Deschideti documentul corespunzător din
mathworld.wolfram.com
c. Trimiteţi tutorelui lucrarea cu atingerea solicitărilor din
a. si b..
52
Teoremele Cauchy
Dacă AB este un arc de curbă plană, dat prin ecuaţiile
parametrice :
x x(t );
y y(t); t a, b
, unde x(t), y(t) sunt
53
Astfel, dacă D este un domeniu triplu conex (fig.1)
atunci sunt necesare şi suficiente două tăieturi T1 şi T2
pentru a obţine un domeniu Δ simplu conex,
Δ D T1...T2 (fig.1)
T2
T1
Γ2
Γ1
Figura 1
f (z)dz 0
Γ
oricare ar fi curba închisă Γ situată în întregime în D.
f (z)dz f (z)dz .
AMB AM B
f ( z )dz f (z)dz .
Γ i 1 Γ
i
55
..............................................................................
..............................................................................
..............................................................................
..............................................................................
6.2. Aplicaţii
dz
3. Calculaţi a) C z a ; C :| z | R | a | .
dz
b) z R z 2 + 1; R 1.
e iz
4. Calculaţi z R z 2 2 dz; R .
5. Calculaţi integralele:
i z
a) ze 2 .
z R z 1
dz; R 1
cosh z
b) �x2 y 2
+ 1;a �b
a 2 b2 z2 +1
dz .
dz
c) C (z 2 + 1) n ; C : x 2 + y 2 - ay 0, a 1.
57
ez
6. I z r z 3 dz .
2i
unde f (z ) e z . Deci, I i .
2
z
cosh
2
C (z + i) 4 dz, C : z + 2i 2.
7. I
2i d3 z 4
Soluţie: I cosh .
3! dz3 2 z -i 24
z100 e iz
8. I C 2
z +1
dz, C : 4x 2 + y 2 4 0.
58
Soluţie:
z100eiz z100eiz
I 2i + 2i 2 sinh .
z+i z i z i z i
Consecinţă: O integrală complexă pe un contur închis este
sau nulă , sau se descompune într-o sumă de n integrale,
“văzute în forme Cauchy”. Calculul se face uşor şi fără erori
dacă în prealabil se reprezintă grafic curba integrală şi
singularităţile funcţiei de sub semnul integralei (numită şi
integrant). O altă formă a teoremei Cauchy este întalnită în
limbajul “reziduurilor” în capitolul următor.
Aplicaţii on-line:
d. Identificaţi pe google documente pdf din acest
capitol.
b. Deschideti documentul corespunzător din
mathworld.wolfram.com
c. Trimiteţi tutorelui lucrarea cu atingerea solicitărilor din
a. si b..
59
7. REZIDUURI ŞI APLICAŢII
1
Coeficientul c 1 al termenului se numeşte reziduul
za
funcţiei f(z) relativ la punctul singular z = a şi se notează
rez(f; a).
Ţinând seama de formula ce dă coeficienţii seriei Laurent
1 f (z)
cn
2i Γa (z a ) n +1 dz avem
1
c 1
2i Γ f (z)dz (2)
Γ
f ( z )dz 2i rez(f , a
k 1
k
) (7)
62
7.2. Integrale cu teorema reziduurilor
P( x )
I. Integralele de forma: Q( x ) dx ,
unde
x2 +1
Exemplu: I
x4 +1
dx
63
b
r O r x
Figura 2
Atunci
z2 +1 R x2 +1 z2 +1
Γ z4 +1
dz
R x 4 + 1
dx +
z4 +1
dz
z2 +1
Funcţia f ( z ) are poli simpli:
z4 +1
+ 2k + 2k
z k cos + i sin , k 0,1,2,3
4 4
din care numai
3 3
z 0 cos + i sin şi z1 cos + i sin
4 4 4 4
se află în interiorul conturului mărginit de curba Γ .
Aplicăm teorema reziduurilor şi avem:
z2 +1 R x +1
2
z2 +1
z 4 + 1dz R x 4 + 1 z 4 + 1dz
dx +
2i[rez ( f , z 0 ) + rez ( f , z1 )]
dacă R z .
64
Dacă trecem la limită în egalitatea precedentă şi ţinem seama
că R lim zf ( z) 0 rezultă:
z
R x2 +1 x 2 +1
lim
R R x4 +1
dx + x 4 +1
dx 2i[ rez(f , z 0 ) + rez(f , z1 )]
.
Calculăm
z2 +1 z(z 2 + 1) i 2
rez(f , z 0 ) lim lim ,
zz0 4z 3 zz0 4z 4
4
z2 +1 z(z 2 + 1) i 2
rez(f , z1 ) lim lim .
z z1 4z 3 z z1 4z 4 4
În consecinţă
x2 +1 i 2 i 2
I
4
x +1
dx 2i
4
2 .
4
2
II. Integralele de forma: I 0 R sin , cos d ,
unde R(u,v) este o funcţie raţională , se pot calcula cu
teorema reziduurilor dacă se face schimbarea de variabilă
z e i . Atunci
sin
1 i
2i
e e i
1 1
z
2i z
cos
1 i
2
1 1
e + e i z +
2 z
65
iar dz ie i d. Când parcurge intervalul 0,2, z
parcurge cercul z 1 o singură dată.
Ca urmare
2
I 0 R sin , cos d
1 1 1 1
x 1 R 2i z z , 2 z + z dz 2i
k
rez( R , z k ).
2 d
Exemplu: I 0 2 sin 2
.
4 zdz 4
� 4 �2 i �rez f , zk
i z 1 z + 6z +1 i
2
k
z
Funcţia are patru poli, iar în cercul z 1 se
z 4 + 6z 2 + 1
află polii z1 3+ 8 , z2 3 + 8.
66
z2
rez f , z1 lim z
lim
z z1 4z 3 + 12 z z z1 4z 4 + 12z 2
2
3+ 8
4 2
4
3+ 8
+ 12 3 + 8
z2
rez f , z 2 lim
z
z z 2 4z 3 + 12z
lim
z z 2 4 z 4 + 12z 2
2
3+ 8
4 2
4
3+ 8
+ 12 3+ 8
Avem
rez f , z1 + rez f , z 2 8
2
2 d 2
şi deci I 0 2 sin 2
8
8
2.
67
Pentru calculul integralei K pe conturul menţionat utilizăm
funcţia: f ( z) F( z ) e iz .
x cos 2 x
Exemplu: I x 2 2x + 5dx .
Pentru calculul integralei I, asociem integrala
x sin 2 x
J x 2 2x + 5dx
xe 2ix
şi împreună cu I avem: K I + iJ
x 2 2x + 5dx .
Evident ReK = I şi ImK = J.
Conturul de integrare Γ AB , unde A(-R, 0) , B(R,
0), iar este semicercul x 2 + y 2 R 2 , y 0 . Calculăm
ze 2iz
integrala
Γ z 2 2z + 5dz care pentru y = 0 se reduce la
K.
Avem
ze 2ix R xe 2ix
Γ z 2 2z + 5dz R z 2 2z + 5dz 2i
k 1
rez(f , z k )
68
ze 2iz
Funcţia f (z ) are polii z1 1 + 2i şi
z 2 2z + 5
z 2 1 2i, din care numai z1 este în interiorul
domeniului mărginit de Γ . În consecinţă
( z z1 ) ze 2iz
rez ( f , z1 ) lim
z � z1 z z z z
1 2
(1 + 2i )e 2i (1+ 2i ) 1 + 2i �
e 4 �
e 2i
4i 4i
Ca urmare
ze 2iz (1 + 2i )e 4 e 2i
z 2 2 z + 5 dz 2i
4i
[cos 2 2 sin 2 + (2 cos 2 + sin 2)]
2e 4
Trecând la limită pentru R obţinem:
� xe 2ix
K � 2 dx
� x 2 x + 5
[cos 2 2sin 2 + i 2 cos 2 + sin 2 ]
2e 4
respectiv
I 4
cos 2 2 sin 2; J 2 cos 2 + sin 2
2e 2e 4
69
7.3. Aplicaţii la reziduuri cu integrale şi formule pentru
reziduuri
Să se arate că:
x4 +1 4
1.
6
x +1
dx
3
x2
2. x 2 + a 2 2 dx 2a
x2
3. 0 x 2 + a 2 2 dx 4a
dx 2a + b
4. x 2 + a 2 x 2 + b 2 2a 3 a + b 2 b
sin 2 z
5. I
C z 2 z 2 + 4iz 3
dz
1
a) C : z 3 ; b) C : z i 1 ; c) C : z + 1
2
70
dz
6. z a z 3 (1 + z 2 ) pentru a < 1, a > 1.
cos n
7. I
b ia cos d
0
1 i sin n
Soluţie: J � d
2 b ia cos
căreia îi ataşăm integrala nulă
1 e in
I + iJ
2
b ia cos
d .
dz z2 +1
Notând e i z, d , cos .
iz 2z
zn
Prin înlocuire obţinem I + iJ az 2 + 2biz + a
x 1
dz .
71
zn
rez ( f ; z1 )
2az + izb z z1
in an
�
n
2i a 2 + b 2 a2 + b2 + b
Rezultă:
i n an
I 2irez f ; z1
2 2
n .
a +b
a 2 + b2 + b
sin x sin nx
8. I
5 4 cos x
dx
e inx sin x
Calculând J + iI iI
5 4 cos x dx.
Punând e ix z
1 z 2 1z n 1dz 1 2irez f ; 1
iI
2
z 1
2z 2
5z + 2 2 2 2 n +1
72
IV. ANALIZA DE CALCUL OPERATIONAL
8. CALCULUL OPERAŢIONAL
(Transformatele Laplace şi Fourier)
73
astfel încât f(t) să fie derivabilă pe fiecare interval x j1 , x j
şi să existe limitele laterale
f ( x j1 + 0), f ( x j 0), f ( x + 0), f ( x 0), j 1, n .
Cazuri concrete.
(d1) Funcţia f ( x ) e bx , cu b real sau complex va avea
creştere exponenţială putând lua a = Reb, M > 1 şi x 0 R
Într-adevăr,
f ( x )e ax e (Re b ) x e ax 1 M ( 1 )
1
L(f ( x ))(p) 0 e bx e px dx 0 e ( p b ) x dx
pb
( 2 )
74
1 1
Deci exp(bx ) şi exp(bx ) ,
pb pb
convergenţa integralei având loc pentru p > Reb dacă p R
şi Rep > Reb dacă p C.
Observaţie.
Nu este greu să vedem că L este liniară. Utilizând
liniaritatea, rezultatul din d1 şi relaţiile lui Euler:
cos t
2
1 i t
e + e i t ,
sin t
1 it
2i
e e it
obţinem:
1 1 1 p
cos t + ( 1 )
2 p i p + i p 2 + 2
1 1 1
sin t ( 2 )
2 p i p + i p 2 + 2
pentru Re p Im .
76
F( p) 0 f ( x ) e px dx
(7)
şi pentru p = 0 obţinem
f (x )
x
dx F(g )dg (11)
0 0
F ( n +1) (g )dg 1 n +1 x
n
f ( x ) e px dx (12)
0 p
77
care pentru p = 0 devine
F
( n +1)
(g )dg 1 n +1 x n f ( x )dx
0 p
(13)
Pentru n = 0, egalitatea (13) devine
f ( x )dx F(q)dq F(q)
0
(14)
0 0
8.2. Aplicaţii.
1. Să se arate că
sin 2 u
u2
du
2
0
78
��sin 2 xu � � 1
L I ( x) L �� 2 du � �2 � L sin 2 xu du
�0 u � 0u
� �
1 1 1 1 �1 p �
� 2 �
L 1 cos 2 xu du 2 �
2 �du
20u 2 0 u �p p + 4u 2 �
�
2 du 1 2u �
� 2 2 2 arc tg .
p 0 4u + p p p 0 2 p2
Deci
L I( x )
2p 2
şi
sin 2 xu
I( x )
2
x
u2
du
0
Trecând la limită pentru x 1 obţinem
sin 2 xu
u2
du
2
.
0
79
Soluţie: Teorema derivării originalului conduce la
L y ( x) p 3 Y ( p ) p 2 y (0) py (0) y (0)
p 3Y ( p ) p 2 p + 1
L y ( x ) p 2 Y ( p ) py (0) y (0)
p 2 Y ( p) p 1
L y ( x) pY ( p ) y (0) pY ( p ) 1
L y ( x ) Y ( p )
10
şi L 5 sin 2 x 2
p +4
Se obţine ecuaţia operaţională
p 3
2p 2 p + 2 Y (p) p 2 + p + 2
2
10
p +4
,
din care
1 1 5 1 1 p 2
Y ( p) + + +
3 p + 1 12 p 2 4 p 2 + 4 p 2 + 4
80
3x x
y( x ) + 4 y( x ) sin sin
2 2
şi
y(0) 1, y(0) 0
din care
1 p 5 p 1 p
Y ( p) +
6 p +1 6 p + 4 2 p2 + 4 2
2 2
Utilizând teorema derivării imaginii se obţine
2 4p
L x sin 2 x 2
p + 4 p + 4 2
2
iar
1 5 1
y( x ) cos x + cos 2 x + x sin 2 x.
6 6 8
81
�
L(af ( x) + bg ( x))( p) �(af ( x) + bg ( x)e px dx
0
� �
a �f ( x )e px dx + b �g ( x )e px dx
0 0
Demonstraţie:
L f ( ax ) ( p ) f ( ax )e px dx (cu ax t)
0
p
1 t 1 p
a 0
f (t )e a
dt
a
L ( f ( x )) .
a
�
L f ( x ) (ap ) � f ( x )e apx
dx (cu ax t)
0
82
1 t 1 x
a 0 f e pt dt L f (p).
a a a
Într-adevăr,
L( f ( x a) H ( x a))( p )
�
�f ( x a ) H ( x a )e px dx
0
�
(cu x-a t) �H(t) f (t )e p (t +a ) dt
0
�
�
e ap �f (t )e pt dt e ap �f (t )e pt dt
0
0
ap
e F ( p ), a 0
Dacă relaţia (3) se foloseşte de la dreapta la stânga,
trebuie ţinut cont că originalul H(x-a) = 0 pentru x < a
şi atunci, fără a mai utiliza funcţia unitate avem:
83
e ap F( p) f ( x a ).
Într-adevăr,
�
L( f ( x + a))( p) ��� f ( x + a)e px dx �
�
� 0 �
�
��� f (t )e p (t a ) dt �
�
� a �
� a
e ap �f ( x )e px dx e ap �f ( x )e px dx.
0 0
Teorema de deplasare
e x f ( x ) F( p + ), C (3)
L e x f ( x ) ( p )
0 f ( x )e p + x dx F( p + )
pF ( p ) f (0+)
Deci f(x) = p F(p) - f(0+) şi pentru n = 1, (5) se verifică.
Pentru n = 2 avem:
f ( x ) f ( x ) pL (f ( x )) f (0+) (ţinând cont de
rezultatul pentru n = 1 obţinem) =
p 2 F( p) pf (0+ ) f (0+ ). Pentru n > 2 iterăm acest
procedeu.
Operaţia de derivare în spaţiul funcţiilor original se transformă
în operaţia de înmulţire cu p în spaţiul funcţiilor imagine,
abstracţie făcându-se de un polinom în p.
85
În ipoteza că xf ( x ), x 2 f(x),..., x n f(x) sunt
funcţii original se poate argumenta posibilitatea derivării sub
integrală în raport cu p în relaţia de definiţie:
F( p) 0 f ( x )e px dx
Obţinem formulele:
�
( x) n f ( x) F ( p),
n
(6)
n 0,1,...
Deci şi acestei operaţii de derivare îi corespunde operaţia de
înmulţire cu “x”.
x
Demonstraţie: Fie g ( x ) 0 f ( t )dt .
0 1
Cum f C ( R ), g C (R + ) şi g(0+) = 0 (din modul
de definire a lui g). g ( x ) pG ( p) g (0+ ) pG ( p)
Deci
L g ( x ) ( p) L f ( x ) ( p) pG ( p) pL
x
0 f ( t )dt (p)
de unde se obţine (7).
f (x)
valabilă când este funcţie original.
x
Demonstraţie:
p F(s)ds f ( x )e sx dx
p 0
f ( x )
f ( x )e sx dx
e sx s p dx
0 0 x s
f ( x ) px f (x)
0 x
e dx L
x
( p).
� a
Din faptul că sin ax , şi
p2 + a2
� p
cos ax cu p > 0
p + a2
2
87
a
ebx sin ax ,
( p b) 2 + a 2
p b
ebx cos ax ,
( p b) 2 + a 2
a 0, b �R.
cu formula:
f (x)
0 x
dx 0 F(q )dq
Într-adevăr,
f (x)
x
p
F(q )dq
se mai scrie
88
f (x)
0 e -px dx 0 F(q )dq
x
Făcând p = 0 obţinem relaţia căutată.
� 1
a) sin x , atunci cu d avem:
p +1
2
sin x dp
0 x
dx
0 2
arctgp
p +1 0 2
b) pentru a > 0, b > 0
e at e bt b
0 t
dt ln
a
Într-adevăr, avem:
e at e bt 1 1
0 t
dt
0
p+a p+b
dp
p+a p a b
ln ln ln .
p+b p 0 b a
89
1
�t e e bt sin mtdt
�
at
0
1.
b a
arctg arctg , a , b, m 0
m m
cos at cos bt b
2.
0 t
dt ln , a , b 0
a
3.
1 ax a
0 x
e sin bxdx arctg , a 0, b 0.
b
90
utilizând tabelul “ f(x) = F(p) “ obţinem soluţia din “spaţiul
original “. Calculul operaţional este calculul care utilizează
transformata Laplace.
Exemple:
91
1. Fie ecuaţia : x 5x + 6 x e t cu
x (0) 1, x (0) 1.
1 t 7
Arătaţi că x ( t ) e 5e 2 t + e 3t
2 2
Soluţie:
2
x�� (0) p 2 X ( p) + p 1
p X ( p) px(0) x�
1
et
p 1
Ecuaţia operaţională este:
2
p 5p + 6 X ( p) p + 6 +
1
p 1
Descompunând în fracţii simple obţinem:
5 7 1
X ( p) + +
p 2 2( p 3) 2( p 1)
De unde:
1 t 7 3t
x(t ) e + e 5e 2 t .
2 2
2. x 4 x + 4 x sin t cu x (0) 1, x (0) 2.
2
Soluţie: p 4p + 4 X ( p) p 2 +
1
1+ p2
92
21 1 3 4p
X ( p) + + +
25( p 2) 5( p 2) 2 25(p 2 + 1) 25( p 2 + 1)
21 2 t 1 2 t 3 4
x(t ) e + te + sin t + cos t
25 5 25 25
Soluţie:
1 1 5 p 1 4
X( p) +
2p 2p 2 2 p 1 2 + 1 p 1 2 + 1
1 t 5 t
x (t ) + e cos t 4e t sin t
2 2 2
2(p + 1) 2
Soluţie: ( p + 2) X (p) 2
+
(p + 1) + 4 (p + 1) 2 + 4
2
X ( p) şi x ( t ) e t sin 2 t
p + 1 2 + 1
5.
x 3y 3x + 3z, y x - t, z -z, x(0) x (0) 0
93
y (0) 0, y(0) -1, z(0) 1, z (0) 0.
Soluţie:
p 2 X 3 Y X + Z , p 2 Y X - Y - 1, p 2 Z -Z + p de
unde:
3( p 1) 3(p - 1) p
X 2 2
, Y 2 2 2
, Z 2
(p + 4)p p (p + 1)(p + 4) p +1
fază în care putem considera problema rezolvată aproape în
totalitate (restul calculelor fiind de rutină şi uzură).
funcţiei f.
Din cele făcute până acum rezultă că pentru a exista
transformata Laplace L (f )(p ) este suficient ca f să aibă
abscisă de R sau şi
convergenţă
transformata are sens pentru p în cazul p R şi
94
Re p în cazul p C. Vom nota cu f , g abscisele
de convergenţă pentru f şi g.
b. Pentru p max f , g , are loc proprietatea de
�
liniaritate af + bg aF ( p ) + bG ( p ).
Aplicaţii on-line:
a. Identificaţi pe google documente pdf din acest
capitol.
b. Deschideti documentul corespunzător din
mathworld.wolfram.com
c. Trimiteţi tutorelui lucrarea cu atingerea solicitărilor din
a. si b..
95
V. ANALIZA STOCHASTICA
9. STRUCTURI PROBABILISTE
Probă. Eveniment
Să considerăm experienţa aruncării unui zar. Este,
evident, vorba de o experienţă aleatoare adică o experienţă al
cărei rezultat variază la întâmplare. (Cuvântul ,,aleator’’ are
sens de intâmplător şi provine de la latinescul alea, care
inseamna zar). Să notăm cu {1} apariţia feţei cu un singur
punct, cu {2} apariţia feţei cu două puncte etc.
În urma unei aruncări obţinem unul din rezultatele
(I) {1}, {2}, {3}, {4}, {5}, {6}.
Alte rezultate nu mai sunt posibile şi unul dintre ele se produce
neapărat. Acestea sunt probele experienţei. În general, probă
96
se numeşte rezultatul unei experienţe aleatoare. În legătură cu
aruncarea zarului ne putem pune o serie de întrebări ale căror
răspunsuri nu le putem cunoaşte decât după apariţia uneia din
probele (I). De pildă, ne putem întreba dacă vom obţine o faţă
cu 4 puncte etc.
Toate aceste situaţii legate de experienţa noastră şi
despre care putem spune cu certitudine că s-au produs sau nu,
după efectuarea experienţei, poartă numele de evenimente .
Pentru precizarea acestor noţiuni să reluăm experienţa
aruncării unui zar şi să presupunem că ne punem întrebarea
dacă vom obţine o faţă cu un număr par de puncte. În acest
caz, experienţa este aruncarea zarului, proba este rezultatul
care se obţine la sfîrşitul experientei, iar evenimentul care ne
interesează este apariţia unui număr par de puncte.
Evenimentul se realizează, dacă se obţine una din
probele {2},{4},{6}, şi nu se realizează în caz contrar. Există
deci, trei probe care realizează evenimentul şi acest eveniment
este perfect determinat de cele trei probe.
Dacă în cadrul aceleiaşi experienţe, ne interesează
apariţia feţei cu un punct, suntem în prezenţa unui eveniment
care poate fi realizat de o singură probă, şi anume proba {1}.
Evenimentul care poate fi realizat de o probă şi numai
de una se numeşte eveniment elementar.
Celelalte evenimente le numim compuse.
97
Alte exemple
1. Să considerăm experienţa constând din aruncarea unui zar
de două ori la rând. Rezultatul experienţei nu va mai putea fi
reprezentat printr-un număr, ci printr-o pereche de numere.
Probele experienţei sunt:
100
Eveniment sigur, Eveniment imposibil
Fiecărei experienţe i se ataşează două evenimente cu
caracter special: evenimentul sigur şi evenimentul
imposibil. Evenimentul sigur este un eveniment care se
realizează cu certitudine la fiecare efectuare a experienţei.
De exemplu, la aruncarea unui zar, apariţia uneia din feţele
1,2,3,4,5,6; este evenimentul sigur al experienţei.
Scoaterea unei bile de culoare albă sau neagră dintr-o
urnă, conţinând numai bile albe şi negre, este de asemenea un
eveniment sigur. Evenimentul imposibil nu se realizează la
nici o efectuare a experienţei, în cazul aruncării unui zar,
apariţia unei alte feţe decît 1,2,3,4,5,6 este un eveniment
imposibil. Extragerea unei bile de altă culoare decît alb sau
negră dintr-o urnă, conţinând numai bile albe şi negre este de
asemenea un eveniment imposibil. Este clar că evenimentul
imposibil constă în nerealizarea evenimentului sigur, sau că
evenimentul sigur constă în nerealizarea evenimentului
imposibil. Evenimentul sigur îl vom nota cu litera E, iar
evenimentul imposibil cu .
Evenimente contrare
Să notăm cu A evenimentul apariţiei uneia din feţele
1,2,3 la aruncarea unui zar şi cu B apariţia uneia din feţele
4,5,6. Se observă că atunci când A nu se realizează, adică nu
obţinem ca rezultat al aruncării una din feţele 1,2,3, se
101
realizează B, adică obţinem una din feţele 4,5,6 şi invers, când
nu se realizează B, se realizează A. În cazul unei extrageri
dintr-o urnă conţinând bile numai albe şi bile negre să notăm cu
A evenimentul extragerii unei bile albe şi cu B evenimentul
extragerii unei bile negre. Şi acum se observa că nerealizarea
lui A este echivalentă cu realizarea lui B, iar nerealizarea lui B
este echivalentă cu realizarea lui A. În ambele cazuri vom
spune că evenimentele A şi B sunt evenimente contrare.
Întotdeauna unui eveniment îi corespunde un eveniment
contrar, a cărui realizare înseamnă prin definiţie, nerealizarea
primului.
Se observă că dacă un eveniment B este contrariul
unui eveniment A, atunci şi A este contrar al lui B. Remarcăm
că evenimentul sigur şi evenimentul imposibil sunt contrare
unul altuia. Evenimentul contrar unui eveniment A îl vom nota
în acest caz cu A sau CA. În acest caz sunt evidente relaţiile:
A A; E ; E.
104
urmărim dacă obţinem sau nu una din probele {2}, {4}, {6}.
Altceva nu ne interesează.
Evenimentul este perfect determinat de mulţimea
formată din aceste trei probe şi deci îl putem identifica.
Vom scrie: A={2, 4, 6}.
Dacă ne interesează apariţia unei feţe cu un număr 4
puncte, urmărim dacă apare una din probele {1}, {2}, {3}, {4} şi
atât. Evenimentul respectiv este determinat de mulţimea de
probe care îl realizează: {1,2,3,4}.
Acum este justificată şi notarea evenimentului
imposibil prin deoarece mulţimea probelor care îl
realizează este vidă. Este justificată de asemenea, notarea
evenimentului contrar lui A prin CA.
Am văzut mai înainte că evenimentul elementar este
realizat de o singură probă. El este perfect determinat de
proba care îl realizează. Am văzut că „A implică B", înseamnă
că realizarea evenimentului A atrage după sine realizarea
evenimentului B, adică printre probele, care realizează pe B se
găsesc toate probele care realizează pe A şi deci pntem scrie
A B.
Operaţii cu evenimente
Câmp de evenimente
Mulţimea tuturor evenimentelor legate de o experienţă
( inclusiv evenimentul sigur şi evenimentul imposibil ) formează
un câmp de evenimente.
Să considerăm o urnă conţinând patru bile numerotate
1,2,3,4. Câmpul de evenimente al experienţei constând într-o
extragere din această urnă este format din: , {1}, {2}, {3},
{4}, {1,2}, {1,3}, {1,4}, {2,3}, {2,4}, {3,4}, {1,2,3}, {1,2,4}, {1,3,4},
{2,3,4}, {1,2,3,4} = E.
106
9.2. Probabilitate, Câmp finit de probabilitate
Frecvenţa
Să considerăm o experienţă şi un eveniment A
corespunzător acestei experienţe. Să repetăm această
experienţă de n ori în condiţii identice, să notăm prin
numărul de realizări ale evenimentului A şi prin (n )
numărul de nerealizări ale lui A.
Numărul fn poartă numele de frecvenţă.
n
De exemplu să aruncăm o monedă de o sută de ori şi să
admitem că faţa conţinând stema a apărut de 51 de ori.
51
Numărul f100 reprezintă frecvenţa apariţiei
100
stemei în aceste 100 de experienţe. Frecvenţa variază de la
experienţă la experienţă. Ea are un caracter empiric
experimental. Numărul
poate varia de la 0 la n inclusiv.
Avem =0, când din n repetări consecutive ale experienţei,
evenimentul A nu s-a realizat niciodată. Dimpotrivă, dacă din n
experienţe consecutive, evenimentul A s-a realizat în toate
experienţele, rezultă n.
În toate celelalte cazuri 0 < < n. De aici rezultă 0 f n 1
, oricare ar fi n.
107
Experienţa arată, că pentru multe fenomene de masă
frecvenţa fn pentru n crescând necontenit, se apropie din ce
în ce mai mult de o anumită valoare. Această proprietate poartă
numele de lege a numerelor mari .
Exemple:
1. Experienţa constă din aruncarea unei monede.
Fie A şi B evenimentele de a ieşi respectiv o faţă sau
cealaltă. Dacă moneda este perfectă, nu avem nici un motiv să
admitem că una din feţe are o şansă mai mare de apariţie decît
alta. Acest lucru se confirmă experimental, prin faptul că
aruncând moneda de un număr mare de ori cele două feţe apar
aproximativ la fel de des. Evenimentele A şi B sunt egal
posibile.
Probabilitate
Exemple:
.i 1, 6, j 1, 6
110
c) Dar probabilitatea ca numărul înscris pe bila extrasă din
prima urnă să fie mai mic decât cel de pe bila extrasă din a
doua urnă ?
111
m
0 1
n
Pentru a demonstra relaţia 2., vom observa că pentru
evenimentul m<n înseamnă că există evenimente elementare
pentru care evenimentul E nu se verifică, ceeea ce este
imposibil‚ întrucât stim că evenimentul E se realizează
întotdeauna.
Relaţia 3. rezultă din faptul că pentru evenimentul Ǿ
avem m = 0. Într-adevăr dacă am avea m > 0 înseamnă că
există cel puţin un caz posibil care realizează pe Ǿ, ceea ce
este absurd.
Pentru a demonstra relaţia 4., notăm prin numărul n
numărul cazurilor egal posibile, prin m numărul cazurilor
favorabile realizării evenimentului A, prin s numărul cazurilor
favorabile realizării evenimentului B.
Deoarece A∩B=Ǿ înseamnă că evenimentele A şi B
sunt incompatibile între ele. Nici unul din cazurile favorabile
realizării evenimentului A nu poate realiza pe B şi invers.
Rezultă că m + s reprezintă numărul cazurilor egal posibile
favorabile evenimentului A U B.
Conform definiţiei probabilităţii, avem:
m s
P ( A) , P ( B) ,
n n
112
m+s
P( A U B)
n
şi deci:
P ( A B) P( A) + P ( B ),
dacă A I B = Ǿ.
Pentru a demonstra relaţia 5., observăm că
A U A E, A I A Ǿ
În baza egalităţii 4., reiese
P ( A) + P ( A ) P ( E ).
Ţinînd seama de 2., obţinem imediat egalitatea 5.
Această ultimă relaţie prezintă deseori interes la
rezolvarea problemelor.
Dacă vrem să calculăm probabilitatea unui eveniment
A observăm că raţionamentul şi calculele sunt mai dificile
decât la calculul probabilităţii evenimentului contrar A
calculăm mai întâi probabilitatea P ( A ) a acestuia din urmă
şi apoi calculăm:
P ( A) P ( A ),
Evenimente independente
Fie A şi B două evenimente, Dacă
P(A∩B)=P(A)xP(B)
evenimentele A şi B sunt, prin definiţie, independente.
Exemple:
113
1. Considerăm experienţa constând din aruncarea a două
zaruri, unul roşu şi celălalt verde. Fie A evenimentul ca zarul
roşu în experienţă să apară cu faţa 1 şi B evenimentul ca zarul
verde să apară cu faţa 5. Sunt evenimentele A şi B
independente ?
Deci
6 1
P( B)
24 4
Pentru A ∩ B, avem două cazuri favorabile
(4 1 2 3), (4 1 3 2).
Rezultă
2 1
P( A B)
24 12
si cum
P( A I B ) �P ( A) �
P( B)
obtinem
1 1 1
� �
12 4 4
Rezultă că evenimentele A şi B nu sunt independente.
116
Astfel evenimentele A, B, C sunt independente dacă avem
toate egalităţile:
P ( A B ) P ( A) P ( B ),
P ( A C ) P ( A) P (C ),
P ( B C ) P ( B ) P (C ),
P ( A B C ) P ( A) P ( B ) P (C ).
Exemple:
1. Urna U conţine şapte bile albe şi trei roşii, iar urna V conţine
trei bile albe. Extragem câte o bilă din fiecare urnă. Care este
probabilitatea ca ambele bile extrase să fie albe ?
117
Fie A evenimentul de a scoate o bilă albă din urna U. Avem 10
cazuri egal posibile dintre care 7 sunt favorabile. Deci:
7
p .
10
Fie evenimentul de a scoate o bilă albă din urna U. Avem tot 10
cazuri egal posibile dintre care numai 3 sunt favorabile. Prin
urmare:
3
q .
10
Să observăm că, fie am scos o bilă albă din U, fie că
am scos o bilă roşie, numărul cazurilor egal posibile din urna V
rămâne tot 10, iar numărul cazurilor egal posibile favorabile tot
3. Evenimentele A şi B sunt independente şi deci putem aplica
regula de înmulţire a probabilităţilor.
Probabilitatea de a scoate doua bile albe este:
7 3 21
10 10 100
Calculul acestei probabilităţi se poate face şi fără aplicarea
teoremei înmulţirii probabilităţilor, însă este mai complicat.
2. Urna U conţine şapte bile albe şi trei roşii, iar urna V trei bile
albe şi şapte roşii. Extragem mai întâi o bilă din urna U şi
indiferent de culoarea ei o punem în urna V. Extragem apoi din
118
urna V o bilă. Care este probabilitatea să obţinem două bile
albe ?
Câmp de probabilitate
Mulţimea tuturor evenimentelor legate de o experienţă
împreună cu probabilităţile respective formează un câmp de
probabilitate.
119
Probabilităţile calculate se referă la evenimente legate
de experienţe având un număr finit de cazuri posibile
(evenimente elementare), deci am avut de-a face cu câmpuri
finite de probabilitate.
Probleme rezolvate
121
Probleme
123
5. O urnă conţine 3 bile albe şi 4 bile negre. Din această urnă
se scoate o bilă şi în locul ei se pune o bilă de cealaltă culoare
şi apoi se face o nouă extragere.
a) Care este probabilitatea ca a doua bilă extrasă să fie
neagră ?
b) Care este probabiltatea să obţinem bile de culori diferite ?
c) Care este probabilitatea ca a doua bilă să fie neagră ştiind
că prima bilă a fost albă ?
1) P( A U B) P( A) + P( B) P( A I B)
Fie n numărul cazurilor egal posibile ale experienţei în
raport cu care A şi B sunt evenimente, m numărul cazurilor
124
favorabile lui A şi s numărul cazurilor favorabile lui B. Să
presupunem că din cele m cazuri favorabile lui A, t sunt
favorabile lui A I B .
Numărul cazurilor favorabile lui AUB este m + s – t (şi nu
m + s) deoarece în acest caz , t cazuri ar fi numărate de 2 ori
şi la A şi la B).
Rezultă:
m s
P( A) ; P( B )
n n
m+ s t
P( A U B) ;
n
t
P( A I B)
n
125
şi deci relaţia
m+ s t m s t
+
n n n n
se poate scrie
P( A U B) P( A) + P( B) P( A I B)
126
Aplicaţie. O urnă conţine 3 bile albe şi 7 bile negre, iar alta
conţine 7 bile albe şi 3 negre. Din fiecare urnă se extrage câte
o bilă.Care este probabilitatea să obţinem cel putin o bilă albă ?
Solutie. Fie A evenimentul extragerii unei bile albe din prima
urnă şi B evenimentul extragerii unei bile albe din a doua urnă.
Avem de calculat probabilitatea evenimentului A B
P ( A B ) P ( A) + P ( B ) P ( A B ).
3 7
P ( A) , P( B)
10 10
Deoarece A şi B sunt independente
3 7 21
P( A B) P( A) P( B )
10 10 100
şi deci
3 7 21 79
P( A B) +
10 10 100 100
totalitatea lor.
Pentru a extrage k bile albe si n- k negre, trebuie
să se realizeze k evenimente Ai si n-k evenimente Ai .
Evenimentul
Ai1 Ai2 ... AiK AiK + 2 ... Ain
se realizează cu probabilitatea
pi1 pi2 ... piK qiK + 2 qiK + 2 ...qin
,
unde litera p apare de k ori cu diferiţi indici , iar q de n-k ori cu
indici care nu apar la p de k ori cu diferiţi indici, iar q de n-k ori
cu indici care nu apar la p.
Se observă uşor că după aceeaşi regulă se
calculează coeficientul lui xk in polinomul
P x p1 x + q1 p 2 x + q 2 ... p n x + q n .
Schema lui Poisson ajută la rezolvarea problemelor in care se
cere probabilitatea realizării de k ori a unui eveniment intr-o
128
experienţă ce consta in efectuarea a n experiente
independente, atunci când cunoaştem probabilitatea realizării
evenimentului in fiecare din cele n experienţe.
adică va fi egala cu C nk p k q n k .
129
Recunoaştem in această expresie termenul general al ridicării
la puterea n a binomului px + q. Pentru acest motiv schema se
mai numeşte binomială.
Deoarece urnele sunt identice, putem considera că
toate extragerile se fac dintr-o singura urnă, bila extrasa
punându-se înapoi in urna după fiecare extragere.
Obţinem astfel schema lui Bernoulli:
Probabilitatea de a scoate k bile albe din n extrageri dintr-o
urnă, punând de fiecare data bila înapoi, este
Pn , k C nk p k q n k ,
unde p este probabilitatea obţinerii unei bile albe dintr-o singură
extragere si q = 1 – p.
Schema lui Bernoulli rezolvă problemele in care se cere
probabilitatea realizării unui eveniment de k ori intr-o serie de n
efectuări a unei experienţe ,atunci când se cunoaşte
probabilitatea evenimentului la o singură efectuare a
experienţei.
Avem
1 1
p ; q ; n 4; k 1
2 2
130
1 3
1 1 1
P4 ,1 C 41 .
2 2 4
Probleme
131
3. O urnă conţine 12 bile numerotate cu 1, 2, …, 12. Se face o
extragere din aceasta urnă. Care este probabilitatea obţinerii
sau a unui număr par ,sau a unui număr mai mic ca 5, sau a
unui pătrat perfect ?
133
9.4. Variabile aleatoare
135
3 in timp ce ia valoarea 1, daca a iesit fata 6. Mai mult, atunci
când una din cele doua variabile ia valoarea 2, cealalta nu
poate lua aceeasi valoare.
Stim ca aruncarea zarului este o experienta care da
nastere la un tip de probabilitate. Am notat multimea
evenimentelor elementare prin E ={1,2,3,4,5,6}. Se vede ca
variabila X este o functie definita pe E, care ia urmatoarele
valori:
X({1})=1 ; X({2})=1; X({3})=2;
X({4})=2; X({5})=3; X({6})=3
La fel :
Y({1})= Y({6})=1 ;
Y({2})=Y({5})=2;
Y({3})=Y({4})= 3.
In general, la orice variabila aleatoare ne intereseaza
probabilitatea ca ea sa ia o anumita valoare. Dar pentru a
putea vorbi de probabilitate, trebuie sa avem in vedere un
camp de probabilitati. Vom spune deci ca o variabila aleatoare
este o functie definita pe multimea evenimentelor elementare
ale unui camp de probabilitate. Egalitatea de forma
x xi
este un eveniment. In exemplul de mai sus , egalitatea X=1
este echivalenta cu evenimentul {1,2}. Intr-adevar aceasta
egalitate se realizeaza, daca si numai daca se realizeaza acest
136
eveniment, cu alte cuvinte jucatorul respectiv capata un punct
daca si numai daca obtine una din fetele 1 sau 2.
Deoarece egalitatile x x1 , x x2 , …. x xn sunt
incompatibile doua câte doua dar si una din ele se realizeaza
neaparat, avem;
P x x1 + P x x2 + ... + P x xn 1 ,
p1 + p2 + ... + pn 1
137
�a+x a + x2 a + xn �
a + X :� 1 �.
� p1 p2 L pn �
138
x x2 xm �y y2 yn �
X : 1 , Y : �1 �,
p1 p2 ... pm �q1 q2 L qn �
z1 z2 zs
Z : ,
r1 r2 ... rs
putem scrie
x + y + z x2 + y 2 + z 2 xm + y n + z s
X + Y + Z : 1 1 1
q111 q112 ... q mns
Cand scriem tabloul de distributie al unei variabile aleatoare e
bine sa avem in vedere ca valorile din primul rând sa fie diferite
de câte doua.
139
unde p ij este probabilitatea realizarii simultane a egalitatilor
X xi , Y y j
Fiind date mai multe variabile aleatoare X,Y, .., V , vom numi
produsul lor variabila XY...V, care ia valoarea xiyi...vk.
De exemplu, fiind date 3 variabile aleatoare
x1 x2 xm �y y2 yn �
X : , Y : �1 �,
p1 p2 ... pm �q1 q2 L qn �
z z2 zs
Z : 1 .
r1 r2 ... rs
Atunci scriem:
x y z x2 y 2 z 2 xm y n z s
XYZ : 1 1 1 ,
p111 p112 ... p mns
unde p ijk este probabilitatea realizarii simultane a relatiilor
X=xi, Y=yj, Z=zk.
140
�x x2 xn �
X : �1 �,
�p1 p2 L pn �
Distributia variabilei Xr este
�x� x2� xn��
: �1
X� �,
�p1 p2 L pn �
142
9.6. Valori medii
Exemplu
a
1
si deci valoarea sa medie va fi egala cu a �
1 = a.
Intr-adevar , fie
�x x2 xn �
X : �1 �,
�p1 p2 L pn �
distributia variabilei X.
Distributia variabilei a+X este
144
�a+x a + x2 a + xn �
A+ X :� 1 �,
� p1 p2 L pn �
Valoarea medie
a + x1 p1 + a + x2 p2 + ... + a + xn pn
a p1 + p2 + ... + pn + p1 x1 + p2 x2 + ... + pn xn .
Fie
x x2 xm �y y2 yn �
X : 1 , Y : �1 �,
p1 p2 ... pm �q1 q2 L qn �
suma variabilelor aleatoare. Ne propunem sa calculam
valoarea medie a variabilei X+Y;
�x + y x2 + y2 xm + yn �
X + Y : �1 1 �
� p11 p12 L pmn �
146
X+Y= p11 (x1+y1) + p12 x1+y2)+…+ (x1+y2)+…+p1n(x1+yn)
+p21(x2+y1) + p22 (x2+y2) + … + p2n(x2+yn) + .... +
y1(p11+p21+…+pm1) +
+ y2(p12+p22+…+pm2) + … + yn(p1n+p2n+…+pmn).
Sa scriem aceasta relatie mai restrans
n
X +Y �x p
i 1
i i1 + pi 2 + ... + pin +
m
+ �y j p j1 + p j 2 + ... + p jn
j 1
147
sau pi= pi1+pi2+…+pin.
Rationand la fel se poate arata ca pj=pj1+pj2+…+pjm.
148
Fiind date doua variabile aleatoare
x1 x2 xm �y y2 yn �
X : , Y : �1 �
p1 p2 ... pm �q1 q2 L qn �
valoarea medie a variabilei
�x y x2 y2 xm yn �
XY : �1 1 �,
�p11 p12 L pmn �
este
M XY p11 x1 y1 + p12 x1 y2 + ... + p1n x1 yn + p21 x2 y1 +
+ p22 x2 y2 + ... + p2 n x2 yn + ... + pm1 xm y1 + pm2 xm y2 + ... pmn xm yn .
In cazul general, aceasta expresie nu se poate simplifica. Daca
insa variabilele aleatoare X si Y sunt independente, se produc
simplificari insemnate. Stim ca in acest caz putem scrie
pij= pi gj
M XY p1 g1 x1 y1 + p1 g 2 x1 y2 + ... + p1 g n x1 yn + p2 g1 x2 y1 +
+ p2 g 2 x2 y2 + ... + p2 g n x2 yn + ... + pm g n xm y1 + pm g 2 xm y2 + ... +
+ pm g n xm yn .
Dând factor p1x1 in prima linie, p2x2 in a doua, ..., pmxm in a m-a
rezulta
M(XY)= p1x1(g1y1+ g2 y2+...+ gnyn)+ p2x2(g1y1+ g2 y2+...+ gnyn)+...+
pmxm(g1y1+ g2 y2+...+ gnyn)
M(XY)= (g1y1+ g2 y2+...+ gnyn) (p1x1+ p2x2+...+ pmxm)
M(XY)= M(X) M(Y)
149
Valoarea medie a unui produs de doua variabile aleatoare
independente este egala cu produsul valorilor medii ale
variabilelor considerate.
Momente
Fiind data o variabila aleatoare X, vom numi
momentul de ordinul k al acestei variabile valoarea medie a
variabilei Xk. Vom scrie Mk(X)= M(Xk).
150
Fiind date o variabila aleatoare X si o constanta , vom numi
abatere de la constanta a variabilei X, variabila
X-. Momentul de ordinul k raportat la constanta al
variabilei X , se defineste ca valoare medie a variabilei
(X-)k.
Pentru .=0, obtinem momentul initial de ordin k, sau, scurt,
momentul de ordinul k al variabilei X pe care l-am definit la
inceputul acestui paragraf. Pentru .=M(X)=m, obtinem
momentul centrat de ordin k al variabilei X.
Variabila X- M(X) se numeste abaterea de la medie a variabilei
X, sau , scurt abaterea variabilei X. Valoarea medie a abaterii
oricarei variabile aleatoare este nula .
Intr-adevar M(X- M(X))=M(X)-M(X)=0
De multe ori, la o variabila aleatoare ne intereseaza cat de mult
se abat valorile variabilei de la valoarea medie.Trebuie sa
stabilim un indicator numeric al imprastierii valorilor variabilei
aleatoare in jurul valorilor medii. Valoarea medie a abaterii de
la medie nu poate caracteriza aceasta imprastiere, deoarece
este nula pentru orice variabila aleatoare. Abaterile diferitelor
valori, avand semne diferite, se compenseaza reciproc. Este
foarte important sa caracterizam imprastierea variabilei X, prin
valoarea medie
a abaterii absolute | X-M(X)| pe care o vom numi abatere
medie.
151
Daca X are distributia
�x x2 xn �
X : �1 �,
�p1 p2 L pn �
atunci distributia abaterii absolute este
�x1 m x2 m xn m �
� �,
� p1 p2 L pn �
unde m=M(X) , iar abaterea medie este
p1|x1-m| + p2|x2-m| +...+ pn |xn-m|.
Se observa ca semnul expresiilor x i-m nu influenteaza
valoarea abaterii medii. Putem lua ca indicator al imprastierii
variabilei X si M[(X-m)2] sau M[(X-m)4] etc. Folosirea abaterii
medii este foarte incomoda in calcul.
Foarte comoda este in schimb folosirea expresiei
M[(X-m)2]
Dispersia
Vom numi dispersie a unei variabile aleatoare X
momentul centrat de ordinul al doilea al acestei variabile.Vom
scrie
2=D2(X)= M[(X-m)2]
unde m=M(X).
Daca X are distributia
152
�x x2 xn �
X : �1 �,
�p1 p2 L pn �
atunci
D2(X)= p1(x1-m)2 + p2(x2-m)2 +...+ pn (xn-m)2=
= p1(x21-2 x1m+m2) + p2(x22-2 x2m+m2) +...+
+pn (x2n-2 xnm+m2)= p1x21+ p2x22 +...+ pnx2n+
+ 2 m( p1x1+ p2x2 +...+ pnxn)+ m2 (p1+ p2+...+ pn)=
=M(X2) – [M(X)]2.
154
D(aX)=aD(X).
Aceste proprietati le are si trebuie sa le aiba orice marime
menita sa masoare gradul de imprastiere. In plus , abaterea
medie patratica are proprietatea ca daca
1=D2(X1), 2=D2(X2) ; ... ; n=D2(Xn)
unde X1, X2 ,…,Xn, sunt variabile aleatoare independente,
atunci
1 D 2 X 1 + X 2 + ... + X n 12 + 22 + ... + n2 .
Probleme
155
�0 1 �
4. Se da variabila aleatoare X :� �. Sa se faca
�0,3 0, 7 �
tabloul pentru Xn (n �N).
�1 0 1 �
5. Se dau variabilele X :� �;
�0,3 0,3 0, 4 �
�1 1 �
Y :� �. Sa se faca tabloul de distributie pentru X si
n
�0,5 0,5 �
Yn (n N).
6. Distributia variabilei X este:
�1 2 3 4�
X :� 2 7 1 1�.
�
�p p �
�
� 4 3 6�
Care este probabilitatea ca X sa ia o valoare mai mica sau
egala cu 3 ?
156
8. Daca variabilele de la problema 5 sunt independente ce
distributie are suma patratelor lor ?
157
�1 1 1 �
�1�
2 2� 3 n n 1 �
X :� �
�1 1 1 �
� L �
�n n n �
158
9.8. Elemente de statistică matematică
159
Exemple:
161
nx
f x ,
n
unde f(x) este frecvenţa relativă a valorii x, nx este frecvenţa
absolută a acestei valori, iar n efectivul total al populaţiei.
Deseori, frecvenţa este dată în procente.
În cazul caracteristicilor cantitative, aceste tabele scot
în evidenţă o corespondenţă între două mulţimi de numere:
mulţimea valorilor caracteristicii şi mulţimea frecvenţelor
corespunzătoare. Se poate remarca analogia cu distribuţia
variabilelor aleatoare.
Dacă distribuţia unei variabile aleatoare o scriem
x1 x2 . . xn
p1 p2 . . pn
unde pe prima linie sunt trecute valorile variabilei,
iar în cea de a doua linie probabilităţile
corespunzătoare acestor valori, tabela
corespunzătoare unei caracteristici cantitative o
scriem ca în tabela următoare.
162
Valoarea Frecvenţa
x1 f1
x2 f2
. .
. .
. .
xn fn
În prima coloană sunt trecute valorile caracteristicii, iar
în cea de a doua frecvenţele corespunzătoare.
De altfel, noţiunile de variabilă aleatoare şi
probabilitate sunt modelele teoretice ale noţiunilor de
caracteristică şi, respectiv frecvenţă, care sunt noţiuni cu
caracter experimental.
De multe ori, în loc de caracteristică se spune
variabilă statistică sau pe scurt variabilă, şi noi vom folosi în
cele ce urmează şi această denumire.
Vom spune că o tabelă de genul precedentei defineşte
distribuţia sau repartiţia statistică a variabilei respective.
Este uşor de observat că suma frecvenţelor relative ale tuturor
valorilor variabilei este 1.
Se numeşte frecvenţă absolută cumulată
crescătoare a unei valori x a variabilei suma frecvenţelor
absolute ale tuturor valorilor variabilei care apar până la x
inclusiv.
163
Se numeşte frecvenţă absolută cumulată
descrescătoare a unei valori x a variabilei suma frecvenţelor
absolute ale tuturor valorilor care apar de la x inclusiv. (în cazul
caracteristicilor cantitative vom considera de aici înainte numai
tabele în care valorile variabilei sunt scrise în ordine
crescătoare).
În mod analog, definim frecvenţa relativă (sau scurt,
frecvenţa) cumulată crescătoare şi frecvenţa cumulată
descrescătoare.
Se numeşte frecvenţă cumulată crescătoare a unei
valori x a variabilei, suma tuturor frecvenţelor valorilor care
apar până la x inclusiv, iar frecvenţa cumulată
descrescătoare suma tuturor frecvenţelor valorilor care apar
de la x inclusiv. Rezultă că frecvenţa cumulată crescătoare este
exprimată şi de raportul dintre frecvenţa absolută cumulată
crescătoare şi efectivul total al populaţiei, iar frecvenţa
cumulată descrescătoare de raportul dintre frecvenţa absolută
cumulată descrescătoare şi efectivul total al populaţiei.
În cele ce urmează, prin frecvenţă cumulată vom
înţelege frecvenţa relativă cumulată crescătoare. Este uşor de
observat că frecvenţa cumulată are ca model teoretic funcţia de
repartiţie. O tabelă statistică poate fi completată cu aceste
frecvenţe.
Serii cronologice
164
Tot în cadrul seriilor statistice sunt incluse şi aşa-
numitele serii cronologice care prezintă evoluţia în timp a
unor mărimi.
În cazul unei tabele corespunzătoare unei serii
cronologice, în prima coloană sunt trecute anumite momente
sau intervale de timp, iar în coloana a doua valorile
corespunzătoare ale mărimii considerate.
a. Reprezentarea în batoane.
166
Această reprezentare se foloseşte mai ales pentru
seriile în care variabila ia un număr mic de valori. De exemplu :
b. Histogramă.
Fiind dată o serie cu clase de valoare de amplitudini
egale obţinem histograma acestei serii, luând pe axa orizontală
167
o succesiune de segmente egale (reprezentând amplitudinea
claselor) şi ridicând, pe fiecare din aceste segmente
considerate ca baze, dreptunghiuri de înălţimi proporţionale cu
frecvenţele (relative sau absolute) ale claselor respective.
25
20
15
10
0
1 2 3 4 5 6 7 8
Figura a
20
10
0
10-20 20-30 30-40 40-50 50-60
Figura b
168
Tabela 1. Distribuţia unor piese după diametrul lor
Mărimea Frecvenţa Frecvenţa
diametrului absolută cumulată
[mm]
10–20 10 10
20–30 15 25
30–40 12 37
40–50 15 52
50–60 8 60
169
e. Poligonul frecvenţelor cumulate (crescătoare) se obţine
unind printr-o linie poligonală punctele (x,y), unde x este
extremitatea dreaptă a intervalului unei clase, iat y frecvenţa
cumulată a clasei respective, la care mai adăugăm punctul
(0,a), unde a este limita inferioară a primei clase.
Figura c
170
Figura d
În mod analog, se defineşte curba frecvenţelor cumulate
descrescătoare. Dacă punctele din figură le unim nu printr-o
linie poligonală, ci printr-o curbă, obţinem curba cumulativă a
seriei considerate.
1. Modul
Modul sau dominanta unei serii statistice se numeşte
valoarea caracteristicii corespunzătoare celei mai mari
frecvenţe, în cazul când valorile caracteristicii sunt date
171
individual şi valoarea centrală a clasei corespunzătoare celei
mai mari frecvenţe, în cazul variabilelor continue, când se dau
clase de variaţie. Această noţiune prezintă interes mai ales în
cazurile când avem o singură dominantă.
2. Mediana
Mediana unei serii este un număr x astfel încât există
tot atâtea unităţi statistice corespunzătoare valorilor <x, ca şi
cele corespunzătoare valorilor >x.
Dacă o caracteristică ia valorile 1, 3, 3, 3, 4, 5, 6, 6, 7,
8, 9 atunci 5 este mediana, deoarece există 5 valori <5 şi 5
valori >5. Dacă avem valorile 1, 3, 3, 3,4, 5, 6, 7, 7, 9 atunci
vom lua ca mediană media aritmetică a numerelor situate la
mijloc (dacă ele au fost scrise în ordinea mărimii). În acest caz,
mediana este 4,5.
Uneori se consideră ca mediană oricare din cele două
numere. Cum se calculează mediana în cazul unei variabile
continue, vom arăta pe un exemplu.
Să considerăm pentru aceasta tabela e. Dacă piesele
ar fi aranjate în ordinea diametrelor lor, noi vrem să calculăm
diametrul celei de-a 30-a.
Diametrul acestei piese este cuprins între 30 şi 40mm.
Clasa 30–40 are frecvenţa absolută 12.
172
Vom presupune că diametrul celor 12 piese
corespunzătoare creşte uniform de la 30 la 40. Deci creşterea
diametrului de la o piesă la următoarea este
40 30
12 .
Pe de altă parte, a 30-a piesă a populaţiei este a 30-25= a 5-a
piesă a clasei (deoarece există 25 de piese cu diametrul <30).
Deci, diametrul celei de-a 30-a piese este
40 30
30 + ( 30 25 ) 34 ,16 mm
12
3. Media aritmetică
Dacă x1 , x2 , ... , xn sunt n valori, se ştie că media lor
x1 + x2 + ... + xn
aritmetică este .
n
Fiind dată distribuţia unei variabile x
Valori Frecvenţe
x1 y1
x2 y2
. .
. .
. .
xn yn
Valoarea medie a variabilei respective este
173
x1 y1 + x2 y2 + ... + xn yn
x
y1 + y2 + ... + yn (1).
Dacă N = y1 + y2 + … + yn este efectivul total al populaţiei,
atunci
y1 y y
x x1 + x2 2 + ... + xn n
N N N
yi
sau, dacă notăm cu f i frecvenţa relativă a valorii xi
N
(i=1, 2, …, n)
x x1 f 1 + x2 f 2 + ... + xn f n
expresia (1) are într-adevăr semnificaţia unei medii aritmetice.
Variabila x ia, după cum reiese din tabelă, de y 1 ori valoarea x1,
de y2 ori valoarea x2 , ş.a.m.d.
Deci pentru a calcula valoarea medie a variabilei, calculăm
media aritmetică a numerelor
x , x ,..., x x , x ,..., x ... x n , x n ,..., x n
1 1 1 2 2 2
y1 ori y 2 ori y n ori
174
Astfel, în exemplul prezentat în tabela 1, se obţin datele din
tabela 2.
Tabela 2
Mărimea Frecvenţa Valoarea xiyi
diametrului absolută yi centrală xi
10-20 10 15 150
20-30 15 25 375
30-40 12 35 420
40-50 15 45 675
50-60 8 55 440
60 2060
2060
x 34 ,3
60
Dacă vrem să lucrăm cu numere mai mici decât cele ce ne sunt
date în tabele, facem următoarele observaţii.
Avem pentru orice i
xi x0 + ( xi x0 ),
xi yi x0 yi + ( xi x0 ) yi
Dând lui i valorile 1, 2, ..., n obţinem n relaţii care adunate
termen cu termen ne da
x1 y1 + x2 y2 + x1 y1 + ... + xn yn x0 ( y1 + y2 + ... + yn ) +
+( x1 x0 ) y1 + ( x2 x0 ) y2 + ... + ( xn x0 ) yn
175
sau împărţind cu N = y1 + y2 + … + yn,
( x1 x0 ) y1 + ( x2 x0 ) y2 + ... + ( xn x0 ) yn
x x0 + .
N
Această relaţie o mai putem scrie
x x0 x x0 .
În exemplul precedent, luând x0 = 35, avem calculele
prezentate în tabela
Tabela 3.
Clase Frecvenţa Valoarea xi x0 xi x0 yi
absolută centrală
yi
10-20 10 15 -20 -200
20-30 15 25 -10 -150
30-40 12 35 0 0
40-50 15 45 10 150
50-60 8 55 20 160
60 -40
176
40
x 35 34 ,3
60
4. Dispersia
Fiind date n valori x1, x2, ..., xn a căror medie este x, se
numeşte dispersia acestor valori mărimea
( x1 x )2 + ( x2 x )2 + ... + ( xn x )2
2
n
Fiind dată seria statistică
177
x1 y1 + x2 y2 + ... + xn yn
unde N = y1 + y2 + ... + yn şi x
N
sunt respectiv efectivul total al populaţiei şi valoarea medie,
dispersia corespunzătoare este
( x1 x )2 y1 + ( x2 x )2 y2 + ... + ( xn x )2 yn
2
N
Mărimea 2 se numeşte abaterea medie pătratică. Ea
178
Dacă am fi înlocuit mărimile xi prin xi-x0 unde x0 este o
constantă, am fi obţinut
( x1 x0 )2 y1 + ( x2 x0 ) 2 y2 + ... + ( xn x0 ) 2 yn
2
N
( x x0 ) 2 (2’)
Unde prin x2 am notat media aritmetică a mărimilor x12, x22, ...,
xn2 cu ponderile y1, y2, ..., yn. Am regăsit o proprietate a
dispersiei unei variabile aleatoare.
Dispersia sau mai bine zis abaterea medie pătratică
indică gradul de împrăştiere a valorilor în jurul valorii medii. O
valoare mică a abaterii indică o pronunţată grupare a valorilor
în jurul mediei aritmetice.
Să calculăm tabela h cu datele necesare calcului
dispersiei 2 . Vom calcula 2 atât direct, cât şi folosind
formula (2’).
Luăm x0=35 şi ştim că x 34,3 . Se obţin datele
din tabela 4.
179
Tabela 4.
Clase
Valoarea centrală xi
yi Frecvenţa
xi x 0 xi x ( xi x 0 )
2 2 2
( xi x ) yi ( xi yxi0 ()xi x )
2
10–20
180
30–40
60 1020 10172,5
0 0
181
Deci
10172,5
2 169 ,5
60
Sau folosind (2’)
1
2 10200 ( 34 ,3 35 )2 169 ,5
60
Se observă că alegând convenabil valoarea lui x0, calculele se
simplifică.
Generalitati
t 0 1 2 3
s 0 0.49 1.96 4.41
Graficul corespunzator este cel din figura h.
182
Daca t este masurat in secunde, iar g in 10 m/s2,
atunci functia de mai sus reprezinta legea caderii corpurilor
in vid. Daca vrem sa verificam experimental aceasta lege,
lasam un corp sa cada si vom masura spatiile parcurse cu S 1,
S2 , S 3, …corespunzatoare momentelor t 1, t 2, t 3, …
Corespondenta intre marimile experimentale t si S defineste o
functie empirica. In general, nu vom obtine la fiecare efectuare
a experientei aceleasi valori pentru S corespunzatoare valorilor
alese pentru t, ele schimbandu-se in raport cu conditiile in care
se efctueaza experienta.
s
t
Figura e
In tabela k sunt trecute datele obtinute experimental.
183
Tabela 6
t (s) 1 2 3 4
S(metri) 0.44 2.1 4.6 7.5
. t
Figura f
184
In fig. i sunt trecute, in raport cu acelasi sistem de axe,
punctele obtinute experimental si graficul functiei teoretice.
Distributiile intalnite mai inainte, ca distributia binominala,
distributia normala etc. sunt distributii teoretice.
Sa presupunem ca se arunca o moneda de 4 ori. Probabilitatea
obtinerii “stemei “ de doua ori si a “banului” de doua ori este
2 2
�1 ��1�
C42 � �� � 0.375 .
�2 ��2�
Acest numar reprezinta pentru noi frecventa teoretica, numarul
in jurul caruia oscileaza frecventele experimentale, daca facem
de un numar mare de ori experienta, in ipoteza ca moneda este
“corecta”. Daca vom efectua 1000 de serii de cate 4 aruncari,
numarul seriilor in care a aparut fiecare fata de 2 ori va fi exact
375, daca moneda este perfect simetrica, va fi un numar
apropiat de 375.
In general, pentru ca o distributie teoretica sa aproximeze o
distributie empirica, trebuie ca frecventele experimentale sa nu
difere mult de cele teoretice. Aceasta face ca si elementele
caracteristice ale celor doua distributii sa ia valori apropiate.
1. AJUSTAREA LINIARA
Tabela 7
X 0 1 2 3 4 5
Y 1.2 1.6 1.7 2.4 2.8 3.8
186
In figura k sunt reprezentate punctele
corespunzatoare tabelei m si dreapta (d).
Figura g
188
AB CD
a= (2)
AC B 2
Se observa ca S rezulta din P, luand
A = Σxi2, B = Σxi, C = n,
D = -Σxiyi, E = -Σyi, F = Σyi2
Se observa ca C = n > 0.
Relatia AC-B2 = nΣ2-(Σxi)2 0 este o inegalitate cunoscuta si
se demonstreaza usor prin inductie.
Rezulta, tinand cont de relatiile (1) si (2) ca S este minim cand
nSxi yi Sxi Syi
a= ,
nSxi 2 (Sxi ) 2
a Σxi+nb-Σyi = 0.
Daca notam :
Σxi = n x ; Σyi = n y
remarcam ca relatia (4) ne asigura ca punctul (x ; y) se gaseste
pe dreapta :
y = ax + b.
189
Sa gasim, folosind metoda celor mai mici patrate, dreapta de
ajustare pentru tabela m. Avem datele din tabela n.
Tabela 8
xi yi xi2 xi yi
0 1.2 0 0
1 1.6 1 1.6
2 1.7 4 3.4
3 2.4 9 7.2
4 2.8 16 11.2
5 3.8 25 19
Total 15 13.5 55 43.4
Avem:
15 13,5
x 2,5 ; x 2, 25 ; n=6
6 6
42, 4 33, 75 8, 65
a 0, 49 .
55 37,5 17,5
Tinand cont de relatia (4) rezulta:
b y ax 2, 25 0, 49 �
2,5 1, 025
Vom arata ca daca inlocuim xi cu xi-x0 si yi cu yi-yo, valoare lui a
nu se schimba. Sa facem aceste schimbari la numaratorul
expresiei care da pe a. Avem:
190
� x x y y n x x y y
i 0 0 i 0 0
�x y x �y y �x + nx y nx �
i i 0 i y+
0 i 0 0
+ nx y + ny x + ny x nx y �x y nx �y.
0 0 0 0 0 i i
Am tinut cont ca nx �x ; ny �y . i i
191
ANEXA
independente.
P A �B
P A | B (2)
P B
Spaţiul de selecţie micşorat este B (evenimentul de
condiţionare).
Toate probabilităţile se referă la întregul spaţiu de
selecţie, P(A) fiind o prescurtare pentru P(A | S), unde S este
întreg spaţiul de selecţie.
încât P A1 ��ǹ
A2 .... An 0, adică probabilitatea
193
X �A1 , X �A2 , X �A3 ,
fiind incompatibile. Urmează că
P X X �A1 + X �A2 + X �A3 şi
ţinând seama de definiţia probabilităţii condoţionate
P X �Ai P Ai P X | Ai , i 1, 2,3 ,
obţinem rezultatul dorit.
(4) poartă numele de formula probabilităţilor totale.
Formula Bayes
194
numite probabilităţile aposteriori (deoarece se calculează
după efectuarea experienţei).
�P A P X | A
(5)
j j
1
195
196
BILETE pentru EXAMEN
Biletul 1
Subiectul 1
a) Calculaţi cu formula sau teorema lui Cauchy:
z
cosh
I
C
2 dz, C : z + 2i 2
z + i
b) Cu teorema reziduurilor, să se calculeze:
+
x cos 2 x
I
x 2 2x + 5dx
Subiectul 2
a) Reprezentaţi printr-o integrală Fourier funcţia neperiodică:
�1 dacă t a
�
f (t ) �
1/ 2 dacă t �a
�0 dacă t a
�
197
e ax e bx
0
x
sin xdx
Subiectul 3
Determinaţi forma canonică şi integraţi ecuaţia:
2u 2u 2u u
+ 2 cos x sin 2 x sin x 0
x 2 xy y 2 y
u
u ( x , y) x4, x.
y sin x y y sin x
care trece prin cercul de ecuaţii z 0, x 2 + y 2 1.
198
Biletul 2
Subiectul 1
a) Calculaţi cu formula sau teorema lui Cauchy:
i z
ze 2
I
z R
z2
dz,
pentru R diferit de 2.
b) Utilizând teorema reziduurilor, calculaţi integrala:
cos n
I
b ia cos d,
0
a 0, b 0
Subiectul 2
a) Reprezentaţi printr-o integrală Fourier funcţia:
0, t ��,�ȥ
� a 0 a,
�
�1, t � a, 0
f (t ) �
�a, t � 0, a
��1, t �a
�
199
b) Să se arate că
sin 2 u
I
0
u2
du
2
Subiectul 3
Determinaţi forma canonică şi integraţi ecuaţia:
l2 x 2 2u
2
x
u
y
2u
0, 0
x y 2
x u
u(x, y) arcsin ,
y 0 l y y 0
200
Biletul 3
Subiectul 1
a) Calculaţi cu formula sau teorema lui Cauchy:
cosh z x2 y2
I
C
z 2
+ 1
dz, C :
a 2
+
b 2
1, b 1, a, b 0
Subiectul 2
a) Reprezentaţi printr-o integrală Fourier funcţia neperiodică:
�1, t �1
�
t , t � 1, 0
f (t ) �
�
�t , t � 0,1
Subiectul 3
201
Determinaţi forma canonică şi integraţi ecuaţia:
2u 2u 2u
4 +8 21 0
x 2 xy y 2
u ( x , y) sin x , u
2y 7x 0 2 y + 3x
202
Biletul 4
Subiectul 1
a) Calculaţi cu formula sau teorema lui Cauchy:
z100 e iz
I
C
z 2
+ 1
dz,
unde C : 4x 2 + y 2 4 0.
b) Utilizând teorema reziduurilor, calculaţi integrala:
sin x sin nx
I
5 cos x
dx
Subiectul 2
a) Rezolvaţi ecuaţia integrală:
I
0 g(u ) cos tudu ( t ),
1 t , t � 0,1
�
�
unde (t ) �
�0, t � 1,
203
Subiectul 3
Determinaţi forma canonică şi integraţi ecuaţia:
2u 2u 2u u
x2 2 xy + y2 + x + y
x 2 x y y 2 x
u
u ( x , y) 1 cos y, 2 y.
x 1 x x 1
204
Biletul 5
Subiectul 1
a) Calculaţi cu formula sau teorema lui Cauchy:
2 +1
ez
C
z3
dz,
unde C : z 2i 1.
b) Utilizând teorema reziduurilor, calculaţi integrala:
cos ax
I
x
0
2
+ b2 2
dx
Subiectul 2
a) Rezolvaţi ecuaţia integrală:
1
f ( t ) cos tx dt x
0
2
+ a2
x 0, a 0
Subiectul 3
205
Determinaţi forma canonică :
2u 2u 2u u u
2 + +2 0.
x 2 xy y 2 x y
206
1. N.Boboc, Funcţii complexe, Editura Didactică şi
Pedagogică, Bucureşti, 1969.
2. R. Cristescu, Matematici superioare, Editura Didactică
şi Pedagogică, Bucureşti, 1956.
3. R. Cristescu, Ordered vector spaces and linear
operators, Editura Academiei and Abacus Press, 1977.
4. R. Cristescu, Analiză funcţională, Editura Didactică şi
Pedagogică, 1979.
5. Gh. Grigore, Analiza numerică, Bucureşti, 1984.
6. A. Haimovici, Ecuaţiile fizicii matematice şi elemente
de calcul variaţional, Editura Didactică şi Pedagogică,
Bucureşti, 1966.
7. O.Onicescu, Gh. Mihoc, Probabilitati si statistica
matematica, Editura Didactică şi Pedagogică, Bucureşti,
1977.
8. O. Mayer, Teoria funcţiilor de o variabilă complexă,
Vol.I,II, Editura Academiei, Bucureşti, 1990.
9. D. Popovici, I.M. Popovici, Analiză matematică –
Structuri fundamentale, Editura Muntenia, 1993.
10. V. Smirnov, Cours de mathématiques supérieurs,
Tome IVF, Editions Moscou, 1975.
11. http://www.mathword.wolfram.com
12. Prob. and Statistics, Calculus and Analysis (CD)
13. http://apollo.imc.ro (CD)
207
14. O. Lipovan, B. Crstici şi alţii, Matematici speciale,
Editura Didactică şi Pedagogică, Bucureşti, 1981.
15. Gh. Marinescu, Spaţii vectoriale normate, Editura
Academiei, Bucureşti, 1956.
16. C. Niculescu, N. Popa, Elemente de teoria spaţiilor
Banach, Editura Academiei, Bucureşti, 1981.
17. D.M. Popovici, I.M. Popovici, Algoritmi pentru metode
numerice, I.M.C., 1994.
18. V. Smirnov, Cours de mathématiques supérieurs,
Tome IVF, Editions Moscou, 1975.
19. H.H. Schaefer, Topological vector spaces, MacMillan,
New-York, 1966.
20. C. Zaman, Econometrie, Pro Democraţia, Bucureşti,
1998.
21. http://www.stieltjes.org/research/proj31.html
22. http://www.stieltjes.org/research/proj32.html
23.http://www.math.mcgill.ca/~malkoun/Escher_Circle_Li
mit_III/Escher_Circle_Limit_III.html
208
Cuprins
209
Tiparul executat în
Tipografia
UNIVERSITĂŢII MARITIME
Constanţa
210