Volterra integral equation solver

I am looking for good (fast and accurate, with recent state of art methods) numerical (not symbolic) Volterra Integral Equation (1st kind especially) solver. A few solvers available on FEX are only based on very naive and ineffective methods.
Thanks in advance for any help.
Motivation note:
I would like to numerically solve the following inverse (Linear Dynamic System) problem:
x'(t) = Ax(t) + Bu(t)
y(t) = Cx(t); x(0) = x0
where A,B,C are known constant matrices and y(t) is known (measured) noisy output discrete signal. The goal is estimation of input signal u.
There exist an analytical form of solution y(t):
y(t) = Integral_0^t (C*exp(A*(t-s))*B*u(s) ds = Integral_0^t (K(s,t)*u(s) ds) [*]
where kernel is K(s,t) = C*exp(A*(t-s))*B, so the eq [*] has the form of Volterra intergral equatioon of first kind.
Question: Is there any other possibility how to effectivelly and robustly estimate uknown input signal u(t)?

8 Comments

Torsten
Torsten on 31 Jan 2023
Edited: Torsten on 31 Jan 2023
Do you know of any (also not written in MATLAB) ? I don't.
I just found this Python solver. But I still trying to find something more "robust and effective".
Torsten
Torsten on 31 Jan 2023
Edited: Torsten on 31 Jan 2023
In the old days, the NAG library was the standard for numerical computations:
I don't know if in times of MATLAB and MATHEMATICA, it can still compete.
Of course, MATHEMATICA is another option:
@Torsten NAG Library ... too far from MATLAB implementation
@Torsten NAG has only one function for Volterra Integral Equation (1st kind), which is based on very specific weakly-singular kernel (Abel) type.
I don't have experience with integral equations, but aren't there more listed here:
?
I'm also not an expert at all on the subject of first kind integral equations, having not touched them for nearly 40 years. But my memory does include the words ill-posed. In turn, that suggests finding an off the shelf solver that will be robust, efficient, etc., will be a difficult thing to hope for. I would contact someone with some expertise in the area.

Sign in to comment.

Answers (0)

Categories

Find more on Mathematics in Help Center and File Exchange

Products

Asked:

on 31 Jan 2023

Commented:

on 3 Feb 2023

Community Treasure Hunt

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

Start Hunting!