レス数が1スレッドの最大レス数(10000件)を超えています。残念ながら投稿することができません。
管理人の独り言２
-
function X = CreateTable(n, c)
for i = 1:n
if i <> 1 then
X(i, i - 1) = -1
end
X(i, i) = c
if i <> n then
X(i, i + 1) = -1
end
end
endfunction
function z = CGAlgorithm(A, b, x_0, n)
r = b - A * x_0
p = r
x = x_0
for k = 1:n
alpha = (r' * p) / (p' * (A * p))
x_new = x + alpha * p
r_new = r - alpha * A * p
bet = -(r_new' * (A * p)) / (p' * (A * p))
p_new = r_new + bet * p
x = x_new
r = r_new
p = p_new
z(k) = norm(A * x - b)
end
endfunction
function z = JacobiAlgorithm(A, b, x_0, n)
x = x_0
for k = 1:n
for i = 1:length(b)
su = 0
for j = 1:length(b)
if i <> j then
su = su + A(i, j) * x(j)
end
end
x_new(i) = (-su + b(i)) / A(i, i)
end
x = x_new
z(k) = norm(A * x - b)
end
endfunction
|
|
掲示板管理者へ連絡
無料レンタル掲示板
