DTMF双音多频

亲，程序界面优美，程序代码简便，值得您下载。新手初学，有做的不好的地方请指教。

tm=[12365;45666;78967;4203568]; % DTMF信号代表的16个数 
N=205;
K=[1820222431343842]; 
f1=[697770852941]; % 行频率向量 
f2=[1209133614771633]; % 列频率向量 
TN=input（‘键入11位电话号码= ‘）; % 输入11位数字 
TNr=0; %接收端电话号码初值为零 
fs=8000;
for l=1:11; 
d=fix（TN/10^（11-l））; 
TN=TN-d*10^（11-l）; 
for p=1:4; 
for q=1:4; 
if tm（pq）==abs（d）; 
    break
end % 检测码相符的列号q 
end 
if tm（pq）==abs（d）;
    break
end % 检测码相符的行号p 
end 
n=0:1023; % 为了发声，加长序列 
x = sin（2*pi*n*f1（p）/fs） + sin（2*pi*n*f2（q）/fs）;% 构成双频信号 
sound（x8000）; % 发出声音（真的就是发出话音信号的声音了）
pause（0.15） 
X=goertzel（x（1:N）K+1）; % 用Goertzel算法计算八点DFT样本 
val = abs（X）; % 列出八点DFT向量 

%subplot（43l）; 
%stem（Kval‘.‘）;
%plot（x）;
%grid;
%xlabel（‘k‘）;
%ylabel（‘|X（k）|‘） % 画出DFT（k）幅度 
%axis（[10 50 0 120]） 

subplot（43l）; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2017-11-27 19:21  DTMF\
     文件       27406  2011-10-19 23:30  DTMF\A.mat
     文件        1275  2017-11-27 19:13  DTMF\dtmf.m
     文件       17942  2017-11-27 19:31  DTMF\DTMFlxy.asv
     文件       45963  2017-11-26 16:40  DTMF\DTMFlxy.fig
     文件       18009  2017-11-27 19:31  DTMF\DTMFlxy.m
     文件         968  2017-11-27 19:06  DTMF\one.asv
     文件         909  2017-11-27 19:17  DTMF\one.m
     文件      208166  2017-11-26 17:07  DTMF\xinhao.m4a