Maple 2023中plot()函数处理复合周期函数异常，求解决方案

问题场景 #

我需要创建一个在0和1之间振荡的周期函数，已实现区间(0,7)内符合需求的振荡函数weekend_function；为实现周期性，编写了修正后的modulo函数（意图解决原生modulo的不连续问题），但将其与原函数复合后，绘图出现以下异常：

• 直接调用plot(fun)时图形显示正确，但修改图像缩放比例后图形损坏；
• 调用plot(fun(t))无法正常工作，仅显示基础线性函数；
• 带参数调用plot也无法正常运行。

相关代码：

sigmoid := t -> 1/(1 + exp(-t));
weekend_function := t -> piecewise(0 <= t and t < 2.5, 0, 
                                  2.5 <= t < 3.5, sigmoid(12*t - 36), 
                                  3.5 <= t < 4.5, 1 - sigmoid(12*t - 48), 
                                  4.5 <= t, 0);
modulo := t -> (floor(t) mod 7) + t - floor(t);
fun := t -> weekend_function(modulo(t));
plot(modulo);
plot(fun);

问题分析 #

1. 自定义modulo函数的表达式结构导致Maple绘图引擎解析歧义：该函数逻辑上等价于非负t下的t mod 7，但floor(t)与mod的组合让Maple难以识别其连续分段特性，缩放时采样逻辑出错。
2. plot(fun(t))未指定变量区间：Maple将t视为符号变量，无法触发数值采样，仅返回表达式的线性框架。
3. 默认采样点数不足：直接plot(fun)依赖自动采样，点数不足以支撑缩放后的细节重建。

解决方案 #

1. 替换为原生mod函数 #

原生t mod 7在非负实数范围内是连续的，完全可以替代自定义modulo，避免解析问题：

fun := t -> weekend_function(t mod 7);

2. 强制指定绘图区间与采样点数 #

明确变量范围并设置足够的采样点，确保绘图引擎捕捉所有分段细节，缩放后图形不会损坏：

plot(fun, t=0..14, numpoints=2000);

3. 正确调用plot(fun(t)) #

必须指定绘图的变量区间，触发数值采样：

plot(fun(t), t=0..14, numpoints=2000);

4. 优化自定义modulo（若需保留） #

将自定义函数改写为更清晰的分段形式，帮助Maple解析：

modulo := t -> piecewise(t >= 0, t - 7*floor(t/7), 
                           t < 0, t - 7*ceil(t/7));

验证代码 #

修改后的完整可运行代码：

sigmoid := t -> 1/(1 + exp(-t));
weekend_function := t -> piecewise(0 <= t and t < 2.5, 0, 
                                  2.5 <= t < 3.5, sigmoid(12*t - 36), 
                                  3.5 <= t < 4.5, 1 - sigmoid(12*t - 48), 
                                  4.5 <= t, 0);
# 使用原生mod实现周期性
fun := t -> weekend_function(t mod 7);

# 正确绘图方式
plot(fun, t=0..14, numpoints=2000);
# 或调用表达式并指定区间
plot(fun(t), t=0..14, numpoints=2000);

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