上传者: 38732315
|
上传时间: 2021-12-13 10:08:17
|
文件大小: 2KB
|
文件类型: -
使用中心差分求解拉普拉斯方程。 脚本 TriDiagMake 创建三对角矩阵 B 和 -I 的 A 三对角矩阵。 然后程序创建具有指定边界条件的向量 b 并求解系统 Au=b。 后处理在最后一部分完成,其中解的向量 u 被制成矩阵 U',然后与 BC 组合以创建矩阵 T。 可以轻松修改以包括 Neumann 和 Robin BC; 此外,Au=b 的求解器是 Matlab 标准的; 它可以被更强大的东西取代(比如并行雅可比方法或快速高斯-赛德尔)。 用于计算方形板上的温度分布或其他类似的潜在问题。 这个想法是创建一个易于理解、不言自明的代码。 希望你喜欢!