如何在非[-1,1]区间（如[-π,π]）重新计算切比雪夫序列常数用于函数逼近？

核心原理 #

切比雪夫多项式的原生区间是[-1,1]，要适配到目标区间[a,b]（比如[-π,π]），只需要做线性变量替换：

• 将目标区间的变量x映射到[-1,1]的变量t：t = 2*(x - (a+b)/2)/(b - a)
• 反向映射（从t转回x）：x = (b-a)*t/2 + (a+b)/2

所有基于[-1,1]的切比雪夫计算（节点、系数、多项式）都可以通过这个映射转换到目标区间。

针对循环代码的具体修改步骤 #

假设你的原代码是计算切比雪夫逼近系数或插值多项式，按以下步骤修改：

1. 定义目标区间参数
   先明确目标区间的上下限：

   a := -Pi;  # 区间左端点
   b := Pi;   # 区间右端点
   

2. 实现映射函数
   定义双向映射方便后续调用：

   # t∈[-1,1] → x∈[a,b]
   x_from_t := t -> (b - a)*t/2 + (a + b)/2;
   # x∈[a,b] → t∈[-1,1]
   t_from_x := x -> 2*(x - (a + b)/2)/(b - a);
   

3. 修改系数计算逻辑
   如果原代码是通过积分计算切比雪夫系数，将原函数f(x)替换为f(x_from_t(t))，保持积分区间仍为[-1,1]，同时保留切比雪夫的权重项：

   n := 5;  # 逼近阶数
   C := Array(0..n);  # 存储系数
   
   for k from 0 to n do
       # 切比雪夫积分权重：k=0时为1/π，其余为2/π
       weight := piecewise(k=0, 1/Pi, 2/Pi);
       # 计算系数：变量替换后积分
       C[k] := weight * evalf(Int( f(x_from_t(t)) * ChebyshevT(k, t) / sqrt(1 - t^2), t=-1..1 ));
   end do;
   

4. 生成目标区间的逼近多项式
   将切比雪夫多项式的变量t用t_from_x(x)替换，得到基于x的多项式：

   # 构造逼近多项式
   approx_poly := sum(C[k] * ChebyshevT(k, t_from_x(x)), k=0..n);
   # 可选：展开为普通多项式形式
   approx_poly_expanded := expand(approx_poly);
   

插值节点的修改（如果用节点插值方式） #

如果你的原代码是基于切比雪夫节点做插值，只需将[-1,1]的节点映射到目标区间：

n := 5;  # 节点数量
# 生成[-1,1]上的切比雪夫节点
t_nodes := [seq( cos( (2*i-1)*Pi/(2*n) ), i=1..n )];
# 映射到[-π,π]区间的节点
x_nodes := map(x_from_t, t_nodes);
# 生成插值多项式
interp_poly := Interpolate(x_nodes, [seq(f(x), x=x_nodes)], x);

内容的提问来源于stack exchange，提问作者jeff_dust