dtmf.asv

来自「基于matlab的DTMF仿真 DTMF为行频和列频（697Hz 770Hz 」· ASV 代码 · 共 27 行

f = [697 770 852 941 1209 1336 1477 1633];
coef = [27980 26956 25701 24219 19073 16325 13085 9315];
%coef = [0.8539    0.8226    0.7843    0.7391    0.5821    0.4982    0.3993    0.2843];
fs = 8000;
N = 125;
k = 8;
lf = 3;
hf = 8;
is = cos(0 : (2*pi*f(lf))/fs : (2*N*pi*f(lf))/fs) + cos(0 : (2*pi*f(hf))/fs : (2*N*pi*f(hf))/fs);
%plot(is);
v(1:k,1:2) = zeros(k,2);
X2(1:k) = zeros(k,1);
for i = 1:N    
    for j = 1:k
        % vk(n) = 2*coef*vk(n-1) - vk(n-2) + x(n)
        tmp = 1*coef(j)*v(j,2)/16384-v(j,1)+is(i),
        v(j,1) = v(j,2),
        v(j,2) = tmp
    end
end
for j = 1:k
    % y(N)*y(N) = vk(N)*vk(N) - 2*coef*vk(N)*vk(N-1) + vk(N-1)*vk(N-1)
    X2(j) = v(j,2)*v(j,2)+v(j,1)*v(j,1)-1*coef(j)*v(j,1)*v(j,2)/16384;
end
X2;
plot(X2,'+');