% script to demonstrate wind Fs = 44100; t = [0:1/Fs:1-1/Fs]; amp = .1; % amplitude wdm = 10; % gating window in milliseconds % % noise % y = amp * randn(1,Fs); % % sine wave % freq = 440; % y = amp * sin ( 2*pi * freq * t ); % % scanner noise % y = wavread('scanner.wav'); % y = y'; % gate sound y1 = wind(wdm,Fs,y); sound(y1,Fs) subplot(2,1,1) plot(t,y1) subplot(2,1,2) specgram(y1,256,Fs)