[yFsbits]=wavread（‘voice.wav‘）;%讀出信號，采樣率和采樣位數。?

sigLength=length（y）;
Y?=?fft（ysigLength）;
%Pyy?=?Y.*?conj（Y）?/?sigLength;%?conj函數，用于求一個復數的復共軛
halflength=floor（sigLength/2）;?%?根據模板大小確定要擴展的邊界
f=Fs*（1:sigLength）/sigLength;
stor=Y（1）;
Y（1）=0;
absY=abs（Y）;
figure（1）
subplot（321）;plot（fabsY）;xlabel（‘Frequency（Hz）‘）;gr