function soundout = upsample( fname, p, fnameo)
%+++++++++++++++++++++++++++++++++++++
%++ Reading digital sound

%-- subsampling factor
%p = 1;
q = 1; 

[sound, fs, bits]=wavread( fname ); 
slength = length(sound);
disp(['Processing ' fname ]);
disp(['            ' num2str( bits) ' bits']);
disp(['            ' num2str( fs ) ' Hz']);
disp(['            ' num2str( slength / fs) ' s']);

soundout = resample(sound,p,q, 30);
newfs = fs*p/q;

disp(['Resampled ']);
disp(['            ' num2str( bits) ' bits']);
disp(['            ' num2str( newfs ) ' Hz']);
disp(['            ' num2str( length( soundout) / newfs) ' s']);

if nargin > 2
  disp(['Writing file ' fnameo]);
  wavwrite( soundout, newfs, bits, fnameo)
end