yading@10: clear all
yading@10: %1.1
yading@10: %% Piano sound
yading@10: 
yading@10: % hfile = '4192__RealRhodesSounds__D4.wav';
yading@10: % [x3, Fs3, nbits3, readinfo3] = wavread(hfile);
yading@10: % %sound(x3, Fs);
yading@10: % y = sinemodel_(x3, hamming(513)', 1024, 100,-35);
yading@10: % wavwrite(y, Fs3,'__RealRhodesSounds_D4.wav');
yading@10: 
yading@10: 
yading@10: %% 2.1 Compute the magnitude values on the to of the main lobe
yading@10: % loc=0;
yading@10: % mag=0;
yading@10: % 
yading@10: % binremainder = round(loc)-loc;
yading@10: % %Blackman-Harris 92dB
yading@10: % lb = [binremainder-4:binremainder+4]'; % main lobe (real value) bins to read
yading@10: % lmag = genbh92lobe(lb,1024)*10.^(mag/20); % lobe magnitudes
yading@10: % figure
yading@10: % plot(lmag);
yading@10: 
yading@10: 
yading@10: %% 2.2 Compute a spectrum of size N from a series of sinusoidal values
yading@10: 
yading@10: % Y = genspecsines([10.3, 20.6, 30.9], [-10, -13, -16], [0, 0, 0], 1024);
yading@10: % 
yading@10: % %inverse fft
yading@10: % y=fftshift(real(ifft(Y)));
yading@10: % plot(y);
yading@10: % 
yading@10: % %the sinusoids 
yading@10: % Y1 = (-10)*exp(1*i);y1=ifft(Y1);
yading@10: % Y2 = (-13)*exp(1*i);y2=fftshift(real(ifft(Y2)));
yading@10: % Y3 = (-16)*exp(1*i);y3=fftshift(real(ifft(Y3)));
yading@10: % %add the sinusoids in the time domain
yading@10: % ytime = y1+y2+y3;
yading@10: % figure
yading@10: % plot(y1);
yading@10: 
yading@10: %% Fugue
yading@10: 
yading@10: % hfile = '22125__Connum__memento_to_Bachs_fugue_II.wav';
yading@10: % [x3, Fs3, nbits3, readinfo3] = wavread(hfile);
yading@10: % %sound(x3, Fs);
yading@10: % y = sinemodel(x3, hamming(2001)', 1024, -60);
yading@10: % wavwrite(y, Fs3,'fugue1.wav');
yading@10: 
yading@10: % %% Soprano
yading@10: % 
yading@10: % hfile = '30084__HerbertBoland__FemalePhrase1.wav';
yading@10: % [x3, Fs3, nbits3, readinfo3] = wavread(hfile);
yading@10: % %sound(x3, Fs);
yading@10: % y = sinemodel(x3, hamming(2001), 1024, -50);
yading@10: % wavwrite(y, Fs3,'female1.wav');
yading@10: 
yading@10: %% Soprano
yading@10: 
yading@10: hfile = '32848__Quilt__Quilty_s_Old_School_Piano.wav';
yading@10: [x3, Fs3, nbits3, readinfo3] = wavread(hfile);
yading@10: %sound(x3, Fs);
yading@10: y = sinemodel_(x3, hamming(2001)', 1024, -50);
yading@10: wavwrite(y, Fs3,'piano1.wav');