DTMF系统：双音多频系统系统：双音多频系统y,fs,nbits=wavread(file)采样值放在向量y中，fs表示采样频率（Hz），nbits表示采样位数sound(x,fs,bits)wavwrite(y, fs, nbits, waveFile)load trainsound(y,Fs)clear load touchtone n = length(y.sig); t = (0:n-1)/y.fs; Fs=y.fs; yy = double(y.sig)/128; plot(t,yy,k) p=abs(fft(yy); f=(0:n-1)*(Fs/n); figure,plot(f,p,g) axis(500 1700 0 600)截取第一段进行分析截取第一段进行分析 %第一个声音第一个声音 n1=find(t0&t0&t0.42); y1=yy(n1); y1=yy(n1); m1=length(y1); m1=length(y1); % %傅里叶分析傅里叶分析 p1=abs(fft(y1); p1=abs(fft(y1); f1=(0:m1-1)*(Fs/m1); f1=(0:m1-1)*(Fs/m1); figure,plot(f1,p1) figure,plot(f1,p1) axis(500 1700 0 600) axis(500 1700 0 600)