Salta al contenuto
Note per Studenti Formulario · Fondamenti di Analisi Matematica 2

FormularioFondamenti di Analisi Matematica 2: definizioni, teoremi e formule delle note, in ordine di capitolo

In questa pagina 14

1. Calcolo infinitesimale per le curve

Curve parametriche, vettore tangente e integrali di funzioni vettoriali

Definizione (prodotto scalare, norma, distanza). Per x‾,y‾∈Rn\underline x, \underline y \in \mathbb{R}^n: x‾⋅y‾=∑i=1nxiyi,∣x‾∣=x‾⋅x‾=x12+⋯+xn2,dist⁡(x‾,y‾)=∣x‾−y‾∣.\underline x \cdot \underline y = \sum_{i=1}^n x_i y_i, \qquad |\underline x| = \sqrt{\underline x \cdot \underline x} = \sqrt{x_1^2 + \dots + x_n^2}, \qquad \operatorname{dist}(\underline x, \underline y) = |\underline x - \underline y|. La norma misura la lunghezza del vettore. Vale ∣λx‾∣=∣λ∣ ∣x‾∣|\lambda \underline x| = |\lambda|\,|\underline x|. Se x‾≠0‾\underline x \neq \underline 0, il versorevettore di lunghezza 1 che indica una direzione e un verso x‾/∣x‾∣\underline x / |\underline x| ha la stessa direzione e verso di x‾\underline x e lunghezza 11.

Esempio. x‾=(1,2,2)\underline x = (1, 2, 2): ∣x‾∣=1+4+4=3|\underline x| = \sqrt{1 + 4 + 4} = 3, versore (1/3,2/3,2/3)(1/3, 2/3, 2/3). Per y‾=(2,−1,0)\underline y = (2, -1, 0): x‾⋅y‾=2−2+0=0\underline x \cdot \underline y = 2 - 2 + 0 = 0, quindi x‾⊥y‾\underline x \perp \underline y.

Teorema (disuguaglianza di Cauchy-Schwarz). ∣x‾⋅y‾∣≤∣x‾∣ ∣y‾∣|\underline x \cdot \underline y| \le |\underline x|\,|\underline y| per ogni x‾,y‾∈Rn\underline x, \underline y \in \mathbb{R}^n.

Teorema (disuguaglianza triangolare) (D). ∣x‾+y‾∣≤∣x‾∣+∣y‾∣|\underline x + \underline y| \le |\underline x| + |\underline y|. Di conseguenza dist⁡(x‾,y‾)≤dist⁡(x‾,z‾)+dist⁡(z‾,y‾)\operatorname{dist}(\underline x, \underline y) \le \operatorname{dist}(\underline x, \underline z) + \operatorname{dist}(\underline z, \underline y).

Esempio. (1,0,2)×(0,1,1)=(0⋅1−2⋅1,  2⋅0−1⋅1,  1⋅1−0)=(−2,−1,1)(1,0,2) \times (0,1,1) = (0\cdot1 - 2\cdot1,\; 2\cdot 0 - 1\cdot 1,\; 1\cdot 1 - 0) = (-2, -1, 1); infatti (−2,−1,1)⋅(1,0,2)=0(-2,-1,1)\cdot(1,0,2) = 0 e (−2,−1,1)⋅(0,1,1)=0(-2,-1,1)\cdot(0,1,1) = 0.

Definizione (limite). Sia t0t_0 un punto di accumulazione di AA e ℓ‾∈Rm\underline\ell \in \mathbb{R}^m. Si scrive lim⁡t→t0r‾(t)=ℓ‾\lim_{t \to t_0}\underline r(t) = \underline \ell se ∀ε>0 ∃δ>0:t∈A, 0<∣t−t0∣<δ  ⟹  ∣r‾(t)−ℓ‾∣<ε,\forall \varepsilon > 0\ \exists \delta > 0:\quad t \in A,\ 0 < |t - t_0| < \delta \implies |\underline r(t) - \underline\ell| < \varepsilon , cioè se ∣r‾(t)−ℓ‾∣→0|\underline r(t) - \underline \ell| \to 0: la distanza tra il punto mobile e ℓ‾\underline\ell tende a zero.

Proposizione (limite per componenti). lim⁡t→t0r‾(t)=ℓ‾  ⟺  lim⁡t→t0rj(t)=ℓj\lim_{t\to t_0}\underline r(t) = \underline\ell \iff \lim_{t\to t_0} r_j(t) = \ell_j per ogni j=1,…,mj = 1, \dots, m.

Esempio. lim⁡t→0(sin⁡tt,  log⁡(1+t2),  et)=(1,0,1)\lim_{t \to 0}\big(\tfrac{\sin t}{t},\; \log(1 + t^2),\; e^t\big) = (1, 0, 1): ogni componente ha il suo limite (Limiti notevoli di funzioniI limiti notevoli per le forme 0/0 (log(1+x)/x, (e^x − 1)/x, ((1+x)^α − 1)/x, sin x/x, (1 − cos x)/x², arctan x/x...) con le dimostrazioni, e la loro forma come asintoticità: log(1+f) ~ f, e^f − 1 ~ f, sin f ~ f, 1 − cos f ~ f²/2 quando f → 0.Limiti notevoli di funzioni →).

Definizione (curva, sostegno). Una curva parametrica (o cammino) è una funzione continua r‾:I→Rm\underline r : I \to \mathbb{R}^m, con II intervallo. Il sostegno è l'immagine Γ=r‾(I)={r‾(t):t∈I}\Gamma = \underline r(I) = \{\underline r(t) : t \in I\}. La curva si indica con γ={r‾,Γ}\gamma = \{\underline r, \Gamma\}; r‾\underline r è una parametrizzazione di Γ\Gamma.

Esempio. r‾(t)=(cos⁡t,sin⁡t)\underline r(t) = (\cos t, \sin t) su [0,2π][0, 2\pi]: sostegno la circonferenza unitaria. La stessa circonferenza è il sostegno di s‾(t)=(cos⁡2t,sin⁡2t)\underline s(t) = (\cos 2t, \sin 2t) su [0,π][0, \pi] (stesso giro in metà tempo) e di (cos⁡t,sin⁡t)(\cos t, \sin t) su [0,4π][0, 4\pi] (due giri). Sono tre curve diverse con lo stesso sostegno, e avranno lunghezze diverse (Lunghezza di una curva e ascissa curvilineaLa lunghezza di una curva è l'estremo superiore delle lunghezze delle poligonali inscritte. Se r è C^1 la curva è rettificabile e ℓ = ∫|r'(t)|dt (grafici: ∫√(1+f'^2); polari: ∫√(f²+f'^2)). La lunghezza non cambia riparametrizzando, anche invertendo il verso. L'ascissa curvilinea s(t)=∫|r'| dà la parametrizzazione a velocità 1: T = dr/ds. Con essa si definiscono curvatura κ=|dT/ds| e normale N; l'accelerazione è a = v' T + κ v² N.Lunghezza di una curva e ascissa curvilinea →).

Definizione (curva chiusa, curva semplice). Sia r‾:[a,b]→Rm\underline r : [a,b] \to \mathbb{R}^m una curva.

  • È chiusa se r‾(a)=r‾(b)\underline r(a) = \underline r(b).
  • È semplice se non passa due volte per lo stesso punto: r‾(t1)=r‾(t2)\underline r(t_1) = \underline r(t_2) con t1≠t2t_1 \neq t_2 è possibile solo per {t1,t2}={a,b}\{t_1, t_2\} = \{a, b\}.

Esempio. (cos⁡t,sin⁡t)(\cos t, \sin t) su [0,2π][0, 2\pi] è chiusa e semplice. Su [0,4π][0, 4\pi] è chiusa ma non semplice (passa due volte per ogni punto). Un "otto" percorso senza staccare la penna è chiuso e non semplice; un segmento è semplice e non chiuso.

Grafico interattivo: Chiusa e semplice (cerchio), semplice non chiusa (spirale di Archimede), chiusa non semplice (lemniscata, 'otto')

Definizione (derivata). r‾\underline r è derivabile in t0t_0 (interno) se esiste finito r‾′(t0)=lim⁡h→0r‾(t0+h)−r‾(t0)h.\underline r'(t_0) = \lim_{h \to 0}\frac{\underline r(t_0 + h) - \underline r(t_0)}{h} . Per il limite per componenti, r‾\underline r è derivabile se e solo se lo sono tutte le componenti, e r‾′(t)=(r1′(t),…,rm′(t))\underline r'(t) = (r_1'(t), \dots, r_m'(t)). Si scrive r‾∈C1\underline r \in C^1 se r‾′\underline r' esiste ed è continua; la derivata seconda è r‾′′=(r1′′,…,rm′′)\underline r'' = (r_1'', \dots, r_m'').

Esempio. r‾(t)=(t2,sin⁡t,e2t)\underline r(t) = (t^2, \sin t, e^{2t}): r‾′(t)=(2t,cos⁡t,2e2t)\underline r'(t) = (2t, \cos t, 2e^{2t}) e r‾′(0)=(0,1,2)\underline r'(0) = (0, 1, 2).

Definizione (curva regolare). r‾:I→Rm\underline r : I \to \mathbb{R}^m è regolare se r‾∈C1(I)\underline r \in C^1(I) e r‾′(t)≠0‾\underline r'(t) \neq \underline 0 per ogni t∈It \in I. Allora è definito il versore tangente T‾(t)=r‾′(t)/∣r‾′(t)∣\underline T(t) = \underline r'(t)/|\underline r'(t)| e la retta tangente in r‾(t0)\underline r(t_0) è x‾(λ)=r‾(t0)+λ r‾′(t0),λ∈R.\underline x(\lambda) = \underline r(t_0) + \lambda\,\underline r'(t_0),\qquad \lambda \in \mathbb{R}. Una curva continua è regolare a tratti se [a,b][a,b] si divide in un numero finito di intervalli su ciascuno dei quali è regolare (le poligonali ne sono l'esempio tipico: nei vertici si hanno punti angolosi).

Esempio 1 (ellisse). r‾(t)=(acos⁡t,bsin⁡t)\underline r(t) = (a\cos t, b\sin t) con a,b>0a,b>0 è chiusa e semplice su [0,2π][0,2\pi], r‾′(t)=(−asin⁡t,bcos⁡t)\underline r'(t) = (-a\sin t, b\cos t) e ∣r‾′(t)∣2=a2sin⁡2t+b2cos⁡2t>0|\underline r'(t)|^2 = a^2\sin^2 t + b^2\cos^2 t > 0 (seno e coseno non si annullano insieme): regolare. Per t0=π/4t_0 = \pi/4: punto (a22,b22)\big(\tfrac{a\sqrt2}{2}, \tfrac{b\sqrt2}{2}\big), direzione (−a22,b22)\big(-\tfrac{a\sqrt2}{2}, \tfrac{b\sqrt2}{2}\big), e la tangente è xa+yb=2\dfrac{x}{a} + \dfrac{y}{b} = \sqrt2. Per t0=π/2t_0 = \pi/2: punto (0,b)(0,b), direzione (−a,0)(-a, 0): la tangente è la retta orizzontale y=by = b.

Grafico interattivo: La cuspide (t², t³), t da −1 a 1: nell'origine r'(0) = 0, curva non regolare

Definizione. Se ogni componente di r‾:[a,b]→Rm\underline r : [a,b] \to \mathbb{R}^m è integrabile (Integrale definito di RiemannL'integrale definito di f limitata su [a,b] è il limite delle somme di Cauchy-Riemann (somme di aree di rettangoli), se esiste finito e uguale per ogni scelta dei punti; equivale a sup delle somme inferiori = inf delle somme superiori. Se f ≥ 0 è l'area del trapezoide, in generale area sopra l'asse meno area sotto. Sono integrabili le funzioni continue, monotone, continue a tratti; la funzione di Dirichlet no.Integrale definito di Riemann →), si pone ∫abr‾(t) dt=(∫abr1(t) dt, …, ∫abrm(t) dt)∈Rm.\int_a^b \underline r(t)\,dt = \Big(\int_a^b r_1(t)\,dt,\ \dots,\ \int_a^b r_m(t)\,dt\Big) \in \mathbb{R}^m .

Esempio. ∫0π(cos⁡t,sin⁡t,2t) dt=(sin⁡π−sin⁡0, −cos⁡π+cos⁡0, π2)=(0,2,π2)\int_0^{\pi}(\cos t, \sin t, 2t)\,dt = (\sin\pi - \sin 0,\ -\cos\pi + \cos 0,\ \pi^2) = (0, 2, \pi^2).

Teorema fondamentale del calcolo per funzioni vettoriali (D). Se r‾∈C1([a,b],Rm)\underline r \in C^1([a,b], \mathbb{R}^m), allora ∫abr‾′(t) dt=r‾(b)−r‾(a)\displaystyle\int_a^b \underline r'(t)\,dt = \underline r(b) - \underline r(a).

Esempio. Per r‾(t)=(t2,t3)\underline r(t) = (t^2, t^3) su [0,2][0,2]: ∫02(2t,3t2) dt=(4,8)=r‾(2)−r‾(0)\int_0^2 (2t, 3t^2)\,dt = (4, 8) = \underline r(2) - \underline r(0). ✓

Teorema (norma dell'integrale). Se r‾\underline r è integrabile, ∣∫abr‾(t) dt∣≤∫ab∣r‾(t)∣ dt\displaystyle\Big|\int_a^b \underline r(t)\,dt\Big| \le \int_a^b |\underline r(t)|\,dt.

Lunghezza di una curva e ascissa curvilinea

Definizione (poligonale inscritta, lunghezza). Sia r‾:[a,b]→Rm\underline r : [a,b] \to \mathbb{R}^m una curva. Per una partizione P={a=t0<t1<⋯<tN=b}\mathcal P = \{a = t_0 < t_1 < \dots < t_N = b\} la poligonale inscritta è quella che congiunge r‾(t0),r‾(t1),…,r‾(tN)\underline r(t_0), \underline r(t_1), \dots, \underline r(t_N), e la sua lunghezza è ℓ(P)=∑i=1N∣r‾(ti)−r‾(ti−1)∣.\ell(\mathcal P) = \sum_{i=1}^N |\underline r(t_i) - \underline r(t_{i-1})| . La curva è rettificabiledotata di lunghezza finita se {ℓ(P)}\{\ell(\mathcal P)\} è limitato superiormente; la lunghezza è ℓ(γ)=sup⁡Pℓ(P)\ell(\gamma) = \sup_{\mathcal P}\ell(\mathcal P).

Esempio. Semicirconferenza di raggio 11, r‾(t)=(cos⁡t,sin⁡t)\underline r(t) = (\cos t, \sin t), t∈[0,π]t\in[0,\pi], con N=2N = 2 intervalli uguali: i punti sono (1,0)(1,0), (0,1)(0,1), (−1,0)(-1,0) e la poligonale vale 2+2=2.83\sqrt2 + \sqrt2 = 2.83, contro la lunghezza vera π=3.14\pi = 3.14. Con N=6N = 6 si ottiene 6⋅2sin⁡(π/12)=3.1066\cdot 2\sin(\pi/12) = 3.106: si avvicina a π\pi dal basso.

Teorema (lunghezza delle curve C1C^1) (D). Se r‾∈C1([a,b],Rm)\underline r \in C^1([a,b], \mathbb{R}^m), allora γ\gamma è rettificabile e ℓ(γ)=∫ab∣r‾′(t)∣ dt.\ell(\gamma) = \int_a^b |\underline r'(t)|\,dt .

Definizione. Se r‾:[a,b]→Rm\underline r : [a,b] \to \mathbb{R}^m è regolare e φ:[c,d]→[a,b]\varphi : [c,d] \to [a,b] è biiettiva e C1C^1, la curva r‾~(u)=r‾(φ(u))\tilde{\underline r}(u) = \underline r(\varphi(u)) è una riparametrizzazione. Ha lo stesso sostegno. Se φ′>0\varphi' > 0 (crescente) le due curve sono equivalenti: stesso punto iniziale e finale, stesso verso. Se φ′<0\varphi' < 0 (decrescente) r‾~\tilde{\underline r} cambia orientazione: percorre il sostegno al contrario.

Esempio. r‾(t)=(Rcos⁡t,Rsin⁡t)\underline r(t) = (R\cos t, R\sin t), t∈[0,2π]t\in[0,2\pi]. Con t=φ(u)=u/Rt = \varphi(u) = u/R, u∈[0,2πR]u\in[0,2\pi R], si ottiene r‾~(u)=(Rcos⁡(u/R),Rsin⁡(u/R))\tilde{\underline r}(u) = (R\cos(u/R), R\sin(u/R)) equivalente, con velocità ∣r‾~′∣=1|\tilde{\underline r}'| = 1 (mentre ∣r‾′∣=R|\underline r'| = R). Con t=−ut = -u, u∈[−2π,0]u\in[-2\pi,0], si ottiene la circonferenza percorsa in senso orario.

Proposizione (la lunghezza non dipende dalla parametrizzazione). Se r‾~\tilde{\underline r} è una riparametrizzazione di r‾\underline r (φ′>0\varphi' > 0 oppure φ′<0\varphi' < 0), allora ℓ(r‾~)=ℓ(r‾)\ell(\tilde{\underline r}) = \ell(\underline r).

Definizione. Sia r‾:[a,b]→Rm\underline r : [a,b]\to\mathbb{R}^m regolare. L'ascissa curvilinea è s(t)=∫at∣r‾′(τ)∣ dτ,s:[a,b]→[0,ℓ(γ)],s(t) = \int_a^t |\underline r'(\tau)|\,d\tau , \qquad s : [a,b]\to[0,\ell(\gamma)] , la lunghezza dell'arco tra r‾(a)\underline r(a) e r‾(t)\underline r(t).

Teorema (parametrizzazione ad arco) (D). La curva r‾~(s)=r‾(t(s))\tilde{\underline r}(s) = \underline r(t(s)), s∈[0,ℓ(γ)]s\in[0,\ell(\gamma)], è una riparametrizzazione equivalente di r‾\underline r e ∣r‾~′(s)∣=1|\tilde{\underline r}'(s)| = 1 per ogni ss.

Esempio (circonferenza). ∣r‾′∣=R|\underline r'| = R, s=Rts = Rt, t=s/Rt = s/R: r‾~(s)=(Rcos⁡sR,Rsin⁡sR)\tilde{\underline r}(s) = (R\cos\frac sR, R\sin\frac sR), s∈[0,2πR]s\in[0,2\pi R].

Definizione. La curvaturamisura di quanto ruota la tangente per ogni unità di lunghezza percorsa è κ=∣dT‾ds∣\kappa = \Big|\dfrac{d\underline T}{ds}\Big| e, dove κ≠0\kappa\ne0, il versore normale principale è N‾=1κdT‾ds\underline N = \dfrac{1}{\kappa}\dfrac{d\underline T}{ds}. Il raggio di curvatura è 1/κ1/\kappa.

Formula (decomposizione dell'accelerazione). a‾=aT T‾+aN N‾\underline a = a_T\,\underline T + a_N\,\underline N, con aT=v′=r‾′⋅r‾′′∣r‾′∣,aN=κv2=∣r‾′×r‾′′∣∣r‾′∣,∣a‾∣2=aT2+aN2.a_T = v' = \frac{\underline r'\cdot\underline r''}{|\underline r'|},\qquad a_N = \kappa v^2 = \frac{|\underline r'\times\underline r''|}{|\underline r'|},\qquad |\underline a|^2 = a_T^2 + a_N^2 . La curvatura si calcola con κ=∣r‾′×r‾′′∣∣r‾′∣3 (R3),κ=∣x′y′′−x′′y′∣(x′2+y′2)3/2 (R2),κ=∣f′′∣(1+f′2)3/2 (grafico).\kappa = \frac{|\underline r'\times\underline r''|}{|\underline r'|^3}\ (\mathbb{R}^3),\qquad \kappa = \frac{|x'y'' - x''y'|}{(x'^2+y'^2)^{3/2}}\ (\mathbb{R}^2),\qquad \kappa = \frac{|f''|}{(1+f'^2)^{3/2}}\ (\text{grafico}).

Integrali curvilinei di prima specie, baricentri e momenti d'inerzia

Definizione (integrale di prima specie). Sia r‾:[a,b]→Rm\underline r : [a,b] \to \mathbb{R}^m una curva C1C^1 con sostegno γ\gamma e sia f:γ→Rf : \gamma \to \mathbb{R} continua. Si pone ∫γf ds=∫abf(r‾(t)) ∣r‾′(t)∣ dt.\int_\gamma f\,ds = \int_a^b f(\underline r(t))\,|\underline r'(t)|\,dt . ds=∣r‾′(t)∣ dtds = |\underline r'(t)|\,dt è l'elemento d'arco.

Esempio (la parete sotto un piano). γ\gamma: semicirconferenza r‾(t)=(cos⁡t,sin⁡t)\underline r(t) = (\cos t, \sin t), t∈[π,2π]t\in[\pi, 2\pi] (la metà inferiore, dove y≤0y \le 0), f(x,y)=−y≥0f(x,y) = -y \ge 0. Qui ∣r‾′∣=1|\underline r'| = 1 e f(r‾(t))=−sin⁡tf(\underline r(t)) = -\sin t: ∫γ(−y) ds=∫π2π(−sin⁡t) dt=[cos⁡t]π2π=1−(−1)=2.\int_\gamma(-y)\,ds = \int_\pi^{2\pi}(-\sin t)\,dt = [\cos t]_\pi^{2\pi} = 1 - (-1) = 2 . È l'area della parete alta 11 al centro e 00 agli estremi, che si appoggia sulla semicirconferenza.

Teorema (D). Sia r‾~(u)=r‾(φ(u))\tilde{\underline r}(u) = \underline r(\varphi(u)), u∈[c,d]u\in[c,d], una riparametrizzazione di r‾:[a,b]→Rm\underline r : [a,b]\to\mathbb{R}^m (φ\varphi biiettiva e C1C^1), sia con φ′>0\varphi' > 0 sia con φ′<0\varphi' < 0 (verso opposto). Allora per ogni ff continua ∫γ~f ds=∫γf ds.\int_{\tilde\gamma}f\,ds = \int_\gamma f\,ds .

Definizione (massa e baricentro). Per una curva γ\gamma con densità μ≥0\mu \ge 0: M=∫γμ ds,xb=1M∫γx μ ds,yb=1M∫γy μ ds,zb=1M∫γz μ ds.M = \int_\gamma\mu\,ds,\qquad x_b = \frac1M\int_\gamma x\,\mu\,ds,\quad y_b = \frac1M\int_\gamma y\,\mu\,ds,\quad z_b = \frac1M\int_\gamma z\,\mu\,ds . Se il filo è omogeneo (μ≡1\mu\equiv1) si ha M=ℓ(γ)M = \ell(\gamma) e si parla di baricentro geometrico (o centroide).

Esempio (quarto di circonferenza). r‾(t)=(Rcos⁡t,Rsin⁡t)\underline r(t) = (R\cos t, R\sin t), t∈[0,π2]t\in[0,\frac\pi2], μ=1\mu = 1: ∣r‾′∣=R|\underline r'| = R, M=πR2M = \frac{\pi R}{2}, xb=2πR∫0π/2Rcos⁡t⋅R dt=2Rπ[sin⁡t]0π/2=2Rπ,yb=2Rπx_b = \frac{2}{\pi R}\int_0^{\pi/2}R\cos t\cdot R\,dt = \frac{2R}{\pi}\big[\sin t\big]_0^{\pi/2} = \frac{2R}{\pi},\qquad y_b = \frac{2R}{\pi} (per la simmetria rispetto alla bisettrice y=xy = x le due coordinate coincidono). Poiché 2/π≈0.6372/\pi\approx0.637, il baricentro (2Rπ,2Rπ)\big(\frac{2R}\pi,\frac{2R}\pi\big) sta a distanza 22πR≈0.90R\frac{2\sqrt2}{\pi}R \approx 0.90R dall'origine: dentro la regione racchiusa dall'arco, non sull'arco.

Definizione. Fissato un asse, sia δ(x,y,z)\delta(x,y,z) la distanza del punto dall'asse. Il momento d'inerzia della curva (densità μ\mu) rispetto a quell'asse è I=∫γδ2 μ ds.I = \int_\gamma\delta^2\,\mu\,ds .

Esempio. Arco di circonferenza r‾(t)=(Rcos⁡t,Rsin⁡t)\underline r(t) = (R\cos t, R\sin t), t∈[0,π2]t\in[0,\frac\pi2], μ=1\mu = 1, rispetto all'asse yy: δ=∣x∣=Rcos⁡t\delta = |x| = R\cos t, ∣r‾′∣=R|\underline r'| = R, I=∫0π/2(Rcos⁡t)2⋅R dt=R3∫0π/2cos⁡2t dt=πR34.I = \int_0^{\pi/2}(R\cos t)^2\cdot R\,dt = R^3\int_0^{\pi/2}\cos^2t\,dt = \frac{\pi R^3}{4}. (Per l'integrale di cos⁡2\cos^2: per parti si ottiene J=∫cos⁡2=[sin⁡tcos⁡t]+∫sin⁡2t=π2−JJ = \int\cos^2 = [\sin t\cos t] + \int\sin^2 t = \frac\pi2 - J, quindi J=π4J = \frac\pi4.) Il filo ha massa M=πR2M = \frac{\pi R}2, quindi I=R22MI = \frac{R^2}{2}M: come se tutta la massa stesse alla distanza R/2R/\sqrt2 dall'asse.

2. Funzioni di più variabili: limiti, continuità e topologia

Funzioni di più variabili, limiti e continuità

Definizione (dominio naturale). Se ff è data da una formula, il dominio naturaleil più grande insieme su cui la formula ha senso è il più grande insieme dove la formula ha senso: radici di indice pari con radicando ≥0\ge 0, logaritmi con argomento >0>0, denominatori ≠0\ne 0, arcsin⁡\arcsin e arccos⁡\arccos con argomento in [−1,1][-1,1].

Esempio. f(x,y)=xy log⁡(x+y−1)f(x,y) = \sqrt{xy}\,\log(x+y-1) richiede xy≥0xy \ge 0 e x+y−1>0x+y-1>0. La prima condizione descrive il primo e il terzo quadrante, bordi inclusi; la seconda il semipiano aperto sopra la retta y=1−xy = 1-x. Nel terzo quadrante x+y≤0x+y\le0, quindi la seconda non vale mai: resta dom⁡f={x≥0, y≥0, x+y>1}.\operatorname{dom}f = \{x\ge0,\ y\ge0,\ x+y>1\}. Gli assi non sono esclusi del tutto: per esempio (2,0)(2,0) è nel dominio (xy=0≥0xy=0\ge0 e 2−1>02-1>0). Invece il segmento da (1,0)(1,0) a (0,1)(0,1) è escluso, perché lì x+y−1=0x+y-1=0 e il logaritmo vuole un argomento strettamente positivo.

Definizione (grafico, insiemi di livello). Il grafico è graf⁡f={(x‾,f(x‾)):x‾∈A}⊆Rn+1\operatorname{graf}f = \{(\underline x, f(\underline x)) : \underline x\in A\}\subseteq\mathbb{R}^{n+1}; per n=2n=2 è una superficie z=f(x,y)z = f(x,y). L'insieme di livello cc è Lc={x‾∈A:f(x‾)=c}L_c = \{\underline x\in A : f(\underline x) = c\}.

Esempio 1. f(x,y)=1−x2−y2f(x,y) = \sqrt{1-x^2-y^2} ha dominio il disco x2+y2≤1x^2+y^2\le1 e valori in [0,1][0,1]. Per c<0c<0 o c>1c>1, Lc=∅L_c = \emptyset; per c=1c = 1, L1={(0,0)}L_1 = \{(0,0)\}; per 0≤c<10\le c<1 si eleva al quadrato (lecito, entrambi i membri ≥0\ge0): x2+y2=1−c2x^2+y^2 = 1-c^2, circonferenza di raggio 1−c2\sqrt{1-c^2} (per c=12c=\frac12, raggio 32\frac{\sqrt3}2). Il grafico è z2=1−x2−y2z^2 = 1-x^2-y^2, z≥0z\ge0: la semisfera superiore.

Grafico interattivo: Insiemi di livello di f = √(1 − x² − y²): circonferenze di raggio √(1 − c²) per c = 0, 0.5, 0.8, e il punto (0,0) per c = 1

Definizione (punto di accumulazionepunto vicino al quale cadono infiniti punti dell'insieme). x‾0∈Rn\underline x_0\in\mathbb{R}^n è di accumulazione per E⊆RnE\subseteq\mathbb{R}^n se ogni intorno di x‾0\underline x_0 contiene un punto di EE diverso da x‾0\underline x_0. Non è detto che x‾0∈E\underline x_0\in E.

Definizione (limite, tre forme equivalenti).

  • Per successioni (Heine). lim⁡x‾→x‾0f=L\lim_{\underline x\to\underline x_0}f = L se per ogni successione x‾k∈X\underline x_k\in X, x‾k≠x‾0\underline x_k\ne\underline x_0, x‾k→x‾0\underline x_k\to\underline x_0 si ha f(x‾k)→Lf(\underline x_k)\to L.
  • Diretta, LL finito. ∀ε>0 ∃δ>0: x‾∈X, 0<∣x‾−x‾0∣<δ⇒∣f(x‾)−L∣<ε\forall\varepsilon>0\ \exists\delta>0:\ \underline x\in X,\ 0<|\underline x - \underline x_0|<\delta\Rightarrow|f(\underline x) - L|<\varepsilon.
  • Diretta, L=+∞L = +\infty. ∀K>0 ∃δ>0: 0<∣x‾−x‾0∣<δ⇒f(x‾)>K\forall K>0\ \exists\delta>0:\ 0<|\underline x - \underline x_0|<\delta\Rightarrow f(\underline x)>K (e simile per −∞-\infty).

Esempio (ε\varepsilon-δ\delta). lim⁡x‾→0‾∣x‾∣=0\lim_{\underline x\to\underline 0}|\underline x| = 0: dato ε\varepsilon, basta δ=ε\delta=\varepsilon perché ∣ ∣x‾∣−0 ∣=∣x‾∣<δ=ε|\,|\underline x|-0\,| = |\underline x|<\delta=\varepsilon.

Teorema (permanenza del segno). Se lim⁡x‾→x‾0f=L>0\lim_{\underline x\to\underline x_0}f = L>0 (anche +∞+\infty), allora esiste un intorno forato di x‾0\underline x_0 in cui f>0f>0.

Definizione. Sia x‾0∈X\underline x_0\in X di accumulazione per XX. ff è continua in x‾0\underline x_0 se lim⁡x‾→x‾0f(x‾)=f(x‾0)\lim_{\underline x\to\underline x_0}f(\underline x) = f(\underline x_0). In un punto isolato ff è continua per definizione. C(X)C(X) è l'insieme delle funzioni continue in ogni punto di XX.

Esempio. f(x,y)=ex+ycos⁡x+2f(x,y)=\dfrac{e^{x+y}}{\cos x + 2} è continua su R2\mathbb{R}^2: x+yx+y è somma di proiezioni, ex+ye^{x+y} è composizione con l'esponenziale; cos⁡x+2≥1>0\cos x+2\ge1>0 non si annulla mai, quindi il quoziente è continuo ovunque.

Proposizione. Se lim⁡x‾→x‾0f=L\lim_{\underline x\to\underline x_0}f = L e γ‾:I→X\underline\gamma : I\to X è una curva continua con γ‾(t0)=x‾0\underline\gamma(t_0) = \underline x_0 e γ‾(t)≠x‾0\underline\gamma(t)\ne\underline x_0 per t≠t0t\ne t_0, allora lim⁡t→t0f(γ‾(t))=L\lim_{t\to t_0}f(\underline\gamma(t)) = L.

Esempio (fascio di rette). f(x,y)=xyx2+y2f(x,y)=\dfrac{xy}{x^2+y^2} in (0,0)(0,0) è del tipo 00\frac00. Sulle rette y=mxy=mx: f(t,mt)=mt2t2(1+m2)=m1+m2f(t,mt)=\dfrac{mt^2}{t^2(1+m^2)} = \dfrac m{1+m^2}, che dipende da mm (12\frac12 per m=1m=1, −12-\frac12 per m=−1m=-1, 00 per m=0m=0). Il limite non esiste.

Teorema (maggiorazione con funzione radiale). Sia L∈RL\in\mathbb{R} e g:(0,+∞)→Rg:(0,+\infty)\to\mathbb{R} con lim⁡ρ→0+g(ρ)=0\lim_{\rho\to0^+}g(\rho)=0. Se esiste r>0r>0 tale che ∣f(x‾)−L∣≤g(∣x‾−x‾0∣)per ogni x‾∈X, 0<∣x‾−x‾0∣<r,|f(\underline x) - L|\le g\big(|\underline x - \underline x_0|\big)\quad\text{per ogni }\underline x\in X,\ 0<|\underline x - \underline x_0|<r, allora lim⁡x‾→x‾0f(x‾)=L\lim_{\underline x\to\underline x_0}f(\underline x) = L.

Esempio. lim⁡(x,y)→(0,0)x2−y∣x‾∣1/2\displaystyle\lim_{(x,y)\to(0,0)}\frac{x^2-y}{|\underline x|^{1/2}}. Sulla parabola y=x2y=x^2 il numeratore è 00: candidato L=0L=0. In polari f~(ρ,θ)=ρ2cos⁡2θ−ρsin⁡θρ1/2=ρ1/2(ρcos⁡2θ−sin⁡θ),∣f~∣≤ρ1/2(ρ+1)=:g(ρ)→0.\tilde f(\rho,\theta)=\frac{\rho^2\cos^2\theta-\rho\sin\theta}{\rho^{1/2}}=\rho^{1/2}\big(\rho\cos^2\theta-\sin\theta\big),\quad|\tilde f|\le\rho^{1/2}(\rho+1)=:g(\rho)\to0 . Il limite vale 00.

Topologia di Rn, teorema di Weierstrass e teorema degli zeri

Definizione.

  • x‾\underline x è interno a EE se ∃r>0: Br(x‾)⊆E\exists r>0:\ B_r(\underline x)\subseteq E. L'insieme dei punti interni è l'interno E˚⊆E\mathring E\subseteq E.
  • x‾\underline x è esterno a EE se è interno al complementare Ec=Rn∖EE^c=\mathbb{R}^n\setminus E: ∃r: Br(x‾)∩E=∅\exists r:\ B_r(\underline x)\cap E=\emptyset.
  • x‾\underline x è di frontiera se non è né interno né esterno, cioè se ogni intorno contiene sia punti di EE sia punti di EcE^c. L'insieme è ∂E\partial E.
  • La chiusura è E‾=E∪∂E\overline E=E\cup\partial E.

Esempio 1. E={(x,y):x≥0, y>0}E=\{(x,y): x\ge0,\ y>0\}. Un punto con x>0,y>0x>0,y>0 ha un disco attorno contenuto in EE: interno. Un punto (0,y0)(0,y_0), y0>0y_0>0, ha a sinistra punti con x<0x<0 fuori da EE: non interno. Quindi E˚={x>0, y>0},∂E={(x,0):x≥0}∪{(0,y):y≥0},E‾={x≥0, y≥0}.\mathring E=\{x>0,\ y>0\},\quad\partial E=\{(x,0):x\ge0\}\cup\{(0,y):y\ge0\},\quad\overline E=\{x\ge0,\ y\ge0\}. Qui (0,1)∈E∩∂E(0,1)\in E\cap\partial E e (1,0)∈E‾∖E(1,0)\in\overline E\setminus E.

Definizione. EE è aperto se E=E˚E=\mathring E, cioè se ogni suo punto è interno. EE è chiuso se EcE^c è aperto.

Teorema. La palla aperta Br(x‾0)B_r(\underline x_0) è un insieme aperto; la palla chiusa {∣y‾−x‾0∣≤r}\{|\underline y-\underline x_0|\le r\} e la sfera {∣y‾−x‾0∣=r}\{|\underline y-\underline x_0|=r\} sono chiuse.

Proposizione. Se ff è continua su Rn\mathbb{R}^n (o su un aperto), allora per ogni cc gli insiemi {f<c}\{f<c\} e {f>c}\{f>c\} sono aperti, mentre {f≤c}\{f\le c\}, {f≥c}\{f\ge c\}, {f=c}\{f=c\} sono chiusi (in Rn\mathbb{R}^n). In generale la controimmagine di un aperto è aperta e quella di un chiuso è chiusa.

Esempio. {0<∣x‾∣<R}={∣x‾∣>0}∩{∣x‾∣<R}\{0<|\underline x|<R\}=\{|\underline x|>0\}\cap\{|\underline x|<R\} è l'intersezione di due aperti: aperto. Non è chiuso (non contiene 0‾\underline 0 né la sfera di raggio RR, che sono di frontiera), è limitato; ∂E={0‾}∪{∣x‾∣=R}\partial E=\{\underline 0\}\cup\{|\underline x|=R\} e E‾={∣x‾∣≤R}\overline E=\{|\underline x|\le R\}. Esempio. {x2+y2<4}∩{y≥x2}\{x^2+y^2<4\}\cap\{y\ge x^2\}: aperto ∩\cap chiuso. Non è né aperto né chiuso (sul bordo del disco non c'è la circonferenza ma c'è la parabola).

Grafico interattivo: E = {0 < |x| < 2} nel piano: disco forato, aperto e limitato; frontiera = origine + circonferenza di raggio 2

Teorema (Weierstrass). Se K⊆RnK\subseteq\mathbb{R}^n è compatto (non vuoto) e f:K→Rf:K\to\mathbb{R} è continua, allora ff ha massimo e minimo assoluti in KK.

Esempio. f(x,y)=x2yf(x,y)=x^2y sul disco x2+y2≤1x^2+y^2\le1: il disco è chiuso (è {g≤1}\{g\le1\} con gg continua) e limitato, ff è un polinomio. Massimo e minimo esistono.

Corollario. Sia C⊆RnC\subseteq\mathbb{R}^n chiuso, f:C→Rf:C\to\mathbb{R} continua. Se lim⁡∣x‾∣→∞, x‾∈Cf(x‾)=+∞\displaystyle\lim_{|\underline x|\to\infty,\ \underline x\in C}f(\underline x)=+\infty, allora ff ha un minimo assoluto in CC. Se il limite è −∞-\infty, ha un massimo assoluto.

Esempio. f(x,y)=x2+y2−xy≥12(x2+y2)f(x,y)=x^2+y^2-xy\ge\frac12(x^2+y^2) (perché xy≤12(x2+y2)xy\le\frac12(x^2+y^2)) →+∞\to+\infty: ha minimo assoluto su R2\mathbb{R}^2 (è 00 nell'origine); non ha massimo, essendo illimitata superiormente.

Definizione. E⊆RnE\subseteq\mathbb{R}^n è connesso per archifatto di un solo pezzo: due punti qualsiasi si uniscono con una curva continua dentro l'insieme se per ogni x‾,y‾∈E\underline x,\underline y\in E esiste una curva continua γ‾:[a,b]→E\underline\gamma:[a,b]\to E con γ‾(a)=x‾\underline\gamma(a)=\underline x, γ‾(b)=y‾\underline\gamma(b)=\underline y (cioè il sostegno sta in EE).

Teorema degli zeri. Sia E⊆RnE\subseteq\mathbb{R}^n connesso, f:E→Rf:E\to\mathbb{R} continua. Se f(x‾)>0f(\underline x)>0 e f(y‾)<0f(\underline y)<0 per due punti di EE, esiste z‾∈E\underline z\in E con f(z‾)=0f(\underline z)=0.

Esempio. f(x,y)=x2−y2−1f(x,y)=x^2-y^2-1. Z={x2−y2=1}Z=\{x^2-y^2=1\} è l'iperbole equilatera con vertici (±1,0)(\pm1,0) e asintoti y=±xy=\pm x, e divide il piano in tre regioni connesse: a sinistra del ramo sinistro R1R_1, la fascia centrale R2R_2, a destra del ramo destro R3R_3. Punti di prova: f(−2,0)=3>0f(-2,0)=3>0, f(0,0)=−1<0f(0,0)=-1<0, f(2,0)=3>0f(2,0)=3>0. Quindi f>0f>0 in R1∪R3R_1\cup R_3, f<0f<0 in R2R_2.

Grafico interattivo: Segno di f = x² − y² − 1: f > 0 nelle due regioni esterne ai rami dell'iperbole, f < 0 tra i rami (contiene l'origine)

3. Calcolo differenziale per funzioni scalari

Derivate parziali, differenziabilità e piano tangente

Definizione (derivata parziale, gradiente). Sia f:X⊆Rn→Rf : X\subseteq\mathbb{R}^n\to\mathbb{R} e x‾0\underline x_0 interno a XX. La derivata parziale rispetto a xix_i è ∂f∂xi(x‾0)=fxi(x‾0)=lim⁡h→0f(x‾0+h e‾i)−f(x‾0)h,\frac{\partial f}{\partial x_i}(\underline x_0) = f_{x_i}(\underline x_0) = \lim_{h\to0}\frac{f(\underline x_0 + h\,\underline e_i) - f(\underline x_0)}{h}, se il limite esiste finito. Se esistono tutte le derivate parziali, ff è derivabile in x‾0\underline x_0 e il gradientevettore che ha come componenti le derivate parziali è il vettore ∇f(x‾0)=(fx1(x‾0),…,fxn(x‾0))\nabla f(\underline x_0) = \big(f_{x_1}(\underline x_0),\dots,f_{x_n}(\underline x_0)\big).

Esempio. f(x,y)=ex2yf(x,y)=e^{x^2y}: fx=2xy ex2yf_x = 2xy\,e^{x^2y}, fy=x2ex2yf_y = x^2e^{x^2y}, quindi ∇f(1,1)=(2e,e)\nabla f(1,1)=(2e,e).

Definizione (differenziabile). ff è differenziabile in x‾0\underline x_0 (interno) se esiste a‾∈Rn\underline a\in\mathbb{R}^n tale che f(x‾0+h‾)=f(x‾0)+a‾⋅h‾+o(∣h‾∣)(h‾→0‾),f(\underline x_0+\underline h) = f(\underline x_0) + \underline a\cdot\underline h + o(|\underline h|)\quad(\underline h\to\underline 0), cioè lim⁡h‾→0‾f(x‾0+h‾)−f(x‾0)−a‾⋅h‾∣h‾∣=0\displaystyle\lim_{\underline h\to\underline 0}\frac{f(\underline x_0+\underline h)-f(\underline x_0)-\underline a\cdot\underline h}{|\underline h|}=0.

Teorema (D). Se ff è differenziabile in x‾0\underline x_0 allora (1) ff è continua in x‾0\underline x_0; (2) esistono tutte le derivate parziali e a‾=∇f(x‾0)\underline a=\nabla f(\underline x_0).

Esempio. f(x,y)=∣x∣ yf(x,y)=\sqrt{|x|}\,y in (0,0)(0,0): f(x,0)=0f(x,0)=0 e f(0,y)=0f(0,y)=0, quindi ∇f(0,0)=(0,0)\nabla f(0,0)=(0,0). L'errore relativo è ∣x∣ yx2+y2=ρ1/2∣cos⁡θ∣sin⁡θ\dfrac{\sqrt{|x|}\,y}{\sqrt{x^2+y^2}} = \rho^{1/2}\sqrt{|\cos\theta|}\sin\theta, in modulo ≤ρ1/2→0\le\rho^{1/2}\to0: ff è differenziabile in (0,0)(0,0) (anche se non lo è in (0,y0)(0,y_0), y0≠0y_0\ne0, dove fxf_x non esiste).

Definizione (iperpiano tangente). Se ff è differenziabile in x‾0\underline x_0, l'iperpiano tangente al grafico in (x‾0,f(x‾0))(\underline x_0,f(\underline x_0)) è z=f(x‾0)+∇f(x‾0)⋅(x‾−x‾0).z=f(\underline x_0)+\nabla f(\underline x_0)\cdot(\underline x-\underline x_0). Per n=1n=1 è la retta tangente; per n=2n=2 è il piano tangente z=f(x0,y0)+fx(x0,y0)(x−x0)+fy(x0,y0)(y−y0),z=f(x_0,y_0)+f_x(x_0,y_0)(x-x_0)+f_y(x_0,y_0)(y-y_0), con vettore normale N‾=(−fx(x0,y0), −fy(x0,y0), 1)\underline N=(-f_x(x_0,y_0),\,-f_y(x_0,y_0),\,1) (dall'equazione −fx x−fy y+z=cost-f_x\,x - f_y\,y + z = \text{cost}).

Esempio. f=ex2yf=e^{x^2y} in (1,1)(1,1), punto del grafico (1,1,e)(1,1,e): f=ef=e, fx=2ef_x=2e, fy=ef_y=e, quindi z=e+2e(x−1)+e(y−1)=2ex+ey−2e,N‾=(−2e,−e,1).z=e+2e(x-1)+e(y-1)=2ex+ey-2e,\qquad\underline N=(-2e,-e,1).

Teorema del differenziale totale. Se in un intorno di x‾0\underline x_0 esistono tutte le derivate parziali di ff e sono continue in x‾0\underline x_0, allora ff è differenziabile in x‾0\underline x_0.

Definizione (C1C^1). Su un aperto XX, ff è di classe C1(X)C^1(X) se è continua, esistono tutte le derivate parziali e sono continue in XX. Per il teorema, f∈C1(X)⇒ff\in C^1(X)\Rightarrow f differenziabile in ogni punto di XX.

Esempio. ex2ye^{x^2y} ha derivate parziali continue su R2\mathbb{R}^2: è C1C^1, quindi differenziabile ovunque, e ha piano tangente in ogni punto.

Teorema (derivazione lungo una curva). Se r‾:I→X\underline r:I\to X è derivabile in tt e ff è differenziabile in r‾(t)\underline r(t), allora (f∘r‾)′(t)=∇f(r‾(t))⋅r‾′(t)=∑i=1nfxi(r‾(t)) ri′(t).(f\circ\underline r)'(t)=\nabla f(\underline r(t))\cdot\underline r'(t)=\sum_{i=1}^n f_{x_i}(\underline r(t))\,r_i'(t).

Teorema (funzione di una variabile composta con gg). Se g:X→Rg:X\to\mathbb{R} è differenziabile in x‾\underline x e φ:R→R\varphi:\mathbb{R}\to\mathbb{R} è derivabile in g(x‾)g(\underline x), allora ∇(φ∘g)(x‾)=φ′(g(x‾)) ∇g(x‾)\nabla(\varphi\circ g)(\underline x)=\varphi'(g(\underline x))\,\nabla g(\underline x).

Esempio. ∇(log⁡(x2+y2))=1x2+y2(2x,2y)\nabla\big(\log(x^2+y^2)\big)=\dfrac{1}{x^2+y^2}(2x,2y). Esempio. r‾(t)=(cos⁡t,sin⁡t)\underline r(t)=(\cos t,\sin t) e f(x,y)=x2yf(x,y)=x^2y: (f∘r‾)′(t)=2xy⋅(−sin⁡t)+x2cos⁡t=−2cos⁡tsin⁡2t+cos⁡3t(f\circ\underline r)'(t)=2xy\cdot(-\sin t)+x^2\cos t=-2\cos t\sin^2t+\cos^3t.

Derivate direzionali, formula del gradiente e curve di livello

Definizione. Sia x‾0\underline x_0 interno al dominio di ff e v‾∈Rn\underline v\in\mathbb{R}^n un versore (∣v‾∣=1|\underline v|=1). La derivata direzionalependenza di f lungo la retta per x0 con direzione v di ff in x‾0\underline x_0 lungo v‾\underline v è Dv‾f(x‾0)=∂f∂v‾(x‾0)=lim⁡t→0f(x‾0+t v‾)−f(x‾0)t=h′(0),h(t)=f(x‾0+tv‾),D_{\underline v}f(\underline x_0)=\frac{\partial f}{\partial\underline v}(\underline x_0)=\lim_{t\to0}\frac{f(\underline x_0+t\,\underline v)-f(\underline x_0)}{t}=h'(0),\qquad h(t)=f(\underline x_0+t\underline v), se il limite esiste finito.

Esempio (dalla definizione). f(x,y)=ex2yf(x,y)=e^{x^2y}, v‾=(a,b)\underline v=(a,b) con a2+b2=1a^2+b^2=1. Con h(t)=e(x0+ta)2(y0+tb)h(t)=e^{(x_0+ta)^2(y_0+tb)} si ha h′(t)=h(t) [2a(x0+ta)(y0+tb)+b(x0+ta)2]h'(t)=h(t)\,[2a(x_0+ta)(y_0+tb)+b(x_0+ta)^2], quindi Dv‾f(x0,y0)=ex02y0 (2a x0y0+b x02)=a fx+b fy.D_{\underline v}f(x_0,y_0)=e^{x_0^2y_0}\,(2a\,x_0y_0+b\,x_0^2)=a\,f_x+b\,f_y . Questa uguaglianza non è un caso.

Teorema (formula del gradiente) (D). Se ff è differenziabile in x‾0\underline x_0, allora per ogni versore v‾\underline v esiste Dv‾f(x‾0)D_{\underline v}f(\underline x_0) e Dv‾f(x‾0)=∇f(x‾0)⋅v‾.D_{\underline v}f(\underline x_0)=\nabla f(\underline x_0)\cdot\underline v .

Esempio. f=ex2yf=e^{x^2y} in (1,1)(1,1), ∇f=(2e,e)\nabla f=(2e,e). Nella direzione v‾=(35,45)\underline v=\big(\frac35,\frac45\big) (norma 11): Dv‾f=6e5+4e5=2eD_{\underline v}f=\frac{6e}{5}+\frac{4e}{5}=2e. Nella direzione (1,1)(1,1) non normalizzata bisogna prima dividere per 2\sqrt2: v‾=12(1,1)\underline v=\frac1{\sqrt2}(1,1) e Dv‾f=3e2D_{\underline v}f=\frac{3e}{\sqrt2}.

Teorema (D). Se ∇f(x‾)≠0‾\nabla f(\underline x)\ne\underline 0:

  • la derivata direzionale è massima per α=0\alpha=0, cioè lungo v‾max⁡=∇f(x‾)∣∇f(x‾)∣\underline v_{\max}=\dfrac{\nabla f(\underline x)}{|\nabla f(\underline x)|}, e vale ∣∇f(x‾)∣|\nabla f(\underline x)|;
  • è minima per α=π\alpha=\pi, cioè lungo v‾min⁡=−∇f(x‾)∣∇f(x‾)∣\underline v_{\min}=-\dfrac{\nabla f(\underline x)}{|\nabla f(\underline x)|}, e vale −∣∇f(x‾)∣-|\nabla f(\underline x)|;
  • è nulla nelle direzioni ortogonali a ∇f(x‾)\nabla f(\underline x) (α=π2\alpha=\frac\pi2). Se ∇f(x‾)=0‾\nabla f(\underline x)=\underline 0, tutte le derivate direzionali sono nulle.

Esempio (la collina). Un'altura ha quota (in metri) z=h(x,y)=400−x2+4y2100z=h(x,y)=400-\dfrac{x^2+4y^2}{100}, con x,yx,y in metri. In (10,5)(10,5): h=400−100+100100=398h=400-\frac{100+100}{100}=398, ∇h=(−x50,−2y25)=(−0.2,−0.4)\nabla h=\big(-\frac{x}{50},-\frac{2y}{25}\big)=(-0.2,-0.4) e ∣∇h∣=0.2≈0.447|\nabla h|=\sqrt{0.2}\approx0.447. La salita più ripida è lungo ∇h∣∇h∣=(−15,−25)\frac{\nabla h}{|\nabla h|}=\big(-\frac1{\sqrt5},-\frac2{\sqrt5}\big) (verso la cima, nell'origine), con pendenza 0.4470.447 m di quota per metro orizzontale (circa 24∘24^\circ); la discesa più ripida è la direzione opposta. Camminando lungo v‾=(25,−15)\underline v=\big(\frac2{\sqrt5},-\frac1{\sqrt5}\big), perpendicolare al gradiente (∇h⋅v‾=−0.4+0.45=0\nabla h\cdot\underline v=\frac{-0.4+0.4}{\sqrt5}=0), si resta, al primo ordine, a quota costante: si segue la curva di livello (la "isoipsa").

Teorema (D). Sia ff differenziabile in (x0,y0)(x_0,y_0) e r‾:I→R2\underline r:I\to\mathbb{R}^2 una curva regolare con r‾(t0)=(x0,y0)\underline r(t_0)=(x_0,y_0) il cui sostegno sta nell'insieme di livello Lc={f=c}L_c=\{f=c\}. Allora ∇f(x0,y0)⋅r‾′(t0)=0\nabla f(x_0,y_0)\cdot\underline r'(t_0)=0: il gradiente è ortogonale al vettore tangente alla curva. Se ∇f(x0,y0)≠0‾\nabla f(x_0,y_0)\ne\underline 0, la retta tangente alla curva di livello in (x0,y0)(x_0,y_0) è fx(x0,y0)(x−x0)+fy(x0,y0)(y−y0)=0.f_x(x_0,y_0)(x-x_0)+f_y(x_0,y_0)(y-y_0)=0 .

Esempio (ellisse). L'ellisse x2a2+y2b2=1\frac{x^2}{a^2}+\frac{y^2}{b^2}=1 è il livello 11 di f=x2a2+y2b2f=\frac{x^2}{a^2}+\frac{y^2}{b^2}, con ∇f=(2xa2,2yb2)\nabla f=\big(\frac{2x}{a^2},\frac{2y}{b^2}\big), mai nullo sull'ellisse. La tangente in (x0,y0)(x_0,y_0) è 2x0a2(x−x0)+2y0b2(y−y0)=0\frac{2x_0}{a^2}(x-x_0)+\frac{2y_0}{b^2}(y-y_0)=0, cioè (dividendo per 22 e usando x02a2+y02b2=1\frac{x_0^2}{a^2}+\frac{y_0^2}{b^2}=1) x0 xa2+y0 yb2=1.\frac{x_0\,x}{a^2}+\frac{y_0\,y}{b^2}=1 . Controllo con la parametrizzazione (acos⁡t,bsin⁡t)(a\cos t,b\sin t) di Curve parametriche, vettore tangente e integrali di funzioni vettorialiUna curva è una funzione continua r : [a,b] → R^m; il sostegno è la sua immagine. Limiti, continuità e derivate si fanno componente per componente. La curva è regolare se r è C^1 e r'(t) ≠ 0 ovunque: allora r'(t) è il vettore tangente e la retta tangente è r(t0) + λ r'(t0). Esempi: circonferenza, segmento, grafico (t, f(t)), curve polari (f(θ)cosθ, f(θ)sinθ) regolari se f' ≠ 0 dove f = 0, elica, cuspide (t²,t³). L'integrale di r si fa per componenti; vale ∫ r' = r(b) − r(a) e ‖∫ r‖ ≤ ∫ ‖r‖.Curve parametriche, vettore tangente e integrali di funzioni vettoriali →: per t=π/4t=\pi/4 si riottiene xa+yb=2\frac xa+\frac yb=\sqrt2.

Derivate seconde, teorema di Schwarz e formula di Taylor del secondo ordine

Teorema (Schwarz). Se fxif_{x_i}, fxjf_{x_j} e le derivate miste fxixjf_{x_ix_j}, fxjxif_{x_jx_i} esistono in un intorno di x‾0\underline x_0 e le due derivate miste sono continue in x‾0\underline x_0, allora fxixj(x‾0)=fxjxi(x‾0)f_{x_ix_j}(\underline x_0)=f_{x_jx_i}(\underline x_0). In particolare, se f∈C2(A)f\in C^2(A) la matrice hessiana è simmetrica in ogni punto.

Definizione (differenziale secondo). Se f∈C2(A)f\in C^2(A), il differenziale secondo in x‾0\underline x_0 è la forma quadratica associata all'hessiana: d2f(x‾0)(h‾)=h‾THf(x‾0)h‾=∑i,j=1nfxixj(x‾0)hihj.d^2f(\underline x_0)(\underline h)=\underline h^TH_f(\underline x_0)\underline h=\sum_{i,j=1}^nf_{x_ix_j}(\underline x_0)h_ih_j . Per n=2n=2: d2f(x0,y0)(h,k)=fxxh2+2fxyhk+fyyk2d^2f(x_0,y_0)(h,k)=f_{xx}h^2+2f_{xy}hk+f_{yy}k^2.

Esempio. f=sin⁡(xy)f=\sin(xy): ∇f=(ycos⁡xy, xcos⁡xy)\nabla f=(y\cos xy,\,x\cos xy); fxx=−y2sin⁡xyf_{xx}=-y^2\sin xy, fxy=fyx=cos⁡xy−xysin⁡xyf_{xy}=f_{yx}=\cos xy-xy\sin xy, fyy=−x2sin⁡xyf_{yy}=-x^2\sin xy. In (0,1)(0,1) (xy=0xy=0): Hf(0,1)=(0110)H_f(0,1)=\begin{pmatrix}0&1\\1&0\end{pmatrix} e d2f(0,1)(h,k)=2hkd^2f(0,1)(h,k)=2hk.

Teorema (Taylor, resto di Peano) (D). Se f∈C2(A)f\in C^2(A) e x‾0∈A\underline x_0\in A, allora per h‾→0‾\underline h\to\underline 0 f(x‾0+h‾)=f(x‾0)+∇f(x‾0)⋅h‾+12 h‾THf(x‾0) h‾+o(∣h‾∣2).f(\underline x_0+\underline h)=f(\underline x_0)+\nabla f(\underline x_0)\cdot\underline h+\frac12\,\underline h^TH_f(\underline x_0)\,\underline h+o(|\underline h|^2).

Esempio 1. f=sin⁡(xy)f=\sin(xy) in (0,1)(0,1): f=0f=0, ∇f=(1,0)\nabla f=(1,0), Hf=(0110)H_f=\begin{pmatrix}0&1\\1&0\end{pmatrix}. Con h=xh=x, k=y−1k=y-1: sin⁡(xy)=0+x+0⋅k+12 [ 2 x(y−1) ]+o(x2+(y−1)2)=x+x(y−1)+o(⋅)=xy+o(x2+(y−1)2),\sin(xy)=0+x+0\cdot k+\tfrac12\,[\,2\,x(y-1)\,]+o(x^2+(y-1)^2)=x+x(y-1)+o(\cdot)=xy+o\big(x^2+(y-1)^2\big), cioè il polinomio di Taylor di secondo grado è proprio xyxy.

Forme quadratiche e criteri di definitezza

Definizione. Sia M∈Rn×nM\in\mathbb{R}^{n\times n} simmetrica. La forma quadraticapolinomio omogeneo di secondo grado nelle variabili associata è q:Rn→Rq:\mathbb{R}^n\to\mathbb{R}, q(h‾)=h‾TMh‾=∑i,j=1nmijhihj.q(\underline h)=\underline h^TM\underline h=\sum_{i,j=1}^nm_{ij}h_ih_j .

Definizione (segno). La forma qq (o la matrice MM) si dice

  • definita positiva se q(h‾)>0q(\underline h)>0 per ogni h‾≠0‾\underline h\ne\underline 0; definita negativa se q(h‾)<0q(\underline h)<0 per ogni h‾≠0‾\underline h\ne\underline0;
  • semidefinita positiva se q(h‾)≥0q(\underline h)\ge0 per ogni h‾\underline h; semidefinita negativa se q(h‾)≤0q(\underline h)\le0 per ogni h‾\underline h;
  • indefinita se esistono h‾1,h‾2\underline h^1,\underline h^2 con q(h‾1)>0>q(h‾2)q(\underline h^1)>0>q(\underline h^2).

Grafico interattivo: Livelli q = c delle forme h² − k² (indefinita: iperboli), h² + k²/4 (definita positiva: ellissi) e h² (semidefinita: rette parallele)

Teorema (criterio per n=2n=2) (D). Sia M=(abbc)M=\begin{pmatrix}a&b\\b&c\end{pmatrix} e det⁡M=ac−b2\det M=ac-b^2.

  • qq è definita positiva   ⟺  det⁡M>0\iff\det M>0 e a>0a>0.
  • qq è definita negativa   ⟺  det⁡M>0\iff\det M>0 e a<0a<0.
  • qq è indefinita   ⟺  det⁡M<0\iff\det M<0.
  • Se det⁡M=0\det M=0, qq è semidefinita (positiva se a>0a>0 o c>0c>0, negativa se a<0a<0 o c<0c<0; se a=c=0a=c=0 allora b=0b=0 e M=0M=0). Se a=0a=0 si può usare cc al posto di aa (stesso enunciato).

Teorema (segno e autovalori) (D).

  • qq semidefinita positiva   ⟺  λi≥0\iff\lambda_i\ge0 per ogni ii; semidefinita negativa   ⟺  λi≤0\iff\lambda_i\le0 per ogni ii.
  • qq definita positiva   ⟺  λi>0\iff\lambda_i>0 per ogni ii; definita negativa   ⟺  λi<0\iff\lambda_i<0 per ogni ii.
  • qq indefinita   ⟺  \iff esistono un autovalore >0>0 e uno <0<0.

Corollario (stima). Se λmin⁡\lambda_{\min} e λmax⁡\lambda_{\max} sono il minimo e il massimo autovalore, per ogni h‾\underline h λmin⁡∣h‾∣2≤q(h‾)≤λmax⁡∣h‾∣2.\lambda_{\min}|\underline h|^2\le q(\underline h)\le\lambda_{\max}|\underline h|^2 .

Esempio. q=3h2+2hk+3k2q=3h^2+2hk+3k^2: M=(3113)M=\begin{pmatrix}3&1\\1&3\end{pmatrix}, det⁡(M−λI)=(3−λ)2−1=0\det(M-\lambda I)=(3-\lambda)^2-1=0, λ=2,4\lambda=2,4; autovettori (1,−1)(1,-1) e (1,1)(1,1). Definita positiva; su ∣h‾∣=1|\underline h|=1 il minimo di qq è 22 (in ±(1,−1)2\pm\frac{(1,-1)}{\sqrt2}) e il massimo 44 (in ±(1,1)2\pm\frac{(1,1)}{\sqrt2}).

Criterio di Sylvester (n=3n=3).

  • MM è definita positiva   ⟺  D1>0\iff D_1>0, D2>0D_2>0, D3>0D_3>0.
  • MM è definita negativa   ⟺  D1<0\iff D_1<0, D2>0D_2>0, D3<0D_3<0 (i segni si alternano).
  • Se det⁡M≠0\det M\ne0 e i minori non rientrano in nessuno dei due schemi, la forma è indefinita (nessun autovalore è nullo e non hanno tutti lo stesso segno).

Esempio 1. M=(210121012)M=\begin{pmatrix}2&1&0\\1&2&1\\0&1&2\end{pmatrix}: D1=2>0D_1=2>0; D2=4−1=3>0D_2=4-1=3>0; D3=2⋅(4−1)−1⋅(2−0)=6−2=4>0D_3=2\cdot(4-1)-1\cdot(2-0)=6-2=4>0. Definita positiva (gli autovalori sono 22, 2±22\pm\sqrt2, tutti >0>0). Esempio 2. M=(120210001)M=\begin{pmatrix}1&2&0\\2&1&0\\0&0&1\end{pmatrix}: D1=1>0D_1=1>0, D2=1−4=−3<0D_2=1-4=-3<0: non definita positiva né negativa; det⁡M=−3<0\det M=-3<0. Indefinita (autovalori 3,−1,13,-1,1). Esempio 3. M=(−2101−2000−1)M=\begin{pmatrix}-2&1&0\\1&-2&0\\0&0&-1\end{pmatrix}: D1=−2<0D_1=-2<0, D2=4−1=3>0D_2=4-1=3>0, D3=3⋅(−1)=−3<0D_3=3\cdot(-1)=-3<0: definita negativa (autovalori −1,−1,−3-1,-1,-3).

4. Ottimizzazione libera e funzioni convesse

Massimi e minimi liberi e natura dei punti critici

Definizione (estremi assoluti). x0\mathbf{x}_0 è un punto di massimo assoluto (o globale) se f(x0)≥f(x)f(\mathbf{x}_0) \ge f(\mathbf{x}) per ogni x∈X\mathbf{x} \in X; di minimo assoluto se f(x0)≤f(x)f(\mathbf{x}_0) \le f(\mathbf{x}) per ogni x∈X\mathbf{x} \in X. Il numero f(x0)f(\mathbf{x}_0) si chiama valore massimo (minimo) e si scrive max⁡Xf\max_X f (min⁡Xf\min_X f).

Definizione (estremi relativi). x0\mathbf{x}_0 è un punto di massimo relativo (o locale) se esiste un intorno sferico UU di x0\mathbf{x}_0 tale che f(x0)≥f(x)f(\mathbf{x}_0) \ge f(\mathbf{x}) per ogni x∈U∩X\mathbf{x} \in U \cap X; per il minimo relativo si rovescia la disuguaglianza. L'estremo è forte (o stretto) se la disuguaglianza è stretta per ogni x≠x0\mathbf{x} \ne \mathbf{x}_0 (in XX per gli assoluti, in U∩XU \cap X per i relativi).

Esempio. f(x,y)=x2+y2f(x,y) = x^2 + y^2 ha in (0,0)(0,0) un minimo assoluto stretto: f(x,y)>0=f(0,0)f(x,y) > 0 = f(0,0) per ogni (x,y)≠(0,0)(x,y) \ne (0,0). Invece f(x,y)=cos⁡x+cos⁡yf(x,y) = \cos x + \cos y ha in (0,0)(0,0) un massimo assoluto (vale 22, e cos⁡≤1\cos \le 1 ovunque) ma ha massimi assoluti anche in (2π,0)(2\pi, 0), (0,2π)(0, 2\pi), ...: l'estremo è assoluto ma non stretto.

Teorema (Fermat). Sia A⊆RnA \subseteq \mathbb{R}^n aperto, f:A→Rf : A \to \mathbb{R}, e sia x0∈A\mathbf{x}_0 \in A un punto di massimo o di minimo relativo. Se ff è derivabile in x0\mathbf{x}_0 (esistono tutte le derivate parziali), allora ∇f(x0)=0,cioeˋ∂f∂xk(x0)=0per k=1,…,n.\nabla f(\mathbf{x}_0) = \mathbf{0}, \qquad \text{cioè} \qquad \frac{\partial f}{\partial x_k}(\mathbf{x}_0) = 0 \quad \text{per } k = 1, \dots, n.

Esempio. f(x,y)=x2+y2−2xf(x,y) = x^2 + y^2 - 2x ha ∇f=(2x−2, 2y)\nabla f = (2x - 2,\ 2y), che si annulla solo in (1,0)(1, 0). Infatti f(x,y)=(x−1)2+y2−1≥−1=f(1,0)f(x,y) = (x-1)^2 + y^2 - 1 \ge -1 = f(1,0): lì c'è il minimo assoluto, e vale −1-1.

Definizione (punto critico). Un punto x0\mathbf{x}_0 in cui ff è derivabile e ∇f(x0)=0\nabla f(\mathbf{x}_0) = \mathbf{0} si chiama punto criticochiamato anche punto stazionario: il gradiente è il vettore nullo (o stazionario) di ff.

Definizione (punto di sella). Un punto critico x0\mathbf{x}_0 che non è né di massimo né di minimo relativo si chiama punto di sella: in ogni intorno di x0\mathbf{x}_0 la differenza f(x)−f(x0)f(\mathbf{x}) - f(\mathbf{x}_0) assume sia valori >0> 0 sia valori <0< 0.

Teorema (condizioni sufficienti del secondo ordine). Sia AA aperto, f∈C2(A)f \in C^2(A), ∇f(x0)=0\nabla f(\mathbf{x}_0) = \mathbf{0} e H=Hf(x0)H = H_f(\mathbf{x}_0).

  1. Se HH è definita positiva (q(h)>0q(\mathbf{h}) > 0 per ogni h≠0\mathbf{h} \ne \mathbf{0}), x0\mathbf{x}_0 è un punto di minimo relativo forte.
  2. Se HH è definita negativa, x0\mathbf{x}_0 è un punto di massimo relativo forte.
  3. Se HH è indefinita (qq assume valori di entrambi i segni), x0\mathbf{x}_0 è un punto di sella.
  4. Se HH è solo semidefinita (o nulla), il criterio non decide.

Esempio. f(x,y)=x2+3y2f(x,y) = x^2 + 3y^2 ha ∇f(0,0)=0\nabla f(0,0) = 0 e H=(2006)H = \begin{pmatrix} 2 & 0 \\ 0 & 6\end{pmatrix}, con autovalori 2,6>02, 6 > 0: definita positiva, quindi minimo relativo forte in (0,0)(0,0).

Corollario (n=2n = 2). Sia ∇f(x0,y0)=(0,0)\nabla f(x_0, y_0) = (0,0) e f∈C2f \in C^2.

  • det⁡H>0\det H > 0 e fxx>0f_{xx} > 0: minimo relativo forte.
  • det⁡H>0\det H > 0 e fxx<0f_{xx} < 0: massimo relativo forte.
  • det⁡H<0\det H < 0: punto di sella.
  • det⁡H=0\det H = 0: il criterio non decide.

Esempio completo. f(x,y)=x3+y3−3xyf(x,y) = x^3 + y^3 - 3xy.

  • Gradiente: ∇f=(3x2−3y, 3y2−3x)\nabla f = (3x^2 - 3y,\ 3y^2 - 3x). Dal sistema y=x2y = x^2 e x=y2x = y^2 si ha x=x4x = x^4, cioè x(x3−1)=0x(x^3 - 1) = 0: x=0x = 0 oppure x=1x = 1. Punti critici (0,0)(0,0) e (1,1)(1,1).
  • Hessiana: H(x,y)=(6x−3−36y)H(x,y) = \begin{pmatrix} 6x & -3 \\ -3 & 6y\end{pmatrix}.
  • In (0,0)(0,0): H=(0−3−30)H = \begin{pmatrix} 0 & -3 \\ -3 & 0\end{pmatrix}, det⁡H=−9<0\det H = -9 < 0: sella, con f(0,0)=0f(0,0) = 0.
  • In (1,1)(1,1): H=(6−3−36)H = \begin{pmatrix} 6 & -3 \\ -3 & 6\end{pmatrix}, det⁡H=27>0\det H = 27 > 0 e fxx=6>0f_{xx} = 6 > 0: minimo relativo forte, con f(1,1)=−1f(1,1) = -1. (Gli autovalori sono 33 e 99.) Il minimo non è assoluto: f(x,0)=x3→−∞f(x, 0) = x^3 \to -\infty per x→−∞x \to -\infty.

Teorema (Weierstrass generalizzato). Sia C⊆RnC \subseteq \mathbb{R}^n chiuso e illimitato e f:C→Rf : C \to \mathbb{R} continua. Se lim⁡∥x∥→+∞x∈Cf(x)=+∞,\lim_{\substack{\|\mathbf{x}\| \to +\infty \\ \mathbf{x} \in C}} f(\mathbf{x}) = +\infty, allora ff ammette minimo assoluto in CC. Se il limite è −∞-\infty, ammette massimo assoluto.

Esempio. f(x,y)=x2+y2f(x,y) = x^2 + y^2 su R2\mathbb{R}^2: tende a +∞+\infty all'infinito e infatti ha minimo assoluto 00; non ha massimo, perché è illimitata superiormente.

Funzioni convesse di più variabili e minimi globali

Definizione (insieme convesso). A⊆RnA \subseteq \mathbb{R}^n è convesso se per ogni x,y∈A\mathbf{x}, \mathbf{y} \in A il segmento [x,y][\mathbf{x}, \mathbf{y}] è contenuto in AA. È strettamente convesso se, per x≠y\mathbf{x} \ne \mathbf{y} in AA, il segmento aperto (senza estremi) sta nella parte interna di AA.

Esempio. Il disco x2+y2≤1x^2 + y^2 \le 1 è convesso; il semidisco {x2+y2≤1, y≥0}\{x^2 + y^2 \le 1,\ y \ge 0\} è convesso ma non strettamente: due punti sul diametro sono uniti da un segmento che sta sul bordo.

Definizione (funzione convessa). Sia A⊆RnA \subseteq \mathbb{R}^n convesso e f:A→Rf : A \to \mathbb{R}. ff è convessa se per ogni x,y∈A\mathbf{x}, \mathbf{y} \in A e ogni t∈[0,1]t \in [0,1] f((1−t)x+ty)≤(1−t) f(x)+t f(y).f\big((1-t)\mathbf{x} + t\mathbf{y}\big) \le (1-t)\,f(\mathbf{x}) + t\,f(\mathbf{y}). È strettamente convessa se per x≠y\mathbf{x} \ne \mathbf{y} e t∈(0,1)t \in (0,1) vale la disuguaglianza stretta. ff è concava se −f-f è convessa (disuguaglianza rovesciata).

Esempio. f(x)=x2f(x) = x^2 è strettamente convessa: per x=0x = 0, y=2y = 2, t=12t = \frac12 si ha f(1)=1<12f(0)+12f(2)=2f(1) = 1 < \frac12 f(0) + \frac12 f(2) = 2.

Teorema (primo ordine). Sia AA aperto e convesso e f:A→Rf : A \to \mathbb{R} differenziabile. Allora ff è convessa se e solo se per ogni x,x0∈A\mathbf{x}, \mathbf{x}_0 \in A f(x)≥f(x0)+∇f(x0)⋅(x−x0).f(\mathbf{x}) \ge f(\mathbf{x}_0) + \nabla f(\mathbf{x}_0)\cdot(\mathbf{x} - \mathbf{x}_0). Inoltre ff è strettamente convessa se e solo se la disuguaglianza è stretta per x≠x0\mathbf{x} \ne \mathbf{x}_0.

Esempio. f(x,y)=x2+y2f(x,y) = x^2 + y^2 in x0=(1,0)\mathbf{x}_0 = (1,0): il piano tangente è z=1+2(x−1)=2x−1z = 1 + 2(x - 1) = 2x - 1, e x2+y2−(2x−1)=(x−1)2+y2≥0x^2 + y^2 - (2x - 1) = (x-1)^2 + y^2 \ge 0: il grafico sta sopra.

Teorema (secondo ordine). Sia AA aperto e convesso e f∈C2(A)f \in C^2(A).

  1. ff è convessa in AA   ⟺  \iff Hf(x)H_f(\mathbf{x}) è semidefinita positiva per ogni x∈A\mathbf{x} \in A.
  2. Se Hf(x)H_f(\mathbf{x}) è definita positiva per ogni x∈A\mathbf{x} \in A, allora ff è strettamente convessa. Il viceversa è falso.

Teorema (minimi globali). Sia AA aperto e convesso, f:A→Rf : A \to \mathbb{R} convessa e differenziabile, e x0∈A\mathbf{x}_0 \in A con ∇f(x0)=0\nabla f(\mathbf{x}_0) = \mathbf{0}. Allora x0\mathbf{x}_0 è un punto di minimo globale. Se ff è strettamente convessa, x0\mathbf{x}_0 è l'unico punto di minimo globale (e un punto critico, se esiste, è unico). Per ff concava: massimo globale.

Esempio. f(x,y,z)=x2+y2+z2−2x−2z−5f(x,y,z) = x^2 + y^2 + z^2 - 2x - 2z - 5. Il gradiente (2x−2, 2y, 2z−2)(2x - 2,\ 2y,\ 2z - 2) si annulla in (1,0,1)(1,0,1); H=2I3H = 2I_3 è definita positiva ovunque. Quindi (1,0,1)(1,0,1) è l'unico minimo globale e min⁡f=1+0+1−2−2−5=−7\min f = 1 + 0 + 1 - 2 - 2 - 5 = -7. Non c'è massimo, perché f→+∞f \to +\infty lungo ogni retta.

5. Funzioni implicite ed estremi vincolati

Funzioni implicite e teorema di Dini nel piano

Definizione (funzione definita implicitamente). Un intervallo aperto II e una funzione g:I→Rg : I \to \mathbb{R} sono definiti implicitamente dall'equazione f(x,y)=0f(x,y) = 0 se f(x,g(x))=0per ogni x∈I.f(x, g(x)) = 0 \quad \text{per ogni } x \in I. Cioè il grafico di gg è contenuto in L0L_0. Analogamente per x=h(y)x = h(y) con f(h(y),y)=0f(h(y), y) = 0.

Teorema (Dini, o delle funzioni implicite). Sia A⊆R2A \subseteq \mathbb{R}^2 aperto, f∈C1(A)f \in C^1(A), (x0,y0)∈A(x_0, y_0) \in A con f(x0,y0)=0efy(x0,y0)≠0.f(x_0, y_0) = 0 \qquad \text{e} \qquad f_y(x_0, y_0) \ne 0. Allora esistono un intervallo aperto I∋x0I \ni x_0 e una unica funzione g:I→Rg : I \to \mathbb{R} (con valori in un opportuno intorno di y0y_0) tali che g(x0)=y0g(x_0) = y_0 e f(x,g(x))=0f(x, g(x)) = 0 per ogni x∈Ix \in I. Inoltre g∈C1(I)g \in C^1(I) e g′(x)=−fx(x,g(x))fy(x,g(x))(x∈I).(FDI)g'(x) = -\frac{f_x(x, g(x))}{f_y(x, g(x))} \qquad (x \in I). \qquad \text{(FDI)} Se f∈Ck(A)f \in C^k(A), allora g∈Ck(I)g \in C^k(I).

Esempio. La circonferenza f=x2+y2−1f = x^2 + y^2 - 1 ha ∇f=(2x,2y)\nabla f = (2x, 2y), nullo solo nell'origine che non sta sulla curva. Quindi tutti i punti sono non critici. Dove y≠0y \ne 0 (fy=2y≠0f_y = 2y \ne 0) si esplicita y=g(x)y = g(x) con g′=−xyg' = -\frac{x}{y}; dove x≠0x \ne 0 si esplicita x=h(y)x = h(y) con h′=−yxh' = -\frac{y}{x}.

Estremi vincolati - metodo diretto e moltiplicatori di Lagrange

Definizione (estremo vincolato). P∈ΓP \in \Gamma è punto di massimo (minimo) relativo vincolato se esiste un intorno UU di PP con f(P)≥f(Q)f(P) \ge f(Q) (≤\le) per ogni Q∈Γ∩UQ \in \Gamma \cap U.

Teorema (moltiplicatori di Lagrange, vincolo in R2\mathbb{R}^2). Sia AA aperto, f,g∈C1(A)f, g \in C^1(A), Γ={g=b}\Gamma = \{g = b\}. Se P∈ΓP \in \Gamma è un punto di estremo relativo vincolato di f∣Γf|_\Gamma e ∇g(P)≠0,\nabla g(P) \ne \mathbf{0}, allora esiste λ∈R\lambda \in \mathbb{R} (il moltiplicatore) tale che ∇f(P)=λ ∇g(P).\nabla f(P) = \lambda\,\nabla g(P).

Grafico interattivo: Vincolo x² + y² = 1 e livelli di f = 2x² + y² − x. Il livello f = 3/4 (verde) tocca la circonferenza in (1/2, ±√3/2): minimi vincolati. Il livello f = 3 (arancio) la tocca in (−1, 0): massimo. Nei punti di tangenza ∇f e ∇g sono paralleli

Definizione (punto critico vincolato). Un punto P∈ΓP \in \Gamma con ∇g(P)≠0\nabla g(P) \ne \mathbf{0} per cui esiste λ\lambda con ∇f(P)=λ∇g(P)\nabla f(P) = \lambda\nabla g(P) si dice punto critico vincolato (o stazionario) di ff su Γ\Gamma.

Esempio (stesso problema di prima). f=2x2+y2−xf = 2x^2 + y^2 - x, g=x2+y2g = x^2 + y^2, b=1b = 1. Il sistema: {4x−1=2λx2y=2λyx2+y2=1\begin{cases} 4x - 1 = 2\lambda x \\ 2y = 2\lambda y \\ x^2 + y^2 = 1\end{cases} La seconda equazione dà y(1−λ)=0y(1 - \lambda) = 0: o y=0y = 0 o λ=1\lambda = 1.

  • y=0y = 0: dal vincolo x=±1x = \pm1; i punti sono (1,0)(1,0) (con λ=32\lambda = \frac32) e (−1,0)(-1,0) (con λ=52\lambda = \frac52).
  • λ=1\lambda = 1: la prima diventa 4x−1=2x4x - 1 = 2x, cioè x=12x = \frac12; dal vincolo y2=34y^2 = \frac34, y=±32y = \pm\frac{\sqrt3}2. Quattro punti, con f(1,0)=1f(1,0) = 1, f(−1,0)=3f(-1,0) = 3, f(12,±32)=34f\left(\frac12, \pm\frac{\sqrt3}2\right) = \frac34: stesso risultato del metodo diretto. La circonferenza è compatta, quindi (Weierstrass) massimo e minimo esistono e sono fra questi candidati: max⁡=3\max = 3, min⁡=34\min = \frac34. Il punto (1,0)(1,0) è un estremo relativo vincolato (massimo locale su Γ\Gamma) ma non assoluto.

6. Funzioni a valori vettoriali e matrice jacobiana

Funzioni a valori vettoriali, matrice jacobiana e trasformazioni di coordinate

Definizione (limite). Se x0\mathbf{x}_0 è punto di accumulazione di AA e L∈Rm\mathbf{L} \in \mathbb{R}^m, lim⁡x→x0f(x)=L  ⟺  lim⁡x→x0∥f(x)−L∥=0.\lim_{\mathbf{x} \to \mathbf{x}_0}\mathbf{f}(\mathbf{x}) = \mathbf{L} \iff \lim_{\mathbf{x} \to \mathbf{x}_0}\|\mathbf{f}(\mathbf{x}) - \mathbf{L}\| = 0.

Proposizione (limite per componenti). f→L\mathbf{f} \to \mathbf{L} se e solo se fi(x)→Lif_i(\mathbf{x}) \to L_i per ogni i=1,…,mi = 1, \dots, m. In particolare f\mathbf{f} è continua se e solo se lo sono tutte le componenti.

Definizione (differenziabilità). f:A→Rm\mathbf{f} : A \to \mathbb{R}^m, A⊆RnA \subseteq \mathbb{R}^n aperto, è differenziabile in x0\mathbf{x}_0 se ogni componente fif_i è differenziabile in x0\mathbf{x}_0, cioè fi(x0+h)−fi(x0)=∇fi(x0)⋅h+o(∥h∥),h→0,i=1,…,m.f_i(\mathbf{x}_0 + \mathbf{h}) - f_i(\mathbf{x}_0) = \nabla f_i(\mathbf{x}_0)\cdot\mathbf{h} + o(\|\mathbf{h}\|), \qquad \mathbf{h} \to \mathbf{0}, \quad i = 1, \dots, m.

Definizione (matrice jacobiana). La matrice jacobiana di f\mathbf{f} in x0\mathbf{x}_0 è la matrice m×nm \times n (mm righe, nn colonne) Df(x0)=Jf(x0)=(∂fi∂xj(x0))i=1,…,mj=1,…,n=(∇f1(x0)T⋮∇fm(x0)T).D\mathbf{f}(\mathbf{x}_0) = J_{\mathbf{f}}(\mathbf{x}_0) = \left(\frac{\partial f_i}{\partial x_j}(\mathbf{x}_0)\right)_{\substack{i = 1,\dots,m \\ j = 1,\dots,n}} = \begin{pmatrix} \nabla f_1(\mathbf{x}_0)^T \\ \vdots \\ \nabla f_m(\mathbf{x}_0)^T\end{pmatrix}. La riga ii è il gradiente della componente fif_i; la colonna jj contiene le derivate di tutte le componenti rispetto a xjx_j.

Corollario (differenziale totale). Se AA è aperto e tutte le componenti sono di classe C1(A)C^1(A), allora f\mathbf{f} è differenziabile in ogni punto di AA (infatti ogni fi∈C1f_i \in C^1 è differenziabile, Derivate parziali, differenziabilità e piano tangenteLa derivata parziale f_{x_i}(x0) è la derivata di t ↦ f(x0 + t e_i); il gradiente ∇f raccoglie le derivate parziali. In più variabili l'esistenza delle derivate parziali NON implica la continuità. f è differenziabile in x0 se f(x0 + h) = f(x0) + a·h + o(|h|); allora f è continua, ha gradiente e a = ∇f(x0) (D). Il grafico ha piano tangente z = f(x0) + ∇f(x0)·(x − x0) e normale (−f_x, −f_y, 1). Teorema del differenziale totale: se le derivate parziali esistono vicino a x0 e sono continue in x0, f è differenziabile; C^1 ⇒ differenziabile ⇒ derivabile e continua, ma i viceversa sono falsi. Regole: Leibniz, catena (f∘r)' = ∇f(r)·r'.Derivate parziali, differenziabilità e piano tangente →).

Esempio. Per f(u,v)=(uv, u+v, u−v)\mathbf{f}(u,v) = (uv,\ u + v,\ u - v): Df(u,v)=(vu111−1)(3×2),Df(1,2)=(21111−1).D\mathbf{f}(u,v) = \begin{pmatrix} v & u \\ 1 & 1 \\ 1 & -1\end{pmatrix} \quad (3 \times 2), \qquad D\mathbf{f}(1,2) = \begin{pmatrix} 2 & 1 \\ 1 & 1 \\ 1 & -1\end{pmatrix}. Vale f(1,2)=(2,3,−1)\mathbf{f}(1,2) = (2, 3, -1) e, per un piccolo spostamento h=(0,1; −0,1)\mathbf{h} = (0{,}1;\ -0{,}1), l'approssimazione lineare è (2,3,−1)+Df(1,2)h=(2,3,−1)+(0,1; 0; 0,2)=(2,1; 3; −0,8)(2,3,-1) + D\mathbf{f}(1,2)\mathbf{h} = (2,3,-1) + (0{,}1;\ 0;\ 0{,}2) = (2{,}1;\ 3;\ -0{,}8), e il valore esatto f(1,1; 1,9)=(2,09; 3; −0,8)\mathbf{f}(1{,}1;\ 1{,}9) = (2{,}09;\ 3;\ -0{,}8) differisce di 0,010{,}01.

Teorema (regola della catena). Siano f:A⊆Rn→B⊆Rm\mathbf{f} : A \subseteq \mathbb{R}^n \to B \subseteq \mathbb{R}^m differenziabile in x0\mathbf{x}_0 e g:B→Rk\mathbf{g} : B \to \mathbb{R}^k differenziabile in f(x0)\mathbf{f}(\mathbf{x}_0). Allora g∘f\mathbf{g}\circ\mathbf{f} è differenziabile in x0\mathbf{x}_0 e D(g∘f)(x0)=Dg(f(x0))⋅Df(x0).D(\mathbf{g}\circ\mathbf{f})(\mathbf{x}_0) = D\mathbf{g}\big(\mathbf{f}(\mathbf{x}_0)\big)\cdot D\mathbf{f}(\mathbf{x}_0).

Esempio (numerico). Con f(u,v)=(uv, u+v, u−v)\mathbf{f}(u,v) = (uv,\ u + v,\ u - v) e g(x,y,z)=(x+yz, x2−z)\mathbf{g}(x,y,z) = (x + yz,\ x^2 - z) (g:R3→R2\mathbf{g} : \mathbb{R}^3 \to \mathbb{R}^2), calcoliamo D(g∘f)(1,2)D(\mathbf{g}\circ\mathbf{f})(1,2). Si ha f(1,2)=(2,3,−1)\mathbf{f}(1,2) = (2,3,-1), Dg(x,y,z)=(1zy2x0−1),Dg(2,3,−1)=(1−1340−1),D\mathbf{g}(x,y,z) = \begin{pmatrix} 1 & z & y \\ 2x & 0 & -1\end{pmatrix}, \qquad D\mathbf{g}(2,3,-1) = \begin{pmatrix} 1 & -1 & 3 \\ 4 & 0 & -1\end{pmatrix}, e quindi D(g∘f)(1,2)=(1−1340−1)(21111−1)=(4−375).D(\mathbf{g}\circ\mathbf{f})(1,2) = \begin{pmatrix} 1 & -1 & 3 \\ 4 & 0 & -1\end{pmatrix}\begin{pmatrix} 2 & 1 \\ 1 & 1 \\ 1 & -1\end{pmatrix} = \begin{pmatrix} 4 & -3 \\ 7 & 5\end{pmatrix}. (Controllo: g∘f=(uv+(u+v)(u−v), u2v2−u+v)\mathbf{g}\circ\mathbf{f} = (uv + (u+v)(u-v),\ u^2v^2 - u + v) e le sue derivate in (1,2)(1,2) sono: ∂u(uv+u2−v2)=v+2u=4\partial_u(uv + u^2 - v^2) = v + 2u = 4, ∂v(uv+u2−v2)=u−2v=−3\partial_v(uv + u^2 - v^2) = u - 2v = -3, ∂u(u2v2−u+v)=2uv2−1=7\partial_u(u^2v^2 - u + v) = 2uv^2 - 1 = 7, ∂v(u2v2−u+v)=2u2v+1=5\partial_v(u^2v^2 - u + v) = 2u^2v + 1 = 5 ✓.)

7. Integrali doppi

Integrali doppi, domini semplici e formule di riduzione

Definizione (integrale doppio su un rettangolo). ff è integrabile (secondo Riemann) su RR se lim⁡n→∞Sn\lim_{n\to\infty}S_n esiste finito e non dipende dalla scelta dei punti PhkP_{hk}. In tal caso ∬Rf(x,y) dx dy=lim⁡n→∞Sn.\iint_R f(x,y)\,dx\,dy = \lim_{n\to\infty}S_n.

Esempio (funzione costante). Se f≡kf \equiv k, ogni Sn=k∑∣Ihk∣=k (b−a)(d−c)S_n = k\sum|I_{hk}| = k\,(b-a)(d-c), quindi ∬Rk dx dy=k ∣R∣\iint_R k\,dx\,dy = k\,|R|: il volume di un parallelepipedo è area di base per altezza.

Teorema (integrabilità). Se ff è continua su un rettangolo chiuso RR, allora ff è integrabile su RR.

Teorema (Fubini sui rettangoli). Se f∈C(R)f \in C(R), R=[a,b]×[c,d]R = [a,b]\times[c,d], ∬Rf dx dy=∫cd(∫abf(x,y) dx)dy=∫ab(∫cdf(x,y) dy)dx.\iint_R f\,dx\,dy = \int_c^d\left(\int_a^b f(x,y)\,dx\right)dy = \int_a^b\left(\int_c^d f(x,y)\,dy\right)dx.

Esempio. I=∬[0,1]×[0,2]11+x+y dx dy\displaystyle I = \iint_{[0,1]\times[0,2]}\frac{1}{1+x+y}\,dx\,dy. La funzione non è definita sulla retta x+y=−1x+y = -1, che non incontra il rettangolo (lì 1+x+y≥11+x+y \ge 1): è continua su RR. Si integra prima in yy: ∫02dy1+x+y=[ln⁡(1+x+y)]y=0y=2=ln⁡(3+x)−ln⁡(1+x).\int_0^2\frac{dy}{1+x+y} = \Big[\ln(1+x+y)\Big]_{y=0}^{y=2} = \ln(3+x) - \ln(1+x). Poi, per parti (∫ln⁡u dx=xln⁡u−∫xudx\int\ln u\,dx = x\ln u - \int\frac{x}{u}dx con u=3+xu = 3 + x o u=1+xu = 1 + x) e usando x3+x=1−33+x\frac{x}{3+x} = 1 - \frac{3}{3+x}, x1+x=1−11+x\frac{x}{1+x} = 1 - \frac{1}{1+x}: I=∫01[ln⁡(3+x)−ln⁡(1+x)] dx=3ln⁡43≈0,863.I = \int_0^1[\ln(3+x) - \ln(1+x)]\,dx = 3\ln\frac43 \approx 0{,}863. (Nell'altro ordine si ottiene lo stesso risultato.)

Definizione. ff è integrabile su Ω\Omega se f~\tilde f è integrabile su RR, e si pone ∬Ωf=∬Rf~\displaystyle\iint_\Omega f = \iint_R\tilde f.

Definizione (dominio yy-semplice, o normale rispetto all'asse xx). Esistono [a,b][a,b] e g1,g2∈C[a,b]g_1, g_2\in C[a,b] con g1≤g2g_1\le g_2 tali che Ω={(x,y):a≤x≤b, g1(x)≤y≤g2(x)}.\Omega = \{(x,y) : a\le x\le b,\ g_1(x)\le y\le g_2(x)\}. Dominio xx-semplice (normale rispetto all'asse yy): Ω={(x,y):c≤y≤d, h1(y)≤x≤h2(y)}\Omega = \{(x,y) : c\le y\le d,\ h_1(y)\le x\le h_2(y)\} con h1≤h2h_1\le h_2 continue. Dominio regolare: unione finita di domini semplici.

Teorema (integrabilità). Se Ω\Omega è un dominio regolare e f∈C(Ω)f \in C(\Omega), allora ff è integrabile su Ω\Omega.

Teorema (formule di riduzione). Sia f∈C(Ω)f\in C(\Omega).

  • Se Ω={a≤x≤b, g1(x)≤y≤g2(x)}\Omega = \{a\le x\le b,\ g_1(x)\le y\le g_2(x)\} è yy-semplice: ∬Ωf dx dy=∫ab(∫g1(x)g2(x)f(x,y) dy)dx.\iint_\Omega f\,dx\,dy = \int_a^b\left(\int_{g_1(x)}^{g_2(x)}f(x,y)\,dy\right)dx.
  • Se Ω={c≤y≤d, h1(y)≤x≤h2(y)}\Omega = \{c\le y\le d,\ h_1(y)\le x\le h_2(y)\} è xx-semplice: ∬Ωf dx dy=∫cd(∫h1(y)h2(y)f(x,y) dx)dy.\iint_\Omega f\,dx\,dy = \int_c^d\left(\int_{h_1(y)}^{h_2(y)}f(x,y)\,dx\right)dy.

Definizione (area). Un insieme limitato Ω\Omega è misurabile (secondo Peano-Jordan) se la funzione costante 11 è integrabile su Ω\Omega, e la sua area (misura bidimensionale) è ∣Ω∣=∬Ω1 dx dy.|\Omega| = \iint_\Omega 1\,dx\,dy.

Teorema. Se Ω\Omega è misurabile, ff è limitata su Ω\Omega e continua tranne che su un insieme di area nulla, allora ff è integrabile.

Esempio 1. f(x,y)=xyx2+y2f(x,y) = \frac{xy}{x^2+y^2} per (x,y)≠(0,0)(x,y)\ne(0,0), f(0,0)=0f(0,0) = 0, sul disco unitario. L'unico punto di discontinuità è l'origine (area nulla). Inoltre ff è limitata: da (∣x∣−∣y∣)2≥0(|x|-|y|)^2\ge0 si ha x2+y2≥2∣xy∣x^2+y^2\ge2|xy| e ∣f∣≤12|f|\le\frac12. Quindi è integrabile.

Definizione (media integrale). Se ∣Ω∣≠0|\Omega|\ne0, la media integrale di ff su Ω\Omega è m=1∣Ω∣∬Ωfm = \frac{1}{|\Omega|}\iint_\Omega f.

Teorema (della media integrale). Se Ω\Omega è un dominio regolare connesso e f∈C(Ω)f\in C(\Omega), esiste (x0,y0)∈Ω(x_0,y_0)\in\Omega con f(x0,y0)=1∣Ω∣∬Ωf dx dy.f(x_0,y_0) = \frac{1}{|\Omega|}\iint_\Omega f\,dx\,dy.

Cambiamento di variabili negli integrali doppi, integrali generalizzati e gaussiana

Definizione (diffeomorfismo). T:D′→DT : D'\to D è un diffeomorfismo (globale) se è biiettiva, di classe C1C^1 e la sua inversa è di classe C1C^1.

Teorema (cambiamento di variabili). Sia D⊆R2D\subseteq\mathbb{R}^2 un dominio regolare, T:D′→DT : D'\to D, T(u,v)=(x(u,v),y(u,v))T(u,v) = (x(u,v), y(u,v)) un diffeomorfismo e f∈C(D)f\in C(D). Allora ∬Df(x,y) dx dy=∬D′f(x(u,v),y(u,v))∣∂(x,y)∂(u,v)∣du dv,∣∂(x,y)∂(u,v)∣=∣det⁡DT(u,v)∣.\iint_Df(x,y)\,dx\,dy = \iint_{D'}f\big(x(u,v),y(u,v)\big)\left|\frac{\partial(x,y)}{\partial(u,v)}\right|du\,dv,\qquad \left|\frac{\partial(x,y)}{\partial(u,v)}\right| = |\det DT(u,v)|.

Esempio (lineare). Area dell'ellisse x2a2+y2b2≤1\frac{x^2}{a^2}+\frac{y^2}{b^2}\le1. Con x=aux = au, y=bvy = bv la matrice è diag⁡(a,b)\operatorname{diag}(a,b), det⁡=ab\det = ab, e l'ellisse corrisponde al disco u2+v2≤1u^2+v^2\le1 di area π\pi. Quindi ∣E∣=ab⋅π=πab|E| = ab\cdot\pi = \pi ab.

Proposizione. Sia DD misurabile e simmetrico rispetto all'asse yy (se (x,y)∈D(x,y)\in D allora (−x,y)∈D(-x,y)\in D) e ff integrabile su DD.

  • Se ff è dispari in xx, cioè f(−x,y)=−f(x,y)f(-x,y) = -f(x,y): ∬Df=0\iint_Df = 0.
  • Se ff è pari in xx: ∬Df=2∬D∩{x≥0}f\iint_Df = 2\iint_{D\cap\{x\ge0\}}f.

Definizione (successione invadente). Una successione di insiemi DjD_j invade R2\mathbb{R}^2 se ogni DjD_j è un dominio regolare, Dj⊆Dj+1D_j\subseteq D_{j+1}, e ⋃jDj=R2\bigcup_jD_j = \mathbb{R}^2. Esempi: i dischi B‾j(0)\overline{B}_j(\mathbf{0}) e i quadrati [−j,j]2[-j,j]^2.

Definizione. ∬R2f dx dy:=lim⁡j→∞∬Djf dx dy.\displaystyle\iint_{\mathbb{R}^2}f\,dx\,dy := \lim_{j\to\infty}\iint_{D_j}f\,dx\,dy.

Teorema. Se f≥0f\ge0 è continua, il valore del limite non dipende dalla successione invadente.

Esempio (potenze di ρ\rho). Sia f(x,y)=1(x2+y2)αf(x,y) = \frac{1}{(x^2+y^2)^\alpha}.

  • Vicino all'origine: ∬ε≤ρ≤1f=2π∫ε1ρ1−2αdρ\iint_{\varepsilon\le\rho\le1}f = 2\pi\int_\varepsilon^1\rho^{1-2\alpha}d\rho. Per ε→0\varepsilon\to0 il limite è finito se e solo se 1−2α>−11-2\alpha>-1, cioè α<1\alpha<1 (per α=1\alpha=1 si ha 2πln⁡1ε→∞2\pi\ln\frac1\varepsilon\to\infty).
  • All'infinito: ∬1≤ρ≤Rf=2π∫1Rρ1−2αdρ\iint_{1\le\rho\le R}f = 2\pi\int_1^R\rho^{1-2\alpha}d\rho ha limite finito per R→∞R\to\infty se e solo se 1−2α<−11-2\alpha<-1, cioè α>1\alpha>1. Quindi nessun α\alpha rende ff integrabile sia vicino all'origine sia all'infinito, ma α<1\alpha<1 lo è vicino a 00 e α>1\alpha>1 vicino a ∞\infty.

Teorema (integrale di Gauss). ∬R2e−x2+y22 dx dy=2π\displaystyle\iint_{\mathbb{R}^2}e^{-\frac{x^2+y^2}2}\,dx\,dy = 2\pi e di conseguenza ∫−∞+∞e−x22 dx=2π\displaystyle\int_{-\infty}^{+\infty}e^{-\frac{x^2}2}\,dx = \sqrt{2\pi}.

Esempio (con un'ellisse). ∬R2e−(x2+4y2)dx dy\iint_{\mathbb{R}^2}e^{-(x^2+4y^2)}dx\,dy. Con x=ux=u, y=v/2y=v/2 (determinante 12\frac12): 12∬e−(u2+v2)du dv=12(∫e−t2dt)2=12π=π2\frac12\iint e^{-(u^2+v^2)}du\,dv = \frac12\left(\int e^{-t^2}dt\right)^2 = \frac12\pi = \frac\pi2.

8. Integrali tripli e solidi di rotazione

Integrali tripli, coordinate cilindriche e sferiche e solidi di rotazione

Definizione (dominio zz-semplice). E={(x,y,z):(x,y)∈D, h1(x,y)≤z≤h2(x,y)}E = \{(x,y,z) : (x,y)\in D,\ h_1(x,y)\le z\le h_2(x,y)\}, con D⊂R2D\subset\mathbb{R}^2 dominio regolare (la proiezione di EE sul piano xyxy) e h1≤h2h_1\le h_2 continue su DD.

Teorema (integrazione per fili). Se f∈C(E)f\in C(E) e EE è zz-semplice, ∭Ef dx dy dz=∬D(∫h1(x,y)h2(x,y)f(x,y,z) dz)dx dy.\iiint_Ef\,dx\,dy\,dz = \iint_D\left(\int_{h_1(x,y)}^{h_2(x,y)}f(x,y,z)\,dz\right)dx\,dy.

Definizione (dominio stratificato rispetto a zz). E={(x,y,z):a≤z≤b, (x,y)∈Ez}E = \{(x,y,z) : a\le z\le b,\ (x,y)\in E_z\}, dove Ez⊂R2E_z\subset\mathbb{R}^2 è la sezione (strato) di EE a quota zz, un dominio regolare.

Teorema (integrazione per strati). Se f∈C(E)f\in C(E), ∭Ef dx dy dz=∫ab(∬Ezf(x,y,z) dx dy)dz.\iiint_Ef\,dx\,dy\,dz = \int_a^b\left(\iint_{E_z}f(x,y,z)\,dx\,dy\right)dz.

Teorema. Se T:E′→ET : E'\to E, (u,v,w)↦(x,y,z)(u,v,w)\mapsto(x,y,z), è un diffeomorfismo tra domini regolari di R3\mathbb{R}^3 e f∈C(E)f\in C(E), ∭Ef dx dy dz=∭E′f(T(u,v,w))∣det⁡DT(u,v,w)∣du dv dw.\iiint_Ef\,dx\,dy\,dz = \iiint_{E'}f(T(u,v,w))\left|\det DT(u,v,w)\right|du\,dv\,dw.

Esempio (ellissoide). L'ellissoide x2a2+y2b2+z2c2≤1\frac{x^2}{a^2}+\frac{y^2}{b^2}+\frac{z^2}{c^2}\le1 è l'immagine della palla unitaria con x=aux=au, y=bvy=bv, z=cwz=cw (determinante abcabc): ∣E∣=abc⋅43π=43πabc|E| = abc\cdot\frac43\pi = \frac43\pi abc.

Teorema (Pappo-Guldino, volume). ∣R(A)∣3=θˉ xB ∣A∣2,\displaystyle|R(A)|_3 = \bar\theta\,x_B\,|A|_2, cioè il volume del solido è l'area di AA per la lunghezza dell'arco percorso dal baricentro (xB,zB)(x_B,z_B) di AA (θˉxB\bar\theta x_B, ossia 2πxB2\pi x_B per una rotazione completa).

Esempio (toro). Il disco di raggio rr con centro a distanza R>rR>r dall'asse, ruotato di 2π2\pi: ∣A∣=πr2|A|=\pi r^2, xB=Rx_B = R (centro del disco): V=πr2⋅2πR=2π2r2RV = \pi r^2\cdot2\pi R = 2\pi^2r^2R.

9. Campi vettoriali, lavoro e forme differenziali

Campi vettoriali, lavoro e forme differenziali

Definizione (campo vettoriale). Sia Ω⊆Rn\Omega\subseteq\mathbb{R}^n un aperto, con n=2n=2 o n=3n=3. Un campo vettoriale è una funzione F:Ω→Rn\mathbf{F}:\Omega\to\mathbb{R}^n, cioè una regola che a ogni punto P∈ΩP\in\Omega associa un vettore F(P)\mathbf{F}(P) "applicato" in PP. In componenti F=(F1,…,Fn)\mathbf{F}=(F_1,\dots,F_n), dove ogni Fi:Ω→RF_i:\Omega\to\mathbb{R} è una funzione scalare. Si dice di classe C1C^1 se lo sono tutte le componenti.

Esempio. F(x,y)=(−y,x)\mathbf{F}(x,y)=(-y,x) in R2\mathbb{R}^2: nel punto (1,0)(1,0) vale (0,1)(0,1), in (0,1)(0,1) vale (−1,0)(-1,0), in (2,0)(2,0) vale (0,2)(0,2). È il campo "di rotazione" attorno all'origine.

Grafico interattivo: Campo di rotazione F = (-y, x): i vettori (disegnati in scala 0,3) sono perpendicolari al raggio e più lunghi lontano dall'origine

Definizione (forma differenziale lineare). Una forma differenziale lineare (o 1-forma) su Ω\Omega è una funzione che a ogni punto associa un covettore: ω=F1 dx+F2 dy+F3 dz\omega=F_1\,dx+F_2\,dy+F_3\,dz, dove dx,dy,dzdx,dy,dz sono i covettori costanti della base duale, dx=(1,0,0)dx=(1,0,0), dy=(0,1,0)dy=(0,1,0), dz=(0,0,1)dz=(0,0,1).

Esempio. ω=3xyz dx+xy2 dy+ezsin⁡x dz\omega=3xyz\,dx+xy^2\,dy+e^z\sin x\,dz corrisponde al campo F=(3xyz, xy2, ezsin⁡x)\mathbf{F}=(3xyz,\ xy^2,\ e^z\sin x); applicata al campo G=(2z, exy, xlog⁡y)\mathbf{G}=(2z,\ e^xy,\ x\log y) dà ω(G)=6xyz2+xexy3+xezsin⁡x log⁡y\omega(\mathbf{G})=6xyz^2+xe^xy^3+xe^z\sin x\,\log y.

Definizione (lavoro). Sia F:Ω→Rn\mathbf{F}:\Omega\to\mathbb{R}^n un campo continuo e γ\gamma una curva regolare a tratti con parametrizzazione r:[a,b]→Ω\mathbf{r}:[a,b]\to\Omega, t↦r(t)t\mapsto\mathbf{r}(t). Il lavoro di F\mathbf{F} lungo γ\gamma (integrale curvilineo di seconda specie, integrale di ω\omega su γ\gamma) è L=∫γω=∫γF⋅dr=∫abF(r(t))⋅r′(t) dt.L=\int_\gamma\omega=\int_\gamma\mathbf{F}\cdot d\mathbf{r}=\int_a^b\mathbf{F}(\mathbf{r}(t))\cdot\mathbf{r}'(t)\,dt . Per una curva regolare a tratti si somma sui pezzi. Se γ\gamma è chiusa il lavoro si chiama anche circuitazione e si scrive ∮γ\oint_\gamma.

Esempio. F=(2xy, x2)\mathbf{F}=(2xy,\ x^2) lungo r(t)=(t,t2)\mathbf{r}(t)=(t,t^2), t∈[0,1]t\in[0,1]. Si ha r′=(1,2t)\mathbf{r}'=(1,2t), F(r(t))=(2t3, t2)\mathbf{F}(\mathbf{r}(t))=(2t^3,\ t^2), quindi F⋅r′=2t3+2t3=4t3\mathbf{F}\cdot\mathbf{r}'=2t^3+2t^3=4t^3 e L=∫014t3 dt=1L=\int_0^14t^3\,dt=1.

Teorema (D, comportamento rispetto al cambio di parametro). Sia r:[a,b]→Ω\mathbf{r}:[a,b]\to\Omega una parametrizzazione regolare e φ:[c,d]→[a,b]\varphi:[c,d]\to[a,b] di classe C1C^1 con φ′≠0\varphi'\ne0, biunivoca; sia r~=r∘φ\tilde{\mathbf{r}}=\mathbf{r}\circ\varphi. Se φ′>0\varphi'>0 allora ∫γ~ω=∫γω\int_{\tilde\gamma}\omega=\int_\gamma\omega; se φ′<0\varphi'<0 allora ∫γ~ω=−∫γω\int_{\tilde\gamma}\omega=-\int_\gamma\omega.

Campi conservativi, irrotazionali e insiemi semplicemente connessi

Definizione (campo conservativo, forma esatta). Sia Ω⊆Rn\Omega\subseteq\mathbb{R}^n aperto e F∈C1(Ω;Rn)\mathbf{F}\in C^1(\Omega;\mathbb{R}^n). F\mathbf{F} si dice conservativo in Ω\Omega se esiste una funzione scalare U:Ω→RU:\Omega\to\mathbb{R} tale che ∇U=Fin Ω,cioeˋ  ∂U∂xi=Fi  ∀i.\nabla U=\mathbf{F}\quad\text{in }\Omega,\qquad\text{cioè }\ \frac{\partial U}{\partial x_i}=F_i\ \ \forall i . UU si chiama potenziale di F\mathbf{F}. Nel linguaggio delle forme, ω=F⋅dr\omega=\mathbf{F}\cdot d\mathbf{r} si dice esatta e UU è una primitiva di ω\omega: dU=ωdU=\omega. Poiché F=∇U\mathbf{F}=\nabla U è C1C^1, UU è di classe C2C^2.

Esempio. U=x2yU=x^2y ha ∇U=(2xy, x2)\nabla U=(2xy,\ x^2), quindi F=(2xy,x2)\mathbf{F}=(2xy,x^2) è conservativo in R2\mathbb{R}^2. Anche U=xyU=xy per (y,x)(y,x), e U=12(x2+y2)U=\tfrac12(x^2+y^2) per il campo radiale (x,y)(x,y).

Teorema (D, teorema fondamentale del calcolo per i campi). Se F=∇U\mathbf{F}=\nabla U in Ω\Omega e γ⊂Ω\gamma\subset\Omega è una curva regolare a tratti con punto iniziale P1=r(a)P_1=\mathbf{r}(a) e punto finale P2=r(b)P_2=\mathbf{r}(b), allora ∫γF⋅dr=U(P2)−U(P1).\int_\gamma\mathbf{F}\cdot d\mathbf{r}=U(P_2)-U(P_1).

Esempio. U=x2yU=x^2y, da (0,0)(0,0) a (1,2)(1,2): lavoro =U(1,2)−U(0,0)=2=U(1,2)-U(0,0)=2, qualunque sia il cammino. Controllo sul segmento (t,2t)(t,2t): F=(4t2,t2)\mathbf{F}=(4t^2,t^2), r′=(1,2)\mathbf{r}'=(1,2), integrando 6t26t^2, integrale 22.

Teorema (D, caratterizzazione dei campi conservativi). Sia Ω\Omega aperto e connesso, F∈C1(Ω;Rn)\mathbf{F}\in C^1(\Omega;\mathbb{R}^n). Sono equivalenti: (a) per ogni coppia di curve regolari a tratti γ1,γ2⊂Ω\gamma_1,\gamma_2\subset\Omega con gli stessi estremi, ∫γ1ω=∫γ2ω\int_{\gamma_1}\omega=\int_{\gamma_2}\omega (indipendenza dal cammino); (b) ∮γω=0\oint_\gamma\omega=0 per ogni curva chiusa regolare a tratti γ⊂Ω\gamma\subset\Omega; (c) F\mathbf{F} è conservativo in Ω\Omega.

Definizione (rotore, campo irrotazionale). Per F=(F1,F2,F3)∈C1(Ω;R3)\mathbf{F}=(F_1,F_2,F_3)\in C^1(\Omega;\mathbb{R}^3) il rotore è rot⁡F=∇∧F=(∂F3∂y−∂F2∂z, ∂F1∂z−∂F3∂x, ∂F2∂x−∂F1∂y),\operatorname{rot}\mathbf{F}=\nabla\wedge\mathbf{F}=\left(\frac{\partial F_3}{\partial y}-\frac{\partial F_2}{\partial z},\ \frac{\partial F_1}{\partial z}-\frac{\partial F_3}{\partial x},\ \frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right), cioè il determinante formale con righe (e1,e2,e3)(\mathbf{e}_1,\mathbf{e}_2,\mathbf{e}_3), (∂x,∂y,∂z)(\partial_x,\partial_y,\partial_z), (F1,F2,F3)(F_1,F_2,F_3) (Prodotto vettoriale e prodotto mistoIl prodotto vettoriale v × w di due vettori di R³ è ortogonale a entrambi, ha norma uguale all'area del parallelogramma che formano e verso dato dalla regola della mano destra; si calcola con un determinante simbolico. Il prodotto misto u · (v × w) = det(u, v, w) è il volume (con segno) del parallelepipedo. Applicazioni: normale a un piano per tre punti, direzione di una retta intersezione di due piani, aree, volumi, distanza tra rette sghembe.Prodotto vettoriale e prodotto misto →). F\mathbf{F} è irrotazionale se rot⁡F=0\operatorname{rot}\mathbf{F}=\mathbf{0} in Ω\Omega; la forma ω\omega è allora chiusa. Nel piano F=(F1,F2)\mathbf{F}=(F_1,F_2) si vede F\mathbf{F} come (F1,F2,0)(F_1,F_2,0) indipendente da zz e rimane solo la terza componente: rot⁡F=∂F2∂x−∂F1∂y(lo scalare rot⁡F⋅e3).\operatorname{rot}\mathbf{F}=\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\quad(\text{lo scalare }\operatorname{rot}\mathbf{F}\cdot\mathbf{e}_3).

Esempio. F=(y+z, x+z, x+y)\mathbf{F}=(y+z,\ x+z,\ x+y): le derivate incrociate sono tutte uguali a 11, quindi rot⁡F=0\operatorname{rot}\mathbf{F}=\mathbf{0} (verificato con sympy).

Proposizione (D, condizione necessaria). Se F∈C1\mathbf{F}\in C^1 è conservativo, allora è irrotazionale (ogni forma esatta è chiusa).

Definizione. Un aperto Ω⊆Rn\Omega\subseteq\mathbb{R}^n è semplicemente connesso se è connesso (per archi) e ogni curva chiusa semplice contenuta in Ω\Omega si può deformare con continuità, restando in Ω\Omega, fino a ridursi a un punto. In pratica nel piano: connesso e senza buchi.

Grafico interattivo: Un disco è semplicemente connesso; la corona circolare no: la circonferenza di raggio 1,5 non si contrae a un punto senza attraversare il buco

Teorema (irrotazionale + semplicemente connesso ⇒\Rightarrow conservativo). Se Ω⊆Rn\Omega\subseteq\mathbb{R}^n è aperto semplicemente connesso e F∈C1(Ω;Rn)\mathbf{F}\in C^1(\Omega;\mathbb{R}^n) è irrotazionale, allora F\mathbf{F} è conservativo in Ω\Omega (ogni forma chiusa è esatta). (Dimostrazione non richiesta al corso; nel piano si deduce dal teorema di Stokes, v. Formule di Gauss-Green e teorema di Stokes nel pianoPer un dominio D del piano con bordo orientato positivamente (D "a sinistra"), ∬ f_y = −∮ f dx e ∬ f_x = ∮ f dy. Se ne ricavano: l'area come ∮ x dy = −∮ y dx = ½∮(x dy − y dx); la formula di Stokes nel piano ∬ (F₂ₓ − F₁ᵧ) dxdy = ∮ F₁dx + F₂dy (flusso del rotore = circuitazione); l'invarianza della circuitazione di un campo irrotazionale tra curve che circondano lo stesso buco.Formule di Gauss-Green e teorema di Stokes nel piano →.)

Corollario (conservatività locale). Se F\mathbf{F} è C1C^1 e irrotazionale in un aperto qualsiasi, ogni punto ha un intorno sferico BB in cui F\mathbf{F} ha un potenziale: la palla è convessa, quindi semplicemente connessa.

10. Formule di Gauss-Green nel piano

Formule di Gauss-Green e teorema di Stokes nel piano

Definizione (orientamento positivo). Si dice che ∂D\partial D è orientato positivamente, e si scrive ∂+D\partial^+D, se è percorso in modo che il dominio resti alla sinistra di chi cammina. Equivalentemente: il versore normale n\mathbf{n} ottenuto ruotando il versore tangente T\mathbf{T} di π/2\pi/2 in senso orario punta verso l'esterno di DD. Se la curva ha parametrizzazione (x(t),y(t))(x(t),y(t)), la normale è n=(y′(t), −x′(t))x′(t)2+y′(t)2.\mathbf{n}=\frac{(y'(t),\,-x'(t))}{\sqrt{x'(t)^2+y'(t)^2}} .

Esempio. Sulla circonferenza antioraria (cos⁡t,sin⁡t)(\cos t,\sin t): T=(−sin⁡t,cos⁡t)\mathbf{T}=(-\sin t,\cos t) e n=(cos⁡t,sin⁡t)\mathbf{n}=(\cos t,\sin t), che punta all'esterno del disco. Dunque per un dominio con un solo bordo (curva chiusa semplice) l'orientamento positivo è quello antiorario.

Grafico interattivo: Corona circolare: bordo esterno antiorario, bordo interno orario (la corona resta a sinistra)

Teorema (D, formule di Gauss-Green). Sia DD un dominio del piano con bordo regolare a tratti orientato positivamente e f∈C1(D)f\in C^1(D) (cioè C1C^1 in un aperto che contiene DD).

  • Se DD è yy-semplice, cioè D={a≤x≤b, φ1(x)≤y≤φ2(x)}D=\{a\le x\le b,\ \varphi_1(x)\le y\le\varphi_2(x)\} con φ1,φ2∈C1([a,b])\varphi_1,\varphi_2\in C^1([a,b]), allora ∬D∂f∂y dx dy=−∮∂+Df dx.(1)\iint_D\frac{\partial f}{\partial y}\,dx\,dy=-\oint_{\partial^+D}f\,dx. \tag{1}
  • Se DD è xx-semplice, cioè D={c≤y≤d, ψ1(y)≤x≤ψ2(y)}D=\{c\le y\le d,\ \psi_1(y)\le x\le\psi_2(y)\}, allora ∬D∂f∂x dx dy=∮∂+Df dy.(2)\iint_D\frac{\partial f}{\partial x}\,dx\,dy=\oint_{\partial^+D}f\,dy. \tag{2}

Esempio. D=[0,1]2D=[0,1]^2, f=y2f=y^2: ∬D2y dx dy=1\iint_D2y\,dx\,dy=1. Il bordo antiorario è formato dai quattro lati; ∮f dx\oint f\,dx è non nulla solo sul lato superiore (y=1y=1, percorso da x=1x=1 a x=0x=0, quindi ∫101 dx=−1\int_1^0 1\,dx=-1) e il lato inferiore (y=0y=0, dà 00). Quindi −∮f dx=−(−1)=1-\oint f\,dx=-(-1)=1.

Corollario (area con un integrale di linea). Se DD ha bordo orientato positivamente, ∣D∣=∬Ddx dy=−∮∂+Dy dx=∮∂+Dx dy=12∮∂+D(x dy−y dx).|D|=\iint_D dx\,dy=-\oint_{\partial^+D}y\,dx=\oint_{\partial^+D}x\,dy=\frac12\oint_{\partial^+D}(x\,dy-y\,dx).

Esempio. Triangolo (0,0),(2,0),(1,1)(0,0),(2,0),(1,1): 12[(0⋅0−2⋅0)+(2⋅1−1⋅0)+(1⋅0−0⋅1)]=1\tfrac12\bigl[(0\cdot0-2\cdot0)+(2\cdot1-1\cdot0)+(1\cdot0-0\cdot1)\bigr]=1 (base 22, altezza 11, area 11).

Grafico interattivo: Regione D sotto un arco di cicloide: |D| = 3π

Teorema (D, formula di Stokes nel piano, o del rotore nel piano). Sia DD un dominio con bordo regolare a pezzi orientato positivamente e F=(F1,F2)∈C1\mathbf{F}=(F_1,F_2)\in C^1 in un aperto che contiene DD. Allora ∬D(∂F2∂x−∂F1∂y)dx dy=∮∂+DF1 dx+F2 dy.\iint_D\left(\frac{\partial F_2}{\partial x}-\frac{\partial F_1}{\partial y}\right)dx\,dy=\oint_{\partial^+D}F_1\,dx+F_2\,dy . Il membro di sinistra è ∬Drot⁡F⋅e3 dx dy\iint_D\operatorname{rot}\mathbf{F}\cdot\mathbf{e}_3\,dx\,dy, quello di destra è la circuitazione ∮∂+DF⋅T ds\oint_{\partial^+D}\mathbf{F}\cdot\mathbf{T}\,ds.

Esempio. F=(−y,x)\mathbf{F}=(-y,x) sul disco unitario: ∂xF2−∂yF1=1−(−1)=2\partial_xF_2-\partial_yF_1=1-(-1)=2, quindi la circuitazione antioraria è 2⋅area=2π2\cdot\text{area}=2\pi. Calcolando direttamente sulla circonferenza (cos⁡t,sin⁡t)(\cos t,\sin t) si ha F⋅r′=sin⁡2t+cos⁡2t=1\mathbf{F}\cdot\mathbf{r}'=\sin^2t+\cos^2t=1 e ∮=∫02π1 dt=2π\oint=\int_0^{2\pi}1\,dt=2\pi: lo stesso valore.

Proposizione (D, invarianza della circuitazione). Siano γ1,γ2\gamma_1,\gamma_2 due curve chiuse semplici regolari a tratti, con γ2\gamma_2 contenuta nella regione interna di γ1\gamma_1, entrambe percorse in senso antiorario. Sia F∈C1\mathbf{F}\in C^1 irrotazionale in un aperto che contiene la regione DD compresa tra le due. Allora ∮γ1F⋅dr=∮γ2F⋅dr\oint_{\gamma_1}\mathbf{F}\cdot d\mathbf{r}=\oint_{\gamma_2}\mathbf{F}\cdot d\mathbf{r}.

11. Superfici e integrali di superficie

Superfici parametriche regolari, area e integrali di superficie

Definizione (superficie parametrica regolare semplice). Sia D⊂R2D\subset\mathbb{R}^2 la chiusura di un aperto connesso e r:D→R3\mathbf{r}:D\to\mathbb{R}^3, r(u,v)=(x(u,v),y(u,v),z(u,v))\mathbf{r}(u,v)=(x(u,v),y(u,v),z(u,v)) continua. La coppia (S,r)(S,\mathbf{r}) con S=r(D)S=\mathbf{r}(D) è una superficie regolare semplice se, indicato con D˚\mathring D l'interno di DD: (i) r\mathbf{r} è iniettiva e di classe C1C^1 in D˚\mathring D; (ii) la matrice jacobiana Jr(u,v)J\mathbf{r}(u,v) (matrice 3×23\times2) ha rango massimo 22 in ogni punto di D˚\mathring D. SS è il sostegno, r\mathbf{r} la parametrizzazione.

Definizione (piano tangente, normale). Nel punto x0=r(u0,v0)\mathbf{x}_0=\mathbf{r}(u_0,v_0), con (u0,v0)(u_0,v_0) interno, il piano tangente è il piano per x0\mathbf{x}_0 generato da ru\mathbf{r}_u e rv\mathbf{r}_v. Il vettore N=ru∧rv\mathbf{N}=\mathbf{r}_u\wedge\mathbf{r}_v è il vettore normale e n=N/∥N∥\mathbf{n}=\mathbf{N}/\|\mathbf{N}\| il versore normale. L'equazione del piano tangente è (ru∧rv)(u0,v0)⋅(x−x0)=0,(\mathbf{r}_u\wedge\mathbf{r}_v)(u_0,v_0)\cdot(\mathbf{x}-\mathbf{x}_0)=0 , cioè a0(x−x0)+b0(y−y0)+c0(z−z0)=0a_0(x-x_0)+b_0(y-y_0)+c_0(z-z_0)=0 con (a0,b0,c0)=N(u0,v0)(a_0,b_0,c_0)=\mathbf{N}(u_0,v_0).

Definizione (area, integrale di superficie). Sia (S,r)(S,\mathbf{r}) una superficie regolare con DD limitato e misurabile. L'area è A(S)=∬D∥ru∧rv∥ du dv=∬D(∂(y,z)∂(u,v))2+(∂(z,x)∂(u,v))2+(∂(x,y)∂(u,v))2 du dv.A(S)=\iint_D\|\mathbf{r}_u\wedge\mathbf{r}_v\|\,du\,dv=\iint_D\sqrt{\left(\tfrac{\partial(y,z)}{\partial(u,v)}\right)^2+\left(\tfrac{\partial(z,x)}{\partial(u,v)}\right)^2+\left(\tfrac{\partial(x,y)}{\partial(u,v)}\right)^2}\,du\,dv . Per f:S→Rf:S\to\mathbb{R} continua l'integrale di ff su SS (integrale di superficie di prima specie) è ∬Sf dσ=∬Df(r(u,v)) ∥ru∧rv∥ du dv,\iint_Sf\,d\sigma=\iint_Df(\mathbf{r}(u,v))\,\|\mathbf{r}_u\wedge\mathbf{r}_v\|\,du\,dv , dove dσ=∥ru∧rv∥ du dvd\sigma=\|\mathbf{r}_u\wedge\mathbf{r}_v\|\,du\,dv è l'elemento d'area; in particolare A(S)=∬SdσA(S)=\iint_Sd\sigma.

Esempio. Sfera, porzione compresa fra i paralleli φ1<φ<φ2\varphi_1<\varphi<\varphi_2 e i meridiani θ1<θ<θ2\theta_1<\theta<\theta_2: A=∫θ1θ2dθ∫φ1φ2R2sin⁡φ dφ=R2(θ2−θ1)(cos⁡φ1−cos⁡φ2).A=\int_{\theta_1}^{\theta_2}d\theta\int_{\varphi_1}^{\varphi_2}R^2\sin\varphi\,d\varphi=R^2(\theta_2-\theta_1)(\cos\varphi_1-\cos\varphi_2). Per φ∈[0,π]\varphi\in[0,\pi], θ∈[0,2π]\theta\in[0,2\pi] si ritrova 4πR24\pi R^2; una calotta attorno al polo nord con colatitudine massima φˉ\bar\varphi ha area 2πR2(1−cos⁡φˉ)2\pi R^2(1-\cos\bar\varphi).

Teorema (D, secondo teorema di Pappo-Guldino). L'area della superficie di rotazione generata da γ\gamma è la lunghezza di γ\gamma per la lunghezza dell'arco di circonferenza percorso dal baricentro (centroide) di γ\gamma: A=θˉ∫γx ds=θˉ xB ℓ(γ),A=\bar\theta\int_\gamma x\,ds=\bar\theta\,x_B\,\ell(\gamma), dove xBx_B è l'ascissa del baricentro della curva e ℓ(γ)\ell(\gamma) la sua lunghezza. In particolare, per una rotazione completa, A=2π∫γx dsA=2\pi\int_\gamma x\,ds.

12. Flusso, divergenza e rotore

Flusso, teorema della divergenza e teorema del rotore

Definizione (flusso). Sia (S,r)(S,\mathbf{r}) una superficie regolare con versore normale n=ru∧rv∥ru∧rv∥\mathbf{n}=\dfrac{\mathbf{r}_u\wedge\mathbf{r}_v}{\|\mathbf{r}_u\wedge\mathbf{r}_v\|} e F:S→R3\mathbf{F}:S\to\mathbb{R}^3 un campo continuo. Il flusso di F\mathbf{F} attraverso SS (nel verso di n\mathbf{n}) è Φ=∬SF⋅n dσ=∬DF(r(u,v))⋅(ru∧rv)(u,v) du dv.\Phi=\iint_S\mathbf{F}\cdot\mathbf{n}\,d\sigma=\iint_D\mathbf{F}(\mathbf{r}(u,v))\cdot(\mathbf{r}_u\wedge\mathbf{r}_v)(u,v)\,du\,dv . Il secondo membro segue dal primo perché n dσ=N∥N∥∥N∥ du dv=N du dv\mathbf{n}\,d\sigma=\dfrac{\mathbf{N}}{\|\mathbf{N}\|}\|\mathbf{N}\|\,du\,dv=\mathbf{N}\,du\,dv: il fattore ∥N∥\|\mathbf{N}\| si semplifica e non servono radici.

Esempio. F=(0,0,z)\mathbf{F}=(0,0,z) attraverso il grafico z=1−x2−y2z=1-x^2-y^2 su x2+y2≤1x^2+y^2\le1, normale verso l'alto: F⋅N=1−x2−y2\mathbf{F}\cdot\mathbf{N}=1-x^2-y^2 e Φ=∬D(1−ρ2)ρ dρ dθ=π2\Phi=\iint_D(1-\rho^2)\rho\,d\rho\,d\theta=\tfrac\pi2.

Definizione. Un dominio regolare di R3\mathbb{R}^3 è la chiusura Ω\Omega di un aperto limitato la cui frontiera ∂Ω\partial\Omega è il sostegno di una superficie regolare (o regolare a pezzi: per esempio i solidi delimitati da un numero finito di grafici, il cubo, il cilindro con le basi). Una superficie con sostegno ∂Ω\partial\Omega si dice chiusa. La normale è uscente (si scrive ne\mathbf{n}_e) se P+tn(P)∉ΩP+t\mathbf{n}(P)\notin\Omega per t>0t>0 piccolo, entrante se P+tn(P)∈Ω˚P+t\mathbf{n}(P)\in\mathring\Omega.

Esempio (campo elettrostatico). E=kqr∥r∥3\mathbf{E}=kq\dfrac{\mathbf{r}}{\|\mathbf{r}\|^3} sulla sfera di raggio RR: ne=r/R\mathbf{n}_e=\mathbf{r}/R e E⋅ne=kqR2\mathbf{E}\cdot\mathbf{n}_e=\dfrac{kq}{R^2} (costante). Quindi Φ=kqR2⋅4πR2=4πkq\Phi=\dfrac{kq}{R^2}\cdot4\pi R^2=4\pi kq, indipendente da RR e proporzionale alla carica: è il teorema di Gauss dell'elettrostatica.

Definizione (divergenza). Per F=(F1,F2,F3)∈C1\mathbf{F}=(F_1,F_2,F_3)\in C^1, div⁡F=∇⋅F=∂F1∂x+∂F2∂y+∂F3∂z(nel piano: ∂xF1+∂yF2).\operatorname{div}\mathbf{F}=\nabla\cdot\mathbf{F}=\frac{\partial F_1}{\partial x}+\frac{\partial F_2}{\partial y}+\frac{\partial F_3}{\partial z}\quad(\text{nel piano: }\partial_xF_1+\partial_yF_2).

Esempio. div⁡(x,y,z2)=1+1+2z\operatorname{div}(x,y,z^2)=1+1+2z. div⁡(x,0,y)=1\operatorname{div}(x,0,y)=1.

Teorema (della divergenza, o di Gauss; D per domini semplici rispetto a tutti gli assi). Sia Ω⊂R3\Omega\subset\mathbb{R}^3 un dominio regolare (o con frontiera regolare a pezzi), ne\mathbf{n}_e la normale uscente su ∂Ω\partial\Omega e F∈C1\mathbf{F}\in C^1 in un aperto che contiene Ω\Omega. Allora ∭Ωdiv⁡F dx dy dz=∬∂ΩF⋅ne dσ.\iiint_\Omega\operatorname{div}\mathbf{F}\,dx\,dy\,dz=\iint_{\partial\Omega}\mathbf{F}\cdot\mathbf{n}_e\,d\sigma . Se si usa la normale entrante ni=−ne\mathbf{n}_i=-\mathbf{n}_e vale ∭Ωdiv⁡F=−∬∂ΩF⋅ni dσ\iiint_\Omega\operatorname{div}\mathbf{F}=-\iint_{\partial\Omega}\mathbf{F}\cdot\mathbf{n}_i\,d\sigma.

Esempio (la sfera di prima). F=(x,0,y)\mathbf{F}=(x,0,y), div⁡F=1\operatorname{div}\mathbf{F}=1: Φ=vol⁡(BR)=43πR3\Phi=\operatorname{vol}(B_R)=\tfrac43\pi R^3, in un passaggio.

Teorema (D, divergenza nel piano). Sia D⊂R2D\subset\mathbb{R}^2 un dominio con frontiera regolare a pezzi, ne\mathbf{n}_e la normale esterna a ∂D\partial D, F∈C1\mathbf{F}\in C^1. Allora ∬Ddiv⁡F dx dy=∮∂DF⋅ne ds\displaystyle\iint_D\operatorname{div}\mathbf{F}\,dx\,dy=\oint_{\partial D}\mathbf{F}\cdot\mathbf{n}_e\,ds.

Esempio. F=(xy,x2)\mathbf{F}=(xy,x^2) uscente dal bordo del quadrato [−12,12]2[-\tfrac12,\tfrac12]^2: div⁡F=y\operatorname{div}\mathbf{F}=y, e ∬y dx dy=0\iint y\,dx\,dy=0 per simmetria; il calcolo lato per lato dà lo stesso.

Definizione (superficie con bordo). Sia D⊂R2D\subset\mathbb{R}^2 la chiusura di un aperto connesso, A⊃DA\supset D un aperto e r:A→R3\mathbf{r}:A\to\mathbb{R}^3 di classe C1C^1, iniettiva in DD con JrJ\mathbf{r} di rango 22 in ogni punto di DD. S=r(D)S=\mathbf{r}(D) è una superficie regolare con bordo, e r(∂D)\mathbf{r}(\partial D) è il bordo. Esempi: un grafico z=f(x,y)z=f(x,y) su DD ha come bordo il grafico di ff ristretto a ∂D\partial D; la semisfera ha come bordo l'equatore; una sfera intera non è una superficie con bordo.

Teorema (del rotore, o di Stokes). Sia (S,r)(S,\mathbf{r}) una superficie regolare con bordo e F∈C1\mathbf{F}\in C^1 in un aperto che contiene SS. Allora ∬Srot⁡F⋅n dσ=∮∂+SF1 dx+F2 dy+F3 dz=∮∂+SF⋅T ds.\iint_S\operatorname{rot}\mathbf{F}\cdot\mathbf{n}\,d\sigma=\oint_{\partial^+S}F_1\,dx+F_2\,dy+F_3\,dz=\oint_{\partial^+S}\mathbf{F}\cdot\mathbf{T}\,ds . Cioè: il flusso del rotore attraverso SS è uguale alla circuitazione di F\mathbf{F} lungo il bordo orientato positivamente.

Esempio (campo V=(−y,x,yez2)V=(-y,x,ye^{z^2}) sulla semisfera). Sia SS la semisfera unitaria z≥0z\ge0 con normale uscente (verso l'alto). Il bordo orientato positivamente è l'equatore (cos⁡t,sin⁡t,0)(\cos t,\sin t,0) antiorario, su cui V=(−sin⁡t,cos⁡t,sin⁡t)V=(-\sin t,\cos t,\sin t) e r′=(−sin⁡t,cos⁡t,0)\mathbf{r}'=(-\sin t,\cos t,0), V⋅r′=1V\cdot\mathbf{r}'=1. Quindi ∬Srot⁡V⋅n dσ=∮V⋅dr=2π\iint_S\operatorname{rot}V\cdot\mathbf{n}\,d\sigma=\oint V\cdot d\mathbf{r}=2\pi. Calcolando direttamente, rot⁡V=(ez2,0,2)\operatorname{rot}V=(e^{z^2},0,2), rot⁡V⋅n=xez2+2z\operatorname{rot}V\cdot\mathbf{n}=xe^{z^2}+2z e ∬S2z dσ=2π\iint_S2z\,d\sigma=2\pi (il termine in xx ha integrale nullo per simmetria): stesso risultato.

Teorema (D, rotore per superfici chiuse). Se S=∂ΩS=\partial\Omega è una superficie chiusa regolare a pezzi e F∈C2\mathbf{F}\in C^2 in un aperto che contiene Ω\Omega, allora ∬Srot⁡F⋅n dσ=0\iint_S\operatorname{rot}\mathbf{F}\cdot\mathbf{n}\,d\sigma=0: il flusso del rotore attraverso una superficie chiusa è nullo.

Definizione. Un campo G∈C1\mathbf{G}\in C^1 in un aperto E⊆R3E\subseteq\mathbb{R}^3 si dice solenoidale se div⁡G=0\operatorname{div}\mathbf{G}=0 in EE. Un campo F\mathbf{F} con G=rot⁡F\mathbf{G}=\operatorname{rot}\mathbf{F} si dice potenziale vettore di G\mathbf{G}.

13. Equazioni differenziali ordinarie: il problema di Cauchy

Problema di Cauchy per sistemi di equazioni differenziali

Definizione (sistema in forma normale). Sia D⊆R×RnD\subseteq\mathbb{R}\times\mathbb{R}^n aperto e f:D→Rn\mathbf{f}:D\to\mathbb{R}^n un campo. Il sistema di nn equazioni del primo ordine in forma normale è y′(t)=f(t,y(t)),cioeˋ{y1′=f1(t,y1,…,yn) ⋮yn′=fn(t,y1,…,yn)\mathbf{y}'(t)=\mathbf{f}(t,\mathbf{y}(t)),\qquad\text{cioè}\qquad\begin{cases}y_1'=f_1(t,y_1,\dots,y_n)\\\ \vdots\\y_n'=f_n(t,y_1,\dots,y_n)\end{cases} Una soluzione è una coppia (I,y)(I,\mathbf{y}): I⊆RI\subseteq\mathbb{R} intervallo, y:I→Rn\mathbf{y}:I\to\mathbb{R}^n derivabile, con (t,y(t))∈D(t,\mathbf{y}(t))\in D e y′(t)=f(t,y(t))\mathbf{y}'(t)=\mathbf{f}(t,\mathbf{y}(t)) per ogni t∈It\in I. L'integrale generale è l'insieme di tutte le soluzioni.

Definizione (problema di Cauchy). Fissato (t0,y0)∈D(t_0,\mathbf{y}_0)\in D, il problema di Cauchy (o ai valori iniziali) è {y′(t)=f(t,y(t))y(t0)=y0.\begin{cases}\mathbf{y}'(t)=\mathbf{f}(t,\mathbf{y}(t))\\ \mathbf{y}(t_0)=\mathbf{y}_0 .\end{cases} Una soluzione è una soluzione (I,y)(I,\mathbf{y}) del sistema con t0∈It_0\in I e y(t0)=y0\mathbf{y}(t_0)=\mathbf{y}_0.

Esempio. y′=yy'=y, y(0)=3y(0)=3: la soluzione è y=3ety=3e^t su R\mathbb{R} (con la costante c=3c=3 scelta dal dato iniziale). Per n=1n=1 il sistema è un'equazione scalare, e si ritrova l'impostazione di Analisi 1.

Teorema (di Peano, esistenza locale). Se f:D→Rn\mathbf{f}:D\to\mathbb{R}^n è continua e (t0,y0)∈D(t_0,\mathbf{y}_0)\in D, esiste δ>0\delta>0 e almeno una soluzione del problema di Cauchy definita su (t0−δ,t0+δ)(t_0-\delta,t_0+\delta).

Teorema (di Cauchy, esistenza e unicità locale). Sia D⊆R×RnD\subseteq\mathbb{R}\times\mathbb{R}^n aperto, (t0,y0)∈D(t_0,\mathbf{y}_0)\in D, f:D→Rn\mathbf{f}:D\to\mathbb{R}^n continua e tale che tutte le derivate parziali ∂fi∂yk\dfrac{\partial f_i}{\partial y_k} (i,k=1,…,ni,k=1,\dots,n) esistano e siano continue in DD (cioè f\mathbf{f} è C1C^1 nelle variabili di stato). Allora esiste δ>0\delta>0 tale che il problema di Cauchy ha una e una sola soluzione definita su (t0−δ,t0+δ)(t_0-\delta,t_0+\delta).

Teorema (esistenza e unicità globale). Sia S=[a,b]×RnS=[a,b]\times\mathbb{R}^n e f:S→Rn\mathbf{f}:S\to\mathbb{R}^n continua, con le derivate ∂fi/∂yk\partial f_i/\partial y_k continue (o localmente lipschitziana in y\mathbf{y}). Se esistono H,K≥0H,K\ge0 tali che ∥f(t,y)∥≤H+K∥y∥∀(t,y)∈S,\|\mathbf{f}(t,\mathbf{y})\|\le H+K\|\mathbf{y}\|\qquad\forall(t,\mathbf{y})\in S , allora per ogni (t0,y0)∈S(t_0,\mathbf{y}_0)\in S il problema di Cauchy ha una e una sola soluzione definita in tutto [a,b][a,b].

14. Sistemi lineari e stabilità degli equilibri

Sistemi lineari ed equazioni del secondo ordine a coefficienti costanti

Definizione. Un sistema lineare del primo ordine è y′(t)=A(t)y(t)+g(t),(S)\mathbf{y}'(t)=A(t)\mathbf{y}(t)+\mathbf{g}(t),\tag{S} con A(t)=(aij(t))A(t)=(a_{ij}(t)) matrice n×nn\times n a coefficienti continui su un intervallo II (anche illimitato) e g∈C(I;Rn)\mathbf{g}\in C(I;\mathbb{R}^n) termine noto. Si dice omogeneo se g≡0\mathbf{g}\equiv\mathbf{0}.

Teorema (D, risolubilità del problema di Cauchy). Per ogni t0∈It_0\in I e y0∈Rn\mathbf{y}_0\in\mathbb{R}^n il problema y′=Ay+g\mathbf{y}'=A\mathbf{y}+\mathbf{g}, y(t0)=y0\mathbf{y}(t_0)=\mathbf{y}_0 ha una e una sola soluzione, definita in tutto II.

Corollario (D, equazioni scalari del secondo ordine). Se b,c,d∈C(I)b,c,d\in C(I), per ogni t0∈It_0\in I e y0,v0∈Ry_0,v_0\in\mathbb{R} il problema y′′+b(t)y′+c(t)y=d(t)y''+b(t)y'+c(t)y=d(t), y(t0)=y0y(t_0)=y_0, y′(t0)=v0y'(t_0)=v_0 ha un'unica soluzione definita su tutto II.

Proposizione (D, principio di sovrapposizione per l'omogeneo). Se y1,…,ym\mathbf{y}_1,\dots,\mathbf{y}_m risolvono y′=Ay\mathbf{y}'=A\mathbf{y}, allora anche ogni combinazione v=∑iαiyi\mathbf{v}=\sum_i\alpha_i\mathbf{y}_i la risolve: le soluzioni formano un sottospazio vettoriale di VV.

Teorema (D, struttura dell'integrale generale dell'omogeneo). L'insieme delle soluzioni di y′=A(t)y\mathbf{y}'=A(t)\mathbf{y} è uno spazio vettoriale di dimensione nn: esistono nn soluzioni linearmente indipendenti y1,…,yn\mathbf{y}_1,\dots,\mathbf{y}_n e ogni soluzione è ∑αiyi\sum\alpha_i\mathbf{y}_i con αi∈R\alpha_i\in\mathbb{R}.

Esempio (n=1n=1). y′=a(t)yy'=a(t)y: lo spazio è di dimensione 11 e l'integrale generale è y=kexp⁡(∫t0ta(s) ds)y=k\exp\left(\int_{t_0}^ta(s)\,ds\right).

Proposizione (D, sovrapposizione per il non omogeneo). Se y~\tilde{\mathbf{y}} è una soluzione particolare di y′=Ay+g\mathbf{y}'=A\mathbf{y}+\mathbf{g}, l'integrale generale è y=z+y~\mathbf{y}=\mathbf{z}+\tilde{\mathbf{y}}, con z\mathbf{z} che descrive l'integrale generale dell'omogeneo associato z′=Az\mathbf{z}'=A\mathbf{z}.

Teorema (D, integrale generale). Siano λ1,2\lambda_{1,2} le radici di λ2+bλ+c=0\lambda^2+b\lambda+c=0.

  1. Se b2>4cb^2>4c (radici reali distinte): y=αeλ1t+βeλ2ty=\alpha e^{\lambda_1t}+\beta e^{\lambda_2t}.
  2. Se b2=4cb^2=4c (radice doppia λ=−b/2\lambda=-b/2): y=(α+βt)eλty=(\alpha+\beta t)e^{\lambda t}.
  3. Se b2<4cb^2<4c (radici complesse μ±iω\mu\pm i\omega, μ=−b/2\mu=-b/2, ω=4c−b2/2\omega=\sqrt{4c-b^2}/2): y=eμt(αcos⁡ωt+βsin⁡ωt)y=e^{\mu t}(\alpha\cos\omega t+\beta\sin\omega t).

Esempio. y′′−y′−12y=0y''-y'-12y=0: radici 1±72=4,−3\frac{1\pm7}2=4,-3, y=αe4t+βe−3ty=\alpha e^{4t}+\beta e^{-3t}. y′′+4y′+5y=0y''+4y'+5y=0: radici −2±i-2\pm i, y=e−2t(αcos⁡t+βsin⁡t)y=e^{-2t}(\alpha\cos t+\beta\sin t).

Stabilità degli equilibri dei sistemi autonomi 2x2

Definizione. Un sistema si dice autonomo se il campo non dipende dal tempo: y′=f(y)\mathbf{y}'=\mathbf{f}(\mathbf{y}). In dimensione 22: {x′=f(x,y)y′=g(x,y)f,g∈C1(D), D⊆R2 aperto.\begin{cases}x'=f(x,y)\\ y'=g(x,y)\end{cases}\qquad f,g\in C^1(D),\ D\subseteq\mathbb{R}^2\text{ aperto}.

Proposizione (D, invarianza per traslazioni temporali). Se r(t)\mathbf{r}(t) risolve y′=f(y)\mathbf{y}'=\mathbf{f}(\mathbf{y}) su (a,b)(a,b), allora per ogni τ∈R\tau\in\mathbb{R} anche s(t)=r(t+τ)\mathbf{s}(t)=\mathbf{r}(t+\tau) è soluzione, su (a−τ,b−τ)(a-\tau,b-\tau); ha la stessa orbita.

Proposizione (D, le orbite non si intersecano). Se due orbite hanno un punto in comune, coincidono.

Definizione. PP è un punto di equilibrio (o stato stazionario) se f(P)=0\mathbf{f}(P)=\mathbf{0}.

Esempio. Per {x′=−2x−23yy′=3x+y\begin{cases}x'=-2x-\frac23y\\y'=3x+y\end{cases} gli equilibri risolvono 3x+y=03x+y=0 (la seconda equazione è anche la prima, a meno di un fattore): un'intera retta.

Definizione (stabilità secondo Lyapunov).

  • PP è stabile se ∀ε>0 ∃δ>0\forall\varepsilon>0\ \exists\delta>0 tale che ∥y0−P∥<δ\|\mathbf{y}_0-P\|<\delta implica che φ(t,y0)\boldsymbol\varphi(t,\mathbf{y}_0) esiste per ogni t≥0t\ge0 e ∥φ(t,y0)−P∥<ε\|\boldsymbol\varphi(t,\mathbf{y}_0)-P\|<\varepsilon per ogni t≥0t\ge0 (le soluzioni che partono vicino restano vicine per sempre).
  • PP è asintoticamente stabile (attrattore) se è stabile e inoltre esiste δ1>0\delta_1>0 tale che ∥y0−P∥<δ1⇒lim⁡t→+∞φ(t,y0)=P\|\mathbf{y}_0-P\|<\delta_1\Rightarrow\lim_{t\to+\infty}\boldsymbol\varphi(t,\mathbf{y}_0)=P. Se vale per ogni dato iniziale è un attrattore globale.
  • PP è instabile se non è stabile.

Teorema (classificazione, caso det⁡A≠0\det A\ne0). Siano λ1,λ2∈C\lambda_1,\lambda_2\in\mathbb{C} gli autovalori di AA. L'origine è

  1. asintoticamente stabile se e solo se Re⁡λ1<0\operatorname{Re}\lambda_1<0 e Re⁡λ2<0\operatorname{Re}\lambda_2<0;
  2. stabile (non asintoticamente) se e solo se Re⁡λ1,Re⁡λ2≤0\operatorname{Re}\lambda_1,\operatorname{Re}\lambda_2\le0 e almeno uno vale 00 (in tal caso gli autovalori sono ±iω\pm i\omega);
  3. instabile se e solo se almeno un autovalore ha Re⁡λ>0\operatorname{Re}\lambda>0.

Proposizione (D). Per AA reale 2×22\times2: l'origine è asintoticamente stabile se e solo se τ<0\tau<0 e d>0d>0.

Esempio. A=(−2α31)A=\begin{pmatrix}-2&\alpha\\3&1\end{pmatrix}: τ=−1<0\tau=-1<0 sempre, d=−2−3α>0  ⟺  α<−23d=-2-3\alpha>0\iff\alpha<-\tfrac23. Per α=−23\alpha=-\tfrac23 è d=0d=0 e c'è una retta di equilibri: non asintoticamente stabile. Quindi l'origine è asintoticamente stabile se e solo se α<−23\alpha<-\tfrac23 (Esercizio - stabilità e problema di Cauchy al variare di un parametro (ricevimento, appello 8 luglio 2025)).