tm=[12365;45666;78967;4203568];?%?DTMF信號代表的16個數(shù)?
N=205;
K=[1820222431343842];?
f1=[697770852941];?%?行頻率向量?
f2=[1209133614771633];?%?列頻率向量?
TN=input（‘鍵入11位電話號碼=?‘）;?%?輸入11位數(shù)字?
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;?%?為了發(fā)聲，加長序列?
x?=?sin（2*pi*n*f1（p）/fs）?+?sin（2*pi*n*f2（q）/fs）;%?構(gòu)成雙頻信號?
sound（x8000）;?%?發(fā)出聲音（真的就是發(fā)出話音信號的聲音了）
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