信号与系统-电话拨号号码识别检测-仅含matlab代码和结果图

【目的】

(1)掌握信号时域和频域分析方法，对系统时域分析基本原理和方法的理解；

(2)掌握利用MATLAB对电话拨号音频进行识别；

(3)培养学生综合利用所学信号处理的理论和技术解决实际问题的能力。

【知识点】

双音多频(DTMF)基本概念、信号时域和频域分析

【背景知识】电话拨号号码识别

【情况描述】现有一个犯罪的场景，嫌疑人在作案后给他的上级打电话，拨号声被现场的录音设备记录了下来，录音文件为dial.wav（见附件。时域波形图如图１所示）。请你根据录音文件，将电话号码分析出来，协助警方破案。

注：此录音文件已经进行了初步处理，每个拨号声0.1s，随即是0.1s的静音（噪声不计算在内），且周围环境的噪声已被去除。

【研讨内容】

1. (*)对录音文件分别进行时域和频域分析，判断在哪个域可以更好地分析出电话号码；
2. (***)利用Matlab语言设计电话拨号识别程序，实现电话号码的检测与识别。

【研讨要求】

(1)编写实现上述研讨内容的MATLAB仿真程序，逐段画出电话拨号音频的时域和频域波形；

(2)要将电话拨号声识别，需要对每个电话号码的两个频率进行判断，说出判断方法；

(3)在研讨过程中，若发现问题，要进行探究，并给出解决问题的方案。

---------------------------------------------------------------------------------------------------------------------------------

问题解决：

[y,fs]=audioread('电话拨号识别/dial.wav');
time=(1:length(y))/fs;
figure;
plot(time,y);
n=length(y)/fs/0.2; %n个拨号键
for i=1:n
    j=(i-1)*fs*0.2+1;
    m=y(j:j+fs*0.2-1);  %0.2秒提取一次
   [M,f]=ctft1(m,fs,length(m));
  subplot(4,3,i);
plot(f,abs(M)/max(abs(M)));
 title(['第',num2str(i),'个号码']);
axis([500,2000,0,1]);
grid;
end

function[X,f]=ctft1(x,fs,n)
X=fftshift(fft(x,n))/fs;
f=-fs/2+(0:n-1)*fs/n;
end