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

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

16.04.2019, 14:21. Показов 9919. Ответов 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
Ответ Создать тему
Новые блоги и статьи
Установка Qt Creator для C и C++: ставим среду, CMake и MinGW без фреймворка Qt
8Observer8 05.04.2026
Среду разработки Qt Creator можно установить без фреймворка Qt. Есть отдельный репозиторий для этой среды: https:/ / github. com/ qt-creator/ qt-creator, где можно скачать установщик, на вкладке Releases:. . .
AkelPad-скрипты, структуры, и немного лирики..
testuser2 05.04.2026
Такая программа, как AkelPad существует уже давно, и также давно существуют скрипты под нее. Тем не менее, прога живет, периодически что-то не спеша дополняется, улучшается. Что меня в первую очередь. . .
Отображение реквизитов в документе по условию и контроль их заполнения
Maks 04.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеСпецтехники", разработанного в конфигурации КА2. Данный документ берёт данные из другого нетипового документа. . .
Фото всей Земли с борта корабля Orion миссии Artemis II
kumehtar 04.04.2026
Это первое подобное фото сделанное человеком за 50 лет. Снимок называют новым вариантом легендарной фотографии «The Blue Marble» 1972 года, сделанной с борта корабля «Аполлон-17». Новое фото. . .
Вывод диалогового окна перед закрытием, если документ не проведён
Maks 04.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: реализовать программный контроль на предмет проведения документа. . .
Программный контроль заполнения реквизитов табличной части документа
Maks 02.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: 1. Реализовать контроль заполнения реквизита. . .
wmic не является внутренней или внешней командой
Maks 02.04.2026
Решение: DISM / Online / Add-Capability / CapabilityName:WMIC~~~~ Отсюда: https:/ / winitpro. ru/ index. php/ 2025/ 02/ 14/ komanda-wmic-ne-naydena/
Программная установка даты и запрет ее изменения
Maks 02.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "СписаниеМатериалов", разработанного в конфигурации КА2. Задача: при создании документов установить период списания автоматически. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru