Две функции на одном графике

@artemka_ss · Регистрация: 01.06.2015

Студворк — интернет-сервис помощи студентам

Здравствуйте!
Есть вот такой вопрос. Необходимо на одном графике построит две функции с разным значением эпсилон. Можно ли на одном графике построить значения вызываемых функций u1 и u2? Или можно как то сделать это внутри "b"?

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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import solve_banded as sb
 
def b(dx=0.1, xmax=1, tmax=1, dt=1e-6, D=1, C=1, ntt=5, nt=1000, title='', beta=1, e=1, Uc=0.5, P=10):
   
    nx = int(xmax / dx + 1)                           # число разбиений по х
    x = np.linspace(0, xmax, nx)
    hx = xmax / nx
                                                    
    t = np.linspace(0, tmax, nt)                      # число разбиений по t
    ht = tmax / nt
                                                            
    a = -D/(hx**2)                                    # Приведенные коэффиценты
    g = a
    
    u = np.zeros((nt,nx),dtype=np.float)              # Создаем массивы для расчета решения системы
    u[0,:] = Uc                                       # Задаем граничные и начальные и условия
    u[:,0] = Uc
    u[:,1] = Uc
    
    ab = np.zeros((3,nx),dtype=np.float)              
    
    r = np.zeros(nx,dtype=np.float)                 
    
    countWhile = 0
    for i in range(1, nt):                          # Цикл по времени,прохождение по слоям
        epsilon = True
        lastResult = u[i-1,:]
        
        while epsilon:
            countWhile = countWhile + 1
            ab[0,:] = a
            ab[1,:] = C / ht + 2 * D / (hx**2) + beta * ((lastResult-Uc) ** 4) - P
            ab[2,:] = g 
            
            r[:] = (C / ht) * lastResult
            
            currentResult = sb((1, 1), ab, r)            # solve_banded преобразует массив ab в a перестановкой элементов: a[i,j] == ab[u + i - j, j]
                                                        # и решает СЛАУ a * x = r  
            dif = lastResult - currentResult
            
            if max(abs(dif)) < e:
                epsilon = False
                
            lastResult = currentResult
 
        u[i,:] = lastResult  
       
    print("Всего слоев: {}. Количество итераций: {}".format(nt, countWhile))
    print("Максимальная температура на последнем слое: {}".format(max(abs(u[nt-1,:]))))
    
    tstep = int(nt / ntt) 
    if tstep == 0:
        tstep = 1
        
    for i in range(1,nt,tstep):
        plt.plot(x, u[i,:], label = 't={0:8.3e}'.format(t[i]))    #отображение результатов
   
    plt.legend(loc = 'best')
    plt.xlabel('X')
    plt.ylabel('U(t, x)')
    plt.show()
    
    return u
 
nt1000 = 1000
print("Исследование 1: сколько итераций необходимо для решения уравнения.")
u1 = b(Uc=0.5, P=10, dx=0.01, dt=0.0001, D=0.1, C=5, nt=nt1000, e=0.001, title= 'dvd')
u2 = b(Uc=0.5, P=10, dx=0.01, dt=0.0001, D=0.1, C=5, nt=nt1000, e=0.000001)
print("Разница температур на последнем слое: {}".format(max(abs(u2[nt1000 - 1,:] - u1[nt1000 - 1,:]))))

Добавлено через 1 час 11 минут
Вопрос решен! Спасибо за внимание)
Вдруг кому понадобится:

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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import solve_banded as sb
 
def b(dx=0.1, xmax=1, tmax=1, dt=1e-6, D=1, C=1, ntt=5, nt=1000, title='', beta=1, e=1, Uc=0.5, P=10, drawGraph = True):
   
    nx = int(xmax / dx + 1)                           # число разбиений по х
    x = np.linspace(0, xmax, nx)
    hx = xmax / nx
                                                    
    t = np.linspace(0, tmax, nt)                      # число разбиений по t
    ht = tmax / nt
                                                            
    a = -D/(hx**2)                                    # Приведенные коэффиценты
    g = a
    
    u = np.zeros((nt,nx),dtype=np.float)              # Создаем массивы для расчета решения системы
    u[0,:] = Uc                                       # Задаем граничные и начальные и условия
    u[:,0] = Uc
    u[:,1] = Uc
    
    ab = np.zeros((3,nx),dtype=np.float)              
    
    r = np.zeros(nx,dtype=np.float)                 
    
    countWhile = 0
    for i in range(1, nt):                          # Цикл по времени,прохождение по слоям
        epsilon = True
        lastResult = u[i-1,:]
        
        while epsilon:
            countWhile = countWhile + 1
            ab[0,:] = a
            ab[1,:] = C / ht + 2 * D / (hx**2) + beta * ((lastResult-Uc) ** 4) - P
            ab[2,:] = g 
            
            r[:] = (C / ht) * lastResult
            
            currentResult = sb((1, 1), ab, r)            # solve_banded преобразует массив ab в a перестановкой элементов: a[i,j] == ab[u + i - j, j]
                                                        # и решает СЛАУ a * x = r  
            dif = lastResult - currentResult
            
            if max(abs(dif)) < e:
                epsilon = False
                
            lastResult = currentResult
 
        u[i,:] = lastResult  
       
    print("Всего слоев: {}. Количество итераций: {}".format(nt, countWhile))
    print("Максимальная температура на последнем слое: {}".format(max(abs(u[nt-1,:]))))
    
    tstep = int(nt / ntt) 
    if tstep == 0:
        tstep = 1
        
    for i in range(1,nt,tstep):
        plt.plot(x, u[i,:], label = 't={0:8.3e}'.format(t[i]))    #отображение результатов
   
    if drawGraph == True:
        plt.legend(loc = 'best')
        plt.xlabel('X')
        plt.ylabel('U(t, x)')
        plt.show()
    
    return u
 
nt1000 = 1000
print("Исследование 1: сколько итераций необходимо для решения уравнения.")
u1 = b(Uc=0.5, P=10, dx=0.01, dt=0.0001, D=0.1, C=5, nt=nt1000, e=0.01, title= 'dvd', drawGraph=False)
u2 = b(Uc=0.5, P=10, dx=0.01, dt=0.0001, D=0.1, C=5, nt=nt1000, e=0.0000001, drawGraph=False)
print("Разница температур на последнем слое: {}".format(max(abs(u2[nt1000 - 1,:] - u1[nt1000 - 1,:]))))
plt.legend(loc = 'best')
plt.xlabel('X')
plt.ylabel('U(t, x)')
plt.show()

Новые блоги и статьи Все статьи Все блоги /
Перемещение выделенных строк ТЧ из одного документа в другой Maks 30.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .	Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .	Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .	Настройка движения справочника по регистру сведений Maks 29.03.2026 Решение ниже реализовано на примере нетипового справочника "ТарифыМобильнойСвязи" разработанного в конфигурации КА2, с целью учета корпоративной мобильной связи в коммерческом предприятии. . . .
Автозаполнение реквизита при выборе элемента справочника Maks 27.03.2026 Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. При выборе "Спецтехники" (Тип Справочник. Спецтехника), заполняется. . .	Сумматор с применением элементов трёх состояний. Hrethgir 26.03.2026 Тут. https:/ / fips. ru/ EGD/ ab3c85c8-836d-4866-871b-c2f0c5d77fbc Первый документ красиво выглядит, но без схемы. Это конечно не даёт никаких плюсов автору, но тем не менее. . . всё может быть. . .	Автозаполнение реквизитов при создании документа Maks 26.03.2026 Программный код из решения ниже размещается в модуле объекта документа, в процедуре "ПриСозданииНаСервере". Алгоритм проверки заполнения реализован для исключения перезаписи значения реквизита,. . .	Команды формы и диалоговое окно Maks 26.03.2026 1. Команда формы "ЗаполнитьЗапчасти". Программный код из решения ниже на примере нетипового документа "ЗаявкаНаРемонтСпецтехники" разработанного в конфигурации КА2. В качестве источника данных. . .