
- •Тема 9. РАБОТА С РАЗРЕЖЕННЫМИ МАТРИЦАМИ
- •Где применяются разреженные матрицы
- •Поля физических величин
- •Математическая модель поля
- •Прямая полевая задача
- •Градиент скалярной функции
- •Лапласиан от скалярного поля
- •Обратная полевая задача
- •уравнения математической физики
- •метод сеток для решения ДУ
- •Получение алгебраической СЛАУ для искомых значений ui
- •aiui 1 biui ciui 1
- •При решении более сложных ДУ методом конечных элементов приходится решать СЛАУ
- •Прямой ход метода Гаусса приводит эту систему к СЛАУ с верхней треугольной матрицей
- •Схема алгоритма
- •Одним из рациональных способов хранения не нулевых элементов матрицы является
- •Эффективность работы с таким массивом определяется введением класса с продуманным набором методов и
- •Тип ячеек массива
- •Методы работы со стеком, которые в дальнейшем понадобятся для реализации метода Гаусса
- •Методы работы со стеком Добавить элемент в начало стека
- •Добавить в упорядоченный по j стек
- •Методы работы со стеком
- •Методы работы со стеком Добавить в упорядоченный по j стек (конец)
- •Методы работы со стеком Читать j – й элемент
- •Методы работы со стеком
- •Методы работы со стеком удалить 1 – й элемент
- •Базовый класс работы с матрицей
- •Методы работы с массивом стеков Создать массив
- •Методы работы с массивом стеков Уничтожить массив
- •Методы работы с массивом стеков Добавить элемент в массив
- •Методы работы с массивом стеков Добавить элемент с упорядочиванием по j
- •Методы работы с массивом стеков Читать элемент без удаления
- •Методы работы с массивом стеков Распечатать матрицу в StringGrid
- •Методы работы со СЛАУ Создать объект СЛАУ
- •Методы работы со СЛАУ Уничтожить объект СЛАУ
- •Методы работы со СЛАУ Добавить элемент в b
- •Методы работы со СЛАУ Распечатать объект СЛАУ
- •Решить СЛАУ с верхней треугольной матрицей
- •Решить СЛАУ с верхней треугольной матрицей
- •Пример работы с классом TSlau
- •Пример работы с классом TSlau
- •Контроль утечки памяти
- •Контрольные вопросы
Тема 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 |