dtmfen.m

Matlab源程序

% 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);   %构成双频信号
    y = awgn(x1,4) ;                                     %加高斯噪声
    gg=y;