Code: Alles auswählen
a_fc = plt.subplot(gs[0,0])
for row in range(30):
f3_1, Pxx_spec = signal.welch(data_FBGS2[:,row], fs, 'flattop', 4096,
scaling='spectrum')
f_new= 1/(f3_1[1:]*3600)
Pxx = Pxx_spec[1:]
a_fc.loglog(f_new, np.sqrt(Pxx),'k-',zorder=3, )
import numpy.polynomial.polynomial as poly
#Regression
xmax,xmin = 23.92,0.04 #Der Bereich in dem Interpoliert werden soll. Wo gebe ich mein Minnimu umd Maximum ein?
logx,logy = np.log(f_new,np.log(Pxx)
coeffs = np.polyfit(logx,logy,deg=1)
poly = np.poly1d(coeffs)
yfit = lambda f_new: np.exp(poly(np.log(f_new)))
plt.loglog(f_new[0:],yfit(f_new),'m-')
Und gibt es einen fit der mir eine Potenzfunktion also sowas wie f^(5/3) ausgiebt bis jetz wird das ja über Polynome gelöst