RILT fmincon

版本 1.0.3 (33.7 MB) 作者: François Liénard
Compute the Regularized Inverse Laplace Transform of an autocorrelation function with the fmincon function. Similar to CONTIN.
0.0
(0)
34.0 次下载
更新时间 2022/11/14

查看许可证

The purpose of these codes is to compute the Regularized Inserve Laplace Transform (RILT) of an autocorrelation function. This is an update of CONTIN (S.W.Provencher, 1982) that was encoded in Fortran.
A similar code has been written by Iari-Gabriel Marino (2022). rilt (https://www.mathworks.com/matlabcentral/fileexchange/6523-rilt).
But this former code was based on fminsearch.
Here, we use the fmincon function that is more suitable and faster.
We use a logarithmic scale but it can be easily changed to a linear one.
'Main_file' is the main file to compute the RILT, and all other files are called from this one.
'Alpha' is a weight on the constraint and you may change it according to your data.
'RILT_fmincon.m' takes around 1s to run with Ns = 200. The complexity cost is more than linear.
'rosenbrock.m' contains the function to minimize.
In 'Data_Folder' you can find a toy signal ('filename.mat') that was built with an AR(1) process.
This algorithm is limited to exponential-like decreasing ACFs but it can be adapted to other kinds of ACFs through the matrix 'A'.
A more general article around this topic is:
'A multiscale time-Laplace method to extract relaxation times from non-stationary dynamic light scattering signals' (2022). https://doi.org/10.1063/5.0088005
This work is provided by François Liénard, Eric Freyssingeas and Pierre Borgnat (ENS de Lyon).

引用格式

François Liénard, Eric Freyssingeas, Pierre Borgnat (2022). RILT fmincon (https://www.mathworks.com/matlabcentral/fileexchange/<...>), MATLAB Central File Exchange. Retrieved October 28, 2022.

MATLAB 版本兼容性
创建方式 R2022a
兼容任何版本
平台兼容性
Windows macOS Linux
标签 添加标签
致谢

参考作品: rilt

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

RILT_fmincon

版本 已发布 发行说明
1.0.3

The fmincon algorithm that is used is 'interior-point' and not 'trust-region-reflective' as it was previously mentioned.
1.0.2

Similar to CONTIN
1.0.1

Similar to CONTIN.
1.0.0