sig2 - tkyswiki @ ウィキ

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

plot(sig2)
xlim([9000 15000]);
sig2a=sig2(10001:10512);

plot(sig2)
xlim([22000 31000]);
sig2i=sig2(25501:26012);

plot(sig2)
xlim([43000 48000]);
sig2u=sig2(45001:45512);

plot(sig2)
xlim([62000 67000]);
sig2e=sig2(64001:64512);

plot(sig2)
xlim([90000 95000]);
sig2o=sig2(93001:93512);

fft_sig2a = fft(sig2a);
log_sig2a = log(abs(fft_sig2a));
ifft_sig2a = ifft(log_sig2a);
zero_sig2a = ifft_sig2a;
zero_sig2a(12:512-12) = 0;
fft2_sig2a = fft(zero_sig2a);
plot(8000/512:8000/512:8000,fft2_sig2a,"r");

fft_sig2i = fft(sig2i);
log_sig2i = log(abs(fft_sig2i));
ifft_sig2i = ifft(log_sig2i);
zero_sig2i = ifft_sig2i;
zero_sig2i(12:512-12) = 0;
fft2_sig2i = fft(zero_sig2i);
plot(8000/512:8000/512:8000,fft2_sig2i,"r");

fft_sig2u = fft(sig2u);
log_sig2u = log10(abs(fft_sig2u));
ifft_sig2u = ifft(log_sig2u);
zero_sig2u = ifft_sig2u;
zero_sig2u(12:512-12) = 0;
fft2_sig2u = fft(zero_sig2u);
plot(8000/512:8000/512:8000,fft2_sig2u,"r");

fft_sig2e = fft(sig2e);
log_sig2e = log10(abs(fft_sig2e));
ifft_sig2e = ifft(log_sig2e);
zero_sig2e = ifft_sig2e;
zero_sig2e(12:512-12) = 0;
fft2_sig2e = fft(zero_sig2e);
plot(8000/512:8000/512:8000,fft2_sig2e,"r");

fft_sig2o = fft(sig2o);
log_sig2o = log10(abs(fft_sig2o));
ifft_sig2o = ifft(log_sig2o);
zero_sig2o = ifft_sig2o;
zero_sig2o(12:512-12) = 0;
fft2_sig2o = fft(zero_sig2o);
plot(8000/512:8000/512:8000,fft2_sig2o,"r");