meta data for this page
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
a_mechanik:kippender_besenstiel:gruppenseiten:gruppe341:start [ 8 January 2021 15:18] – [Dokumentation von Fehlern] alexandersteding | a_mechanik:kippender_besenstiel:gruppenseiten:gruppe341:start [22 January 2021 13:18] (current) – [Computerprogramm] alexandersteding | ||
---|---|---|---|
Line 11: | Line 11: | ||
{{: | {{: | ||
- | Gemessen wurde die Länge der Seite h. Die Länge der Seite L ist hier jeweils die Länge des Besenstieles. Über die Cosinusbeziehung | + | Gemessen wurde die Länge der Seite h. Die Länge der Seite L ist hier jeweils die Länge des Besenstiels. Über die Cosinusbeziehung |
\begin{equation} | \begin{equation} | ||
Cos=\frac{Ankathete}{Hypotenuse} | Cos=\frac{Ankathete}{Hypotenuse} | ||
Line 19: | Line 19: | ||
Cos=\frac{h}{l} | Cos=\frac{h}{l} | ||
\end{equation} | \end{equation} | ||
- | kann der Winkel α berechnet werden. Gemäß dem Wechselwinkelsatzt | + | kann der Winkel α berechnet werden. Gemäß dem Wechselwinkelsatz |
===== Versuchsdurchführung ===== | ===== Versuchsdurchführung ===== | ||
An dieser Stelle soll kurz auf die Versuchsdurchführung eingegangen werden. Zunächst wurde der Versuch ohne Berücksichtigung der Luftreibung durchgeführt. Dazu wurden zuerst die beiden Stäbe mittels einem Maßband vermessen.{{: | An dieser Stelle soll kurz auf die Versuchsdurchführung eingegangen werden. Zunächst wurde der Versuch ohne Berücksichtigung der Luftreibung durchgeführt. Dazu wurden zuerst die beiden Stäbe mittels einem Maßband vermessen.{{: | ||
- | Der Messfehler der Länge wurde mit 50% der kleinsten messbaren Größe, in diesem Fall 50% von 1mm ensprächend 0,5mm angenommen. Um den Nachbarschaftsfrieden nicht zusehr | + | Der Messfehler der Länge wurde mit 50% der kleinsten messbaren Größe, in diesem Fall 50% von 1mm ensprächend 0,5mm angenommen. Um den Nachbarschaftsfrieden nicht zu sehr überzustrapazieren, wurde der Küchenboden mit einigen Handtüchern abgefedert. Daraus resultierende Fehler wurden in der Auswertung betrachtet. Wie oben beschrieben wurde im zweiten Schritt der Winkel durch Messen |
{{: | {{: | ||
Line 30: | Line 30: | ||
Im nächsten Schritt wurde mit einem Geräusch die akustische Stoppuhr in Phyphox gestartet während gleichzeitig der Besen fallen gelassen wurde. Dabei wurde darauf geachtet, dass dieser fallengelassen wird ohne diesem Anschwung oder sonst eine Beschleunigung zusätzlich zugeben. Die akustische Stoppuhr stoppt die Messung sobald der Besen den Boden ereicht. Das folgende Video zeigt den Versuchsablauf: | Im nächsten Schritt wurde mit einem Geräusch die akustische Stoppuhr in Phyphox gestartet während gleichzeitig der Besen fallen gelassen wurde. Dabei wurde darauf geachtet, dass dieser fallengelassen wird ohne diesem Anschwung oder sonst eine Beschleunigung zusätzlich zugeben. Die akustische Stoppuhr stoppt die Messung sobald der Besen den Boden ereicht. Das folgende Video zeigt den Versuchsablauf: | ||
- | {{ : | + | {{ : |
Auf diese Weise wurden zunächst für beide Stiele für insgesammt fünf Winkel φ (5°, 20°, 35°, 50°, 65°) jeweils fünf mal die Fallzeit gemessen. | Auf diese Weise wurden zunächst für beide Stiele für insgesammt fünf Winkel φ (5°, 20°, 35°, 50°, 65°) jeweils fünf mal die Fallzeit gemessen. | ||
==== Messung des Luftwiederstandes ==== | ==== Messung des Luftwiederstandes ==== | ||
- | Um den Einfluss der Reibung durch Luftwiederstand | + | Um den Einfluss der Reibung durch Luftwiderstand |
{{: | {{: | ||
Line 53: | Line 53: | ||
===== Messung der Fallzeit ===== | ===== Messung der Fallzeit ===== | ||
- | Hier wurden die gemessenen Fallzeiten für den Stiel 1 mit Länge 1,33 m ± 0,5 mm und Stiel 2 mit 74,9 m ± 0,5 mm aufgetragen. Weiterhin wurden die gemessenen Längen von h zur Berechnung des Winkels, sowie der Mittelwert und berechneter | + | Hier wurden die gemessenen Fallzeiten für den Stiel 1/ Stiel mit Pappe mit Länge 1,33 m ± 0,5 mm und Stiel 2 mit 74,9 m ± 0,5 mm aufgetragen. Weiterhin wurden die gemessenen Längen von h zur Berechnung des Winkels, sowie der Mittelwert und berechneter |
=== 5° === | === 5° === | ||
Line 68: | Line 68: | ||
| ||| | | | ||| | | ||
^ Mittelwert | ^ Mittelwert | ||
- | ^ Standartfehler | + | ^ Standardfehler |
| |||| | | |||| | ||
^ Länge h [cm] | 132,30 ±0, | ^ Länge h [cm] | 132,30 ±0, | ||
Line 87: | Line 87: | ||
| ||| | | | ||| | | ||
^ Mittelwert | ^ Mittelwert | ||
- | ^ Standartfehler | + | ^ Standardfehler |
| |||| | | |||| | ||
^ Länge h [cm] | 124,90\\ ±0, | ^ Länge h [cm] | 124,90\\ ±0, | ||
Line 104: | Line 104: | ||
| | | | ||
^ Mittelwert | ^ Mittelwert | ||
- | ^ Standartfehler | + | ^ Standardfehler |
| |||| | | |||| | ||
^ Länge h [cm] | 108, | ^ Länge h [cm] | 108, | ||
Line 122: | Line 122: | ||
| ||| | | | ||| | | ||
^ Mittelwert | ^ Mittelwert | ||
- | ^ Standartfehler | + | ^ Standardfehler |
| |||| | | |||| | ||
^ Länge h [cm] | 85, | ^ Länge h [cm] | 85, | ||
Line 139: | Line 139: | ||
| | | | ||
^ Mittelwert | ^ Mittelwert | ||
- | ^ Standartfehler | + | ^ Standardfehler |
| |||| | | |||| | ||
^ Länge h [cm] | 56, | ^ Länge h [cm] | 56, | ||
+ | === Unsicherheit der Winkel === | ||
+ | Für die berechneten Winkel wurden die Fehler gemäß Auswertung berechnet: | ||
+ | ^ ^ Stab 1 und Stab mit Pappe ^ | ||
+ | | Winkel in rad | Fehler | ||
+ | | 0,09 | 0,04 | | ||
+ | | 0,347 | 0,011 | | ||
+ | | 0,609 | 0,005 | | ||
+ | | 0,871 | 0,003 | | ||
+ | | 1, | ||
+ | |||
+ | ^ ^ Stab 2 ^ | ||
+ | | Winkel in rad | Fehler | ||
+ | | 0,085 | 0,011 | | ||
+ | | 0, | ||
+ | | 0, | ||
+ | | 0, | ||
+ | | 1, | ||
====== Computerprogramm ====== | ====== Computerprogramm ====== | ||
- | Für die numerische Berechnung der Bewegunsgleichung wurde der folgende Code in Python geschrieben. | + | Für die numerische Berechnung der Bewegunsgleichung wurde der folgende Code in Python geschrieben. In der überarbeiteten Fassung berechnet er für alle Winkel im Intervall [0,05;1,5] Radiant die Fallzeiten. Die Daten werden automatisch als Plot und formatierte Textdatei zur weiteren Verwendung(z.B. QTIPlot) gespeichert. |
<code python [enable_line_numbers=" | <code python [enable_line_numbers=" | ||
import math | import math | ||
- | # Voreinstellungen | + | import numpy as np |
- | Anfangsbeschleunigung= 9.81 # | + | import matplotlib.pyplot as plt |
- | Anfangsgeschwindigkeit = 0 # Geschwindigkeit zum Zeitpunkt t=0 in m/s | + | ############## |
- | Meta=[] # Liste in der die Berechneten Fallzeiten gespeichert werden | + | Erdbeschleunigung= 9.81 # Erdbeschleunigung |
+ | Meta=[] # Liste in der die Berechneten Fallzeiten gespeichert werden | ||
+ | winkelliste=[] | ||
+ | Berechnungswinkel=0.05 | ||
+ | while Berechnungswinkel <= 1.51: | ||
+ | winkelliste.append(format(Berechnungswinkel,' | ||
+ | Berechnungswinkel += 0.01 | ||
+ | Winkelliste=list(map(float, | ||
- | # | + | ############## |
print(" | print(" | ||
Stablänge = float(input(" | Stablänge = float(input(" | ||
print(" | print(" | ||
Zeitschritt =float(input(" | Zeitschritt =float(input(" | ||
- | Tau= math.sqrt((2*Stablänge)/ | + | ############### |
+ | Tau= math.sqrt((2*Stablänge)/ | ||
+ | for Winkel in Winkelliste: | ||
+ | listewes=[] | ||
+ | Momentanwinkel=Winkel | ||
+ | Momentangeschwindigkeit=0 | ||
+ | while Momentanwinkel < | ||
+ | Beschleunigung= math.sin(Momentanwinkel)/ | ||
+ | Geschwindigkeit= Momentangeschwindigkeit + Zeitschritt*Beschleunigung | ||
+ | neuerwinkel= Momentanwinkel + Zeitschritt*Geschwindigkeit | ||
+ | listewes.append(neuerwinkel) | ||
+ | Momentanwinkel=neuerwinkel | ||
+ | Momentangeschwindigkeit=Geschwindigkeit | ||
- | print(" | + | k=format((len(listewes))*Zeitschritt, |
- | wahl= input(" | + | |
- | wahlf= float(wahl) | + | |
- | if wahlf< | + | ############### |
- | print(" | + | data = np.array([np.array(Winkelliste).astype(np.float) , np.array(Meta).astype(np.float) ]).T #Umwandeln der Winkel und Fallzeiten in ein Array und Transponieren der Daten |
- | Winkeleins | + | datafile_path |
- | | + | with open(datafile_path, |
- | else : | + | |
- | print("Bitte Winkel in Radiant eingeben" | + | |
- | | + | ############### |
- | | + | xstart= 0 |
- | | + | xend= 1.6 |
- | Winkelvier | + | xticks= 0.05 |
- | Winkelfünf = input(" | + | ystart= 0 |
- | Anfangswinkellistenamen=[Winkeleins,Winkelzwei, Winkeldrei, Winkelvier, Winkelfünf] | + | yend= 1.5 |
- | + | plt.figure(figsize=(20,10), dpi=100) | |
- | Anfangswinkelliste= list(map(float, | + | plt.grid(True, |
- | + | plt.axis([xstart,xend, | |
- | #Berechnungen | + | plt.title(' |
- | for Anfangswinkel in Anfangswinkelliste: | + | plt.xticks(np.arange(xstart, |
- | Startbeschleuni= math.sin(Anfangswinkel)/ | + | a=plt.scatter(Winkelliste, |
- | | + | plt.legend([a], |
- | | + | plt.xlabel(' |
- | | + | plt.ylabel(' |
- | gest= Startgeschwendi | + | plt.savefig('graph'+'_'+str(Stablänge)+'_'+str(Zeitschritt)+' |
- | | + | plt.close |
- | | + | |
- | listewes=[] | + | |
- | b=Zeitschritt | + | |
- | + | ||
- | while Winkelt | + | |
- | | + | |
- | listebes.append(bes) | + | |
- | ges= gest + b*bes # | + | |
- | listeges.append(ges) | + | |
- | wes= Winkelt | + | |
- | | + | |
- | | + | |
- | gest=ges | + | |
- | + | ||
- | k=format((len(listewes))*Zeitschritt, '.3g') # Berechnung der Fallzeit t und Vermeidung des Floating Point Errors durch Signifikante Stellen | + | |
- | Meta.append(k) #Einhängen der Fallzeiten in die Liste Meta | + | |
- | + | ||
- | for item, | + | |
- | print('Die Fallzeit für den Startwinkel',nextitem,'Radiant beträgt', item, 'Sekunden') | + | |
</ | </ | ||
===== Überprüfung des Code===== | ===== Überprüfung des Code===== | ||
Um den Code zu verifizieren wird die Fallzeit bei einer Stablänge von 1.45 m bei einem Startwinkel von 0.25 rad und einem Zeitschritt von 0.01 s berechnet. Diese beträgt 0.8 Sekunden. Der Vergleich mit der Grafik aus der Aufgabenstellung bestätigt die Berechnung. | Um den Code zu verifizieren wird die Fallzeit bei einer Stablänge von 1.45 m bei einem Startwinkel von 0.25 rad und einem Zeitschritt von 0.01 s berechnet. Diese beträgt 0.8 Sekunden. Der Vergleich mit der Grafik aus der Aufgabenstellung bestätigt die Berechnung. | ||
+ | |||
+ | {{ : | ||
+ | |||
+ | |||
====== Dokumentation von Fehlern ====== | ====== Dokumentation von Fehlern ====== | ||
- | Ein häufig aufgetauchter | + | Ein häufig aufgetauchter |
{{ : | {{ : | ||
- | Die dadurch zusätzlich entstehenden Schallwellen | + | Die dadurch zusätzlich entstehenden Schallwellen |
- | ===== Syntax und Funktionen im Wiki ===== | + | |
- | Hier noch Links zu | + | |
- | * den [[doku> | + | |
- | * [[: | + | |
- | * [[: | + | |