我如何在不影响我的（x-a）术语的情况下获得“a”处的函数的导数任何建议？答案

【问题标题】：how do i get the derivatives of the functions at "a" with out compromising my (x-a) terms any suggestions?我如何在不影响我的（x-a）术语的情况下获得“a”处的函数的导数任何建议？
【发布时间】：2014-06-04 15:32:03
【问题描述】：

我需要这个函数来输出一个以 a 为中心的 P n 次泰勒多项式这是我拥有的代码，它工作得几乎完美，但我不确定如何在不影响我的（x-a）术语的情况下在“a”处获得函数的导数...... 从我看到的 diff(f,k) 计算 f 的第 k 个导数，但不能插入在一个。看起来我的代码需要另一个 matlab 函数来完成任何建议？

function [P] = mytaylor(f,a,n) 
  f = sym(f);
  syms x;
  terms = 0;
  for k = 0:n
    fk = diff(f,k);
    terms = terms + fk*(x-a)^(k)/factorial(k);   
  end
  P = terms;
end

【问题讨论】：

标签： matlab function taylor-series

【解决方案1】：

您需要使用subs 函数在a 处评估您的符号表达式fk。

【讨论】：

应该进入for循环还是之后？
您必须先取导数，然后才能对其进行评估。由于k 会随着您的循环而变化（因此您的派生fk 会发生变化），因此您需要在循环中执行此操作。我只想将一行更改为terms = terms + subs(fk,x,a)*(x-a)^(k)/factorial(k);