《数字信号处理实验报告--电话拨号音合成与识别》由会员分享，可在线阅读，更多相关《数字信号处理实验报告--电话拨号音合成与识别（11页珍藏版）》请在金锄头文库上搜索。
1、. . . 理工大学实验报告学院（系）： 专业： 班级： 姓 名： 学号： 组： _ 实验时间： 实验室： 实验台： 指导教师签字： 成绩： 实验二 拨号音合成与识别一、实验结果与分析数字拨号音频谱分析（代码与频谱）：. . n=1:410;fs=8192;d1=sin(2*pi*697/fs*n)+sin(2*pi*1209/fs*n);D1=fft(d1);t1=(0:length(d1)-1)*fs/length(d1)-fs/2;figure(1);plot(t1,fftshift(abs(D1);title(按键1的频谱);d1=sin(2*pi*697/fs*n)+sin(2*pi
2、*1336/fs*n);D1=fft(d1);figure(2);plot(t1,fftshift(abs(D1);title(按键2的频谱);d1=sin(2*pi*697/fs*n)+sin(2*pi*1477/fs*n);D1=fft(d1);figure(3);plot(t1,fftshift(abs(D1);title(按键3的频谱);d1=sin(2*pi*770/fs*n)+sin(2*pi*1209/fs*n);D1=fft(d1);figure(4);plot(t1,fftshift(abs(D1);title(按键4的频谱);d1=sin(2*pi*770/fs*n)+si
3、n(2*pi*1336/fs*n);D1=fft(d1);figure(5);plot(t1,fftshift(abs(D1);title(按键5的频谱);d1=sin(2*pi*770/fs*n)+sin(2*pi*1477/fs*n);D1=fft(d1);figure(6);plot(t1,fftshift(abs(D1);title(按键6的频谱);d1=sin(2*pi*852/fs*n)+sin(2*pi*1209/fs*n);D1=fft(d1);figure(7);plot(t1,fftshift(abs(D1);title(按键7的频谱);d1=sin(2*pi*852/fs
4、*n)+sin(2*pi*1336/fs*n);D1=fft(d1);figure(8);plot(t1,fftshift(abs(D1);title(按键8的频谱);d1=sin(2*pi*852/fs*n)+sin(2*pi*1477/fs*n);D1=fft(d1);figure(9);plot(t1,fftshift(abs(D1);title(按键9的频谱);d1=sin(2*pi*941/fs*n)+sin(2*pi*1209/fs*n);D1=fft(d1);figure(10);plot(t1,fftshift(abs(D1);title(按键*的频谱);d1=sin(2*pi
5、*941/fs*n)+sin(2*pi*1336/fs*n);D1=fft(d1);figure(11);plot(t1,fftshift(abs(D1);title(按键0的频谱);d1=sin(2*pi*941/fs*n)+sin(2*pi*1477/fs*n);D1=fft(d1);figure(12);plot(t1,fftshift(abs(D1);title(按键#的频谱);频谱图如下： 图形拨号面板的制作：首先是框架的搭建：利用callback在.m文件中对各个模块进行代码的编辑（下面列出主要部分）：1、数字09的显示与发声（以数字7为例）：% - Executes on but
6、ton press in pushbutton7.function pushbutton7_Callback(hObject, eventdata, handles)% hObject handle to pushbutton7 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)n0=strcat(get(handles.numshow,string),7); % 获取数字
7、set(handles.numshow,string,n0); % 显示 n=1:handles.DTMFnum; % 每个数字 410 个采样点表示d0=sin(2*pi*852/handles.fs*n)+sin(2*pi*1209/handles.fs*n); % 对应行频列频叠加 space=zeros(1,handles.DTMFnum); %400 个 0 模拟静音信号 phone=handles.NUM,d0; handles.NUM=phone,space; % 存储连续的拨号音信号 guidata(hObject, handles);wavplay(d0,8192);2、删除
8、键*的代码：% - Executes on button press in pushbutton11.function pushbutton11_Callback(hObject, eventdata, handles)% hObject handle to pushbutton11 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)n=1:1000; num=get(han
9、dles.numshow,string); l=length(num); n11=strrep(num,num,num(1:l-1); d11=sin(0.7217*n)+sin(0.9273*n); set(handles.numshow,string,n11); L=length(handles.NUM); handles.NUM=handles.NUM(1:L-820); guidata(hObject, handles);wavplay(d11,8192);3、确认键#的代码：% - Executes on button press in pushbutton12.function p
10、ushbutton12_Callback(hObject, eventdata, handles)% hObject handle to pushbutton12 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)n0=strcat(get(handles.numshow,string),#); % 获取数字 set(handles.numshow,string,n0); %
11、 显示 n=1:handles.DTMFnum; % 每个数字 410 个采样点表示d0=sin(2*pi*941/handles.fs*n)+sin(2*pi*1477/handles.fs*n); % 对应行频列频叠加 guidata(hObject, handles);wavplay(d0,8192);4、清空键Reset的代码：% - Executes on button press in btFW.function btFW_Callback(hObject, eventdata, handles)% hObject handle to btFW (see GCBO)% eventd
12、ata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)handles.NUM=;set(handles.numshow,string,); % 显示 set(handles.numRec,string,); % 显示 guidata(hObject, handles);5、显示键的代码：% - Executes on button press in btRec.function btRec_Callback(hOb
13、ject, eventdata, handles)% hObject handle to btRec (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)L=length(handles.NUM); n=L/handles.DTMFsum; number=; for i=1:n j=(i-1)*handles.DTMFsum+1; d=handles.NUM(j:j+(handles.DTMFnum-1); % 截取出每个数字 f=fft(d,8192); % 以 N=2048 作 FFT 变换 a=abs(f); p=a.*a/handles.fs; % 计算功率谱% p=a.*a/10000; % 计算功率谱 num(1)=find(p(1:1000)=max(p(1:1000); % 找行频 num(2)=1000+find(p(1000:1700)=max(p(1000:1700); % 找列频 if (num(1) 730)
