Mathesis Wiki

Benutzer-Werkzeuge

Webseiten-Werkzeuge

Sie befinden sich hier: start » ws1920 » arbeitstermin_06
Zuletzt angesehen: • astrophysik • zwischenziele • partikelsimulation • n-body-simulation • protokolle • 29.05.2019 • protokolle • schrifterkennung • boids • hier_geht_es_zu_unserem_code_mit_convolutional_layern
ws1920:arbeitstermin_06

Unterschiede

Hier werden die Unterschiede zwischen zwei Versionen gezeigt.

Link zu dieser Vergleichsansicht

Beide Seiten der vorigen Revision Vorhergehende Überarbeitung
Nächste Überarbeitung 		Vorhergehende Überarbeitung
ws1920:arbeitstermin_06 [2020/02/09 20:55]
Zetraeder 		ws1920:arbeitstermin_06 [2020/03/25 16:28] (aktuell)
Zetraeder
Zeile 4: Zeile 4:
 **16.01.2020** **16.01.2020**
  
-[[ws1920:update_00|<<]] | [[ws1920:​scutoids|Home]] | [[ws1920:​update_00|>>​]]+[[ws1920:arbeitstermin_05|<<]] | [[ws1920:​scutoids|Home]] | [[ws1920:​cd_06|Code]] | [[ws1920:​update_00|>>​]]
  
-Wir haben uns heute darum gekümmert, dass die verschiedene Fehler noch etwas seltener auftauchen als zuvor (anpassung des Epsilonbereiches). Dann haben wir uns mit den Möglichkeiten auseinandergesetzt,​ dass tatsächliche Voronoi zu generieren. 
- 
-{{:​ws1920:​rechtewinkelbild.jpg?​400|}} ​ 
  
 Abb. Generationsstand 17:02 Uhr Abb. Generationsstand 17:02 Uhr
  
-Außerdem haben wir eine Methode geschrieben,​ die den Winkel zwischen zwei verschiedenen Vektoren ausgibt. So können wir herausfinden,​ welche beiden Linien links und rechts von einer Verbindungslinie sind (wichtig für die Auswahl der Voronoi-Linien,​ bei denen die Schnittpunkte gefunden werden sollen).+{{:​ws1920:​rechtewinkelbild.jpg?175 |}} 
  
-<code python>+Wir haben uns heute darum gekümmert, dass die verschiedene Fehler noch etwas seltener auftauchen als zuvor (anpassung des Epsilonbereiches). Dann haben wir uns mit den Möglichkeiten auseinandergesetzt,​ dass tatsächliche Voronoi zu generieren.
  
-# -*- coding: utf-8 -*- +Außerdem haben wir eine Methode geschrieben, die den Winkel zwischen zwei verschiedenen ​Vektoren ​ausgibt. So können wir herausfinden,​ welche beiden Linien links und rechts von einer Verbindungslinie sind (wichtig für die Auswahl ​der Voronoi-Linien, ​bei denen die Schnittpunkte gefunden ​werden sollen).
-"""​ +
-Created on Thu Dec 19 15:31:18 2019 +
- +
-@author: Lukas +
-"""​ +
- +
-import numpy as np +
-import math +
-import turtle +
-import random as rn +
- +
-turtle.speed(0) +
-turtle.ht() +
- +
-#   ​VARIABLEN +
- +
-anzahlKerne = 15 +
-max_x = 300 +
-max_y = 300 +
- +
- +
- +
-#   ​KLASSEN +
- +
-class Gerade(object):​ +
-    ov = (0,0)    #​Ortsvektor +
-    rv = (0,0)    #​Richtungsvektor +
-    l  = 0;          #Laenge +
- +
-    kerne = (0,0) +
-     +
-    def toString():​ +
-        print("​OV : "​+Gerade.ov+"​ RV : "​+Gerade.rv+"​ Kerne "​+Gerade.kerne) +
- +
-#   ​METHODEN ​    +
-     +
-def generiereKernListe(n,​maxx,​maxy):​ +
-    #   ​generiert ​eine Liste von Punkten im R^2 +
-    liste=[] +
-     +
-    for c in range(n): +
-        x = rn.randint(-maxx,maxx)#​rn.uniform(0,​maxx) +
-        y = rn.randint(-maxy,​maxy)#​rn.uniform(0,​maxy) +
-         +
-        liste.append((x,​y)) +
-         +
-    return liste +
- +
-def zeichneLinie(start,​ende):​ +
-    #   ​zeichnet eine Linie zwischen start und ende mithilfe einer turtle +
-    turtle.penup() +
-    turtle.goto(start) +
-    turtle.pendown() +
-    turtle.goto(ende) +
-     +
-def richtungLinie(start,​ende):​ +
-    #   gibt den Richtungsvektor von start nach ende aus +
-    richtung ​ = np.subtract(ende,​start) +
-    return richtung +
- +
-def schnittTest(g1,​g2):​ +
-    #   ​testetm ob g1 und g2 sich schneiden +
-    return np.linalg.solve([[g2.rv[0],​-g1.rv[0]],​[g2.rv[1],​-g1.rv[1]]],​[[g1.ov[0]-g2.ov[0]],​[g1.ov[1]-g2.ov[1]]]) +
- +
-def vektorLänge(v2):​ +
-    #   ​errechnet ​die Laenge eines Vektors +
-    v1 = np.array([0,​0]) +
-    return math.sqrt((abs(v1[0]-v2[0])**2)+(abs(v1[1]-v2[1])**2)) +
- +
-def vektorNorm(v):​ +
-    #   ​normiert einen Vektor +
-    return(1/​math.sqrt(v[0]**2+v[1]**2)*v) +
-     +
-def punktAufLinie(p,​o,​r):​ +
-    #   liegt punkt p auf der geraden g = o + m*r?     +
-    if(np.array_equal(p,​o)):​ return False +
-    elif(np.array_equal(r,​[0,​0])):​ return False +
-    elif(r[1] == 0):  +
-        a = (p[0]-o[0])/​r[0] +
-        if (o[1]+a*r[1] == p[1]): +
-            if(a<1 and a>0): return True +
-            else: return False +
-        else: return False +
-    else:  +
-        a = (p[1]-o[1])/​r[1] +
-        if (o[0]+a*r[0] == p[0]): +
-            if(a<1 and a >0): return True +
-            else: return False +
-        else: return False +
- +
-def rechtW(v):​ +
-    #   ​findet einen Vektoren im rechten Winkel zu dem gegebenen Vektor +
-    v_neu = np.array((-v[1],​v[0])) +
-    return v_neu +
- +
-def zwischenWinkel(a,​b):​ +
-    #   ​gibt ​den Winkel zwischen zwei Vektoren ​aus +
-    w = (a[0]*b[0]+a[1]*b[1])/​(math.sqrt(a[0]**2+a[1]**2)*math.sqrt(b[0]**2+b[1]**2)) +
-    w = math.degrees(math.acos(w)) +
-    return w +
- +
-#   ​MAIN +
-     +
-#   ​genriert ​die Liste der Kerne +
-ker = generiereKernListe(anzahlKerne,​max_x,​max_y) +
- +
-print("​Kerne : ",​ker) +
- +
- +
-listeGeraden = [] +
- +
-#   ​generiert alle Linien ​zwischen den Kernen +
-for i in range(0,anzahlKerne):​ +
-    for k in range(i,​anzahlKerne):​ +
-        if (i!=k): +
-             +
-            eineGerade ​    = Gerade() +
-            eineGerade.ov ​ = np.array(ker[i]) +
-            eineGerade.rv ​ = np.subtract(ker[k],​ker[i]) +
-            eineGerade.l ​  = vektorLänge(eineGerade.rv) +
-             +
-            # welche Kerne verbindet ​die Linie? +
-            eineGerade.kerne = (i,k) +
-     +
-            listeGeraden.append(eineGerade) +
- +
- +
-#   ​zeichnet alle Linien zwischen allen Punkten +
-turtle.color("​cyan"​) +
-for i in range(0, len(listeGeraden)):​ +
-    zeichneLinie(listeGeraden[i].ov,​listeGeraden[i].ov+listeGeraden[i].rv) +
- +
-schnitt = 0 +
-l = len(listeGeraden) +
- +
-#testet, welche geraden gelöscht ​werden sollen +
- +
-i = 0 +
-k = 0 +
- +
-epsilon = 0.0001 #​*anzahlKerne +
- +
-while (i < l): +
-    k = i+1 +
-    while (k < l): +
-        z = [0,0] +
-        if(i!=k): +
-            #for k in range(i+1, l): +
-            try: +
-                #   ​schneiden sich die Linien? +
-                z = schnittTest(listeGeraden[i],​listeGeraden[k]) +
-                 +
-                if (z[0][0]<​1-epsilon and z[0][0]>​0+epsilon and z[1][0]<​1-epsilon and z[1][0]>​0+epsilon):​ +
-                    #​print("​\n\n ",​i,"​ OV ",​listeGeraden[i].ov," ",​i,"​ RV ",​listeGeraden[i].rv,"​ ",​k,"​ OV ",​listeGeraden[k].ov,"​ ",​k,"​ RV ",​listeGeraden[k].rv) +
-                    #​print("​vergleiche ",​i,"​ und ", k) +
-                    #​print("​Schnitt\n",​z[0][0],​z[1][0]) +
-                    schnitt=schnitt+1 +
-                    #del listeGeraden[i] ​  #​ACHTUNG:​ EIGENTLICH SOLL DIE LÄNGERE GELÖSCHT WERDEN +
-                     +
-                    #​print("​laenge a : ",​listeGeraden[i].l,"​ laenge b : ",​listeGeraden[k].l,"​\n",​10*"​ - ") +
-                    if (listeGeraden[i].l < listeGeraden[k].l):​ +
-                        del listeGeraden[k] +
-                        i=0 +
-                    else:  +
-                        del listeGeraden[i] +
-                        k=0 +
-            except: +
-                #​print("#​ end of this iteration"​) +
-                break +
-        k=k+1 +
-    i=i+1 +
-         +
-#for i in range(0, l): +
-#    for k in range(i+1, l): +
-        +
- +
- +
-listeRGeraden = [] +
- +
-for i in range(0, len(listeGeraden)):​ +
-    neueGerade = Gerade() +
-    neueGerade.ov = listeGeraden[i].ov+(.5*listeGeraden[i].rv) +
-    neueGerade.rv = rechtW(listeGeraden[i].rv) +
- +
-    listeRGeraden.append(neueGerade) +
- +
- +
- +
-print("​Schnitte : ", schnitt) +
- +
-for i in range(0, len(listeGeraden)):​ +
-    turtle.color("​black"​) +
-    turtle.pensize(2) +
-    zeichneLinie(listeGeraden[i].ov,​listeGeraden[i].ov+listeGeraden[i].rv) +
-    print(listeGeraden[i].kerne) +
-     +
-    turtle.color("​red"​) +
-    turtle.pensize(1) +
-    zeichneLinie(listeRGeraden[i].ov,​listeRGeraden[i].ov+listeRGeraden[i].rv) +
- +
- +
-turtle.Screen().exitonclick() +
-#​turtle.done() +
- +
- +
- +
-# Gute Kernlisten:​ +
- +
-# beste : [(52, 199), (101, 243), (-227, -48), (5, 213), (26, -294), (197, 174), (57, -275), (-61, -63), (206, 6), (-112, 62), (26, 105), (247, 121), (287, 161), (152, 272), (209, -94)] +
-# beste : [(-92, 125), (-172, -127), (282, 26), (43, 114), (52, -6), (-155, -212), (-187, 135), (125, -30), (-78, -80), (8, -165), (-292, -8), (140, -7), (-252, -254), (-47, -19), (292, -225)] +
- +
-# 15 : [(-214, -223), (65, 86), (113, -226), (227, -199), (8, -5), (94, 193), (-218, -280), (108, -30), (206, 217), (-21, -292), (215, 112), (-289, -86), (-127, 50), (-257, 298), (-194, 254)] +
-# 15 : [(-158, 54), (181, -2), (-55, -230), (158, 192), (106, -216), (150, -88), (156, -245), (-22, 261), (67, 2), (-47, -27), (81, -28), (81, -252), (290, -54), (282, 286), (-32, -188)] +
-# 15 : [(166, 27), (-152, 245), (-152, 79), (-64, 298), (259, -222), (98, 177), (38, 227), (-40, -201), (55, 68), (-291, 272), (60, -232), (143, -282), (-86, 49), (53, 270), (160, -177)] +
-# 25 : [(-98, 152), (8, 295), (-153, -130), (-272, 86), (-92, 57), (157, 255), (-121, 134), (-42, -17), (-57, -238), (-290, -26), (-233, 8), (-271, 160), (27, 43), (-32, 292), (108, -297), (-41, 207), (67, -111), (-209, -224), (172, -288), (255, -136), (105, -229), (166, 42), (-145, -288), (33, 193), (77, 78)] +
-#  7 : [(76, -55), (-297, -51), (162, -110), (-50, -9), (125, 18), (-21, 35), (-259, 175), (-194, 23), (-243, -147), (181, -288)] +
- +
-# 30 : [(-168, -165), (94, -217), (-10, 185), (-90, 224), (-140, 246), (189, 28), (297, -126), (-211, 290), (71, -6), (-298, -79), (-78, 51), (277, -97), (135, 56), (47, 107), (182, -138), (258, -241), (-292, 3), (257, -188), (227, -261), (90, -2), (-228, 19), (153, -44), (-169, -146), (225, 33), (78, 166), (142, 238), (0, -105), (217, 83), (211, -195), (-142, -79)] +
-# 30 : [(253, 63), (-202, 25), (117, -194), (271, 247), (-84, -96), (145, 215), (-81, -145), (157, -297), (-157, 80), (217, 164), (175, -253), (-190, -48), (-122, 193), (253, -268), (-211, -81), (-121, 64), (-249, -25), (244, 298), (-42, 98), (-144, 218), (42, -56), (-141, 118), (125, 63), (-119, -21), (-79, 238), (21, 272), (100, -107), (147, 25), (-170, -54), (260, 267)] +
- +
-</​code>​+
  
 ---- ----
  
ws1920/arbeitstermin_06.1581278110.txt.gz · Zuletzt geändert: 2020/02/09 20:55 von Zetraeder

Seiten-Werkzeuge

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki