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 件のコメント

Torsten
Torsten 2023 年 1 月 31 日
編集済み: Torsten 2023 年 1 月 31 日
Do you know of any (also not written in MATLAB) ? I don't.
Michal
Michal 2023 年 1 月 31 日
I just found this Python solver. But I still trying to find something more "robust and effective".
Torsten
Torsten 2023 年 1 月 31 日
編集済み: Torsten 2023 年 1 月 31 日
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:
Michal
Michal 2023 年 1 月 31 日
@Torsten NAG Library ... too far from MATLAB implementation
Torsten
Torsten 2023 年 1 月 31 日
編集済み: Torsten 2023 年 1 月 31 日
Michal
Michal 2023 年 1 月 31 日
@Torsten NAG has only one function for Volterra Integral Equation (1st kind), which is based on very specific weakly-singular kernel (Abel) type.
Torsten
Torsten 2023 年 1 月 31 日
I don't have experience with integral equations, but aren't there more listed here:
?
John D'Errico
John D'Errico 2023 年 2 月 3 日
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.

サインインしてコメントする。

回答 (0 件)

カテゴリ

ヘルプ センター および File ExchangeMathematics についてさらに検索

製品

質問済み:

2023 年 1 月 31 日

コメント済み:

2023 年 2 月 3 日

Community Treasure Hunt

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

Start Hunting!

Translated by