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

@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 01.04.2026 Алгоритм из решения ниже реализован на нетиповом документе "ВыдачаОборудованияНаСпецтехнику" разработанного в конфигурации КА2, в дополнении к предыдущему решению. На форме документа создается. . .	К слову об оптимизации kumehtar 01.04.2026 Вспоминаю начало 2000-х, университет, когда я писал на Delphi. Тогда среди программистов на форумах активно обсуждали аккуратную работу с памятью: нужно было следить за переменными, вовремя. . .	Идея фильтра интернета (сервер = слой+фильтр). Hrethgir 31.03.2026 Суть идеи заключается в том, чтобы запустить свой сервер, о чём я если честно мечтал давно и давно приобрёл книгу как это сделать. Но не было причин его запускать. Очумелые учёные напечатали на. . .	Модель здравосоХранения 6. ESG-повестка и устойчивое развитие; углублённый анализ кадрового бренда anaschu 31.03.2026 В прикрепленном документе раздумья о том, как можно поменять модель в будущем
10 пpимет, которые всегда сбываются Maks 31.03.2026 1. Чтобы, наконец, пришла маршрутка, надо закурить. Если сигарета последняя, маршрутка придет еще до второй затяжки даже вопреки расписанию. 2. Нaдоели зима и снег? Не надо переезжать. Достаточно. . .	Перемещение выделенных строк ТЧ из одного документа в другой Maks 31.03.2026 Реализация из решения ниже выполнена на примере нетипового документа "ВыдачаОборудованияНаСпецтехнику" с единственной табличной частью "ОборудованиеИКомплектующие" разработанного в конфигурации КА2. . . .	Functional First Web Framework Suave DevAlt 30.03.2026 Sauve. IO Апнулись до NET10. Из зависимостей один пакет, работает одинаково хорошо как в режиме проекта так и в интерактивном режиме. из сложностей - чисто функциональный подход. Решил. . .	Автоматическое создание документа при проведении другого документа Maks 29.03.2026 Реализация из решения ниже выполнена на нетиповых документах, разработанных в конфигурации КА2. Есть нетиповой документ "ЗаявкаНаРемонтСпецтехники" и нетиповой документ "ПланированиеСпецтехники". В. . .

Опции темы