// At first we create 0.5 seconds of sound parameters. t=soundsec(0.5); // Then we generate the sound. s=sin(440*t)+sin(220*t)/2+sin(880*t)/2; [nr,nc]=size(t); s(nc/2:nc)=sin(330*t(nc/2:nc)); // We can easily make a Fourier analysis of it. ;analyze(s);
// Save the file in WAV format. // we renormalize s in order to check that save+load is invariant s=s-sum(s)/prod(size(s)); s=s/max(abs(s)); savewave("test.wav",s); // Load it back. s1=loadwave("test.wav"); if maxi(abs(s1-s)) < 1.e-4;end // Now we can make a complete picture of the sound. ;mapsound(s);