Here它说如果没有连接端点,矩阵应该更简单。我认为通过仅改变矩阵我可以实现线性迭代。那么如何将以下原始代码更改为线性版本以打开活动轮廓?
以下是原始代码
% Populating the penta diagonal matrix
A = zeros(m,m);
b = [(2*alpha + 6 *beta) -(alpha + 4*beta) beta];
brow = zeros(1,m);
brow(1,1:3) = brow(1,1:3) + b;
brow(1,m-1:m) = brow(1,m-1:m) + [beta -(alpha + 4*beta)]; % populating a template row
for i=1:m
A(i,:) = brow;
brow = circshift(brow',1)'; % Template row being rotated to egenrate different rows in pentadiagonal matrix
end
[L U] = lu(A + gamma .* eye(m,m));
Ainv = inv(U) * inv(L); % Computing Ainv using LU factorization
end
这段代码循环迭代,但我想以线性方式迭代。这意味着未连接起点和终点。 如果有人能提供任何见解,我将不胜感激。