Практикум j практическое примщенше численных методов

виде двумерного массива. Для примера рассмотрим представление решения (рис. 3.19), которое строиться па равномерной цилиндрической сетке.
1


Рис. 3.19 Полярная сетка

.8 1.6
1.4
1.2
1.0

0. 
   8 0.6 0.4
   

import matplotlib.pyplot as pit
import numpy as np
ro = np.linspace(l., 2., 21)
phi = np.linspace(np.pi/6, np.pi/3, 21)
Ro, Phi = np.meshgrid(ro, phi)
X = Ro*np.cos(Phi)
Y = Ro*np.sin(Phi)
Z = np.exp(-Ro)
pit.contourf(X, Y, Z, 20)
pit.axis(’equal’)
pit. showO
Для данных па произвольном множестве узлов задача визуализации в пакете Matplotlib решается следующим образом. Сначала по заданному множеству узлов строиться триангуляция Делоне и данные интерполируются на струк­турированную сетку, после чего и рисуется график. Задача интерполяции решается с использованием функции griddataO из модуля mlab. Другая воз­
можность связана с использованием расширения Natgrid (mpl.toolkits.natgrid инсталлируется отдельно от пакета Matplotlib).
Пусть х,y,z есть одномерные массивы координат х, у и заданных значений интерполируемой функции Прямоугольная сетка задастся одномерными или двумерными массивами X и Y как обычно. Двумерный массив значений па этой повой сетке определяется функцией griddata(х, 2/, У). Пример