Salta al contenuto
Note per Studenti Esercizio - Trasformazione quadratica, determinante jacobiano e inversa (lezione 19)

Esercizio - Trasformazione quadratica, determinante jacobiano e inversa (lezione 19)

In questa pagina 3

Teoria: Funzioni a valori vettoriali, matrice jacobiana e trasformazioni di coordinateUna funzione f : A ⊆ Rⁿ → Rᵐ ha m componenti scalari f₁,…,f_m; il limite si fa componente per componente. È differenziabile in x₀ se f(x₀+h) − f(x₀) = Df(x₀)h + o(‖h‖), dove la matrice jacobiana Df(x₀) (m righe, n colonne) ha elemento (i,j) = ∂f_i/∂x_j: le righe sono i gradienti delle componenti. Se le componenti sono C¹ la f è differenziabile. Casi particolari: m = 1 (riga = gradiente), n = 1 (colonna = vettore tangente), f affine (Df = A). Regola della catena: D(g∘f)(x₀) = Dg(f(x₀))·Df(x₀). Trasformazioni di coordinate: polari (det = ρ), cilindriche (det = ρ), sferiche (det = ρ² sin φ); |det Df| è il fattore di variazione di aree o volumi.Funzioni a valori vettoriali, matrice jacobiana e trasformazioni di coordinate → e Cambiamento di variabili negli integrali doppi, integrali generalizzati e gaussianaCambiando coordinate (x,y) = T(u,v), con T diffeomorfismo C¹ tra D' e D, ∬D f(x,y)dxdy = ∬{D'} f(T(u,v)) |det DT(u,v)| dudv: il fattore |det DT| è la dilatazione locale delle aree (per le trasformazioni lineari è |det M|, dimostrato col prodotto vettoriale). Polari: x = x₀ + ρcosθ, y = y₀ + ρsinθ, dxdy = ρ dρdθ; area di un settore di corona ½(R₂² − R₁²)Δθ; disco x² + y² ≤ 2ax: ρ ≤ 2a cosθ. Coordinate ellittiche x = aρcosθ, y = bρsinθ, dxdy = abρ dρdθ (area dell'ellisse πab). Simmetrie: D simmetrico rispetto all'asse y, f dispari in x ⇒ ∬ = 0 (D), f pari ⇒ ∬ = 2∬_{x≥0}. Integrali generalizzati di f ≥ 0 continua su R²: limite su una successione invadente, indipendente dalla successione; gaussiana (D): ∬ e^{−(x²+y²)/2} = 2π, quindi ∫_R e^{−x²/2}dx = √(2π).Cambiamento di variabili negli integrali doppi, integrali generalizzati e gaussiana →. Fonte: esercizio nello stile delle lezioni 17 e 19 di Fondamenti di Analisi Matematica 2, Ing. Aerospaziale UniPD.

Testo

Sia T(u,v)=(u2−v2, 2uv)\mathbf{T}(u,v) = (u^2 - v^2,\ 2uv), da R2\mathbb{R}^2 a R2\mathbb{R}^2.

  • (a) Calcolare DTD\mathbf{T} e det⁡DT\det D\mathbf{T}; dire dove det⁡DT=0\det D\mathbf{T} = 0.
  • (b) Calcolare il fattore di dilatazione delle aree in (1,0)(1,0) e in (0,1)(0,1).
  • (c) Mostrare che T\mathbf{T} non è iniettiva su R2\mathbb{R}^2 ma lo è sul semipiano {v>0}\{v > 0\}.
  • (d) Calcolare l'area di T(E)\mathbf{T}(E), dove E={u2+v2≤1, v>0}E = \{u^2 + v^2 \le 1,\ v > 0\}.
  • (e) Calcolare det⁡DT−1\det D\mathbf{T}^{-1} nel punto (x,y)=T(u,v)(x,y) = \mathbf{T}(u,v) con v>0v > 0.

Soluzione

(a) Jacobiana

DT(u,v)=(2u−2v2v2u),det⁡DT=4u2+4v2=4(u2+v2).D\mathbf{T}(u,v) = \begin{pmatrix} 2u & -2v \\ 2v & 2u\end{pmatrix}, \qquad \det D\mathbf{T} = 4u^2 + 4v^2 = 4(u^2 + v^2). Il determinante si annulla solo in (0,0)(0,0).

(b) Fattore di dilatazione delle aree

È ∣det⁡DT∣|\det D\mathbf{T}|: in (1,0)(1,0) vale 44 e in (0,1)(0,1) vale 44. Un quadratino di lato piccolo ε\varepsilon vicino a (1,0)(1,0) viene trasformato in un quadratino di area circa 4ε24\varepsilon^2 (infatti DT(1,0)=2ID\mathbf{T}(1,0) = 2I è una dilatazione di fattore 22 nelle due direzioni, con area moltiplicata per 44).

(c) Iniettività

Si osserva che, con z=u+ivz = u + iv, T\mathbf{T} è il quadrato complesso z2=(u2−v2)+i 2uvz^2 = (u^2 - v^2) + i\,2uv. Se T(u,v)=T(u′,v′)\mathbf{T}(u,v) = \mathbf{T}(u',v') allora z2=z′2z^2 = z'^2, cioè (z−z′)(z+z′)=0(z - z')(z + z') = 0: z′=zz' = z oppure z′=−zz' = -z. Quindi T(u,v)=T(−u,−v)\mathbf{T}(u,v) = \mathbf{T}(-u,-v): per esempio T(1,0)=T(−1,0)=(1,0)\mathbf{T}(1,0) = \mathbf{T}(-1,0) = (1,0), e T\mathbf{T} non è iniettiva su R2\mathbb{R}^2. Sul semipiano {v>0}\{v > 0\} i due punti (u,v)(u,v) e (−u,−v)(-u,-v) non possono essere entrambi presenti (uno ha v>0v > 0, l'altro v<0v < 0): T\mathbf{T} è iniettiva.

(In polari ρ=∣z∣\rho = |z|, θ=arg⁡z\theta = \arg z, la trasformazione manda (ρ,θ)(\rho,\theta) in (ρ2,2θ)(\rho^2, 2\theta): raddoppia gli angoli, perciò un semipiano θ∈(0,π)\theta \in (0,\pi) viene mandato su tutto il piano meno il semiasse positivo delle xx.)

(d) Area dell'immagine

EE è il semidisco superiore aperto. Per il teorema del cambiamento di variabili (con T\mathbf{T} iniettiva, C1C^1 e det⁡DT≠0\det D\mathbf{T} \ne 0 su EE) l'area di T(E)\mathbf{T}(E) è ∣T(E)∣=∬E∣det⁡DT∣ du dv=∬E4(u2+v2) du dv.|\mathbf{T}(E)| = \iint_E |\det D\mathbf{T}|\,du\,dv = \iint_E 4(u^2 + v^2)\,du\,dv. In polari u=ρcos⁡θu = \rho\cos\theta, v=ρsin⁡θv = \rho\sin\theta con ρ∈(0,1)\rho \in (0,1) e θ∈(0,π)\theta \in (0,\pi) (qui il fattore ρ\rho del cambio di variabili si aggiunge): ∣T(E)∣=∫0πdθ∫014ρ2⋅ρ dρ=π⋅4⋅14=π.|\mathbf{T}(E)| = \int_0^\pi d\theta\int_0^1 4\rho^2\cdot\rho\,d\rho = \pi\cdot 4\cdot\frac14 = \pi. È coerente con la geometria: T(E)\mathbf{T}(E) è il disco unitario privato del segmento [0,1]×{0}[0,1]\times\{0\} (che ha area nulla), di area π\pi.

(e) Determinante dell'inversa

Poiché det⁡DT−1(x,y)=1det⁡DT(u,v)\det D\mathbf{T}^{-1}(x,y) = \frac{1}{\det D\mathbf{T}(u,v)} con (x,y)=T(u,v)(x,y) = \mathbf{T}(u,v), bisogna esprimere u2+v2u^2 + v^2 con x,yx,y. Da x2+y2=(u2−v2)2+4u2v2=(u2+v2)2x^2 + y^2 = (u^2 - v^2)^2 + 4u^2v^2 = (u^2 + v^2)^2 si ha u2+v2=x2+y2u^2 + v^2 = \sqrt{x^2 + y^2}. Quindi det⁡DT−1(x,y)=14x2+y2.\det D\mathbf{T}^{-1}(x,y) = \frac{1}{4\sqrt{x^2 + y^2}}.

Errori tipici

  • Dimenticare il valore assoluto del determinante nel fattore di dilatazione delle aree (qui è comunque positivo).
  • Applicare il cambio di variabili senza controllare l'iniettività: su tutto il disco {u2+v2≤1}\{u^2 + v^2 \le 1\} l'integrale di ∣det⁡DT∣|\det D\mathbf{T}| vale 2π2\pi, il doppio dell'area dell'immagine, perché ogni punto dell'immagine viene contato due volte.

Versione ripasso

Testo. T(u,v)=(u2−v2,2uv)\mathbf{T}(u,v) = (u^2 - v^2, 2uv): jacobiana, fattori di dilatazione, iniettività, area di T(E)\mathbf{T}(E) per EE semidisco superiore, det⁡DT−1\det D\mathbf{T}^{-1}.

Soluzione.

  • DT=(2u−2v2v2u)D\mathbf{T} = \begin{pmatrix} 2u & -2v \\ 2v & 2u\end{pmatrix}, det⁡=4(u2+v2)\det = 4(u^2 + v^2), nullo solo in (0,0)(0,0). In (1,0)(1,0) e (0,1)(0,1): fattore 44.
  • T=z2\mathbf{T} = z^2: T(u,v)=T(−u,−v)\mathbf{T}(u,v) = \mathbf{T}(-u,-v), non iniettiva su R2\mathbb{R}^2; iniettiva su v>0v > 0.
  • ∣T(E)∣=∬E4(u2+v2)=∫0π∫014ρ3 dρ dθ=π|\mathbf{T}(E)| = \iint_E 4(u^2 + v^2) = \int_0^\pi\int_0^1 4\rho^3\,d\rho\,d\theta = \pi (disco meno un segmento).
  • x2+y2=(u2+v2)2x^2 + y^2 = (u^2 + v^2)^2 ⇒\Rightarrow det⁡DT−1(x,y)=14x2+y2\det D\mathbf{T}^{-1}(x,y) = \frac{1}{4\sqrt{x^2 + y^2}}.

Errori: cambio di variabili senza iniettività (sul disco intero si otterrebbe 2π2\pi).

Passaggi. T=z2\mathbf{T} = z^2: (z−z′)(z+z′)=0(z - z')(z + z') = 0. In polari T(ρ,θ)=(ρ2,2θ)\mathbf{T}(\rho,\theta) = (\rho^2, 2\theta), il semipiano v>0v > 0 va sul piano privato del semiasse positivo. ∬E4(u2+v2)=∫0π∫014ρ3dρ dθ=π\iint_E4(u^2 + v^2) = \int_0^\pi\int_0^14\rho^3d\rho\,d\theta = \pi. x2+y2=(u2+v2)2x^2 + y^2 = (u^2 + v^2)^2.

Teoria collegata