Also für Matalab gibt es da eine vorgefertigte Funktion:
https://www.mathworks.com/help/signal/ref/pmtm.html

Hi,

ich versuche von einem Signal das Hintergrundrauschen zu fitten. Dabei handelt es sich um rotes rauschen. Diesen Fit benötige ich um dann confidenc levels zu erzeugen.

Für den Fit kann amn Robust autoregessive fitting methode (AR.1) verwenden aber isch weiß nicht wie!

Es kommt noch dazu das ...

Hi zusammen,
ich hab ein Zeitsignal was keine dieskreten abstände hab und möchte und möchte eine Frequenz analyse mit multi-taper methode durchführen (mit 4 tapern). Ist das in Python möglich. Also die Lücken können auch mal größer sein. Ich hab es bis jetzt immer mit lombscargle gemacht
Hier ein ...

Hi Tobis,
danke! Was hälst du von dieser Lösung:

from scipy import stats
import numpy as np

arr = Ca_Fe

alpha = 0.05 # significance level = 5%
df = 8 # degress of freedom = 8
t = stats.t.ppf(1 - alpha/2, df) # t-critical value for 95%
s = np.std(arr, ddof=1) # sample standard deviation
n = len ...

Hallo,
ich war lange nict mehr hier. Ich hoffe ihr könnt mir weiter helfen.

Ich habe ein Signal welches ich mit

f, spec = LombScargle(depth,Ca).autopower(minimum_frequency=min_freq, maximum_frequency=max_freq)


in sein Frequenzspektrum überführt habe. Die Abstände der Daten sind irregulär. Nun ...

Hallo Forum,
Ich habe Probleme mit if statements...
ich habe Koordinaten
x1 y1 x2 y2
2 3 4 5
6 4 3 2

Und jetz will ich schauen ob:
y1 > y2
wenn ja ist x1, y1 der Startpunkt und x2,y2 der Endpunkt
und wenn nein dann ist x2 ,y2 der Startpunkt und x1, x2 der Endpunkt.

Kann mir jemand Helfen??

Grüße

Hi, ich hab eine Frage...
Und zwar habe ich ein Array mit x,y Koordinaten. Diese stell ich dar. Jetzt zeichne ich zwei schräge parralel zuainander liegende Geraden. Und jetzt will ich alle Punkte ermitteln die zwischen den beiden Geraden liegen. In folgendem Beispiel sollten die Outlier nicht ...

Hi,
ich hab mal ne Frage. Also ich hab eine Tabelle in der ich vorne Jahreszahlen stehen hab also von 1 bis 1000 in 1 Jahr Abstand und da hinter einen Wert.
Jetzt würde ich gerne FFT durchfüren um zu schauen irgendwelche Sachen wiederholen.
Das würde ich wie folgt machen:

N1 = 600
# sample spacing ...

He Ho,

ich hab ein Problem...und zwar will ich bei Matplotlib eine zweite x-achse einfügen die aber in der x-Richtung leicht verschoben ist. Also zweite achse ist kein Problem die in y-Achse mit
xaxis3.spines["top"].set_position(("axes", +1.05))
zu verschieben auch nicht aber sie in x-Richtung zu ...

Also weiß nicht genau. Eigentlich nicht das schränkt mich wieder ein oder. Ich will es in eine csv Datei schreiben und von Excel format weg kommen.
Ist das Excel-Datum: datetime.datetime(2016, 2, 26, 14, 3, 16)

Also Daten Einlesen transponieren und Datum umschreiben...

Also jetzt hab ichs noch mal geändert. Und zwar will ich es ja die Daten nur transponieren und das Datum in einem anderen Format ausgeben.
Also sieht jetzt folgender maßen aus:

from openpyxl import load_workbook
import numpy as np

wb = load_workbook(filename = 'Moxa_raman_8_05_1520160226140316 ...

Ok, sehr gut. Jetzt hab ich alles eingeladen. Jetzt möchte ich es noch richtig in einer csv Datei ausgeben.
Also sieht jetzt folgender maßen aus:
from openpyxl import load_workbook
import numpy as np

wb = load_workbook(filename = 'Moxa_raman_8_05_1520160226140316.xlsx')

ws = wb.active

ws = wb ...

Hi,

kann mir jemand helfen ich versuche ein xlsx-sheet einzulesen. und zwar hat dieses folgende Form:

Trace 0
Date 26. Feb. 16
Time 14:03:16
Points distance temperature
m °C
0 0,00 17,52
1 0,50 17,44

Dieses Sheet würde ich mir dann gerne als csv Datei ausgeben lassen und sollte dann ...

O.K. bin jetzt schon ein stück weiter ich hab nur noch ein problem mit den Leeren zeilen

import numpy as np


input = 'Dip_Cor_Moxa_raman_8_05_1520160122122520.csv'

output = 'Dip_Cor_Moxa_raman_8_05_1520160122122520_korr.csv'


with open (input, 'r') as in_file,open (output, 'w ...

Ja und vielleicht sollte ich das mit dem Befehl if machen also so was hier inder Arte:

file = open('*.csv','r')

for line in file:
if line < 112: continue #ersten 112 Zeilen werden weg gelassen
if line == "":
break #bei leerer zeile aufhören

Jetzt könnte ich ihn noch sagen beginn bei der ...

Also so bekomm ich schon mal Werte:
with open('Dip_Cor_Moxa_raman_8_05_1520160131000534.csv','r') as temp_file:

text = temp_file.read()

#trace = text.split('[Trace.*')

sec=re.findall(r'Time.Second.*',text,re.MULTILINE)
minu=re.findall(r'Time.Minute.*',text,re.MULTILINE)
hour=re.findall ...

Ich würde halt alles gern in Blöcke teilen nach dem "[Trace." zum beispiel und jeden block einzeln nach minuten Sekunden und so weiter durchsuchen. Zum Beispiel:
with open('Dip_Cor_Moxa_raman_8_05_1520160131000534.csv','r') as temp_file:

text = temp_file.read()

trace = text.split('[Trace ...

Hi,
und zwar hab ich eine sehr lange txt Datei die folgender maßen aufgebaut ist:

[Trace.8389]
0;0;26.092227935791;-6.77717351913452;-0.122081476822495
1;0.5;25.7060890197754;-6.78669500350952;-0.12099138610065
2;1;25.8580360412598;-6.78311920166016;-0.12584039842089
3;1.5;26.1734962463379;-6 ...

Hi,
und zwar habe ich ein Problem mit dem gridden von xyz-Daten. Ich habe einen sehr große txt-Datei mit xyz Werten. Wenn ich nur einen kleinen Ausschnitt der Datei gridde geht das Problem los aber bei allen Daten bekomme ich eine Fehlermeldung. Mein Code:
from matplotlib import cm
import numpy as ...

Hallo,
Ich versuche das gmtpy mit Anaconda zu installieren. Nach folgender Anleitung:
http://www.mondooceano.it/tutorial/gmtpy.html

Ich habe GMT installiert und gmtpy mit

python setup.py install

Und habe in der Setup.py die GMT Versionsnummer geändert.
Jetzt weiß ich nicht wie ich GMT und gmtpy ...