clf;
fs=100000;
yy=wavread('D:/shuju/file301.wav');
yj=resample(yy,1,10);
y=yj(1:900000);%取90s的信號
subplot(411);
plot(y);%時域圖
axis([0 900000 0 1.5]);
yf=abs(fft(y));
n=length(y);
f=linspace(0,fs/10,n);
subplot(412);
plot(f,yf);%頻域圖
axis([0 700 0 2000]);
y0=y(1:100000);
y1=y(50001:150000);y2=y(100001:200000);
y3=y(150001:250000);y4=y(200001:300000);
y5=y(250001:350000);y6=y(300001:400000);
y7=y(350001:450000);y8=y(400001:500000);
y9=y(450001:550000);y10=y(500001:600000);
y11=y(550001:650000);y12=y(600001:700000);
y13=y(650001:750000);y14=y(700001:800000);
y15=y(750001:850000);y16=y(800001:900000);
yf0=abs(fft(y0));yf1=abs(fft(y1));yf2=abs(fft(y2));yf3=abs(fft(y3));
yf4=abs(fft(y4));yf5=abs(fft(y5));yf6=abs(fft(y6));yf7=abs(fft(y7));
yf8=abs(fft(y8));yf9=abs(fft(y9));yf10=abs(fft(y10));yf11=abs(fft(y11));
yf12=abs(fft(y12));yf13=abs(fft(y13));yf14=abs(fft(y14));yf15=abs(fft(y15));
yf16=abs(fft(y16));
yz=(yf0+yf1+yf2+yf3+yf4+yf5+yf6+yf7+yf8+yf9+yf10+yf11+yf12+yf13+yf14+yf15+yf16)/17;
m=length(yf3);
fn=linspace(0,fs/10,m);
subplot(413);
plot(fn,yz);
axis([0 200 0 400]);
subplot(414);
plot(fn,yz);
axis([10 95 0 400]);