Le 03/02/2023 à 21:25, Samuel Gougeon a écrit :
Le 03/02/2023 à 20:57, Samuel Gougeon a écrit :
.../...
Here is a draft proposal:
1) build the (let's say row) vector A = (dPHI/dt) of sampled data at
sampled values t
2) build the row vector B = f(-t) of sample data at t values
3) build the matrix
Le 03/02/2023 à 20:57, Samuel Gougeon a écrit :
.../...
Here is a draft proposal:
1) build the (let's say row) vector A = (dPHI/dt) of sampled data at
sampled values t
2) build the row vector B = f(-t) of sample data at t values
3) build the matrix C of (padded) A and the matrix D of (shiffted
Le 03/02/2023 à 11:24, Heinz Nabielek a écrit :
On 03.02.2023, at 11:13, Stéphane Mottelet wrote:
Thanks for the code.
Just a remark on the notations, you should write :
F(T) = Int_{0}^{T} PHI(t) . f(T-t) . dt
i.e. not F(t) since t is mute.
However, you should pay attention to the delay
I can't figure out weither the problem resides in your model or in the
computations, but when I had (many years ago) to use discrete convolution to
approximate continuous convolution, I noticed that the best way to obtain
coherent results is to use the composed midpoint rule to approximate the
I can't figure out weither the problem resides in your model or in the
computations, but when I had (many years ago) to use discrete convolution to
approximate continuous convolution, I noticed that the best way to obtain
coherent results is to use the composed midpoint rule to approximate the
This is my latest code version: the 'convoluted secondary failure fraction' is
nicely below primary failure, but seems too low.
Heinz
m=2; // Weibull modulus in mechanism #1
k=1E-7; // corrosion rate(s-1) in mechanism #1
n=500; // number of time steps in hourly intervals
On 03.02.2023, at 11:13, Stéphane Mottelet wrote:
>
> Thanks for the code.
>
> Just a remark on the notations, you should write :
>
> F(T) = Int_{0}^{T} PHI(t) . f(T-t) . dt
>
> i.e. not F(t) since t is mute.
>
> However, you should pay attention to the delay notion associated with
>
Thanks for the code.
Just a remark on the notations, you should write :
F(T) = Int_{0}^{T} PHI(t) . f(T-t) . dt
i.e. not F(t) since t is mute.
However, you should pay attention to the delay notion associated with
convolution and the relationships between discrete convolution and
continuous