
% plot replica time series
% mbp 8/96, Universidade do Algarve

clear all

% set up source spectrum

fmin = 380;
fmax = 420;
T = 17.5;
fs = 1200;

[ s, time ] = lfm( fmin, fmax, T, fs );

s = s / sqrt( mean( s.^ 2 ) );  % scale so that RMS power is unity
s = s / 2;    % knock it down by 6 dB

% hanning weighting to reduce ripple
%s = s .* ( hanning( length( s ) ).^(1/4) )';

sout = zeros( 30 * fs, 1 );
sout( 1: length( s ) ) = s;

wavwrite( sout, fs, 'lfm_in' )
