Sep21-11, 07:11 PM
I'm having a problem with backward finite difference method. I'm dealing with this problem using Mathematica software (and I'm going to paste some code below, maybe you don't exactly know the syntax, but it's pretty easy I think)..
I need to use the backward finite difference method with the Grunwald-Letnikov fractional derivative definition (probably it's an uncommon topic..).
Long story short, G-L definition for the fractional derivative is this(two links here):
and can be usually used for backw.FD in order to approximate some fractional differential equation.
This is my problem.
I have a first order mass-spring differetial equation and and I have to write a code in Mathematica to approximate it (with backwFD method) and then compare it with the fractional differential equation case. (ps: the same problem appears in the second order equation (mass-spring-damper) approximation, but let's just limit this on the first order).
The problem is that, for the ordinary case I have no problem in using FD method, instead I have a big problem in using FD method for the fractional case: I'll write some code to be more precise.
first order - ordinary
first order - fractional
everything as before, but with this FD cicle
The first code (ordinary case) works flawlessly, but the second doesn't: the curve is completely different depending on the "deltah" imposed (it gives no problems when deltah=1 , but smaller values are generally used in approximations...and that's where this code gives problems).
Now, i've tried quite everything with this code but to no avail...the problem is for sure given by deltah but i cannot find a way.
It would be awesome if you could give me some hints on this. I think you could give me some hints also if you don't know Mathematica.
PS: Since it's my first post here (by the way...this forum is awesome and I say Hi to everybody :) ), please excuse me if I'm posting it in the wrong section or if I'm going against some netiquette.
Thanks a lot in advance.
|Register to reply|
|Finite Difference Discretization of a Fourth Order Partial Differential Term||Differential Equations||4|
|Difference between Forward and Backward Fourier Transforms?||Calculus||1|
|Numerical differentiation using forward, backward and central finite difference||Differential Equations||2|
|System of ODE Boundary Value Problem with 2nd Order Backward Difference||Differential Equations||1|
|Definition of Fractional Occupation||Classical Physics||0|