对信号stft变换，绘制出stft频谱图，谱峰图，用sinusoid重建信号以后的图MATLAB

以下是一个简单的例子，展示如何进行STFT变换、绘制STFT频谱图和谱峰图、以及如何使用sinusoid重建信号。在这个例子中，我们使用MATLAB的信号处理工具箱来实现。

% 生成一个测试信号
fs = 1000;
t = 0:1/fs:1-1/fs;
f = [50 100 200 300];
A = [1 0.8 0.5 0.3];
x = sum(A.*sin(2*pi*f.'*t));

% 设置STFT参数
window = 128;
noverlap = 64;
nfft = 256;

% 计算STFT
[S,F,T,P] = spectrogram(x, window, noverlap, nfft, fs);

% 绘制STFT频谱图
figure;
imagesc(T,F,10*log10(abs(S)));
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('STFT Spectrogram');

% 计算谱峰图
[max_val, max_idx] = max(abs(S));
f_peak = F(max_idx);

% 绘制谱峰图
figure;
plot(T, f_peak);
xlabel('Time (s)');
ylabel('Peak Frequency (Hz)');
title('Spectral Peak');

% 重建信号
x_recon = zeros(size(x));
for i = 1:length(f_peak)
    x_recon = x_recon + sin(2*pi*f_peak(i)*t);
end

% 绘制原始信号和重建信号
figure;
plot(t, x, 'b');
hold on;
plot(t, x_recon, 'r--');
xlabel('Time (s)');
ylabel('Amplitude');
title('Sinusoid Reconstruction');
legend('Original Signal', 'Reconstructed Signal');

运行该脚本后，会生成三个图形窗口。第一个窗口显示的是STFT频谱图，第二个窗口显示的是谱峰图，第三个窗口显示的是原始信号和使用sinusoid重建信号的对比。

需要注意的是，这只是一个简单的例子，实际应用中需要根据具体的需求选择合适的STFT参数，并进行适当的信号处理和重建。