Code covered by the BSD License  

Highlights from
Swept-Sine Analysis

4.0

4.0 | 1 rating Rate this file 29 Downloads (last 30 days) File Size: 3.89 KB File ID: #29187
image thumbnail

Swept-Sine Analysis

by Oygo

 

30 Oct 2010 (Updated 02 Dec 2010)

Acquire the linear impulse response of a device under test

| Watch this File

File Information
Description

This method, known as swept-sine, log-chirp, swept-tone, or sine sweep uses an increasing frequency sinusoid to stimulate an acoustic device under test. The response is processed to compress the response to an impulse response. This method is robust against noise, time variance, and weak nonlinearities.

% SYNTHSWEEP Synthesize a logarithmic sine sweep.
% [sweep invsweepfft sweepRate] = synthSweep(T,FS,f1,f2)
% generates a logarithmic sine sweep that starts at frequency f1 (Hz),
% stops at frequency f2 (Hz) and duration T (sec) at sample rate FS (Hz).

% EXTRACTIR Extract impulse response from swept-sine response.
% [irLin, irNonLin] = extractIR(sweep_response, invsweepfft)
% Extracts the impulse response from the swept-sine response. Use
% synthSweep.m first to create the stimulus; then pass it through the
% device under test; finally, take the response and process it with the
% inverse swept-sine to produce the linear impulse response and
% non-linear simplified Volterra diagonals. The location of each
% non-linear order can be calculated with the sweepRate - this will be
% implemented as a future revision.

MATLAB release MATLAB 7.8 (R2009a)
Tags for This File  
Everyone's Tags
chirp(3), impulse response(2), measurement(2), signal processing(2), sweptsine
Tags I've Applied
Add New Tags Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
01 Jul 2012 Ben

You may want to check this. The frequency spectrum of a logarithmic swept sine should decrease with frequency at a rate of 3 dB/oct. The sweep generated by this code do not follow that trend.

07 Dec 2010 Dongwon

This tool is very usefull and is what I was looking for.
One thing that might be hopefully added is to compute and provide a fft plot of the swept sine analysis, albeit the impulse response might be enough.

Thanks.

Thomas

Updates
02 Nov 2010

Added function call and help file information.

02 Dec 2010

Added grpdelay2phase.m to the submission

Contact us