【信号提取】用于在单个帧（段）上进行特征提取（Matlab实现）

原创已于 2024-10-18 18:55:41 修改 · 526 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#matlab #人工智能 #开发语言

于 2024-10-18 18:49:35 首次发布

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

💥1 概述

用于在单个帧（段）上进行特征提取。这个过程涉及从特定的帧或片段中分析并提取特定的特征或属性。它可以应用于各种领域，如图像处理、视频分析或信号处理。通过从单个帧中提取特征，可以深入了解该特定时刻或片段的内容或特性。这有助于完成诸如对象识别、模式检测或分类等任务。提取的特征可以用于进一步处理或作为机器学习算法的输入以执行更复杂的任务。

📚2 运行结果

主函数部分代码：

clear, clc, close all

%% get a section of the sound file
[x, fs] = audioread('DR2_MJAR0_SI2247.wav');    % load an audio file
x = x(:, 1);                                    % get the first channel
x = x/max(abs(x));                              % normalize the signal
N = length(x);                                  % signal length
t = (0:N-1)/fs;                                 % time vector

%% signal framing
frlen = round(20e-3*fs);                        % frame length
hop = round(frlen/2);                           % hop size
[FRM, tfrm] = framing(x, frlen, hop, fs);       % signal framing

%% determine the Short-time energy
STE = sum(FRM.^2);

%% determine the Short-time zero-crossing rate
STZCR = sum(abs(diff(FRM > 0)))/size(FRM, 1);

%% plot the results
% plot the signal waveform
figure(1)
subplot(3, 1, 1)
plot(t, x)
grid on
xlim([0 max(t)])
ylim([-1.1*max(abs(x)) 1.1*max(abs(x))])
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14)
xlabel('Time, s')
ylabel('Amplitude')
title('The signal in the time domain')

% plot the STE
subplot(3, 1, 2)
plot(tfrm, STE)
grid on
xlim([0 max(tfrm)])
ylim([0 1.1*max(STE)])
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14)
xlabel('Time, s')
ylabel('Energy')
title('Short-time Energy')

% plot the STZCR
subplot(3, 1, 3)
plot(tfrm, STZCR)
grid on
xlim([0 max(tfrm)])
ylim([0 1.1*max(STZCR)])
set(gca, 'FontName', 'Times New Roman', 'FontSize', 14)
xlabel('Time, s')
ylabel('ZCR')
title('Short-time ZCR')

%% mark the signal
% Note: (i) the activity flag is rised when there is high energy frame
% in the signal stream; (ii) the consonant flag is rised when there are
% high ZCR activity and high energy in the signal stream.
AF = STE > 0.1*median(STE);
CF = AF & STZCR>2*median(STZCR);

subplot(3, 1, 1)
hold on
plot(tfrm, AF, 'r', 'LineWidth', 1.5)
plot(tfrm, CF, 'g', 'LineWidth', 1.5)
legend('signal', 'activity flag', 'consonant flag', ...
       'Location', 'southeast')