写一段注释清晰的matlab代码，采用Tchebychev 截断级数的办法求出f（x）的近似最佳致逼近三次多项式

好的，以下是注释清晰的matlab代码：

% 定义函数
f = @(x) x.^2 + x.*sin(x) + cos(2*x);

% 定义区间及步长
a = -2;  % 区间左端点
b = 5;   % 区间右端点
n = 1000;  % 步长
h = (b-a)/n;

% 求出数据点
x = a:h:b;  % 等距节点
y = f(x);

% 计算Tchebychev点
n = 4;  % Tchebychev点的个数
t = cos((2*(1:n)-1)*pi/(2*n));  % Tchebychev点

% 将Tchebychev点转换为实际数据点
x_t = (b+a)/2 + (b-a)/2*t;  % Tchebychev点映射为实际数据点

% 计算Tchebychev系数
c = zeros(1, n);  % Tchebychev系数
for k = 1:n
    c(k) = sum(y.*Tchebychev(k-1, (2*x_t-a-b)/(b-a)))/n;  % 计算Tchebychev系数
end

% 计算近似最佳平方逼近三次多项式
p = @(x) c(1)/2 + c(2)*Tchebychev(1, (2*x-a-b)/(b-a)) + ...
    c(3)*Tchebychev(2, (2*x-a-b)/(b-a)) + c(4)*Tchebychev(3, (2*x-a-b)/(b-a));

% 作图比较
x_plot = linspace(a, b, 1000);
y_plot = f(x_plot);
p_plot = p(x_plot);

plot(x_plot, y_plot, 'b-', x_plot, p_plot, 'r--');
xlabel('x');
ylabel('y');
title('f(x)及其三次多项式逼近');
legend('原函数', '三次多项式');

这段代码实现了采用Tchebychev截断级数的方法求出f（x）的近似最佳致逼近三次多项式。代码中通过注释清晰地解释了每一步的操作，方便理解和修改。其中函数f（x）可以根据需要自行更改。