Lab 3
ADVANCED DSP
201501107
3) Speech
clc; clear all; close all;
[wav1,f]=audioread('speech1.wav');
[wav2,fs]=audioread('speech2.wav');
wave1=fft2(wav1);
wave2=fft2(wav2);
maxx=max(length(wave1),length(wave2));
wave1(length(wave1)+1:maxx,1)=0;
wave2(length(wave2)+1:maxx,1)=0;
phase1=angle(wave1);
phase2=angle(wave2);
mag1=abs(wave1);
mag2=abs(wave2);
rewave1=mag1.*exp(phase2*1j);
rewave2=mag2.*exp(phase1*1j);
rwave1=real(ifft2(rewave1));
rwave2=real(ifft2(rewave2));
subplot 411
title('Original wave1');
plot(wave1);
subplot 412
title('Original wave2');
plot(wave2);
subplot 413
title('Mag of 1 and Phase of 2');
plot(rwave1);
subplot 414
title('Mag of 2 and Phase of 1');
plot(rwave2);
audiowrite('M_1_P_2.wav',rwave1,fs);
audiowrite('M_2_P_1.wav',rwave2,f);