Trying to make an Adams-Bashforth method with Richardson error estimate

3 次查看（过去 30 天）

显示更早的评论

ce 2022-12-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1875187-trying-to-make-an-adams-bashforth-method-with-richardson-error-estimate

编辑： Mayank Sengar 2023-1-12

在 MATLAB Online 中打开

%This program solves the initial value problem
%   y' = f(x,y), x0 <= x<= b, y(x0)=y0
%Initializing vaiables 
 %f'(x,y)= 
 f = @(x,y) cos(y).^2;              %derivative in question
 g = @(x) atan(x);                  %this is the actual solution
 x0 = 0;                             %initial value of x
 x_end = 10;                         %end of approximation
 h = 0.1;                           %size of decimal place (0.1,0.001,etc
 y0=0;                               %initial value of y
 
  n = fix((x_end-x0)/h)+1;
 x = linspace(x0,x_end,n);
 y = zeros(n,1);
 y(1) = y0;
 f1 = f(x(1),y(1));
 y(2) = y(1)+h*f1;
 
 %need to add error
 for i = 3:n 
     f2 = f(x(i-1),y(i-1));
     y(i) = y(i-1)+h*(3*f2-f1)/2;
     f1 = f2;
    
 fprintf('%5.4f  %11.8f\n', x(i), y(i));
 plot(x(i),y(i),'b.'); grid on; 
 fplot(g,[x0,x_end]);
 xlabel('x values'); ylabel('y values');
 hold on;
 end
 
 

I'm not sure how I would add the Richardson error to this code. I see the formula in my textbook, but don't understand how I would make it work.

. Like I don't really know what that means. I understand the AB method for solving DefEqs, but not ther errors

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Torsten 2022-12-10

编辑：Torsten 2022-12-10

I'm confident that after reading this article

https://en.wikipedia.org/wiki/Richardson_extrapolation

you will know how Richardson extrapolation works.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Mayank Sengar 2023-1-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1875187-trying-to-make-an-adams-bashforth-method-with-richardson-error-estimate#answer_1147005

编辑：Mayank Sengar 2023-1-12

Go through https://en.wikipedia.org/wiki/Richardson_extrapolation

for understand Richardson expoitation, here is the pseudocode for it:

https://ece.uwaterloo.ca/~dwharder/NumericalAnalysis/12Differentiation/richardson/matlab.html#:~:text=Matlab%20may%20be%20used%20as,f(x)%20is%20defined.

There is also a generalized Richardson extrapolation routine in file exchange: https://www.mathworks.com/matlabcentral/fileexchange/24388-richardson-extrapolation

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

类别

MATLAB Mathematics Numerical Integration and Differential Equations Ordinary Differential Equations

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Translated by