经典课程指导：《数字处理系统课程设计与实践》指导(2)DTMF之四 自适应检测

算法工匠

1月前浏览239

关注

继续讲解！循序渐进的讲解需要有听众，可惜自己线下课堂的学生只有一两位能跟着学！那还有写的必要吗？当然有！网上很多人想学的。数字信号处理技术是可以倚靠吃饭的家伙，学精了，还愁没有好工作吗？

面对学生们的懒惰，只能不断的劝说，效果甚微。一学期结束后，老师的心态也需要不断的调整。我自己做了几年的老师后，发现真的很难有成就感。经常感慨：教书不可为不用心，但听我的学生又有几人？不管学生们的表现如何，我自己是一直有危机感的。至少直到现在，我还在坚持每天写代码和看书。我会坚持下去，因为这也成为了我的生活方式。

关于流信号和段信号处理方式的差异已经录成视频课程（数字信号处理仿真），在仿真秀网站上可以看到！本文接近四千字，仿真程序内容较多，非专业人士莫入！

流处理！

如何把goertzel算法应用于流信号处理？需要大家好好思考了！现在请跟着老师的代码一步一步的学下去！引导式教学需要大家动脑筋，这样掌握了就是自己的知识！先看程序的仿真结果！

大家看的懂仿真结果吗？再来看具体代码！

% Generate DTMF tones

Fs  = 8000;       % Sampling frequency 8 kHz

N   = 800;        % Tones of 100 ms

t   = (0:N-1)/Fs; % 800 samples at Fs

pit = 2*pi*t;

tones = zeros( N, numel(symbol));

for i = 1 : numel(symbol)

    switch (symbol{i})

        case '1'

            tones(:,i) = sum(sin(f(:,1)*pit))'; % 合并两个音频信号

        case '2'

            tones(:,i) = sum(sin(f(:,2)*pit))';

        case '3'

            tones(:,i) = sum(sin(f(:,3)*pit))';

        case '4'

            tones(:,i) = sum(sin(f(:,4)*pit))';

        case '5'

            tones(:,i) = sum(sin(f(:,5)*pit))';

        case '6'

            tones(:,i) = sum(sin(f(:,6)*pit))';

        case '7'

            tones(:,i) = sum(sin(f(:,7)*pit))';

        case '8'

            tones(:,i) = sum(sin(f(:,8)*pit))';

        case '9'

            tones(:,i) = sum(sin(f(:,9)*pit))';

        case '*'

            tones(:,i) = sum(sin(f(:,10)*pit))';

        case '0'

            tones(:,i) = sum(sin(f(:,11)*pit))';

        case '#'

            tones(:,i) = sum(sin(f(:,12)*pit))';

    end

end

for i = 1:12

    p = audioplayer(tones(:,i),Fs,16);

    % 在默认音频设备上播放该音频对象

    play(p)

    % 暂停执行0.5秒

    pause(0.5)

end

aaa = 1;

%% 用 Goertzel 算法估计 DTMF 音调

Nt = 205;

original_f = [lfg(:); hfg(:)];  % Original frequencies

k = round(original_f/Fs*Nt);  % Indices of the DFT

estim_f = round(k*Fs/Nt);

% Frequencies at which the DFT is estimated

tones = tones(1:205,:); 

% 截断音调，只保留205个采样或25.6毫秒

toneChoice = 12    % Select tone

tone = tones(:,toneChoice);

% Estimate DFT using Goertzel

% 自编函数

ydft = func_goertzel_v1(tone,k);

% 注意参数的差异！！！

figure

subplot(4,1,1),plot(abs(ydft(1,:)))

subplot(4,1,2),plot(abs(ydft(2,:)))

subplot(4,1,3),plot(abs(ydft(3,:)))

subplot(4,1,4),plot(abs(ydft(4,:)))

figure

subplot(3,1,1),plot(abs(ydft(5,:)))

subplot(3,1,2),plot(abs(ydft(6,:)))

subplot(3,1,3),plot(abs(ydft(7,:)))

继续深入下去！

修改信号源进行检测！

但同时要考虑到算法和数据的长度有关联！改动Nt的值会发现 Goertzel 算法的结果也在变化！

% Generate DTMF tones

Fs  = 8000;       % Sampling frequency 8 kHz

N   = 800;        % Tones of 100 ms

t   = (0:N-1)/Fs; % 800 samples at Fs

pit = 2*pi*t;

tones = zeros( N, numel(symbol));

for i = 1 : numel(symbol)

    switch (symbol{i})

        case '1'

            tones(:,i) = sum(sin(f(:,1)*pit))'; % 合并两个音频信号

        case '2'

            tones(:,i) = sum(sin(f(:,2)*pit))';

        case '3'

            tones(:,i) = sum(sin(f(:,3)*pit))';

        case '4'

            tones(:,i) = sum(sin(f(:,4)*pit))';

        case '5'

            tones(:,i) = sum(sin(f(:,5)*pit))';

        case '6'

            tones(:,i) = sum(sin(f(:,6)*pit))';

        case '7'

            tones(:,i) = sum(sin(f(:,7)*pit))';

        case '8'

            tones(:,i) = sum(sin(f(:,8)*pit))';

        case '9'

            tones(:,i) = sum(sin(f(:,9)*pit))';

        case '*'

            tones(:,i) = sum(sin(f(:,10)*pit))';

        case '0'

            tones(:,i) = sum(sin(f(:,11)*pit))';

        case '#'

            tones(:,i) = sum(sin(f(:,12)*pit))';

    end

end

来源：通信工程师专辑