Broken-Powelaw-Plot
Verfasst: Donnerstag 5. Dezember 2013, 16:17
Hallo zusammen,
ich habe ein Problem mit dem Plot einer zusammengesetzten Funktion. Genauer: Mit einer "Broken-Powerlaw", also einem Potenzgesetz, welches ab einem
bestimmten x die Potenz ändert. Einzelt lassen sich die Funktionen wunderbar plotten, nur zusammen wird mir keine Funktion angezeigt - einen Fehler spuckt Python jedoch auch nicht aus. Ich glaube das der Fehler in der Schleife liegt, dass das so einfach nicht funktioniert. Mit fällt aber auch keine andere Lösung ein und im Internet habe ich bis jetzt auch nicht gefunden.
Vielen Dank für die Hilfe!
ich habe ein Problem mit dem Plot einer zusammengesetzten Funktion. Genauer: Mit einer "Broken-Powerlaw", also einem Potenzgesetz, welches ab einem
bestimmten x die Potenz ändert. Einzelt lassen sich die Funktionen wunderbar plotten, nur zusammen wird mir keine Funktion angezeigt - einen Fehler spuckt Python jedoch auch nicht aus. Ich glaube das der Fehler in der Schleife liegt, dass das so einfach nicht funktioniert. Mit fällt aber auch keine andere Lösung ein und im Internet habe ich bis jetzt auch nicht gefunden.
Vielen Dank für die Hilfe!
Code: Alles auswählen
import matplotlib as plt
import matplotlib.pyplot as pyplot
from pylab import *
import math
import numpy as np
a = -1.74956
b = -2.24956
br = 0.800358
norm = 1.77567*(10**(-4))
for x in np.arange(0.0001, 10):
if x < br:
y = norm*(x**a)
else:
y = norm*(br**(a-b))*(x**b)
xlabel("Energy [keV]")
pyplot.xscale('log')
ylabel("F [mJy]")
pyplot.yscale('log')
matplotlib.pyplot.autoscale(enable=True, axis='both', tight=None)
plot(x,y)
show()
