Программирование МКЭ в МATLAB. Даутов P.З. - 16 стр.

UptoLike

Составители: 

Глава 2
ПОСТРОЕНИЕ СЕТОК В MATLAB
§ 1. Cоздание и хранение разреженных матриц
В MatLab использование разреженных и плотных (полных) мат-
риц фактически не различается, за исключением того, что
разреженная матрица должна быть описана (создана) функцией
sparse;
имена некоторых функций для разреженных матриц могут быть
другими.
Имеется три способа формирования разреженной матрицы.
1) Создание разреженной матрицы из полной. Например,
B =[1 3 0
0 4 0
2 0 5 ] ; % В i s f u l l matrix
A = s p a r s e (B) % А i s s p a r s e
A = ( 1 , 1 ) 1
( 3 , 1 ) 2
( 1 , 2 ) 3
( 2 , 2 ) 4
( 3 , 3 ) 5
Для изображения разреженной матрицы А используется его так на-
зываемое координатное представление, когда указываются индексы
и значения ненулевых элементов, которые перечисляются последо-
вательно по столбцам. Таким образом А представляется на экране
тремя массивами одинаковой длины (i, j) a.
2) Элементное наполнение матрицы. Предыдущую матрицу можно
ввести также следующим образом.
A=s p ar se ( 3 , 3 ) ; % a l l e l e m e nts o f A e q u a l 0 .
A( 1 ,1)=1 ; A( 3 ,1)=2; A( 1 ,2 )=3; A(2 ,2 ) =4; A(3 , 3 ) =5;
3) Основной способ. В этом случае разреженная матрица получается
из ее координатного представления. Например, следующие команды
также создают матрицу А из предыдущих примеров.
                                                      Глава 2
                 ПОСТРОЕНИЕ СЕТОК В MATLAB


             § 1. Cоздание и хранение разреженных матриц

   В MatLab использование разреженных и плотных (полных) мат-
риц фактически не различается, за исключением того, что
    • разреженная матрица должна быть описана (создана) функцией
      sparse;
    • имена некоторых функций для разреженных матриц могут быть
      другими.
    Имеется три способа формирования разреженной матрицы.
1) Создание разреженной матрицы из полной. Например,
B =[1 3 0
    0 4 0
    2 0 5];         % В i s f u l l matrix
A = s p a r s e (B) % А i s s p a r s e

A = (1 ,1)                   1
    (3 ,1)                   2
    (1 ,2)                   3
    (2 ,2)                   4
    (3 ,3)                   5

Для изображения разреженной матрицы А используется его так на-
зываемое координатное представление, когда указываются индексы
и значения ненулевых элементов, которые перечисляются последо-
вательно по столбцам. Таким образом А представляется на экране
тремя массивами одинаковой длины (i, j) a.
2) Элементное наполнение матрицы. Предыдущую матрицу можно
ввести также следующим образом.
A=s p a r s e ( 3 , 3 ) ; % a l l e l e m e n t s o f A e q u a l 0 .
A( 1 , 1 ) = 1 ; A( 3 , 1 ) = 2 ; A( 1 , 2 ) = 3 ; A( 2 , 2 ) = 4 ; A( 3 , 3 ) = 5 ;

3) Основной способ. В этом случае разреженная матрица получается
из ее координатного представления. Например, следующие команды
также создают матрицу А из предыдущих примеров.