matlab频谱感知代码

以下是一个简单的 MATLAB 频谱感知代码示例：

% 读取音频文件
[y, Fs] = audioread('audio_file.wav');

% 设置参数
N = 1024;                  % FFT 点数
n_overlap = N/2;           % 帧重叠长度
n_frames = floor(length(y)/(N-n_overlap))-1;   % 帧数
freq_range = [1000 4000];  % 感兴趣的频率范围

% 分帧并计算每一帧的频谱
spectrogram_data = zeros(N, n_frames);
for i = 1:n_frames
    frame_data = y((i-1)*(N-n_overlap)+1:i*(N-n_overlap)+N);
    spectrogram_data(:,i) = abs(fft(frame_data, N));
end

% 计算每一帧中感兴趣的频率范围内的能量
energy_data = zeros(n_frames, 1);
for i = 1:n_frames
    freq_start = round(freq_range(1)/(Fs/N));
    freq_end = round(freq_range(2)/(Fs/N));
    energy_data(i) = sum(spectrogram_data(freq_start:freq_end,i).^2);
end

% 判断是否有信号
threshold = 0.1;  % 设定阈值
if max(energy_data) > threshold
    disp('检测到信号');
else
    disp('未检测到信号');
end

这个代码的基本思路是将音频文件分成一帧帧的数据，并计算每一帧的频谱。然后，它将每一帧中感兴趣的频率范围内的能量求和，以估计信号的强度。最后，代码使用一个预设的阈值来判断是否有信号。如果能量超过阈值，则认为检测到信号。