[已解决]添加代码，对代码中的信号y添加合适的窗函数，进行信号修正

臧雁德

添加代码，使代码中的信号y添加合适的窗函数，进行信号修正。

1. import numpy as np
   
2. import matplotlib.pyplot as plt
   
3. from matplotlib.pylab import mpl
   
4. 
   
5. mpl.rcParams['font.sans-serif'] = ['SimHei']  # 显示中文
   
6. mpl.rcParams['axes.unicode_minus'] = False  # 显示负号
   
7. 
   
8. # 生成正弦信号 #
   
9. fs = 1000  # 定义采样率
   
10. t = np.linspace(0,1,fs)  # 定义时间序列
    
11. y = np.sin(2*np.pi*20*t) + 2*np.sin(2*np.pi*40*t)  # 定义正弦信号
    
12. 
    
13. # 进行FFT分析 #
    
14. Y = np.fft.fft(y)  # FFT变换
    
15. normalization_y = 2*np.abs(Y)/fs  # 归一化,能量集中处理
    
16. freqs = np.fft.fftfreq(len(y), 1/fs)  # 自动生成变换后的频率序列
    
17. 
    
18. # 绘制输入波形与频谱 #
    
19. fig,axes = plt.subplots()  # 添加画布
    
20. 
    
21. plt.subplot(211)
    
22. plt.plot(t[0:100],y[0:100])  # 取前一百采样频率点
    
23. plt.title('输入波形(部分)')
    
24. plt.xlabel('时间(s)')
    
25. plt.ylabel('振幅')
    
26. 
    
27. plt.subplot(212)
    
28. plt.plot(freqs[range(int(fs/2))], normalization_y[range(int(fs/2))])  # 显示单边频率
    
29. plt.title('单边频谱')
    
30. plt.xlabel('频率 (Hz)')
    
31. plt.ylabel('振幅')
    
32. plt.xlim(0,100)
    
33. 
    
34. fig.tight_layout()  #自动调整布局，使标题不干涉
    
35. plt.show()

复制代码

最佳答案

月排行榜 / 总排行榜

Richard0824

2023-10-13 20:17:13

window = np.array([0.5 - 0.5 * np.cos(2 * np.pi * n / (1000 - 1)) for n in range(1000)])
y = window * y  666

跳转到最佳答案楼层

Richard0824

window = np.array([0.5 - 0.5 * np.cos(2 * np.pi * n / (1000 - 1)) for n in range(1000)])
y = window * y  666

臧雁德

window = np.array([0.5 - 0.5 * np.cos(2 * np.pi * n / (1000 - 1)) for n in range(1000)])
y = window * y