sig1 - tkyswiki @ ウィキ

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

plot(sig1)
xlim([5000 10000]);
sig1a=sig1(7001:7512);

plot(sig1)
xlim([20000 25000]);
sig1i=sig1(21501:22012);

plot(sig1)
xlim([35000 40000]);
sig1u=sig1(36001:36512);

plot(sig1)
xlim([47000 53000]);
sig1e=sig1(48501:49012);

plot(sig1)
xlim([60000 65000]);
sig1o=sig1(62001:62512);

fft_sig1a = fft(sig1a);
log_sig1a = log10(abs(fft_sig1a));
ifft_sig1a = ifft(log_sig1a);
zero_sig1a = ifft_sig1a;
zero_sig1a(12:512-12) = 0;
fft2_sig1a = fft(zero_sig1a);
plot(8000/512:8000/512:8000,fft2_sig1a,"r");

fft_sig1i = fft(sig1i);
log_sig1i = log10(abs(fft_sig1i));
ifft_sig1i = ifft(log_sig1i);
zero_sig1i = ifft_sig1i;
zero_sig1i(12:512-12) = 0;
fft2_sig1i = fft(zero_sig1i);
plot(8000/512:8000/512:8000,fft2_sig1i,"r");

fft_sig1u = fft(sig1u);
log_sig1u = log10(abs(fft_sig1u));
ifft_sig1u = ifft(log_sig1u);
zero_sig1u = ifft_sig1u;
zero_sig1u(12:512-12) = 0;
fft2_sig1u = fft(zero_sig1u);
plot(8000/512:8000/512:8000,fft2_sig1u,"r");

fft_sig1e = fft(sig1e);
log_sig1e = log10(abs(fft_sig1e));
ifft_sig1e = ifft(log_sig1e);
zero_sig1e = ifft_sig1e;
zero_sig1e(12:512-12) = 0;
fft2_sig1e = fft(zero_sig1e);
plot(8000/512:8000/512:8000,fft2_sig1e,"r");

fft_sig1o = fft(sig1o);
log_sig1o = log10(abs(fft_sig1o));
ifft_sig1o = ifft(log_sig1o);
zero_sig1o = ifft_sig1o;
zero_sig1o(12:512-12) = 0;
fft2_sig1o = fft(zero_sig1o);
plot(8000/512:8000/512:8000,fft2_sig1o,"r");