MATLAB计算积分时出现‘Matrix dimensions must agree’错误求助排查

你遇到的「Matrix dimensions must agree」错误，核心原因是你给a1传入了一个q向量，但integral函数默认是逐标量处理输入的——当integral迭代计算z的积分值时，z是单个标量，而q是你传入的向量，这时候z.^2 + q.^2就会因为标量和向量维度不匹配而报错。

下面给你两种可行的解决方案：

方案1：用arrayfun逐个处理q向量元素 #

这种方法最直接，让a1保持原来的标量输入逻辑，用arrayfun把a1应用到q向量的每一个元素上：

% 定义原函数
a = @(z,q)(z.^2 + q.^2);
% 保持原积分逻辑，接受标量q
a1 = @(q) integral(@(z) a(z,q), 1, 10);

% 生成q的向量，并逐个计算a1的值
q_vec = linspace(0, 3e8);
a1_results = arrayfun(a1, q_vec);

% 绘图
plot(q_vec, a1_results);

方案2：开启integral的向量化支持（更高效） #

如果你想让a1直接支持向量输入，可以给integral加上'ArrayValued', true参数，告诉MATLAB被积分的函数会返回数组，这样就能一次性处理整个q向量：

% 定义原函数
a = @(z,q)(z.^2 + q.^2);
% 修改a1，开启ArrayValued选项，支持q为向量
a1 = @(q) integral(@(z) a(z,q), 1, 10, 'ArrayValued', true);

% 直接传入q向量绘图
q_vec = linspace(0, 3e8);
plot(q_vec, a1(q_vec));

额外验证：解析解对照 #

其实这个积分可以用解析方法计算，对z从1到10积分z²+q²的结果是：
$$\int_{1}^{10} (z^2 + q^2) dz = \left[\frac{z^3}{3} + q^2 z\right]_{1}^{10} = 333 + 9q^2$$
你可以用这个解析结果和代码运行的曲线对比，确保结果正确。

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