追赶法的matlab实现,只需要输入A和d
在方程Ax=d中,只需要输入A和d即可,如果不是三对角阵就会报错。
下面是代码
function [x] = chase(A,d)
%求解线性方程组 Ax=d, 其中 A是三对角阵,如果输入的不是三对角矩阵就会报错
%A是系数矩阵
%d是方程组的右端向量
if ~isequal(tril(A,-1)-diag(diag(tril(A,-1),-1),-1),zeros(size(A...
原创
2019-07-24 11:42:30 ·
8636 阅读 ·
1 评论