[yFsbits]=wavread（‘music.wav‘）;%讀出信號，采樣率和采樣位數(shù)。?
sigLength=length（y）;
Y?=?fft（ysigLength）;
%Pyy?=?Y.*?conj（Y）?/?sigLength;%?conj函數(shù)，用于求一個復(fù)數(shù)的復(fù)共軛
halflength=floor（sigLength/2）;?%?根據(jù)模板大小確定要擴(kuò)展的邊界
f=Fs*（1:sigLength）/sigLength;
stor=Y（1）;
Y（1）=0;
absY=abs（Y）;
figure（1）
subplot（221）;plot（fabsY）;xlabel（‘Frequency（Hz）‘）;grid?on
axis（[0f（end）/20500]）;
t=（0:sigLength-1）/Fs;
subplot（2