MATLAB中热方程三重求和实现异常结果的排查问询

问题描述 #

我得到热方程的通解为三重求和形式，在MATLAB中实现后结果不符合预期：当x=y=z=5时，代码输出值为403.4368，但我预期结果最大值不超过80。请排查是代码问题还是求和公式本身的问题。

我的MATLAB代码如下：

% Define the bounds for the summation
n = 50; % Upper limit for the first index
m = 50; % Upper limit for the second index
p = 50; % Upper limit for the third index
to = 350;
tc = 80;
Lx = 10;
Ly = 10;
Lz = 10;
k = 1;

x=5;
y=5;
z=5;

% Initialize the sum variable
sum_val = 0;

% Perform the triple summation
for a = 1:n
    alpha = (2*a-1)*pi/Lx;
    for b = 1:m
        beta = (2*b-1)*pi/Ly;
        for c = 1:p
            gamma = (2*c-1)*pi/Lz;
            lambda = (alpha^2 + beta^2 + gamma^2);
            % Compute the term to be added to the sum
            term = (sin(alpha*x) * sin(beta*y) * sin(gamma*z) * exp(-lambda*k*0)) / ...
                ((2*a-1) + (2*b-1) + (2*c-1));

            
            % Add the term to the sum
            sum_val = sum_val + term;
        end
    end
end

u = to + sum_val*(64*(to-tc)/pi^3);

% Display the result
disp(['The triple summation value is: ' num2str(u)]);

问题排查与修正 #

核心错误出在求和项的分母：

• 代码中错误地将分母写成了三个奇数的和 (2*a-1) + (2*b-1) + (2*c-1)，但热方程分离变量得到的通解中，分母应为三个奇数的乘积 (2*a-1)*(2*b-1)*(2*c-1)。
• 这个错误会导致每一项的数值被大幅放大（比如a=b=c=1时，分母从1变为3，项的数值是正确值的3倍；索引增大时，偏差会更显著），最终求和结果远超合理范围。

另外可简化的细节：

• exp(-lambda*k*0) 等于1，可直接从代码中移除，不影响计算结果。

修正后的关键代码片段：

term = (sin(alpha*x) * sin(beta*y) * sin(gamma*z)) / ...
    ((2*a-1)*(2*b-1)*(2*c-1));

修正后运行代码，当x=y=z=5时，结果会符合预期（具体数值取决于求和项的截断次数n/m/p）。

补充说明 #

若你的预期是最大值不超过80，需确认初始与边界条件的设定逻辑：当前代码中to=350是初始温度，tc=80是边界温度，t=0时中心(x=y=z=5)温度应等于初始温度350；当t趋近于无穷大时，温度会趋近于边界温度80。如果是计算稳态（t→∞）结果，需将exp(-lambda*k*0)中的0替换为足够大的时间值，此时指数项趋近于0，sum_val趋近于0，u会趋近于80，匹配你的预期。

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