creat_dtmf.m

来自「基于MATLAB的双音多频信号DTMF的生成。」· M 代码 · 共 47 行

%**************************************************************
%Name：digital oscillation
%Function：   
%Input：
%Output：
%Creater： huiqian.cai
%Create date：2008/1/26
%**************************************************************
clear all
fd=697;               %fd振荡器频率，fs采样频率（即中断输出频率）
fs=8000;
A=2*cos(2*pi*fd/fs);   %A,B,C作为系数值,可以预先存入寄存器 
B=-1;
C=2*sin(2*pi*fd/fs);
y(1)=0;                %初始值
y(2)=C;
for n=3:100,
y(n)=A*y(n-1)+B*y(n-2);
end;
y1=y;
subplot(3,1,1);
plot(y);
xlabel('697HZ');
 

fd=1209;               %fd振荡器频率，fs采样频率（即中断输出频率）
fs=8000;
A=2*cos(2*pi*fd/fs);   %A,B,C作为系数值,可以预先存入寄存器 
B=-1;
C=2*sin(2*pi*fd/fs);
y(1)=0;                %初始值
y(2)=C;
for n=3:100,
y(n)=A*y(n-1)+B*y(n-2);
end;
y2=y;
subplot(3,1,2);
plot(y);
xlabel('1209HZ');
y3=y1+y2;
subplot(3,1,3);
plot(y3);
xlabel('dial"1"');