On Mar 9, 1:50 pm,
spellu...@fb04373.mathematik.tu-darmstadt.de (Peter
Spellucci) wrote:
In article <1173452229.704156.190...@p10g2000cwp.googlegroups.com>, "Freddy" <zfreddy...@gmail.com> writes:
Hi Again,
I was just wondering as Freddy wrote that y should be written as
y(x,t,a) actually the only reason that I said "a" is a matrix is
because I have 4 parameters in my function. a = [a1, a2, a3, a4]
so could that be written as y = y(x,t,a1,a2,a3,a4).
and if so how accurate is this representation
[dy(x,t,a1,a2,a3,a4)/da1] = [y(x,t,a1+eps,a2,a3,a4)-y(x,t,a1-
eps,a2,a3,a4)]/(2.eps)
this would be a finite difference replacement, end assuming everything
else is smooth, with an error of O(eps^2)
but I cannot see what this should help. if you are able to
integrate the equation for dy/dt, then you will also be able
to integrate the variational equation for dy/da , provided you are
able to write down the Jacobian of f with respect to a, but again:
the x on the left hand side, where it is gone?
or have you simply a misprint and y does not depend on x?
in that case the integral on the right would turn over in a factor for y(t)
and anything would become pretty easy. but as you wrote it it is
not well defined, even if oe is willing to read this as an integro differential
equation (then again a variable would be missing, this time under the integral
clearly something is wrong
hth
peter
is that possible or I just deserve to be shot for writing
that..lol...and broken many rules.
Cheers,
-F-
Freddy
I was going to respond to your earlier message, and almost changed my
mind after all the repeats.
One way of getting what you want (if you can do it) is to also solve the
differential equation
d(dy/da)/dt = \int_{xmin}^{xmax} [ (\partial f(x:a)/\partial a) y(x,t) +
f(x:a) dy/da ] dx
with initial conditions dy/da = 0.
Note I have mainly followed your notation, but y more properly should be
written y(x,t,a). If y is a vector then dy/da is a matrix, else dy/da is
a vector. I'm assuming that "a" can be thought of as a vector rather
than as a matrix as you describe it.
Regards,
Fred
with initial conditions (dy/da) = 0.
Actually I'm trying to solve the population balance equation for
flocculation.
so I'm trying to solve for y (number density of particles of a certain
diameter "x" at a time "t")
and the RHS does take into account the diameter, the number density,
and various other functions like the flocculation frequency for
example. and it's on the RHS that we have the model constants that I'm
trying to match to experimental data.
I'm really still not sure what might be missing under the integral
though as you mentioned. But hope that gives an idea about the problem
I'm dealing with.
Thanks a lot for the help you've been providing guys.
Freddy