Skript
Skript Skript
Es seien f und g verallgemeinert differenzierbar, also g ∈ H 1 (C) und f ∈ H 1 (R n ) für den Sobolevraum H 1 . Die übliche Norm in H 1 enthält die 2–Norm der Funktion und ihrer Ableitungen, es gilt also ∑ ||f|| H 1 = ||(D α f)|| 2 2 |α|≤1, α∈N n = ∑ || ̂D α f|| 2 2 |α|≤1 = ∑ ||ξ α 2 ̂f|| 2 = = |α|≤1 ∫ (1 + ξ1 2 + ξ2 2 + . . . + ξn)| 2 ̂f(ξ)| 2 dξ R ∫ n (1 + ||ξ|| 2 2)| ̂f(ξ)| 2 dξ. R n Allgemein definieren wir für β ∈ R ≥0 und f mit (1 + ||ξ|| 2 ) β | ̂f(ξ)| 2 integrierbar ∫ ||f|| 2 H = (1 + ||ξ|| 2 ) β | ̂f(ξ)| 2 dξ β R n und ||g|| 2 H β = ∫ S n−1 ∫ R (1 + σ 2 ) β |ĝ(θ, σ)| 2 dσ dθ. Für ganzzahliges β ist die so definierte Norm äquivalent zur üblichen Norm in H β . Satz 4.2.31. (Schlechtgestelltheit der Radontransformation) Sei K ∈ R n kompakt. Für alle α ≥ 0 gibt es c, c ′ (abhängig von K und α), so dass c||f|| H α (R n ) ≤ ||Rf|| H α+(n−1)/2 (C) ≤ c′ ||f|| H α (R n ) ∀f ∈ C ∞ 0 (K). Insbesondere gilt für g = Rf ||R −1 g|| L 2 (R n ) ≤ 1 c ||g|| H (n−1)/2 (C) ≤ c′ c ||R−1 g|| L 2 (R n ) R −1 ist also genau dann stetig, falls die Norm im Datenraum Ableitungen der Ordnung (n − 1)/2 enthält, also ist die Aufgabe der Inversion von R schlecht gestellt von der Ordnung (n − 1)/2. □ 63
Beweis: Sei ohne Einschränkung α = 0. Dann gilt ∫ ∫ ||Rf|| 2 = (1 + σ 2 ) (n−1)/2 | ̂Rf(θ, σ)| 2 dσ dθ H (n−1)/2 S n−1 R = (2π) n−1 ∫ S n−1 = 2(2π) n−1 ∫ R n ∫ (1 + σ 2 ) (n−1)/2 | ̂f(σθ)| 2 dσ dθ R (1 + ||ξ|| 2 ) (n−1)/2 ||ξ|| n−1 | ̂f(ξ)| 2 dξ. Der Bruch ist immer größer als 1, also gilt insbesondere ||Rf|| 2 ≥ 2(2π) n−1 ||f|| 2 H (n−1)/2 L 2 (R n ) . Zum Beweis der zweiten Ungleichung: Für ||ξ|| ≥ 1 ist der Bruch durch 2 (n−1)/2 nach oben beschränkt. Für ||ξ|| ≤ 1 existiert das Integral über den Bruch. Wir ziehen aus dem Integral das Supremum von ̂f(ξ) 2 heraus und beachten, dass ∫ ∫ | ̂f(ξ)| = | e −ixξ f(x) dx| ≤ |f(x)| dx ≤ |K| ||f|| L 2 (K) K K mit Cauchy–Schwartz. □ Bemerkung: Insbesondere ist die wichtigste Anwendung der Radon–Transformation für die Computertomographie schlecht gestellt von der Ordnung 1/2 und damit deutlich besser gestellt als die Aufgabe der Berechnung der ersten Ableitung, die schlecht gestellt ist von der Ordnung 1 (nach Definition des Grades der Schlechtgestelltheit). Bemerkung: Die Schlechtgestelltheit ist bereits in der Inversionsfomel 4.2.27 sichtbar. Die Fehler in der Datenfunktion ĝ(θ, σ) werden im Zwischenergebnis mit |σ| (n−1) verstärkt, für großes |σ| erhalten wir also auch beliebig hohe Rekonstruktionsfehler bei kleinem Datenfehler. Allerdings wird durch die anschließende Integration noch einmal geglättet, so dass sich der Grad der Schlechtgestelltheit insgesamt auf die Hälfte reduziert. 4.3 Implementation Im letzten Abschnitt haben wir nachgewiesen, dass die Inversion der Radon– Transformation schlecht gestellt ist von der Ordnung 1/2. Zur Implementation der Algorithmen müssen wir also insbesondere regularisieren. Wir werden dies durch Projektion auf stabile Unterräume tun, wir suchen die Lösung also in Unterräumen, die von Singulärvektoren mit großen zugehörigen Singulärwerten 64
- Seite 13 und 14: Wir entwickeln die gesuchte Funktio
- Seite 15 und 16: Kapitel 2 Grundlagen Wie wir in den
- Seite 17 und 18: • A ist stetig bezüglich L 2 (Ω
- Seite 19 und 20: Definition 2.1.3 (kleinste Quadrate
- Seite 21 und 22: Bemerkung 2.1.12. karl • A ∗ Ax
- Seite 23 und 24: Da A kompakt ist, besitzt (y n ) n
- Seite 25 und 26: Theorem 2.2.10. Sei A ∈ L(X, Y )
- Seite 27 und 28: Durch Koeffizientenvergleich folgt:
- Seite 29 und 30: Kapitel 3 Regularisierungstheorie I
- Seite 31 und 32: Beweis. Sei y ∈ D(T † ) beliebi
- Seite 33 und 34: und weiterhin gilt. lim δC α(δ)
- Seite 35 und 36: Normalengleichungen. Außerdem ist
- Seite 37 und 38: Durch die Unstetigkeit an 0 wird di
- Seite 39 und 40: Kapitel 4 Inverse Probleme der Inte
- Seite 41 und 42: zu θ(ϕ k ) durchgeführt, daher h
- Seite 43 und 44: 4.2 Funktionalanalytische Grundlage
- Seite 45 und 46: ̂f a (ξ) = (2π) −n/2 ∫ R n =
- Seite 47 und 48: Satz 4.2.5. (Faltungssatz) Seien f,
- Seite 49 und 50: aber f ∉ L 2 (R 2 ). Die Fouriert
- Seite 51 und 52: Wir definieren also allgemein die A
- Seite 53 und 54: keinen Sinn, den das hintere Integr
- Seite 55 und 56: S(C) definiert durch (Rf)(θ, s) =
- Seite 57 und 58: = (2π) 1/2 ̂f(ξ) □ Satz 4.2.20
- Seite 59 und 60: ist ˜f(x) zumindest schon mal posi
- Seite 61 und 62: Theorem 4.2.27. Sei 0 ≤ α < n. D
- Seite 63: Cormack eine alternative Formel her
- Seite 67 und 68: g k = (2π/h) −n ∫ g(ξ)e i(kh)
- Seite 69 und 70: Man kann sich fragen, was die Bandb
- Seite 71 und 72: Beweis: Sei w ∈ S(R n ) eine Test
- Seite 73 und 74: c k = 1 (2π) 1/2 ∫ ||y|| |x| (De
- Seite 75 und 76: 4.4 Was fehlt? Fanbeam. Interlaced
- Seite 77 und 78: die sogenannte Landweber iteration.
- Seite 79 und 80: Primal Adjungiert messungen unbekan
- Seite 81 und 82: Nichtlineare Probleme Wir diskutier
- Seite 83 und 84: We finally mention that any suitabl
- Seite 85 und 86: 5.5 Konjugierte Gradienten (Frank)
- Seite 87 und 88: Satz 5.5.2. Gesucht sei die Lösung
- Seite 89 und 90: 3. Sei m j = 1, ω = 1, C k = (A k
- Seite 91 und 92: linearisiert. Seien also nun die R
- Seite 93 und 94: entsprechend statistisch optimale R
- Seite 95 und 96: L(f) heißt Likelihoodfunktion. Man
- Seite 97 und 98: und das ist die erweiterte Form der
- Seite 99 und 100: Literaturverzeichnis [1] Ȧ. Björc
Es seien f und g verallgemeinert differenzierbar, also g ∈ H 1 (C) und f ∈<br />
H 1 (R n ) für den Sobolevraum H 1 . Die übliche Norm in H 1 enthält die 2–Norm<br />
der Funktion und ihrer Ableitungen, es gilt also<br />
∑<br />
||f|| H 1 =<br />
||(D α f)|| 2 2<br />
|α|≤1, α∈N n<br />
= ∑<br />
|| ̂D α f|| 2 2<br />
|α|≤1<br />
= ∑<br />
||ξ α 2 ̂f|| 2<br />
=<br />
=<br />
|α|≤1<br />
∫<br />
(1 + ξ1 2 + ξ2 2 + . . . + ξn)| 2 ̂f(ξ)| 2 dξ<br />
R<br />
∫<br />
n (1 + ||ξ|| 2 2)| ̂f(ξ)| 2 dξ.<br />
R n<br />
Allgemein definieren wir für β ∈ R ≥0 und f mit (1 + ||ξ|| 2 ) β | ̂f(ξ)| 2 integrierbar<br />
∫<br />
||f|| 2 H<br />
= (1 + ||ξ|| 2 ) β | ̂f(ξ)| 2 dξ<br />
β<br />
R n<br />
und<br />
||g|| 2 H β =<br />
∫<br />
S n−1<br />
∫<br />
R<br />
(1 + σ 2 ) β |ĝ(θ, σ)| 2 dσ dθ.<br />
Für ganzzahliges β ist die so definierte Norm äquivalent zur üblichen Norm in<br />
H β .<br />
Satz 4.2.31. (Schlechtgestelltheit der Radontransformation)<br />
Sei K ∈ R n kompakt. Für alle α ≥ 0 gibt es c, c ′ (abhängig von K und α), so<br />
dass<br />
c||f|| H α (R n ) ≤ ||Rf|| H α+(n−1)/2 (C) ≤ c′ ||f|| H α (R n ) ∀f ∈ C ∞ 0 (K).<br />
Insbesondere gilt für g = Rf<br />
||R −1 g|| L 2 (R n ) ≤ 1 c ||g|| H (n−1)/2 (C) ≤ c′<br />
c ||R−1 g|| L 2 (R n )<br />
R −1 ist also genau dann stetig, falls die Norm im Datenraum Ableitungen der<br />
Ordnung (n − 1)/2 enthält, also ist die Aufgabe der Inversion von R schlecht<br />
gestellt von der Ordnung (n − 1)/2.<br />
□<br />
63