Дифференциация тканей головы Исследуемая цилиндрическая область заменяется на цилиндрическую сетку (она может быть неравномерной по переменной r) и задача. - презентация

1

2

3 Дифференциация тканей головы

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

5 Упорядочение точек цилиндрической сетки происходит, как показано на рисунке. После заполнения самого нижнего слоя нумерация переносится на второй слой и т. д. вплоть до исчерпания всех слоев. Упорядочение точек отдельного слоя цилиндра

6 Такие виды матриц получаются после естественного упорядочения неизвестных в дискретизированном уравнении Пуассона с переменными коэффициентами в цилиндрических координатах Матрица коэффициентов проводимости

7 Проблемы решения линейных систем Решение линейных систем с разреженными матрицами порядка больше прямыми методами (как-то метод Гаусса) неэффективно. В таких случаях применяются итерационные методы с различного рода переобуславливателями. Основное внимание сконцентрировано на методе бисопряженных градиентов с построением переобуславливателя Фурье на каждой итерации, что хотя и увеличивает вычислительную сложность одной итерации, но существенно сокращает время поиска решения системы.

8

9 Сравнение методов с использованием переобуславливателей BiCG с переобуславливателем Фурье