% clear
% clc
% close all
% % 生成随机时程信号
% fs=100;
% signalLength = fs*60*2;
% time = 1/fs:1/fs:signalLength/fs;
% randomSignal = 2*sin(2*pi*0.5*time)+3*cos(2*pi*1*time)+randn(1, signalLength);
function [frequencyPP]=funct_peak(signal,Hz)
% 生成随机时程信号
fs=Hz;
randomSignal=signal;
signalLength = length(randomSignal);
time = 1/fs:1/fs:signalLength/fs;
%功率谱密度图
pwel.window=1024*9;
pwel.noverlap=pwel.window/2
pxx,f\]=pwelch(randomSignal,pwel.window,pwel.noverlap,\[\],fs); plot(f,10\*log10(pxx)) xlim(\[0 2\]) % 寻找峰值点 threshold = 20; % 定义阈值，即前后20个点的纵坐标差值的阈值 \[pks, locs\] = findpeaks(10\*log10(pxx), 'MinPeakHeight', threshold); % 标识峰值点 hold on scatter(f(locs), pks, 'r', 'filled') % 在每个峰值点的位置添加文本标签 for i = 1:length(locs) text(f(locs(i)), pks(i), sprintf('%.2f', f(locs(i))), 'VerticalAlignment', 'bottom', 'HorizontalAlignment', 'center'); frequencyPP(i)=round(f(locs(i)),3);, end hold off