import os
import numpy as np
from scipy.fftpack import fft,ifft
import matplotlib.pyplot as plt
import xlwt
def fft0_7(x,y):#傅里叶变换
#采样点选择1400个,因为设置的信号频率分量最高为600Hz,根据采样定理知采样频率要大于信号频率2倍,所以这里设置采样频率为1400Hz(即一秒内有1400个采样点)
# x=np.linspace(0,1,401)
#
# #设置需要采样的信号,频率分量有180,390和600
# y=7*np.sin(2*np.pi*50*x) + 1.5*np.sin(2*np.pi*100*x)+5.1*np.sin(2*np.pi*150*x)+1
#
yy=fft(y) #快速傅里叶变换
# yreal = yy.real # 获取实数部分
# yimag = yy.imag # 获取虚数部分
yf=abs(yy) # 取模
yf1=yf/((len(x)/2)) #归一化处理
yf1[0]=yf1[0]/2
yf2 = yf1[range(int(len(x)/2))] #由于对称性,只取一半区间
xf = np.arange(<
傅里叶变换并把结果保存在excel文件中
最新推荐文章于 2023-09-03 23:21:07 发布