Matlab中使用int和integral函数求解多变量积分问题求助

报错原因 #

你之前的方案报错是因为int是符号积分函数，要求输入为符号表达式，直接传入数值匿名函数会触发类型不匹配，符号转数值的额外步骤也容易引入兼容问题，本场景下完全不需要引入符号计算逻辑，直接用纯数值积分即可。

基础实现方案 #

% 定义rho_delta相关参数
p = 3;
delta = 0.1;
A_p = 1/integral(@(s) exp(-s.^2),-p,p);
rho_delta = @(t) A_p*exp(-t.^2/delta^2)/delta;

% 替换为你实际的被乘函数f(s)，此处以f(s)=s为示例
f = @(s) s;
% 替换为你需要的积分上限x
x = 1;

% 构造被积函数，注意使用点运算保证向量化兼容
integrand = @(s) rho_delta(x - s) .* f(s);
% 直接调用integral完成数值积分
result = integral(integrand, 0, x);
disp("积分结果：" + result);

进阶批量计算 #

如果需要对多个x值计算对应积分结果，可结合arrayfun实现：

% 定义x的取值序列
x_list = 0:0.1:2;
% 批量计算所有x对应的积分值
result_list = arrayfun(@(curr_x) integral(@(s) rho_delta(curr_x - s).*f(s), 0, curr_x), x_list);

% 可选：绘制积分结果随x变化的曲线
plot(x_list, result_list, 'LineWidth', 1.5);
xlabel('积分上限x');
ylabel('积分结果');
grid on;

注意事项 #

• 构造被积函数时必须使用点运算符（.*/.^/./），保证integral传入向量化参数时运算正常，否则会触发维度不匹配错误。
• 若被积函数存在奇点、强振荡等特性，可通过增加integral的精度参数优化结果：integral(integrand, 0, x, 'RelTol', 1e-8, 'AbsTol', 1e-10)。

可选符号积分方案 #

如果确实需要用符号积分做验证，可按如下方式实现，注意该方案对于复杂f(s)可能无法得到闭式解，运算速度远低于数值积分：

syms t s x sym_f(s)
p = sym(3);
delta = sym(0.1);
A_p = 1/int(exp(-s^2),s,-p,p);
rho_delta(t) = A_p*exp(-t^2/delta^2)/delta;

% 定义符号形式的f(s)，此处以f(s)=s为例
sym_f(s) = s;
% 计算符号积分
sym_result = int(rho_delta(x - s)*sym_f(s), s, 0, x);
% 代入x的取值转成数值结果
num_result = double(subs(sym_result, x, 1));
disp(num_result);

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