快速傅里叶变换-(精品课件).pptVIP

DTMF系统：双音多频系统 [y,fs,nbits]=wavread(file) 采样值放在向量y中， fs表示采样频率（Hz）， nbits表示采样位数 sound(x,fs,bits) wavwrite(y, fs, nbits, waveFile) load train sound(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(t0t0.42); y1=yy(n1); m1=length(y1); %傅里叶分析 p1=abs(fft(y1)); f1=(0:m1-1)*(Fs/m1); figure,plot(f1,p1) axis([500 1700 0 600])