Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2012-03-19, 21:47   #1
gda1990
初级会员
 
gda1990 的头像
 
注册日期: 2010-06-08
年龄: 34
帖子: 5
声望力: 0
gda1990 正向着好的方向发展
默认 语谱图分析!!!!

请问语谱图是如何看出共振峰和基频,谐波分量的?请好心人举例说明(看了半天“横条竖条”都没看懂T_T)
窄带与宽带语谱图是如何设置呢(如下面的程序部分该改哪个参数才能改为宽带语谱图)
若有一段8khz的语音信号,维持时间约6s,那么窗长应当设置为多少?(分别对应于宽窄带)
给出程序(Beijing.wav是语音“北京我来了”)的窄带语谱图
clear all;
[x,sr]=wavread('Beijing.wav'); %sr为采样频率
if (size(x,1)>size(x,2)) %size(x,1)为x的行数,size(x,2)为x的列数
x=x';
end
s=length(x);
w=round(44*sr/1000); %窗长,取离44*sr/100最近的整数
n=w; %fft的点数
ov=w/2; %50%的重叠
h=w-ov;
specgram('Beijing.wav',n,sr,w,w/2);
% win=hanning(n)'; %哈宁窗
上传的图像
文件类型: jpg F1XW1OX87O`LY1$LI635$NE.jpg (88.1 KB, 3 次查看)
__________________
If you fall, dust it off, don’t let up
gda1990 当前离线   回复时引用此帖
旧 2012-03-19, 21:48   #2
gda1990
初级会员
 
gda1990 的头像
 
注册日期: 2010-06-08
年龄: 34
帖子: 5
声望力: 0
gda1990 正向着好的方向发展
默认 回复: 语谱图分析!!!!

win=hamming(n)'; %哈宁窗
c=1;
ncols=1+fix((s-n)/h); %fix函数是将(s-n)/h的小数舎去
d=zeros((1+n/2),ncols);
for b=0:hs-n)
u=win.*x((b+1)b+n));
t=fft(u);
d(:,c)=t(11+n/2))';
c=c+1;
end
__________________
If you fall, dust it off, don’t let up
gda1990 当前离线   回复时引用此帖
旧 2012-03-19, 21:49   #3
gda1990
初级会员
 
gda1990 的头像
 
注册日期: 2010-06-08
年龄: 34
帖子: 5
声望力: 0
gda1990 正向着好的方向发展
默认 回复: 语谱图分析!!!!

tt=[0:hs-n)]/sr;
ff=[0n/2)]*sr/n;
imagesc(tt/1000,ff/1000,20*log10(abs(d)));

axis xy
xlabel('时间/s');
ylabel('频率/kHz');
__________________
If you fall, dust it off, don’t let up
gda1990 当前离线   回复时引用此帖
旧 2012-03-19, 21:51   #4
gda1990
初级会员
 
gda1990 的头像
 
注册日期: 2010-06-08
年龄: 34
帖子: 5
声望力: 0
gda1990 正向着好的方向发展
默认 回复: 语谱图分析!!!!

被系统弄成表情的是: (
__________________
If you fall, dust it off, don’t let up
gda1990 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 12:41


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.