Hello all,
Suppose I am solving an IVP using one of MATLAB solvers, and the
solution of such problem is expected to oscillate with time. My
interst is to find the averge or mean value of one of the solution
components over time. One straightforward way to do this is to solve
the problem over large time interval in order to get rid of the effect
of the distrubance the will result from the initial condition, and
then use MATLAB function mean to find the averge value.
This technique, however, seems for me ineffiecent since there should
be no need to repeat the simulation over and over as long as the
solution is perodic and everything would be repeated. So, I am sharing
my thoughts with you to see if people around have ideas or tricks to
do this.
One thing I am thinking of is SOMEHOW to measure the wave-length of
the solution during the simulation and to stop when t reaches 10*wave-
length, for example. But, the diffuculty is that I am doing a
continuation simulation & the problem depends on some parameters, like
Reynolds number, that will change each time we solve the problem and
accordingly the wave length will be different.
I would appreciate your ideas. Thanks for any help in advance,
Tai