Привет, я новичок в Pyhton, и мне нужна помощь с этим решением. Я пытаюсь закодировать уравнение преобразования Йоуковского для построения профиля. Мне дается круг с радиусом 1,5 центра (-0,5 (ч), 0,0 (к)), и лямбда равна 1. При построении круга это дает желаемый результат, однако нанесение кода аэродинамического профиля дает мне ошибку, и результат графика близок к профилю, но все же странно. Пожалуйста, кто-нибудь, помогите мне с решением.
Изображения кодов моего ноутбука Jupyter:
![my jupyter notebook code image](https://i.stack.imgur.com/0vqd2.jpg)
Пожалуйста, проверьте эти изображения для лучшего понимания.
h = -0.5
k = 0.0
r = 1.5 #radius
lam = 1
center = h + 1j*k
x = r*numpy.cos(theta) + h
y = r*numpy.sin(theta) + k
theta = numpy.linspace(0.0,2*math.pi,100)
alpha = numpy.linspace(0.0,2*math.pi,100)
def circle(alpha,r):
Z = r*numpy.exp(1j*(alpha))
return Z
Z = circle(alpha,r)
pyplot.plot(x,y)
pyplot.grid(True)
pyplot.gca().set_aspect('equal')
def eigen(Z,center,alpha,lam):
z = (Z-center)*(-1j*(alpha)) + (lam**2)*(1j*(alpha))/(Z-center)
return z
e = eigen(Z,center,alpha,lam)
a = ((Z-center) +1/(Z-center))*numpy.cos(alpha)
b = (1/(Z-center) -(Z-center))*numpy.sin(alpha)
pyplot.plot(a,b)
pyplot.gca().set_aspect('equal')
pyplot.grid(True)