无需Symbolic Toolbox构造可积分多项式函数的技术问询

核心问题分析 #

你定义的匿名函数无法被integral执行，大概率是向量化处理不严谨或n的取值未匹配系数数组长度导致的。integral要求被积函数接收向量输入时，返回同维度的向量/数组，你的写法可能在维度匹配或转置处理上存在隐性问题。

最优解决方案：使用内置polyval函数 #

Matlab的polyval是专门为多项式求值设计的函数，天生支持向量化输入，完美适配integral，且代码简洁不易出错。

假设系数数组C=[c_n, c_{n-1}, ..., c_0]对应多项式c_n x^n + c_{n-1}x^{n-1} + ... + c_0，直接用polyval构造函数：

g = @(x) polyval(C, x);

此时调用integral(g, min_val, max_val)完全正常。如果要计算g(x)*f(x)的积分（f(x)为通用函数），只需将两个函数点乘（确保f(x)也支持向量化）：

% 示例：C对应2x²+3x+1，f(x)=exp(-x)，积分区间[0,1]
C = [2, 3, 1];
f = @(x) exp(-x);
integrand = @(x) polyval(C, x) .* f(x);
result = integral(integrand, 0, 1);

手动实现多项式求值的修正方案 #

如果不想用polyval，可以修正你的原写法，确保维度匹配和向量化正确：

1. 确保n的取值为length(C)-1（系数数组长度减1就是最高次幂）
2. 使用非共轭转置C.'代替transpose(C)（避免复数场景下的意外共轭）

n = length(C) - 1;
g = @(x) x.^(n:-1:0) * C.';

此时x为向量时，x.^(n:-1:0)生成m×(n+1)矩阵（m是x的长度），与C.'（列向量）相乘后得到m×1列向量，符合integral的输入输出要求。

关键注意点 #

• 所有参与运算的函数（包括你的f(x)）必须支持向量化：即输入向量x时，返回同长度的向量，不能只处理单个标量。
• 优先使用polyval，它是Matlab优化过的内置函数，能处理各种边界情况（比如x=0时的高次幂计算），代码可读性和可靠性更高。

内容的提问来源于stack exchange，提问作者José Pedro Ferreira