С Новым годом! Форум программистов, компьютерный форум, киберфорум
Python: Научные вычисления
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.83/47: Рейтинг темы: голосов - 47, средняя оценка - 4.83
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 19

Построить прямые в Matplotlib и конформно отобразить их

16.04.2019, 14:21. Показов 9820. Ответов 9

Студворк — интернет-сервис помощи студентам
Здравствуйте!!!
Нужна ваша помощь! Дано уравнение y=a*x, где x=[-1, 1], a=[0,1]. Нужно построить несколько прямых в этих пределах и конформно отобразить их на этом же рисунке с помощью формулы - complex(x, y)**2.

Вот мои попытки построить прямые. Вопрос. Можно ли как нибудь задать 'a' так же как 'x', что бы не прописывать его от руки?
Или можно просто упростить всю эту конструкцию?
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
%matplotlib inline
import matplotlib.pyplot as plt
import math, cmath
 
# построим прямые функции y = a*x, x=[-1,1], a=[0,1]
 
 
x = [-1, -0.5, -0.2, 0.4, 0.7, 1]
y = [xx*0.1 for xx in x]
y2 = [xx*0.3 for xx in x]
y3 = [xx*0.5 for xx in x]
y4 = [xx*0.7 for xx in x]
y5 = [xx*0.9 for xx in x]
y6 = [xx*1 for xx in x];
 
plt.plot(x,y)
plt.plot(x,y2)
plt.plot(x,y3)
plt.plot(x,y4)
plt.plot(x,y5)
plt.plot(x,y6);
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
16.04.2019, 14:21
Ответы с готовыми решениями:

Отобразить конформно область
Отобразить конформно на верхнюю полуплоскость область {{z\in C: \alpha < arg z <\pi - \alpha }} \{ {z\in C: z = iy, y \geq 1}}, где 0...

Конформно отобразить область на верхнюю полуплоскость
Отобразить область {z: |z|<1, z не лежит на отрезке } на верхнюю полуплоскость :help:

Как конформно отобразить верхний полушар единичного радиуса в нижний (с помощью зеркального отображения)
Здравствуйте. Помогите, пожалуйста, конформно отобразить верхний полушар единичного радиуса в нижний (с помощью зеркального отображения,...

9
 Аватар для Wolflind
127 / 107 / 31
Регистрация: 17.03.2016
Сообщений: 445
Записей в блоге: 1
16.04.2019, 14:38
Лучший ответ Сообщение было отмечено artemka_ss как решение

Решение

Цитата Сообщение от artemka_ss Посмотреть сообщение
Можно ли как нибудь задать 'a' так же как 'x'
да

Python
1
2
3
y = []
for a in xrange(1,10):
 y.append([xx*a/10 for xx in x])
Добавлено через 3 минуты
либо просто
Python
1
plt.plot(x,[xx*a/10 for xx in x])
в место
Python
1
y.append([xx*a/10 for xx in x])
1
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 19
16.04.2019, 15:03  [ТС]
Wolflind, можешь расписать прямо на моем примере? А то я не совсем понимаю. Я подставил то, что ты предложил сразу после того как задал 'Х' и мне выдает ошибку: NameError: name 'xrange' is not defined
0
1303 / 843 / 409
Регистрация: 12.03.2018
Сообщений: 2,305
16.04.2019, 15:13
Лучший ответ Сообщение было отмечено artemka_ss как решение

Решение

попробуйте написать range
1
16.04.2019, 15:24

Не по теме:

в Python2 есть range и xrange они отличаются принципом работы
xrange Python2 тоже самое что range Python3

0
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 19
16.04.2019, 15:49  [ТС]
ioprst, Спасибо! Работает!

Добавлено через 13 минут
А как мне теперь преобразовать эти прямые по формуле - complex(x,y)**2 ?
На паре приводили вот такой пример, но это просто какая то жуть, хоть и красиво:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
%matplotlib inline
import matplotlib.pyplot as plt
import math, cmath
 
def conform(f, n = 10, m = 100,
           figsize=(6,6), dpi=300):
    '''
    Функция, рисующая комформное отображение с помощью f(z)
    n- число разбиений единичного квадрата
    m - число точек на линии
    figsize - размер рисунка в дюймах 
    dpi - пространственное разрешения в пикселях на дюйм
    '''
    fig = plt.figure(figsize = figsize, dpi = dpi)
    plt.title('конформное отображение')
    
    # рисуб линии по горизонтали
    for i in range (n+1):
        yy = i/n
        # формирую массивы
        x = [j/m for j in range(m+1)]
        y = [yy for j in range(m+1)]
        # массив комплексных чисел
        z = [ complex(x[j], y[j]) for j in range(m+1)]
        # отображение
        w = [f(z[j]) for j in range(m+1)]
        # делаю массивы с действительной и мнимой частью комплексного числа
        u = [w[j].real for j in range(m+1)]
        v = [w[j].imag for j in range(m+1)]
        plt.plot(u,v, 'r-')
        
# вертикальные линии
    for i in range (n+1):
        xx = i/n
        w = [f(complex(xx, j/m)) for j in range(m+1)]
        u = [w[j].real for j in range(m+1)]
        v = [w[j].imag for j in range(m+1)]
        #print(u,v)
        plt.plot(u,v, 'b-')
    
conform(f= lambda z:z**2, m=10) # вызываю функцию
0
16.04.2019, 19:26

Не по теме:

Цитата Сообщение от artemka_ss Посмотреть сообщение
хоть и красиво:
скорей просто
Цитата Сообщение от artemka_ss Посмотреть сообщение
жуть
как минимум оптимизации 0
зачем столько for как минимум

Python
1
2
3
4
5
6
7
8
9
       x = [j/m for j in range(m+1)]
        y = [yy for j in range(m+1)]
        # массив комплексных чисел
        z = [ complex(x[j], y[j]) for j in range(m+1)]
        # отображение
        w = [f(z[j]) for j in range(m+1)]
        # делаю массивы с действительной и мнимой частью комплексного числа
        u = [w[j].real for j in range(m+1)]
        v = [w[j].imag for j in range(m+1)]
можно 1 циклом сделать, а не 6 будит. Итого в 6 раз быстрей
и вообще желательно разбить на 2 функции а не все в одной
также если Python2 то в место range нужно использовать xrange

0
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 19
17.04.2019, 12:05  [ТС]
Wolflind, а как мне объединить то что ты прислал с этим:

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
%matplotlib inline
import matplotlib.pyplot as plt
import math, cmath
 
# построим прямые функции y = a*x, x=[-1,1], a=[0,1]
 
# размер холста
fig = plt.figure(figsize = (6, 6), dpi = 96)
 
# оси
plt.plot([-1,1], [0,0], color='black', lw=2) # ось x
plt.plot([0,0], [-1,1], color='black', lw=2) # ось y
 
x = [-1, -0.5, -0.2, 0.4, 0.7, 1]
y = []
for a in range(1, 10):
    plt.plot(x,[xx*a/10 for xx in x], 'r-')
0
 Аватар для Wolflind
127 / 107 / 31
Регистрация: 17.03.2016
Сообщений: 445
Записей в блоге: 1
17.04.2019, 12:15
Цитата Сообщение от artemka_ss Посмотреть сообщение
а как мне объединить
что значит объединить?
0
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 19
17.04.2019, 12:33  [ТС]
Цитата Сообщение от Wolflind Посмотреть сообщение
Сообщение от artemka_ss
а как мне объединить
что значит объединить?
Мне нужно теперь эти прямые конформно отобразить на этом же рисунке с помощью формулы complex(x, y)**2.

Вот тут например мы рисовали единичный квадрат и потом его там же отражали:
Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
def conform(f, n):
    '''
    Конформное отображение с помощью функции f(z)
    Сторона квадрата делится на n частей
    '''
    # вверх - вправо - вниз - влево
    zeros = [0]*(n+1)
    ones = [1]*(n+1)
    seq1 = [i/n for i in range(n+1)]
    seq2 = [1-i/n for i in range(n+1)]
    
    xx = zeros + seq1 + ones + seq2
    yy = seq1 + ones + seq2 + zeros
    
    plt.plot(xx, yy, 'b-')
    
    
    # конформное отображение
    m = len(xx)
    # комплексные числа из xx+jyy
    zz = [complex(xx[i], yy[i]) for i in range(m)]
    # отображение
    ww = [f(z) for z in zz]
    uu = [w.real for w in ww]
    vv = [w.imag for w in ww]
    plt.plot(uu, vv, 'r-')
    
 
plt.figure(figsize=(4,4))    
conform(lambda z:z**2, 10)
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
17.04.2019, 12:33
Помогаю со студенческими работами здесь

Отобразить график используя библиотеку matplotlib
Всем привет! Нужно нарисовать графики, используя библиотеку matplotlib, графики изображены на картинках. Заранее благодарю за помощь

Matplotlib и PyQt5 Как отобразить график ?
Привет всем! Прошу помощи! Есть график который выводится по средствам matplotlib'a в Jupytere. fig = plt.figure(figsize=(20, 10)) ...

Нужно корректно отобразить функцию Гаусса [matplotlib]
import numpy as np import scipy.stats as stats import pylab as pl h = sorted() fit = stats.norm.pdf(h, np.mean(h),...

Python matplotlib отобразить график функции в заданной точке собственной формы
Язык Python.Есть график функции. Не важно какой. Например такой 3D график. При компиляции построенный график открывается в другом окне,...

Построить график в matplotlib
Построить график единичного прямоугольного импульса шириной 300 мс, с центром в момент времени 4 с. С помощью matplotlib


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
10
Ответ Создать тему
Новые блоги и статьи
Восстановить юзерскрипты Greasemonkey из бэкапа браузера
damix 15.01.2026
Если восстановить из бэкапа профиль Firefox после переустановки винды, то список юзерскриптов в Greasemonkey будет пустым. Но восстановить их можно так. Для этого понадобится консольная утилита. . .
Изучаю kubernetes
lagorue 13.01.2026
А пригодятся-ли мне знания kubernetes в России?
Сукцессия микоризы: основная теория в виде двух уравнений.
anaschu 11.01.2026
https:/ / rutube. ru/ video/ 7a537f578d808e67a3c6fd818a44a5c4/
WordPad для Windows 11
Jel 10.01.2026
WordPad для Windows 11 — это приложение, которое восстанавливает классический текстовый редактор WordPad в операционной системе Windows 11. После того как Microsoft исключила WordPad из. . .
Classic Notepad for Windows 11
Jel 10.01.2026
Old Classic Notepad for Windows 11 Приложение для Windows 11, позволяющее пользователям вернуть классическую версию текстового редактора «Блокнот» из Windows 10. Программа предоставляет более. . .
Почему дизайн решает?
Neotwalker 09.01.2026
В современном мире, где конкуренция за внимание потребителя достигла пика, дизайн становится мощным инструментом для успеха бренда. Это не просто красивый внешний вид продукта или сайта — это. . .
Модель микоризы: классовый агентный подход 3
anaschu 06.01.2026
aa0a7f55b50dd51c5ec569d2d10c54f6/ O1rJuneU_ls https:/ / vkvideo. ru/ video-115721503_456239114
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR
ФедосеевПавел 06.01.2026
Owen Logic: О недопустимости использования связки «аналоговый ПИД» + RegKZR ВВЕДЕНИЕ Введу сокращения: аналоговый ПИД — ПИД регулятор с управляющим выходом в виде числа в диапазоне от 0% до. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru