有限差分法边界值问题
4视图(30天)
显示旧的评论
嗨,我迷失在这个边值问题需要使用有限差分法:
我附上我所做的和提供的提示我们的教授,我也尝试过不同的方法从其他matlab和youtube的问题,但它总是一些错误运行。
下面的代码我已经到目前为止,它的工作原理,但它看起来不正确,谁能帮助好吗?我将非常感谢!
% %初始化参数
clc,清晰
T_B = 600;%边界条件
T_E = 350;%边界条件
T_Inf = 300;
M = 4;
L = 2.0;
%离散化
num_nodes = 51;
num_segments = num_nodes-1;
h = L / num_segments;
% x = T_B: - h: T_E;
xj = 0: h: L;
Q = @ (T) (4 * (300 - T));
%建立矩阵
对角线= [1 / h;2 / h ^ 2 * 1 (num_segments-1 1);1 / h);
一个=诊断接头(对角线);
为我= 2:num_segments + 1
(张,我)= 1 / h ^ 2;(我,张)= 1 / h ^ 2;% % % % % % % % % % % % % %的迹象
结束
(1、2)= 0;(num_segments + 1, num_segments) = 0;
b = [T_B / h;问(xj (2: end-1));T_E / h];% % % % % % % % % % % % % % % Q (xj (2 _end-1))不只是问
T = \ b;
情节(xj, T)网格