功率谱密度

clear
clc
close all

% 生成间谐时程信号
fs = 1000; % 采样频率为 1000 Hz
t = (0:1/fs:300)'; % 时间从0到300秒，步长为采样间隔
signal = sin(2*pi*0.1*t) + sin(2*pi*0.1*t) +rand(length(t),1);

% 使用pwelch计算功率谱密度函数
[pxx, f] = pwelch(signal, [], [], [], fs);

% 绘制信号波形
figure;
subplot(2,1,1);
plot(t, signal);
title('间谐时程信号');
xlabel('时间 (秒)');
ylabel('幅值');

% 绘制功率谱密度函数
subplot(2,1,2);
plot(f, 10*log10(pxx));
title('功率谱密度函数');
xlabel('频率 (Hz)');
ylabel('功率谱密度 (dB/Hz)');
xlim([0 1])

% 标记峰值点
[peak_values, peak_indices] = findpeaks(10*log10(pxx), 'MinPeakHeight', 10); % 这里设置了一个峰值的最小高度
hold on;
scatter(f(peak_indices), peak_values, 'r', 'filled');
text(f(peak_indices), peak_values, num2str(round(f(peak_indices), 2)), 'VerticalAlignment', 'bottom', 'HorizontalAlignment', 'right');

hold off;