cholupdate

Оцените 1 обновление факторизации Холесского

Синтаксис

R1 = cholupdate(R,x)
R1 = cholupdate(R,x,'+')
R1 = cholupdate(R,x,'-')
[R1,p] = cholupdate(R,x,'-')

Описание

R1 = cholupdate(R,x) где R = chol(A) исходная факторизация Холесского A, возвращает верхний треугольный Фактор Холесского A + x*x', где x вектор-столбец соответствующей длины. cholupdate использование только диагональный и верхний треугольник R. Более низкий треугольник R проигнорирован.

R1 = cholupdate(R,x,'+') совпадает с R1 = cholupdate(R,x).

R1 = cholupdate(R,x,'-') возвращает Фактор Холесского A - x*x'. Сообщение об ошибке сообщает, когда R не является допустимый Фактор Холесского или когда downdated матрица не положительна определенный и так не имеет факторизации Холесского.

[R1,p] = cholupdate(R,x,'-') не возвратит сообщение об ошибке. Если p 0R1 Фактор Холесского A - x*x'. Если p больше 0R1 Фактор Холесского исходного A. Если p 1, cholupdate отказавший, потому что downdated матрица не положительна определенный. Если p 2, cholupdate отказавший, потому что верхний треугольник R не был допустимый Фактор Холесского.

Примеры

A = pascal(4)
A =

     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    20

R = chol(A)
R =

     1     1     1     1
     0     1     2     3
     0     0     1     3
     0     0     0     1
x = [0 0 0 1]';

Это называется рангом одним обновлением A начиная с rank(x*x') 1:

A + x*x' 
ans =
     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    21

Вместо того, чтобы вычислить Фактор Холесского с R1 = chol(A + x*x'), мы можем использовать cholupdate:

R1 = cholupdate(R,x)
R1 =
    1.0000    1.0000    1.0000    1.0000
         0    1.0000    2.0000    3.0000
         0         0    1.0000    3.0000
         0         0         0    1.4142

Затем уничтожьте положительную определенность (и на самом деле сделайте матрицу сингулярной) путем вычитания 1 от последнего элемента A. downdated матрица:

A - x*x'
ans =
 
     1     1     1     1
     1     2     3     4
     1     3     6    10
     1     4    10    19

Сравнение chol с cholupdate:

R1 = chol(A-x*x')
Error using chol
Matrix must be positive definite.
R1 = cholupdate(R,x,'-')
Error using cholupdate
Downdated matrix must be positive definite.

Однако вычитание 0.5 от последнего элемента A производит положительную определенную матрицу, и мы можем использовать cholupdate вычислить его Фактор Холесского:

x = [0 0 0 1/sqrt(2)]';
R1 = cholupdate(R,x,'-') 
R1 =
    1.0000    1.0000    1.0000    1.0000
         0    1.0000    2.0000    3.0000
         0         0    1.0000    3.0000
         0         0         0    0.7071

Советы

cholupdate работает только на полные матрицы.

Расширенные возможности

Смотрите также

|

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте