Kramersovy-Kronigovy relace

Z ωικι.matfyz.cz
Verze z 11. 6. 2012, 18:03, kterou vytvořil 178.248.252.203 (diskuse)

(rozdíl) ← Starší verze | zobrazit aktuální verzi (rozdíl) | Novější verze → (rozdíl)
Přejít na: navigace, hledání

Kramersovy–Kronigovy relace umožňují spočítat reálnou část odezvy lineárního pasivního systému, známe-li imaginární části odezvy při všech frekvencích (nebo naopak určit imaginární část ze znalosti části reálné). Při analýze optických konstant hrají důležitou roli a jsou hojně využívány, protože platí např. pro elektrickou vodivost σ (vystupující v ohmově zákoně j(ω)=σ(ω)E(ω) ). Abychom mohli Kramers–Kronigovu analýzu provést, musí funkce odezvy α(ω)=α1(ω)+iα2(ω) splňovat:

  1. Póly α(ω) jsou všechny nad reálnou osou
  2. Při integraci přes nekonečně velkou polokružnici v horní polorovině komplexní roviny, je integrál z α(ω)/ω roven nule
  3. Pro $ \omega $$ \in $$ \mathbb{R} $ je α1(ω) sudá a α2(ω) lichá

Potom platí:

$ \alpha $1(ω) = 2/π VP ∫0 s α2(s) / (s2 - ω2)ds.

a

α2(ω) = -2/π VP ∫0 ω α1(s) / (s2 - ω2) ds = - 2ω/π VP ∫0 α1(s) / (s2 - ω2)ds.

P značí hlavní hodnotu integrálu.

Důkaz[editovat | editovat zdroj]

Budeme integrovat po křivce znázorněné v http://en.wikipedia.org/wiki/Kramers–Kronig_relations#Derivation.

Protože uvnitř oblasti ohraničené integrační křivkou nejsou žádné póly, je integrál po této křivce z funkce α(s)/(s-ω) roven nule.

Protože funkce odezvy α klesá u nekonečen dostatečně rychle, je integrál přes velkou polokružnici roven nule.

Zbývá tedy rovnost -∞ω + ∫polokružnice kolem ω + ∫ω = 0. Sem dosadíme vyjádření Reziduové věty nebo Cauchyovy věty (možná se ta věta jmenuje po někom úplně jiném). Ta věta, kterou chci využít, říká, že když integruji po části kružnice kolem pólu, tak výsledek je polokružnice kolem ω f(s) ds = πi Res(f).

Pozor! Právě zapsaný výsledek platí, pokud integruji v kladném směru, což je PROTI směru hodinových ručiček. Když si představíte průběh funkce eix, tak pro x stoupající máte "kladný směr obíhání" - a ten je proti směru hodinových ručiček.

Vzpomeneme na výpočet reziduí: http://cs.wikipedia.org/wiki/Reziduum#V.C3.BDpo.C4.8Det_rezidu.C3.AD . Naše zkoumaná funkce, kterou integrujeme, je α(s)/(s-ω), proto reziduum je α(ω). Máme tedy rovnici

-∞ω - πiα(ω) + ∫ω = 0

Nakonec tedy:

α(ω) = 1/πi VP ∫-∞ α(s)/(s-ω)ds.

Tuto rovnost si napíšeme zvlášť pro reálnou a pro imaginární část:

α1(ω) = Re{1/πi VP ∫-∞1(s) + iα2(s)]/(s-ω)ds} = 1/πi VP ∫-∞2(s)/(s-ω)ds = 1/π VP ∫-∞ α2(s)/(s-ω)ds.
α2(ω) = Im{1/πi VP ∫-∞1(s) + iα2(s)]/(s-ω)ds} = Im{-i/π VP ∫-∞1(s) + iα2(s)]/(s-ω)ds} = -1/π Re{VP ∫-∞1(s) + iα2(s)]/(s-ω)ds} = -1/π VP ∫-∞ α1(s)/(s-ω)ds.

Nakonec je nutné integrály přes záporné frekvence převést na integrály přes kladné frekvence. Při tom využijeme předpokladů (α1 sudá a α2 lichá). Převádíme tedy (t=-s, dt = -ds):

-∞0 α2(s)/(s-ω)ds = ∫0 α2(-t)/(-t-ω)-dt = ∫02(t)/(-t-ω)dt = ∫0 α2(t)/(t+ω)dt a
-∞0 α1(s)/(s-ω)ds = ∫0 α1(-t)/(-t-ω)-dt = ∫0 α1(t)/(-t-ω)dt = -∫0 α1(t)/(t+ω)dt.

Nyní mám tedy:

α1(ω) = 1/π [∫0 α2(t)/(t+ω)dt + VP ∫0 α2(s)/(s-ω)ds] a
α2(ω) = -1/π [-∫0 α1(t)/(t+ω)dt' + VP ∫-∞ α1(s)/(s-ω)ds].

Protože všechny integrály jsou od nuly do nekonečna, sloučím je. Budu tedy vždy integrovat součet dvou zlomků - dám je na společného jmenovatele:

  • při výpočtu α1(ω) budu integrovat α2(s)/(s+ω) + α2(s)/(s-ω) = [(s-ω)α2(s) + (s+ω)α2(s)] / (s+ω)(s-ω) = 2sα2(s) / (s2 - ω2),
  • při výpočtu α2(ω) budu integrovat 1(s)/(s+ω) + α1(s)/(s-ω) = [-(s-ω)α1(s) + (s+ω)α1(s)] / (s+ω)(s-ω) = 2ωα1(s) / (s2 - ω2).

Využití pro odrazivost[editovat | editovat zdroj]

Použijeme Kramersovy-Kronigovy relace pro amplitudovou odrazivost r(ω). Protože intenzitní odrazivostní koeficient je R(ω)=r(ω)r(ω)*, dostáváme

ln(r(ω)) = ln[R1/2(ω)] +iθ(ω) (zde je krásně zřetelné rozdělení ne reálnou a imaginární část).

Nyní můžeme vypočítat fázi amplitudové odrazivosti:

θ(ω) = - ω/π VP ∫0 ln[R(s)] / (s2 - ω2) ds = -ω/π VP ∫0 ln[R(s)] × 1/(s2 - ω2) ds.

Využijeme integrování per partes: logaritmus zderivujeme, podíl 1/(s22) zintegrujeme a otočíme znaménko před integrálem. Integrál podílu je

∫1 / (s2 - ω2) ds = ∫ 1/2ω / (s-ω) - 1/2ω / (s+ω) ds = 1/2ω [ln(s-ω) - ln(s+ω)] = 1/2ω ln[(s-ω)/(s+ω)].

Dostáváme tedy:

θ(ω) = + ω/π VP ∫0 dln[R(s)]/ds × {1/2ω ln[(s-ω)/(s+ω)]} ds = 1/2π VP ∫0 dln[R(s)]/ds ln[(s-ω)/(s+ω)] ds.

Spektrální oblasti, ve kterých je odrazivost konstantní, k integrálu nepřispívají. Navíc ani neřispívají příliš oblasti s≫ω a s≪ω, protože tam je hodnota funkce ln[(s-ω)/(s+ω)] malá.