📖 Nauka z symulacjami · #33

Analiza zespolona w fizyce i inżynierii

Analiza zespolona — nauka o funkcjach zmiennej zespolonej różniczkowalnych w sensie zespolonym — to jedno z najpotężniejszych narzędzi fizyki teoretycznej. Funkcja, która jest różniczkowalna raz w sensie zespolonym, jest automatycznie nieskończenie wiele razy różniczkowalna, daje się przedstawić jako szereg potęgowy i podlega globalnym ograniczeniom niewidocznym w analizie rzeczywistej. Od przepływu płynów 2D i elektrostatyki po całki konturowe definiujące propagatory kwantowe i metody spektralne przetwarzania sygnałów — analiza zespolona jest niezastąpiona.

Liczby zespolone i płaszczyzna Arganda

Liczba zespolona z = x + iy to punkt na płaszczyźnie Arganda, z częścią rzeczywistą Re(z) = x i częścią urojoną Im(z) = y. Moduł |z| = √(x² + y²) i argument arg(z) = arctan(y/x) dają postać biegunową z = r e, gdzie wzór Eulera

Wzór Eulera — najpiękniejsze równanie w matematyce e = cos θ + i sin θ
e + 1 = 0   (tożsamość Eulera)
Dowód: szereg Taylora e = 1 + iθ − θ²/2! − iθ³/3! + … = cos θ + i sin θ

Mnożenie w postaci biegunowej: z₁ z₂ = r₁r₂ ei(θ₁+θ₂) — moduły się mnożą, argumenty się dodają. Twierdzenie de Moivre'a: (e)ⁿ = einθ wyraża sin(nθ) i cos(nθ) przez potęgi sin/cos. n-te pierwiastki z jedności to zk = e2πik/n, rozmieszczone równomiernie na okręgu jednostkowym.

Funkcje analityczne i równania Cauchy'ego-Riemanna

Funkcja f(z) = u(x,y) + iv(x,y) jest różniczkowalna w sensie zespolonym w punkcie z₀, jeśli granica (f(z₀+h) − f(z₀))/h istnieje dla h → 0 z dowolnego kierunku na płaszczyźnie zespolonej. Nakłada to silne ograniczenie — równania Cauchy'ego-Riemanna:

Równania Cauchy'ego-Riemanna ∂u/∂x = ∂v/∂y     i     ∂u/∂y = −∂v/∂x

Konsekwencja: zarówno u, jak i v spełniają równanie Laplace'a
∇²u = 0     i     ∇²v = 0
u i v nazywane są funkcjami harmonicznie sprzężonymi

Funkcję różniczkowalną w każdym punkcie pewnego obszaru nazywamy tam analityczną (lub holomorficzną). Funkcje analityczne mają niezwykłe własności:

Popularne funkcje analityczne

Wielomiany, ez, sin z, cos z, ln z (z cięciem gałęziowym) są analityczne poza swoimi punktami osobliwymi. Kluczowe tożsamości:

sin z = (eiz − e−iz)/(2i)  |  cos z = (eiz + e−iz)/2
sinh z = (ez − e−z)/2  |  sin(iz) = i sinh z
ln z = ln|z| + i arg(z) (wielowartościowa; cięcie gałęziowe na ujemnej osi rzeczywistej)

Całki konturowe i twierdzenie Cauchy'ego

Funkcje zespolone całkujemy wzdłuż dróg (konturów) na płaszczyźnie zespolonej. Twierdzenie Cauchy'ego: jeśli f(z) jest analityczna w obszarze jednospójnym, to ∮C f(z) dz = 0 dla dowolnego konturu zamkniętego C w tym obszarze. To zespolony odpowiednik niezależności całki od drogi dla pól zachowawczych.

Wzór całkowy Cauchy'ego f(z₀) = (1/2πi) ∮C f(z)/(z − z₀) dz   (C obejmuje z₀)

f⁽ⁿ⁾(z₀) = (n!/2πi) ∮C f(z)/(z − z₀)ⁿ⁺¹ dz

Konsekwencja: wartości analitycznej f wewnątrz C są w pełni wyznaczone przez wartości brzegowe

Wzór całkowy Cauchy'ego ma nietrywialną moc: znajomość f na brzegu koła wyznacza ją wszędzie wewnątrz. Dlatego funkcje analityczne są tak sztywne w porównaniu z gładkimi funkcjami rzeczywistymi. Wzór daje też za darmo wszystkie pochodne, dowodząc, że analityczna ⟹ C∞.

Osobliwości, szeregi Laurenta i twierdzenie o residuach

Tam, gdzie funkcje przestają być analityczne, mają punkty osobliwe. W pierścieniu r₁ < |z − z₀| < r₂ wokół punktu osobliwego f rozwija się w szereg Laurenta:

Szereg Laurenta i residuum f(z) = ∑n=−∞ aₙ (z − z₀)ⁿ

Residuum: Res(f, z₀) = a₋₁ (współczynnik przy (z−z₀)⁻¹)

Typy osobliwości:
 • Usuwalna: część główna (wyrazy z n<0) = 0, osobliwość można usunąć
 • Biegun rzędu m: część główna skończona, aₙ = 0 dla n < −m
 • Istotna: nieskończona część główna (np. e1/z w 0)

Twierdzenie o residuach to koń pociągowy stosowanej analizy zespolonej:

Twierdzenie o residuachC f(z) dz = 2πi ∑ Res(f, zₖ)   (suma po biegunach wewnątrz C)

Także: Res(f, z₀) = limz→z₀ [(z−z₀) f(z)] dla bieguna prostego
Dla bieguna rzędu m: Res = (1/(m−1)!) lim dm−1/dzm−1 [(z−z₀)ᵐ f(z)]

Dzięki sprytnemu doborowi konturów (zwykle duży półokrąg zamknięty w górnej lub dolnej półpłaszczyźnie) twierdzenie o residuach pozwala obliczyć całki rzeczywiste, które opierają się metodom elementarnym. Klasyczne przykłady:

∫₋∞ dx/(1+x²) = π   (jeden biegun prosty w i)
∫₀ sin(x)/x dx = π/2   (lemat Jordana, kontur „dziurka od klucza")
∫₀ dθ/(a + b cos θ) = 2π/√(a²−b²)   (kontur okręgu jednostkowego)

Odwzorowania konforemne

Funkcja analityczna w = f(z) z f′(z) ≠ 0 jest odwzorowaniem konforemnym: zachowuje kąty między krzywymi. To czyni odwzorowania konforemne kluczowym narzędziem rozwiązywania dwuwymiarowych zagadnień brzegowych: odwzorować skomplikowany obszar na prosty (koło, półpłaszczyznę), rozwiązać tam, odwzorować z powrotem.

Przekształcenia Möbiusa

Rodzina w = (az+b)/(cz+d), ad−bc ≠ 0, przekształca okręgi i proste w okręgi i proste. Każde przekształcenie Möbiusa jest złożeniem przesunięć, obrotów, skalowań i inwersji w = 1/z. Trzy punkty wyznaczają jednoznacznie przekształcenie Möbiusa; dwustosunek jest zachowywany.

Transformacja Żukowskiego

Odwzorowanie Żukowskiego w = z + R²/z zamienia okrąg na płaszczyźnie z w profil lotniczy na płaszczyźnie w. W połączeniu z cyrkulacją Γ (warunek Kutty) daje przepływ generujący siłę nośną wokół skrzydła:

Profil Żukowskiego i twierdzenie Kutty-Żukowskiego o sile nośnej w = z + c²/z (odwzorowanie Żukowskiego, c = R dla profilu symetrycznego)
Potencjał zespolony: F(z) = U(z + R²/z) + (iΓ/2π) ln(z/R)
Siła nośna na jednostkę rozpiętości: L = ρ U Γ (twierdzenie Kutty-Żukowskiego)
Γ = 4πUR sin α (dla cienkiego symetrycznego profilu pod kątem natarcia α)

Odwzorowanie Schwarza-Christoffela

Wzór Schwarza-Christoffela odwzorowuje górną półpłaszczyznę (lub koło jednostkowe) na wnętrze wielokąta o zadanych kątach wierzchołkowych: dw/dz = K ∏ (z − xₖ)(αₖ/π − 1), gdzie x₁<x₂<…<xₙ to rzeczywiste przeciwobrazy wierzchołków wielokąta, a αₖ to kąty wewnętrzne. Leży to u podstaw obliczania pól elektrycznych wewnątrz falowodów i wokół ostrych krawędzi przewodników.

Zastosowania fizyczne

Bezwirowy przepływ płynu 2D

Dla nieściśliwego, bezwirowego przepływu 2D pole prędkości wyprowadza się z potencjału prędkości φ: v = ∇φ, gdzie ∇²φ = 0. Ponieważ φ spełnia równanie Laplace'a, jest ona częścią rzeczywistą analitycznego potencjału zespolonego F(z) = φ(x,y) + iψ(x,y), gdzie ψ to funkcja prądu (∇φ · ∇ψ = 0: linie prądu ⊥ ekwipotencjały). Prędkość zespolona to dF/dz = vₓ − ivᵧ. Odwzorowania konforemne przekształcają wtedy proste przepływy (jednostajny, wir, źródło) w dowolne geometrie — podstawa klasycznej aerodynamiki sprzed ery CFD.

Elektrostatyka 2D i metoda obrazów

Potencjał elektrostatyczny V spełnia ∇²V = −ρ/ε₀. W obszarach wolnych od ładunku ∇²V = 0, więc V jest harmoniczny. Metoda obrazów umieszcza fikcyjne ładunki poza obszarem, aby spełnić warunki brzegowe: ładunek q w pobliżu uziemionej płaszczyzny wprowadza ładunek lustrzany −q w punkcie odbicia. Odwzorowania konforemne rozszerzają to na geometrie cylindryczne, sferyczne i bardziej złożone geometrie przewodników.

Impedancja obwodów prądu zmiennego

Dla sygnałów sinusoidalnych o częstości kołowej ω napięcia i prądy są reprezentowane jako zespolone fazory V = V₀ eiωt. Impedancje są zespolone: ZR = R, ZL = iωL, ZC = 1/(iωC). Impedancja zespolona Z(ω) koduje jednocześnie moduł (|Z|) i kąt fazowy (arg Z). Rezonans zachodzi, gdy Im Z = 0; pochłaniana moc wynosi Re(V I*)/2. Dlatego liczby zespolone wykłada się na każdym kierunku elektrotechniki.

Transformaty Fouriera i Laplace'a, mechanika kwantowa

Transformata Fouriera F(ω) = ∫ f(t) e−iωt dt jest zespoloną funkcją częstotliwości. Bieguny F(ω) na zespolonej płaszczyźnie ω wyznaczają zachowanie w dziedzinie czasu: bieguny w dolnej półpłaszczyźnie → wykładniczo tłumione oscylacje (stabilność). Transformata Laplace'a F(s) = ∫₀ f(t) e−st dt stosowana jest dla sygnałów przyczynowych; jej odwrotność to całka konturowa Bromwicha zamknięta w lewej półpłaszczyźnie. W mechanice kwantowej całka po trajektoriach Feynmana i propagator K(x,t; x₀,t₀) są w naturalny sposób wyrażane za pomocą całek konturowych, a bieguny macierzy S na zespolonej płaszczyźnie pędu odpowiadają stanom związanym i rezonansom.

Przedłużenie analityczne i funkcja dzeta Riemanna

Jeśli dwie funkcje analityczne pokrywają się na krzywej lub zbiorze otwartym, pokrywają się wszędzie na wspólnym obszarze analityczności — to przedłużenie analityczne. Pozwala ono rozszerzyć funkcję poza jej pierwotny obszar: funkcja dzeta Riemanna ζ(s) = ∑ n−s jest zbieżna tylko dla Re(s) > 1, ale przedłużenie analityczne rozszerza ją na całą ℂ z wyjątkiem s = 1. Hipoteza Riemanna (wszystkie nietrywialne zera leżą na Re(s) = 1/2) jest równoważna najostrzejszym znanym oszacowaniom funkcji liczącej liczby pierwsze π(x) i pozostaje jednym z Problemów Milenijnych.

Funkcja dzeta Riemanna i równanie funkcyjne ζ(s) = ∑n=1 n−s = ∏p pierwsze (1 − p−s)⁻¹
Równanie funkcyjne: ζ(s) = 2ˢ πˢ⁻¹ sin(πs/2) Γ(1−s) ζ(1−s)
Zera trywialne: s = −2, −4, −6, …
Zera nietrywialne: s = ½ + iγₙ (hipoteza Riemanna: wszystkie mają Re(s) = ½)

Najważniejsze wnioski

  • Różniczkowalność zespolona (analityczność) jest znacznie bardziej restrykcyjna niż różniczkowalność rzeczywista: analityczna ⟹ C∞, szereg potęgowy, ograniczenia globalne.
  • Równania Cauchy'ego-Riemanna łączą różniczkowalność zespoloną z analizą harmoniczną; część rzeczywista i urojona funkcji analitycznych spełniają równanie Laplace'a.
  • Twierdzenie o residuach zamienia zamknięte całki konturowe w sumy residuów, co pozwala systematycznie obliczać rzeczywiste całki niewłaściwe.
  • Odwzorowania konforemne zachowują kąty i przekształcają skomplikowane brzegi w proste, rozwiązując zagadnienia z równaniem Laplace'a w przepływie płynów 2D, elektrostatyce i przewodnictwie cieplnym.
  • Zespolone fazory, transformaty Fouriera/Laplace'a i propagatory kwantowe wykorzystują algebraiczną i geometryczną moc ℂ.
  • Przedłużenie analityczne nadaje sens funkcjom takim jak ζ(s) poza ich pierwotnym obszarem; hipoteza Riemanna łączy zera ζ z głęboką strukturą liczb pierwszych.