matlab求微分方程组数值解,科学网-用MATLAB编写预估-校正法程序求分数阶常微分方程组数值解-王福昌的博文...

主函数比较

% 根据文献自己编写

q = 0.5; %分数阶阶数

fdefun = @(t,y) [2/gamma(3-q)*t^(2-q)-1/gamma(2-q)*t^(1-q)-y+t^2-t]; %一元微分方程

y0 =0;%初值y0为列向量

h = 2^(-6);%步长

tspan = [0,2];

[t,y] = fdewfc(q,fdefun,tspan,y0,h);

figure(1)

plot(t,y(1,1:end)) ;

xlabel('t'); ylabel('y(t)');

hold on

plot(t,t.^2-t+0.5,'r-.')%解析解

%plot(t,-2./(t+1),'r')

%%%%%%与意大利 Roberto Garrappa, University of Bari, Italy 结果比较

t0=tspan(1); tfinal = tspan(2);

[t, y_fde12] = fde12(q,fdefun,t0,tfinal,y0,h);

plot(t,y_fde12+1,'k.')

Jfdefun = @(t,y)[-1];%Jacobi阵 此处退化为导数

[t, y_flmm2] =