如何使用scipy绘制python中的fft，并修改频率范围，使其显示中心的两个峰值频率？

import numpy as np
import matplotlib.pyplot as plt
from scipy.fft import fft, fftfreq

# 创建时间序列
t = np.linspace(0, 1, 1000, endpoint=False)
# 创建输入信号，例如正弦波
input_signal = np.sin(2*np.pi*50*t) + np.sin(2*np.pi*120*t)

# 对输入信号进行FFT变换
fft_output = fft(input_signal)

# 计算频率轴
freq_axis = fftfreq(len(t), 1/len(t))

# 绘制FFT结果
plt.plot(freq_axis, np.abs(fft_output))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude')
plt.grid(True)
plt.show()

# 修改频率范围
plt.xlim(0, 200)  # 设置频率范围为0到200Hz

import numpy as np
import matplotlib.pyplot as plt
from scipy.fft import fft, fftfreq

# 创建时间序列
t = np.linspace(0, 1, 1000, endpoint=False)

# 创建输入信号，例如正弦波
input_signal = np.sin(2*np.pi*50*t) + np.sin(2*np.pi*120*t)

# 对输入信号进行FFT变换
fft_output = fft(input_signal)

# 计算频率轴
freq_axis = fftfreq(len(t), 1/len(t))

# 绘制FFT结果
plt.plot(freq_axis, np.abs(fft_output))
plt.xlabel('Frequency (Hz)')
plt.ylabel('Amplitude')
plt.grid(True)

# 修改频率范围
plt.xlim(0, 200)  # 设置频率范围为0到200Hz

plt.show()