Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Презентации 2часть / Лекция_31_Разреженные матрицы.ppt
Скачиваний:
19
Добавлен:
11.05.2015
Размер:
482.3 Кб
Скачать

Тема 9. РАБОТА С РАЗРЕЖЕННЫМИ МАТРИЦАМИ

Где применяются разреженные матрицы

Базовый класс для работы с разреженными матрицами на основе массива из стеков

Реализация метода Гаусса решения СЛАУ с разреженной матрицей

07/02/19

1

Где применяются разреженные матрицы

Матрицы, содержащие большое число нулевых элементов называются разреженными

Решение многих прикладных задач (методы теории графов, численные методы решения дифференциальных уравнений методом конечных элементов, методы обработки изображений, методы криптографии) связано с обработкой матриц большой размерности n имеющих малое число ненулевых элементов

Использование специальных способов размещения в памяти разреженных матриц и разработка специальных алгоритмов работы с ними позволяет во многих случаях существенно снизить затраты компьютерных ресурсов и уменьшить время вычислений.

07/02/19

2

Поля физических величин

Любое физическое явление или процесс представляет собой распределение и изменение каких-либо физических величин (скалярных, векторных, тензорных) в некоторой области пространства и во времени.

Распределение некоторой величины в пространстве и во времени получило название поле.

Такие поля окружают нас -температурное поле, электромагнитное поле, поле скоростей, поле вероятности.

Основной задачей теоретической физики, а так же многих технических приложений является исследование полей физических величин, (полевые задачи).

Для их исследования разработано огромное программное обеспечение.

07/02/19

3

Математическая модель поля

Математической моделью поля является функция нескольких

переменных, обычно F (x, y, z,t)

T(x, y, z), E (x, y, z,t), H (x, y, z,t), V (x, y, z,t), (x, y, z,t), ij (x, y, z)

Поля бывают

скалярные

векторные

Тензорные (т.е.описывается несколькими векторами).

стационарные

Нестационарные

При исследовании полей выделяют две задачи – прямую и обратную

07/02/19

4

Прямая полевая задача

Задано поле F ( x, y, z,t) ; требуется установить характер этого поля, например быстроту его изменения от точки к точке.

Математическая теория поля занимается изучением дифференциальных и интегральных свойств различных полей.

Здесь для векторной

u ux x0 uy y0 uz z0

 

и скалярной

(x, y, z)

функций введены операторы дифференцирования:

r

,

r

r

r

,

2

divu

v

rotu

, v

 

07/02/19

5

Градиент скалярной функции

grad

r

 

r

 

r

x x0

y y0

z z0

•Характеризует направление наибольшего возрастания функции в каждой точке поля.

Ротор векторного поля

r

 

 

 

 

u

 

 

 

 

 

 

 

 

u

 

 

 

r

 

 

 

 

 

y

 

 

uz

 

 

y

 

ux

 

uz

 

 

r

 

ux

 

r

 

 

 

 

rot u

 

 

 

 

x

 

 

y

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

0

 

 

y

 

z

0

 

 

x

0

 

 

 

y

 

 

 

 

 

 

 

z

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дивергенция векторного поля

r

 

ux

 

uy

 

uz

 

divu

x

 

z

 

y

6

07/02/19

 

 

 

 

 

Лапласиан от скалярного поля

2 div ( )

2

 

2

 

2

x2

y2

z2

Характеризует распределение источников скалярного поля

Например имеется распределение температуры вдоль плоской пластины T(x,y).

Градиент температуры указывает направление максимального

распространения тепла в данной точке q T (x, y)

Тогда

 

div( T ) 2T qv ( x, y)

 

 

дает распределение мощности источников тепла

qv (x, y)

Зная распределение источников можно получить распределение

 

температуры из уравнения:

2 T

 

 

 

 

2 T

 

2 T

 

 

qv (x,T )

 

 

 

2

 

2

q (x, y)

x

2

 

 

 

 

 

 

v

 

 

 

 

 

x

 

 

y

 

 

 

 

 

7

 

07/02/19

 

 

 

 

 

 

 

 

Обратная полевая задача

Известны условия, в которых находится физический объект требуется найти распределение в пространстве некоторой физической величины, т.е. конкретного вида математического поля.

Чаще всего задача нахождения поля, удовлетворяющего требуемым условиям, приводит к решению краевой задачи для дифференциального (или интегрального) уравнения.

07/02/19

8

уравнения математической физики

Методы составления и, главное, решения уравнений

такого рода изучаются в разделе математической физики – теория дифференциальных уравнений в

частных производных и теория интегральных уравнений. Эти уравнения исторически получили название «уравнения математической физики».

Совокупность теории поля и теории дифференциальных уравнений в частных

производных образует так называемую

классическую математическую физику.

Основной метод решения - проекционно-сеточный метод, который получил название метод конечных

элементов

07/02/19

9

метод сеток для решения ДУ

Будем иллюстрировать реализацию метода сеток на решении простейшей одномерной краевой задачи для ДУ

 

 

 

 

 

 

 

 

 

u

f

 

x

 

;

u(0) ;

u(b) .

 

 

 

 

 

 

 

 

x x

 

 

 

 

 

 

 

u

Искомое

 

 

решение

 

 

 

 

 

 

 

0

b

x

 

 

 

07/02/19

 

u

 

Решение в виде

 

ui

таблицы

 

 

 

 

 

 

 

 

 

 

0

xi

b

x

 

 

 

 

 

10