Function 求变迹函数

Function 求变迹函数,function,signals,Function,Signals,你好,我想在我的图上画一个切趾。但它不起作用。由于我不太擅长编程,我希望有人能帮助我找到正确的切趾函数。我读过,Welch函数应该在这里工作。我只是无法在代码中实现它 import numpy as np from scipy.fftpack import fft,fftfreq,ifft import matplotlib.pyplot as plt from scipy.signal import find_peaks from scipy import signal dataset =

你好,我想在我的图上画一个切趾。但它不起作用。由于我不太擅长编程,我希望有人能帮助我找到正确的切趾函数。我读过,Welch函数应该在这里工作。我只是无法在代码中实现它

import numpy as np
from scipy.fftpack import fft,fftfreq,ifft
import matplotlib.pyplot as plt
from scipy.signal import find_peaks
from scipy import signal

dataset = np.genfromtxt(fname='Halogenlampe1.txt',skip_header=1)
# x
t = dataset[:,0]
# y
sig = dataset[:,1]
npts=len(t)

iFFT = abs(ifft(sig))
freqs = fftfreq(npts, t[1]-t[0])

plt.subplot(211)
plt.ylabel('Intensität')
plt.xlabel('Spannung [V]')
plt.plot(t[:npts], sig[:npts])
plt.subplot(212)
plt.plot(freqs,iFFT)
f, Pxx_den = signal.welch(t, sig)#, nperseg=1024)
plt.semilogy(f, Pxx_den)
plt.xlabel('Frequenz [1/V]')
plt.ylabel('Amplitude')
plt.xlim([0,5])
plt.ylim([0,0.75])
peaks,_= find_peaks(iFFT,height=0.5)
plt.plot(freqs[peaks],iFFT[peaks], 'x')
#plt.xticks(np.arange(0, 5.5, step=0.5))
plt.grid(True)
plt.savefig('Kalilbrierug_Laser', bbox_inches='tight', dpi=400)
plt.show()enter code here