Есть массивы с точками:
X[0.0, 0.0, 0.0, 0.0],Y[0.0, 0.3333, 0.6667, 1.0],
kX[0.0, 0.0, 0.0],kY[0.1667, 0.5, 0.8333]
Нужно сформировать двумерные массивы R2, U и V размера(3,4), элементы которых считаются по формулам:
r2 = math.sqrt(math.pow((kX[k]-X[i]), 2) + math.pow((kY[k]-Y[i] ) , 2))
u = 1 / (2 * math.pi) * ((kY[k]-Y[i]) / math.pow(r2[k][i], 2))
v = - 1 / (2 * math.pi) * ((kX[k]-X[i]) / math.pow(r2[k][i], 2))
Поскольку я новичек, меня особенно интересует подробное описание следующего кода:
###псевдокод####
r2 = math.sqrt(math.pow((kX[k]-X[i]), 2) + math.pow((kY[k]-Y[i] ) , 2))
### реализация на питоне ####
r2 = [[(cp.y - p.y)**2 + (cp.x - p.x)**2 for p in points]
for cp in control_points]
где используется collections.namedtuple
и вообще как обращаться к элементам массивов без использования циклов, желательно на этом примере.
Случай с точками двух массивов понятен:
u = [[1 / (2 * math.pi)*(kx - x) for x in X] for kx in kX]
v = [[1 / (2 * math.pi)*(ky - y) for y in Y] for ky in kY]
А как быть, если нужно:
u = [[1 / (2 * math.pi)*(kx - x) for x in X] for kx in kX] / r2[k][i]
v = [[1 / (2 * math.pi)*(ky - y) for y in Y] for ky in kY] / r2[k][i]
scipy.spatial.distance
модуль, который различные функции по вычислению пространственных метрик предоставляет, чтобы руками их заново не реализовывать.¶ r2k это ряд—не номера в ряду, это сам ряд в матрице r2. Попробуйте,for r2k in r2: print(r2k)
. Аналогично, r2ki это не номер, это сам элемент. В моём ответе индексы вообще не используются. Если вы все элементы по порядку просматриваете, то нет необходимости по индексам обращаться. По практикуйтесь с...