inverses Kurvenintegral

titan99_

Also die Bezierkurve berechne ich eigentlich 2 oder 3-Dimensional, und habe dazu noch 2 Parallele Kurven (http://en.wikipedia.org/wiki/Parallel_curve).
Die Parallelen Kurven sind aber ziemlich unregelmässig, vorallem bei grosser Krümmung. Und ich dachte es sei sowieso besser konstante Abstände zu haben, weil daraus nachher Flächen werden.

Jester

.filmor schrieb:

Wofür brauchst du das eigentlich? Denn damit bekommst du nur auf der Kurve äquidistante Punkte, im Raum sind sie das ja nicht.

Wenn man zum Beispiel eine Kamera entlang der Kurve fahren lassen möchte, ist das genau das richtige, weil die Kamera dann konstante Geschwindigkeit hat.

Jester

Mathematikker schrieb:

Wenn eine Abbildung/Funktion f: IR² -> IR mehr als eine Dimension aus dem Definitionsbereich nutzt (was die Notation IR² eigentlich suggeriert), dann ist sie nicht invertierbar

Nein. R^n und R sind gleichmächtig, also gibt es eine Bijektion, diese ist insbesondere invertierbar.

Allerdings gibt es keine stetige Abbildung, deren Umkehrung ebenfalls stetig ist.

Jester schrieb:

Mathematikker schrieb:

Wenn eine Abbildung/Funktion f: IR² -> IR mehr als eine Dimension aus dem Definitionsbereich nutzt (was die Notation IR² eigentlich suggeriert), dann ist sie nicht invertierbar

Nein. R^n und R sind gleichmächtig, also gibt es eine Bijektion, diese ist insbesondere invertierbar.

Allerdings gibt es keine stetige Abbildung, deren Umkehrung ebenfalls stetig ist.

Diese Erkenntnis bringt uns nicht weiter, da sie praktisch keine Relevanz hat. Sonst müsste man sich ja nicht mit dem IR^n abquälen, wenn man sinnvoll auf IR arbeiten könnte und das Ergebnis wieder nach IR^n übertragen könnte.
Kurz: sie sind nicht isomorph

Aber du hast Recht, da ich die Gleichmächtigkeit nicht beachtet habe, ich habe schon ein paar Stufen praktischer gedacht

.filmor

Aber „isomorph“ ist auch etwas zu stark. Hier relevant ist ja, wie Jester schon sagte, dass sie nicht homöomorph sind. Jeder Isomorphismus ist auch ein Homöomorphismus (sogar ein Diffeomorphismus denke ich) aber nicht andersherum.

titan99 schrieb:

Also die Bezierkurve berechne ich eigentlich 2 oder 3-Dimensional, und habe dazu noch 2 Parallele Kurven (http://en.wikipedia.org/wiki/Parallel_curve).
Die Parallelen Kurven sind aber ziemlich unregelmässig, vorallem bei grosser Krümmung. Und ich dachte es sei sowieso besser Konstante Abstände zu haben, weil daraus nachher Flächen werden.

Warum? Mit dem normalen Verfahren hast du doch mehr Punkte im Bereich starker Krümmung wenn ich mich nicht täusche, das erscheint mir schon recht praktisch für weitere Berechnungen. Berechnest du die Fläche die von der Kurve und ihrer parallelen Kurve umrandet wird? Das ließe sich doch sicher irgendwie abkürzen …

pasti

.filmor schrieb:

Aber „isomorph“ ist auch etwas zu stark. Hier relevant ist ja, wie Jester schon sagte, dass sie nicht homöomorph sind. Jeder Isomorphismus ist auch ein Homöomorphismus (sogar ein Diffeomorphismus denke ich) aber nicht andersherum.

Wann eine Abbildung ein Isomophismus bzw. ein (Homo)Morphismus ist wird dadurch festgelegt, welche Struktur man den zugrundeliegenden Mengen zuordnet.

Ein Isomorphismus ist eine bijjektive Abbildung, welche in beide Richtugen strukturerhaltend ist.

Man kann IR als Menge von Punkten ohne jede Struktur auffassen, dann bezeichnet man jede bijektive Abbildung als Isomorphismus.

Man kann IR auch als Gruppe auffassen, dann sind die bijektiven Gruppenisomorphismen die Isomorphismen, denn in diesem Fall kann man zeigen, dass die Umkehrung eines bjiektiven Gruppenhomomorphismus wieder ein Gruppenhomorphismus ist.

Dasselbe gilt für Ringe und Ringisomorphismen.

Wenn man IR als topologischen Raum auffasst, sind die stetigen Abbildungen die Morphismen. Es gibt jedoch bijektive stetige Funktionen, wo die Umkehrabbildung nicht stetig ist. Also ist nich jede solche Funktion ein Isomophismus. Man muss also zusätlich noch fordern, dass die Umkehrfunktion auch stetig ist.

Link: http://de.wikipedia.org/wiki/Kategorientheorie

.filmor

Kann man sich nicht darauf einigen, dass wenn man nichts dazu sagt man einen Vektorraumisomorphismus meint? :p
Hast natürlich trotzdem recht…

Wollen wir uns nicht darauf einigen, dass es hier ein praktisch zu lösendes Problem gibt und die Kategorientheorie hier fehl am Platz ist?

pasti

Mathematikker schrieb:

Wollen wir uns nicht darauf einigen, dass es hier ein praktisch zu lösendes Problem gibt und die Kategorientheorie hier fehl am Platz ist?

Pha... und du nennst dich Matimatiker.

Das ursprüngliche Problem wurde von Jester schon etwa im 3. Post gelöst.

Walli

Etwa im 3. Post klingt aber auch mehr nach Ingenieur als nach Mathematiker .

In der derivierten Kategorie von perversen C^2(R)-Modulgarben wird das Problem mit der Fouriertransformation trivial.

pasti schrieb:

Mathematikker schrieb:

Wollen wir uns nicht darauf einigen, dass es hier ein praktisch zu lösendes Problem gibt und die Kategorientheorie hier fehl am Platz ist?

Pha... und du nennst dich Matimatiker.

Das ursprüngliche Problem wurde von Jester schon etwa im 3. Post gelöst.

Du hast recht, wir sollten einfach die Existenz einer Lösung zeigen, mehr braucht man ja nicht

pasti

Mathematikker schrieb:

Du hast recht, wir sollten einfach die Existenz einer Lösung zeigen, mehr braucht man ja nicht

Woot? Man kann die Parametrisierung immer expizit als Integral darstellen. In diesem Fall wird es wohl sogar noch ausreichen dieses dann numerisch zu Lösen.

http://de.wikibooks.org/wiki/Diffgeo:_Kurventheorie:_Parameterisierung_nach_Bogenlänge

titan99_

.filmor schrieb:

Warum? Mit dem normalen Verfahren hast du doch mehr Punkte im Bereich starker Krümmung wenn ich mich nicht täusche, das erscheint mir schon recht praktisch für weitere Berechnungen. Berechnest du die Fläche die von der Kurve und ihrer parallelen Kurve umrandet wird? Das ließe sich doch sicher irgendwie abkürzen …

Die Flächen berechne ich aus den parallelen Kurven. Eine der parallelen Kurven hat hat dann aber im Bereich starker Krümmung weniger Punkte als im Bereich schwacher Krümmung.
Die Punkte der Bezier-Kurven berechne ich aus einem Polygon, ich weiss auch nicht wie ich das abkürzen soll.

Kann man eine Funktion numerisch invertieren??

Also irgendwie pack ich es nicht mit dem parametrisieren nach der Bogenlänge. Kann mir vielleicht jemand am Beispiel der Normalparabel y = x^2 erklären wie das geht??

Integrieren tu ich sonst mit dem Wolfram-Integrator http://integrals.wolfram.com/index.jsp

.filmor

Wie immer wenn's numerisch wird hilft wohl das Newtonverfahren.

Du hast ja s(t) gegeben (in Integralform) und möchtest t sd. s(t) = n·∆s ist, also eine Nullstelle von f(t) = s(t) - n·∆s.
Dann bekommst du das Verfahren t_n+1 = t_n - (s(t_n) - n·∆s)/s'(t_n).

s'(t) kennst du (|B'(t)|), vielleicht kann man durch entsprechendes Einsetzen von s(t) da auch noch was vereinfachen, aber zur Not nimmst du da halt da auch schon nen Näherungswert, gibt ja genug gute Quadraturformeln.

Die Funktion hat genau eine Nullstelle, weil s streng monoton steigt und bijektiv ist, ich bin da ganz zuversichtlich, dass das dann auch konvergiert

Damit bekommst du dein t und kannst entsprechend einsetzen.

Aber geht das nicht auch noch irgendwie geschickter? Zum Beispiel indem du die Fläche von Dreiecken nimmst? Die passen auch viel besser in so ne Kurve rein …

titan99_

Also wenn ich die Bezier-Kurve 3-Dimensional berechne, sind die Viereckigen Flächen gekrümmt, kann aber die ja interpolieren.

Wenn ich das richtig verstehe, kann ich die Funktionsgleichnung y(x)=x mit dem Newtonverfahren nach x auflösen und so die Funktion invertieren?

titan99_

Kenner des Problems schrieb:

In der derivierten Kategorie von perversen C^2(R)-Modulgarben wird das Problem mit der Fouriertransformation trivial.

Kann man eine Fouriertransformation trivial invertieren???

.filmor

Ja, kann man (z.B. indem du noch dreimal fouriertransformierst)

Das Newtonverfahren berechnet dir erstmal nur eine Nullstelle einer Funktion (die natürlich noch ein paar Sachen erfüllen muss, aber die sind hier gegeben). Wenn du aber als die betrachtete Funktion f(x) - y nimmst, dann bekommst du natürlich für x = f^-1(y) eine Nullstelle, ist f invertierbar ist das sogar die einzige. Wenn das Verfahren also konvergiert, bekommst du auch das richtige Inverse.

titan99_

Merci für den Tipp, es dauert vermutlich eine Weile bis ich das mit dem Newtonverfahren begreife.

titan99_

Falls noch jemand Zeit hat.

Das Kurvenintegral ist das Integral vom Betrag der Ableitung des Vektors in dem die Funktion steht (x(t), y(t) und z(t) sind die Kompenenten). Bei der natürlichen Parametrisierung (Oder Parametrisierung nach Bogenlänge) http://www.iadm.uni-stuttgart.de/LstAnaMPhy/Weidl/analysis/vorlesung-analysis/node236.html steht in der letzten Zeile dass die Ableitung von r dem Tangentenvektor http://de.wikipedia.org/wiki/Frenetsche_Formeln entspricht, wenn ich das richtig verstehe. Bloss ist der Tangentenvektor R^x und r sollte eigentlich R^1 sein, wie wird dann der Tangentenvektor zu einem Skalar? Sonst steht im zweitletzten Term, dass r S^-1 entspricht, und das ist ja tatsächlich das inverse Kurvenintegral.

Eigentlich soll das Programm aus einem von Zufallspunkten gesuchtem konkaven Polygon, entlang der Polygonlinien eine Fläche bilden, die Ecken sollen aber abgerundet sein. Ich habe es schon mit kubischen Splines versucht, hatte dann aber bei spitzen Winkeln schlaufen und entschied mich dann für Bézier. Vielleicht hat ja sonst noch jemand eine Idee wie man das tun könnte.

Habe Link zu *.exe entfernt und Links zu Screenshots geändert