sig3 - tkyswiki @ ウィキ

sig1 = wavread("oc1.wav");
sig2 = wavread("oc2.wav");
sig3 = wavread("oc3.wav");

plot(sig3);
xlim([5000 11000]);
sig3a=sig3(7001:7512);

plot(sig3);
xlim([23000 28000]);
sig3i=sig3(25001:25512);

plot(sig3);
xlim([42000 48000]);
sig3u=sig3(44001:44512);

plot(sig3);
xlim([70000 77000]);
sig3e=sig3(72501:73012);

plot(sig3);
xlim([90000 100000]);
sig3o=sig3(94001:94512);

fft_sig3a = fft(sig3a);
log_sig3a = log(abs(fft_sig3a));
ifft_sig3a = ifft(log_sig3a);
zero_sig3a = ifft_sig3a;
zero_sig3a(12:512-12) = 0;
fft2_sig3a = fft(zero_sig3a);
plot(8000/512:8000/512:8000,fft2_sig3a,"r");

fft_sig3i = fft(sig3i);
log_sig3i = log(abs(fft_sig3i));
ifft_sig3i = ifft(log_sig3i);
zero_sig3i = ifft_sig3i;
zero_sig3i(12:512-12) = 0;
fft2_sig3i = fft(zero_sig3i);
plot(8000/512:8000/512:8000,fft2_sig3i,"r");

fft_sig3u = fft(sig3u);
log_sig3u = log(abs(fft_sig3u));
ifft_sig3u = ifft(log_sig3u);
zero_sig3u = ifft_sig3u;
zero_sig3u(12:512-12) = 0;
fft2_sig3u = fft(zero_sig3u);
plot(8000/512:8000/512:8000,fft2_sig3u,"r");

fft_sig3e = fft(sig3e);
log_sig3e = log(abs(fft_sig3e));
ifft_sig3e = ifft(log_sig3e);
zero_sig3e = ifft_sig3e;
zero_sig3e(12:512-12) = 0;
fft2_sig3e = fft(zero_sig3e);
plot(8000/512:8000/512:8000,fft2_sig3e,"r");

fft_sig3o = fft(sig3o);
log_sig3o = log(abs(fft_sig3o));
ifft_sig3o = ifft(log_sig3o);
zero_sig3o = ifft_sig3o;
zero_sig3o(12:512-12) = 0;
fft2_sig3o = fft(zero_sig3o);
plot(8000/512:8000/512:8000,fft2_sig3o,"r");