Matlab识别拨号音,大家帮我看下这段用哪个FFT做根据拨号音识别电话号码的...

本帖最后由 想回到那一天 于 2018-4-30 20:29 编辑

%function number=duandianjiance() %函数返回的就是读到的拨号值

%读一段语音信号 这里要填写语音信号的绝对路径

clc

clear

[x, Fs]=audioread('matlab1.wav');

x = double(x);

x = x / max(abs(x));%幅度归一化到[-1,1]

fp1=[500 1700]; fp2=[300 2100];

[n1,Wn1]=ellipord(fp1*2/Fs,fp2*2/Fs,0.5,40);

[b1,a1]=ellip(n1,0.5,40,Wn1);

Output=filter(b1,a1,x);

x=Output;

y=zeros(1,length(x)); %声明y为去掉空白音之后保存新的信号

%常数设置

FrameLen = 240;%帧长

FrameInc = 80;%分帧未重叠部分

amp1 = 30;  %短时能量阈值

amp2 = 30; %即设定能量的两个阈值。

zcr1 = 20;%过零率阈值

zcr2 = 5; %过零率的两个阈值

maxsilence = 8; %用无声的长度来判断语音是否结束

minlen  = 40;   %判断是语音的最小长度

status  = 0;%记录语音段的状态

count   = 0; %语音序列的长度

silence = 0; %无声的长度

%计算过零率

tmp1  = enframe(x(1:end-1), FrameLen, FrameInc);

tmp2  = enframe(x(2:end)  , FrameLen, FrameInc);

signs = (tmp1.*tmp2)<0;

diffs = (tmp1 -tmp2)>0.02;

zcr   = sum(signs.*diffs, 2);

%计算短时能量

amp = sum(abs(enframe(filter([1 -0.9375], 1, x), FrameLen, FrameInc)), 2);

%调整能量门限

amp1 = min(amp1, max(amp)/2);

amp2 = min(amp2, max(amp)/2.5);

%开始端点检测

x1 = 0;

x2 = 0;

length_now=1;

for n=1:length(zcr)

goto = 0;

switch status

case {0,1}                   % 0 = 静音, 1 = 可能开始

if amp(n) > amp1          % 确信进入语音段

x1 = max(n-count-1,1);

status  = 2;

silence = 0;

count   = count + 1;

elseif amp(n) > amp2 | ... % 可能处于语音段

zcr(n) > zcr2

status = 1;

count  = count + 1;

else                       % 静音状态

status  = 0;

count   = 0;

end

case 2,                       % 2 = 语音段

if amp(n) > amp2 | ...   % 保持在语音段

zcr(n) > zcr2

count = count + 1;

else                       % 语音将结束

silence = silence+1;

if silence < maxsilence % 静音还不够长，尚未结束

count  = count + 1;

elseif count < minlen   % 语音长度太短，认为是噪声

status  = 0;

silence = 0;

count   = 0;

else                    % 语音结束

status  = 3;

end

end

case 3,

count = count-silence/2;

x2 = x1 + count -1;%记录语音段结束点 x1为起始点，x2为结束点

%将每一个语音信号长度限制到6900方便后面识别

if ((x2*FrameInc-x1*FrameInc)<6900)

x3=[x(x1*FrameInc:x2*FrameInc)',zeros(1,(6900-x2*FrameInc+x1*FrameInc))]';

else

x3=x(x1*FrameInc:(x1*FrameInc+6900));

end

y(length_now:(length_now+length(x3)-1))=x3;%将每次得到的一个拨号音添加到y序列里面

length_now=length_now+length(x3); %记录非空白语音长度

status=0;

count=0;

end

end

%这些都是简单的小时

y1=y(1:length_now);

%subplot(211)

figure(1)

plot(x) %显示原语音信号的时域谱

%subplot(212)

figure(2)

plot(y1) %显示端点检测滤波之后的时域谱

%wavplay(y1,Fs);

number=shibie(y1,Fs);%调用识别函数进行识别

其中的shibei和enframe子函数是有的，程序没什么语法错误，但是检测不了电话号码出来，是为啥？

2018-4-30 20:29 上传

点击文件名下载附件

53.48 KB, 下载次数: 22