Как построить профиль с помощью преобразования Йоуковски? - PullRequest
0 голосов
/ 05 ноября 2018

Привет, я новичок в Pyhton, и мне нужна помощь с этим решением. Я пытаюсь закодировать уравнение преобразования Йоуковского для построения профиля. Мне дается круг с радиусом 1,5 центра (-0,5 (ч), 0,0 (к)), и лямбда равна 1. При построении круга это дает желаемый результат, однако нанесение кода аэродинамического профиля дает мне ошибку, и результат графика близок к профилю, но все же странно. Пожалуйста, кто-нибудь, помогите мне с решением.

Изображения кодов моего ноутбука Jupyter:
my jupyter notebook code image
my jupyter notebook code image

Пожалуйста, проверьте эти изображения для лучшего понимания.

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)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...