[yFsbits]=wavread（‘music.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（221）;plot（fabsY）;xlabel（‘Frequency（Hz）‘）;grid?on
axis（[0f（end）/20500]）;
t=（0:sigLength-1）/Fs;
subplot（2