clear all

plp = imread('wuh.jpg');
%imshow(plp)

fs = 44100
L = length( plp(:,1,1))
W = length( plp(1,:,1))
frq = zeros(1, L * W);
amp = zeros(1, L * W);
t = 1:L*W;
t = t/fs;

tp = rgb2hsv(plp);

for i = 1:L-1
    frq(W*i+1:W*(i+1)) = plp(i,:,1); %value = frequency
    %amp(W*i+1:W*(i+1)) = plp(i,:,2); %saturation = amplitude
    
end

frq = frq+440;
amp = 1;

wav = amp .* sin( frq .* t);
soundsc(wav(1:882000),fs)
audiowrite('Wuh1.wav',wav,fs)

Wuh1

clear all

plp = imread('wuh.jpg');
%imshow(plp)

fs = 44100
L = length( plp(:,1,1))
W = length( plp(1,:,1))
frq = zeros(1, L * W);
amp = zeros(1, L * W);
t = 1:L*W;
t = t/fs;

tp = rgb2hsv(plp);

for i = 1:L-1
    frq(W*i+1:W*(i+1)) = plp(i,:,3);%value = frequency
    %amp(W*i+1:W*(i+1)) = plp(i,:,2); %saturation = amplitude
    
end


frq = (440 + frq*300*2*pi;
amp = 1;

wav = amp .* sin( frq .* t);
soundsc(wav(1:882000),fs)
audiowrite('Wuh1.wav',wav,fs)
%soundsc(wav,fs)

Wuh2

clear all

[j, fs] = audioread('Eskim.mp3');
ender = 99

jj = fft(j(:,1));
    aj = abs(jj);
    ij = imag(jj);
    rj = real(jj);
    mi = max(abs(ij));
    wav = zeros(length(j),1);
    step = floor(length(j)/(ender+1))
    fiij = ij > 0;
    fiij = fiij - (1-fiij);
    firj = rj > 0;
    firj = firj - (1-firj);
    
for i = 0:ender
    fij = abs(ij) > i/ender*mi;
    frj = abs(rj) > i/ender*mi;
    ii = fiij.*fij*i/ender*mi + (1-fij).*ij;
    rr = firj.*frj*i/ender*mi + (1-frj).*rj;
    tp = rr + ii;
    ok = ifft(tp);
    wav(i*step+1:(i+1)*step) = ok(i*step+1:(i+1)*step);
end
%plot(amp)
t = 1:length(tp);
wav = real(10*wav);
mx = max(wav(44100:length(wav)-44100))
fi = abs(wav) < mx;
wav = wav.*fi;
wav = wav/max(wav);
audiowrite('Eskim5.wav',wav,fs)
soundsc(wav,fs)

whatwhatwhat!