帮我写一段matlab代码 实现对1khz频率正弦波进行采样率1mhz采样点数1024的快速傅里叶变换并绘图

以下是实现对1kHz频率正弦波进行1MHz采样率、1024采样点数的快速傅里叶变换并绘图的MATLAB代码：

% 定义采样率和采样点数
Fs = 1e6; % 采样率为1MHz
N = 1024; % 采样点数为1024

% 生成1kHz频率正弦波
f = 1000; % 正弦波频率为1kHz
t = (0:N-1) / Fs; % 生成时间序列
x = sin(2*pi*f*t); % 生成正弦波信号

% 进行快速傅里叶变换
X = fft(x);

% 计算频率向量
fvec = (0:N-1) * Fs / N;

% 绘制频谱图
plot(fvec/1000, 20*log10(abs(X)/N)); % 频率轴以kHz为单位，以dB为单位绘制幅度谱
xlabel('频率 (kHz)');
ylabel('幅度 (dB)');
title('1kHz正弦波的频谱图');

这段代码首先定义了采样率和采样点数，然后生成了1kHz频率的正弦波信号。接着，使用MATLAB自带的快速傅里叶变换函数fft()对正弦波信号进行了快速傅里叶变换。然后计算了频率向量，用于绘制幅度谱的横坐标。最后，使用MATLAB的plot()函数绘制了1kHz正弦波的频谱图。