(44)MATLAB读取语音信号进行频谱分析

文章目录


前言

语音信号是我们最常见的一种信号,本文使用MATLAB读取一段语音信号画出其波形,然后使用FFT变换给出其频谱,对其频谱进行分析。


一、MATLAB代码

读取语音数据并得出频谱的代码如下:

c 复制代码
%% 读取语音数据
load('speech_8ksps.mat', 'speech_8k')     % load variable data

L = length(speech_8k);
fs = 8000;                                % 语音信号采样率8ksps
t = (0:1/fs:(length(speech_8k)-1)/fs);    % 语音信号的时间刻度

figure()
plot(t,speech_8k)
grid on
title('一段语音信号(fs=8kSPS)')
xlabel('Time(s)')
ylabel('Amplitude')

%% 频域(FFT)结果解析
N = length(speech_8k);
Y = fft(speech_8k,N);                     % 信号的傅里叶变换

% 信号的双边幅度频谱
P2 = abs(Y/L);

% 信号的双边幅度频谱
figure()
plot(0:(fs/N):(fs-fs/N),P2(1:N),'LineWidth',1)
grid on
title(['语音的双边幅度频谱(fs=',num2str(fs),' samples/s)'])
xlabel('频率/Hz')
ylabel('幅度|Y(n)|')

%% 信号的单边幅度频谱
% 计算信号的单边幅度频谱
P1 = P2(1:N/2+1);
P1(2:end-1) = 2*P1(2:end-1);

figure()
plot(0:(fs/N):(fs/2),P1(1:N/2+1),'LineWidth',1)
grid on
title(['语音的单边幅度频谱(fs=',num2str(fs),' samples/s)'])
xlabel('频率(Hz)  (单边谱的频率范围0Hz到fs/2)')              % 频率范围0Hz到fs/2
ylabel('频谱幅度|Y|')

二、仿真结果画图

语音信号的时域波形图:

语音信号的双边幅度谱:

语音信号的单边幅度谱:

三、频谱分析

语音信号我们一般取300Hz-3000Hz的频率范围。一般情况下,男子的语音频谱处于相对较低的范围,如上面图中所示,处于300Hz-1500Hz的频率范围。女子的语音频谱处于相对较高的范围。


相关推荐
cnbestec2 分钟前
协作机器人UR7e与UR12e:轻量化设计与高负载能力助力“小而美”智造升级
人工智能·机器人·协作机器人·ur协作机器人·ur7e·ur12e
zskj_zhyl5 分钟前
毫米波雷达守护银发安全:七彩喜跌倒检测仪重构居家养老防线
人工智能·安全·重构
秋说1 小时前
【PTA数据结构 | C语言版】一元多项式求导
c语言·数据结构·算法
gaosushexiangji1 小时前
利用sCMOS科学相机测量激光散射强度
大数据·人工智能·数码相机·计算机视觉
Maybyy1 小时前
力扣61.旋转链表
算法·leetcode·链表
ai小鬼头2 小时前
AIStarter新版重磅来袭!永久订阅限时福利抢先看
人工智能·开源·github
说私域3 小时前
从品牌附庸到自我表达:定制开发开源AI智能名片S2B2C商城小程序赋能下的营销变革
人工智能·小程序
卡卡卡卡罗特3 小时前
每日mysql
数据结构·算法
飞哥数智坊3 小时前
新版定价不够用,Cursor如何退回旧版定价
人工智能·cursor
12点一刻4 小时前
搭建自动化工作流:探寻解放双手的有效方案(2)
运维·人工智能·自动化·deepseek