% code sub-programme編碼子程序
function dtmfen=DTMFen(d)
global n;
global gg;
symbol=abs(d);                        %求它的ASCII碼   
tm=[49,50,51,65;52,53,54,66;55,56,57,67;42,48,35,68];  %16個ASCII碼
for p=1:4;
    for q=1:4;
        if tm(p,q)==abs(d);break,end             %檢測碼相符的列號
    end
         if tm(p,q)==abs(d);break,end            %檢測碼相符的行號
    end
    f1=[697,770,852,941];                        %行頻率向量
    f2=[1209,1336,1477,1633];                    %列頻率向量
    n=0:2040;
    x1=sin(2*pi*n*f1(p)/8000)+sin(2*pi*n*f2(q)/8000);   %構(gòu)成雙頻信號
    y = awgn(x1,4) ;                                     %加高斯噪聲
    gg=y;