每卡普托实现分数阶导数的定义

17日视图(30天)

显示旧的评论

莫斯科Martinelli 2022年10月21日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1832283-implements-a-fractional-order-derivative-per-caputo-s-definition

评论道: 深2023年5月13日

答:接受 Torsten

我有下面的代码实现了一个每卡普托分数阶导数的定义。我把它从一个稍微老书,它看起来像代码不是最新的,MATLAB抛出一个错误说10号线变量“dy”变化大小与每个迭代。我想通过一个方法来修正这个错误或代码更新,因为我需要计算分数阶导数的函数。

                          函数dy =卡普托(t0, f, gam)
                         
                          m =地板(gam);一个= gam-m;dy = 0;
                         
                          如果gam > 0,信谊t;fd = diff (f t m + 1);其他的;一个= gam;结束
                         
                          为t1 = t0(2:结束)
                         
                          如果gam > 0
                         
                          f = @ (x)潜艇(fd t x) / (t1-x)。^γ(1 a);
                         
                          其他的
                         
                          f = @ (x)潜艇(f t x) / (t1-x)。^ (1 a) /γ(一个);
                         
                          结束
                         
                          dy = [dy;积分(f, 0, t1)];
                         
                          结束
                         
                          结束
                         
                          > > sym t;f =罪(3 * t + 1);t0 = 0:0.01:5;dy =卡普托(t0, f, 0.3);情节(t0) dy)