matlab用代码写泰勒函数（syms）

$$f(x)=f(x_0)+f^{\prime }(x_0)(x-x_0)+\frac{f^{\prime }(x_0)}{2!}(x-x_0{)}^2+\cdots +\frac{f^{\prime }(x_0)}{n!}(x-x_0{)}^n+r_n(x)$$

syms x;  % 声明 x 为符号变量

% 定义待展开的函数和展开的点
f = sin(x);
a = 0;

% 指定展开的阶数
order = 5;

% 初始化展开项的和为 0
taylor_expansion = sym(0);

% 逐项计算展开并累加
for n = 0:order
    % 计算展开项的系数
    coefficient = subs(diff(f, x, n), x, a) / factorial(n);
    
    % 计算展开项并累加到总和中
    term = coefficient * (x - a)^n;
    taylor_expansion = taylor_expansion + term;
end

% 打印结果
 disp(taylor_expansion)
% 
% function result = factorial(n)
%     result = 1;
%     for i = 1:n
%         result = result * i;
%     end
% end

注：
1.在代码中，coefficient = subs(diff(f, x, n), x, a) / factorial(n); 的作用是计算在泰勒展开中的每个展开项的系数。

具体来说，diff(f, x, n) 表示对函数 f 在变量 x 上进行 n 次求导。这里的 diff 函数是用于计算导数的符号函数，可以传入三个参数：函数 f、变量 x 和求导的次数 n。

接着，subs(diff(f, x, n), x, a) 则表示将导数表达式 diff(f, x, n) 中的变量 x 替换为值 a。这里的 subs 函数用于符号代入，可以传入三个参数：表达式、要替换的变量和替换的值。

最后，factorial(n) 表示计算阶乘的值，即将整数 n 的阶乘求出。这里的 factorial 是用于计算阶乘的符号函数，可以传入一个参数：要计算阶乘的数。

最终，通过将导数求值、变量替换和阶乘计算相结合，coefficient 变量会获得在泰勒展开中的每个展开项的系数。