Maxima泰勒级数展开中生成不必要虚数单位%i的问题求助

问题复现 #

运行以下Maxima代码：

rr:sqrt(1-a*cos(θ)^2);
taylor(rr, θ, 0, 2);

得到的结果包含不必要的虚数单位%i：

sqrt(a-1)*%i-(sqrt(a-1)*%i*a*θ^2)/(2*a-2)+...

而实际需要的实数域结果为：

sqrt(1-a)-(sqrt(1-a)*a*θ^2)/(2*a-2)+...

已尝试assume(a<1)但未生效。

解决方案 #

1. 先设置实数域再添加假设 #

Maxima默认可能未严格遵循实数域规则，先指定实数域再添加a<1的假设，重新计算即可得到实数形式结果：

domain: real;
assume(a < 1);
rr:sqrt(1-a*cos(θ)^2);
taylor(rr, θ, 0, 2);

2. 提取泰勒展开结果的实部 #

如果上述方法未生效，直接对展开结果使用realpart()函数，可自动消除虚数单位并转化为等价实数表达式：

assume(a < 1);
rr:sqrt(1-a*cos(θ)^2);
realpart(taylor(rr, θ, 0, 2));

3. 用radcan()化简结果 #

结合a<1的假设，对泰勒展开结果应用radcan()函数，可自动化简掉虚数单位：

assume(a < 1);
rr:sqrt(1-a*cos(θ)^2);
radcan(taylor(rr, θ, 0, 2));

原理说明 #

Maxima默认优先处理复数域的根式分支，即使添加了a<1的假设，也需要明确指定实数域或通过化简函数触发实数规则。设置domain: real;会强制Maxima在实数域内运算，从根源避免虚数单位的引入。

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