1
$\begingroup$

I know this question mayn't be new, but here is my problem: I have to solve this integral equation (numerically, of course, except from very few special cases):

$$\phi_{\nu}(t) = 1 - q\int_0^t \frac{d\psi_{\nu}}{d t'} \phi_{\nu}(t - t')\ dt'$$

Where of course I know what $\psi$ is. I can choose $q = 1$.

I tried to take a look at some past answers, but I did not find what I was looking for. The problem may be that this is a convolution integral.

$\endgroup$
4
  • $\begingroup$ I recently looked into a similar matrix equation. The proposed solution is that if your model allows it to change the lower boundary of the integral to be $t-c$ with $c$ a constant, then you can discretize the integral using some formula (e.g. Trapezoidal rule) and solve the obtained delay equation. $\endgroup$ Commented Jun 10, 2017 at 13:40
  • $\begingroup$ Perhaps, this question belongs in Mathematics instead of Mathematica.SE, $\endgroup$
    – bbgodfrey
    Commented Jul 12, 2017 at 4:38
  • $\begingroup$ @bbgodfrey Since I wanted to have a mathematica code to solve that, I don't see how it should belong to Math SE, people over there won't provide me for a code in that sense... $\endgroup$
    – Enrico M.
    Commented Jul 12, 2017 at 12:44
  • $\begingroup$ Mathematica certainly can reproduce the calculation by @CraigTracy, but no further progress can be made without an expression for g[s]. Given that expression, it may be possible to use InverseLaplaceTransform to obtain the desired answer. Am alternative approach is to discretize your integral equation to form a matrix equation, which can be inverted to obtain a numerical answer. $\endgroup$
    – bbgodfrey
    Commented Jul 12, 2017 at 13:01

1 Answer 1

3
$\begingroup$

This problem can be solved using the Laplace transform of a function $f$: $$ \hat{f}(s):=\int_0^\infty e^{-s t} f(t)\, dt. $$ Your equation can be written as $$ f(t)=1-q\int_0^t f(t-t') g(t')\, dt'. $$ Taking the Laplace transform of both sides and using the convolution property of the Laplace transform one gets $$ \hat{f}(s) = 1/s - q\, \hat{f}(s)\,\hat{g}(s). $$ Solving this for $\hat{f}$ gives $$ \hat{f}(s)=(1/s) \left(1+q\,\hat{g}(s)\right)^{-1}.$$ Now use the inverse transform to find $f(t)$.

$\endgroup$
0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.