Матричная лаборатория MatLab

         

PCG stopped after the maximum



Пример 1

» pcg(A.B)
Warning: PCG stopped after the maximum 4 iterations 
without converging to the desired tolerance le-006
The iterate returned (number 4) has relative residual 0.46 
> In C:\MATI_AB\toolbox\matlab\sparfun\pcg.m at line 347 
ans =


1.7006
1.2870
-2.0535
8.2912
В данном случае решение к успеху не привело, поскольку матрица А —несимметрическая. Новая функция rrrinres не требует, чтобы матрица А была положительно определенной. Достаточно, чтобы она была квадратной и симметрической. В отличие от peg минимизируется не относительная невязка, а абсолютная. Но и эта функция не может решить наш пример:
» minres(A.В.1е-6.1000000)
minres stopped at iteration 1000000 without converging 
to the desired tolerance le-006 
because the maximum number of iterations was reached. 
The iterate returned (number 1000000) has relative residual 0.011 
ans =
1.9669
3.7757
3.0789
1.9367
В MATLAB 6 появилась еще одна новая функция symmlq, которая использует LQ-алгоритм итерационного метода сопряженных градиентов и также не требует, чтобы ее входной аргумент — квадратная симметрическая матрица — была положительно определенной. Эта функция тоже не может решить наш пример, так как наша матрица А — квадратная, но не симметрическая.

Содержание раздела