load_sound.sci

signal procesing toolbox

function [x,fs] = load_sound(name, n0, options)

// load_sound - load a sound from a file.
//
//   [x,fs] = load_sound(name, n0, options);
//
//   Load from a .wav (windows), .au (mac) or .asc (Wavelab)
//   sound file.
//
//   Some cool built in examples include
//   spectrogram.wav, tiger.au, bell.wav
//
// Copyright (c) 2006 Gabriel Peyr巈

options.null = [];
if argn(2)<2
    n0 = [];
end

name = lower(name);

rep = getoptions(options, 'rep', './toolbox_signal/');

// find extension 
I = strfind(name, '.');
if isempty(I)
    ext = [];

        if exist( strcat([rep name '.wav']) )==2
            ext = 'wav';
        end
        if exist( strcat([rep name '.asc']))==2
            ext = 'asc';
        end
        if exist( strcat([rep name '.au']))==2
            ext = 'au';
        end

    if isempty(ext)
        error('Unable to determine extension');
    end
else
    I = I(end);
    ext = name(I+1:end);
    name = name(1:I-1);
end


fs = 22000/2; // sampling rate
if strcmp(ext, 'asc')
        fid = fopen(strcat([rep name]), 'r');
        if fid<0
            error(['Unknown file ' name '.' ext]);
        end
        x = fscanf(fid,'%g');
        fclose(fid);
elseif strcmp(ext,  'wav')
        // MS files
        [x,fs,nbits] = wavread(strcat([rep name '.' ext]));
elseif strcmp(ext, '.au')
        // Sun files
        [x,fs,nbits] = auread(strcat([rep name '.' ext]));
else
    x = load_signal(name,n);
end

if size(x,1)<size(x,2)
    x = x';
end
if size(x,2)>1
    // mono signal
    x = x(:,1);
end
x = x(:);

// specific cropping to remov blank

if strcmp(name, 'bird')
//    sel = [1:9400 6500:9000 16500:18700 26000:28500];
    sel = [1:6000 12500:15000 22500:24000 32500:3400];
    x(sel) = [];
end
if strcmp(name, 'spectrogram')
    x = x(5800:end);
end
if strcmp(name, 'aeio')
    x = x(15000:end);
end
if strcmp(name, 'acha')
    x = x(14000:end);
end
if strcmp(name,'aherohasfallen')
    x = x(5000:2^14+5000);
end
if strcmp(name,'sentence')
    x = x(3300:end);
end
// specific cropping
if strcmp(name, 'tiger')
    x = x(2000:end);
end
if strcmp(name, 'drums')
    x = x(4800:12500);
end

n = size(x,1);
subsampling = getoptions(options, 'subsampling', 1);
if subsampling~=1    
    // sub-sampling 1/10
    fs = fs*subsampling;
    t = linspace(0,1,n);
    ti = linspace(0,1,round(n*subsampling));
    x = interp1( t,x,ti );     
    x = x(:);
end

if not(isempty(n0)) & n0<n
    x = x(1:n0);
end

// rescale to [-1,1]
x = x/max(abs(x));

endfunction