向前,向后和中央的差异

1.043视图(30天)

显示旧的评论

迈克阁楼 2015年4月27日

2
链接

这个问题直接联系

https://it.mathworks.com/matlabcentral/answers/213823-forward-backward-and-central-differences

评论道: Stephen Owino Omondi2022年10月1日

答:接受穆罕默德Abouali

嗨请我试图区分这个函数:e ^ y (x) = (- x) *罪(3倍),使用前进、后退和中央差异利用从x = 0到x = 101点4。和情节估计和实际功能衍生品。这是我的代码:

f = @ (x) exp (- x) *罪(3 * x);%的实际函数的导数fprime = @ (x) exp (- x) *罪(3 * x) + 3 * exp (- x) * cos (3 * x);

                         %步长:
                        
                         h = 0.04;

%向前差分dfdx_forward = (f (2 + h) - f (2)) / h Error_forward = fprime (2) -dfdx_forward %错误

                         % bacward区别
                        
                         dfdx_backward = (f (2) - f (2 h)) /小时
                        
                         Error_backward = fprime -dfdx_backward (2)%的错误
                        
                         %中心差分
                        
                         dfdx_central = (f (2 + h) - f (2 h)) / (2 * h)
                        
                         Error_central = fprime -dfdx_central (2)%的错误

请让我知道如果这是正确的,我错误的地方

0评论
显示1年长的评论藏1年长的评论

登录置评。

在回答这个问题。

接受的答案

穆罕默德Abouali 2015年4月27日

17
链接

直接链接到这个答案

https://it.mathworks.com/matlabcentral/answers/213823-forward-backward-and-central-differences answer_176911

                             有趣= @ (x) exp (- x)。*罪(3 * x);
                            
                             dFun = @ (x) exp (- x)。*罪(3 * x) + 3 * exp (- x)。* cos (3 * x);
                            
                             x = linspace (0, 4101);
                            
                             F (x) =乐趣;
                            
                             h = x - x (1) (2);
                            
                             xCentral = x (2: end-1);
                            
                             dFCenteral = (F(3:结束)- F (1: end-2)) / (2 * h);
                            
                             xForward = x (1: end-1);
                            
                             dFForward = (F(2:结束)- F (1: end-1)) / h;
                            
                             xBackward = x(2:结束);
                            
                             dFBackward = (F(2:结束)- F (1: end-1)) / h;
                            
                             情节(x, dFun (x));
                            
                             持有在
                            
                             情节(xCentral dFCenteral,“r”)
                            
                             情节(xForward dFForward,“k”);
                            
                             情节(xBackward dFBackward,‘g’);
                            
                             传奇(“分析”,“中央”,“前进”,“落后”)