exp.sweep and impulse response

[sweep,invsweep]=expsweep(T,f1,f2,silence,fs)
generate sine sweep with exponential frequency dependent energy decay
over time and time/frequency inverse for impulse response calculation
inputs:
T = sweep duration in seconds
f1 = start frequency in Hz
f2 = end frequency in Hz
silence = silence before and after sweep in seconds (default: 0)
fs = sampling frequency (default: 44100)

[h,h_dist]=sweepIR(rec,Nimp,T,f1,f2,offset,fs)
calculate impulse response and distortion from exponential sweep recording
inputs:
rec = mono sweep recording
T = sweep duration in seconds
f1 = start frequency in Hz
f2 = end frequency in Hz
offset = offset length before impulse response in samples
outputs:
h = impulse response
h_dist = impulse responses of the first 6 nonlinear distortion harmonics