No BSD License  

Highlights from
Sinc Resample

4.33333

4.3 | 3 ratings Rate this file 15 Downloads (last 30 days) File Size: 7.32 KB File ID: #12268

Sinc Resample

by Malcolm Lidierth

 

14 Sep 2006 (Updated 22 Feb 2008)

SincResample resamples signals to arbitrary lengths/frequencies and/or timeshifts them

| Watch this File

File Information
Description

If the input is a bandlimited signal containing components at DC - 0.5Fs where Fs is the sample rate, SINCRESAMPLE will return the signal that would have been seen with a higher sampling rate (and the same filter settings). The algorithm depends on the fact that the bandlimited signal is completely defined when sampling at Fs.

SINCRESAMPLE returns the data convolved with a set of time-shifted windowed sinc functions, one for each of the samples [1..size(x,1)] in the input signal. End-effects are reduced by adding a reflected and mirrored copy of the data at the beginning and end of each column before resampling.
Note that points at the end of Y (and/or beginning if SHIFT is negative) will be extrapolated beyond the boundaries of X.
 
SINCRESAMPLE is a generalization of the method discussed in:
T.Blanche & N.V.Swindale (2006)
Nyquist interpolation improves neuron yield in multiunit recordings
Journal of Neuroscience Methods 155, 81-91
http://dx.doi.org/doi:10.1016/j.jneumeth.2005.12.031

MATLAB release MATLAB 7.2 (R2006a)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (6)
17 Apr 2007 Mike eMike  
13 Nov 2007 cho yong-soo

I am korean user that is interested in this file

06 Mar 2008 karan vvv

ok

08 Jun 2008 James Carraeo

Useful, thanks.

05 Jul 2008 Ted Pavlic

The description of this script makes it sound like an improved version fo the INTERPFT function (which just calls fft, pads, and then calls ifft).

I put together a modified INTERPFTW that allows you to resample aliases from different Nyquist windows. In other words, instead of padding the FFT purely at its ends, it splits the padding so that some is at DC and some is at the ends (i.e., it uses a bandpass filter rather than a low-pass filter).

http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=20524&objectType=file

I imagine that you could make a similar modification to your script.

08 Jul 2010 kk KKsingh

It will be much better if it is posted with some example

Please login to add a comment or rating.
Updates
22 Feb 2008

Improved end-effect handling. Lanczos mode added.

Tag Activity for this File
Tag Applied By Date/Time
dsp Malcolm Lidierth 22 Oct 2008 08:39:18
resample interp lanczos interpolation Malcolm Lidierth 22 Oct 2008 08:39:18

Contact us at files@mathworks.com