I cannot seem to get a second line for the forward finite difference

Question

Michael Corbett 2019-9-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/481960-i-cannot-seem-to-get-a-second-line-for-the-forward-finite-difference

评论： Walter Roberson 2019-9-25

clc
clear
close all
h = [-2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 0 .25 .5 .75 1 1.25 1.5 1.75 2];
x = [-2 -1.75 -1.5 -1.25 -1 -.75 -.5 -.25 0 .25 .5 .75 1 1.25 1.5 1.75 2];
y = @(x) x.^3 - 2*x + 4;   %start function
yd = @(x) 3*x.^2 - 2;      %first deriv
ydd = @(x) 6*x;            %second deriv
for i = 1:2
    for j = 1:17
        xi = x(i);
        hj = h(j);
        
        d1 = yd(xi);
        d2 = ydd(xi);
        
        %forward diff
        f1 = (y(xi+hj) - y(xi))/hj;
        f2 = (y(xi+2*hj) - 2*y(xi+hj))/(hj.^2);
        
        %backward diff
        b1 = (y(xi) - y(xi-hj))/hj;
        b2 = (y(xi)-2*y(xi-hj)+y(xi-2*hj))./(h.^2);
        
        %centered diff
        c1 = (y(xi+hj) - y(xi-hj))/(2*hj);
        c2 = (y(xi+hj) - 2*y(xi)+y(xi-hj))/(hj.^2);
        
        plot(x,yd(x));
        hold on
        plot(x,d2,'r');
        hold off
        legend('Analytic','Forward','Backward','Center')
    end
end