Функция не видит переменные

@sqrt322 · Регистрация: 03.10.2018

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

Функция не видит переменные, какую бы не вставил вместо c0.Как это можно пофиксить?Фотография 1
И еще такой вопрос, правильно ли я сделал систему?То есть мне нужно одной функцией посчитать x2 и подставить во вторую.Если нет,то подскажите как надо, пожалуйста.Фотография 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
import numpy as np
import scipy.integrate as integrate
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
 
#define universal variables
c0 = 0.03
c1 = 10
c3 = 3
m0=1
m1=2
m=m0/m1
#just a little extra, quite unimportant
def f(x):
    f = ((-m)*x)+(1/2)*((m0+m1)/m1)*(abs(x+1.0)-abs(x-1.0))
    return f
#the actual function calculating
def x(G):
    return np.array([(-c1/c3)*f(G[1]-G[0]),
                     (-1/c3)*(f(G[1]-G[0])+G[2]),
                     c3*G[1]])
def dH_dt(H, t=0):
    return np.array([(-c1/c3)*f(H[1]-H[0]),
                     (-1/c3)*(f(H[1]-H[0])+H[2])+(с0/с3)*(G[1]-H[1]),
                     c3*H[1]])
#computational time steps
t = np.linspace(0, 500, 4000)
#x, y and z initial conditions
G0 = [0.7, 0.7, 0.7]
H0 = [0.7, 0.7, 0.7]
H, infodict = integrate.odeint(dH_dt, H0, t, full_output=True)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot(H[:,0], H[:,1], H[:,2])
plt.show()

@Welemir1 · 24.04.2020, 13:48

Сообщение от sqrt322

def y(x):
x1=(c1/с3)*f(x2-x1),
x2=(-1/с3)*(f(x2-x1)+x3),
x3=c3*x2
print('x равно', x2)
return x2

тоже вызывает вопросы:
1) зачем есть аргумент х, ты его нигде внутри функции не используешь
2) зачем мы рассчитываем х3 - нигде потом не используем
3) х1 и х2 в своем рассчете используют друг друга, это вообще по идее не должно работать никак

попробуй просто функцию у(игрек) вызвать с каим то аргументом -разве она работает?

@sqrt322 · 24.04.2020, 13:52 **[ТС]**

Welemir1,Вот если кракто,то система имеет вид(картинка 1),то есть мы считаем в первой системе х2 с точкой и подставляем во вторую систему.Зачем мы считаем х3 мне тоже не особо ясно,но так видимо надо.Я немного переделал код, по моей ужасной логике-так правильнее.Но возникает ошибка, которую я писал выше.

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
import numpy as np
import scipy.integrate as integrate
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
 
#define universal variables
с0 = 0.03
c1 = 10
c3 = 3
m0=1
m1=2
m=m0/m1
x=1
x2=2
x3=3
x1=1
 
#just a little extra, quite unimportant
def f(x):
    f = ((-m)*x)+(1/2)*((m0+m1)/m1)*(abs(x+1.0)-abs(x-1.0))
    return f
#the actual function calculating
def y(x1,x2,x3):
    x11=(c1/c3)*f(x2-x1),
    x22=(-1/c3)*(f(x2-x1)+x3),
    x33=c3*x2
    print('x равно', x22)
    return x22
 
 
def dH_dt(H, t=0):
    return np.array([(-c1/c3)*f(H[1]-H[0]),
                     (-1/c3)*(f(H[1]-H[0])+H[2])+(с0/c3)*(y(x1,x2,x3)-H[1]),
                     c3*H[1]])
#computational time steps
t = np.linspace(0, 500, 4000)
#x, y and z initial conditions
G0 = [0.7, 0.7, 0.7]
H0 = [0.7, 0.7, 0.7]
H, infodict = integrate.odeint(dH_dt, H0, t, full_output=True)
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.plot(H[:,0], H[:,1], H[:,2])
plt.show()

@Garry Galler · 24.04.2020, 18:33

Сообщение от sqrt322

Но возникает ошибка, которую я писал выше

Все элементы массива numpy должны быть одного типа. Иначе они нафиг не нужны, ибо становятся типа object.
Что у вас в массиве мне неизвестно и как его кастануть к одному типу не скажу.
Попробуйте его вывести принтом и посмотреть что у вас там. Если там инты вперемешку с флоатами и чем-то еще - приводите все к одному типу.
И еще массив numpy не должен быть рваным, то иметь разную размерность по числу элементов в строках.

То есть вот та кнормально:

Python
1
2
>>> np.array([(0,1.3),(0, 1.2)]).dtype
dtype('float64')

Вот так плохо:

Python
1
2
3
>>> np.array([(0,1.3,0),(0, 1.2)]).dtype
dtype('O')
>>>

P.S. Перестаньте давать переменным и функциям математич. имена.
Программирование это не мелом по доске черкать.
Если вы хотите, чтобы ваш код прочитали - давайте человекочитаемые имена переменным.
И делайте, пожалуйста, пробелы между переменными и операторами.

Добавлено через 27 минут
В общем вы написали такой страшный говнокод, что я понял где ошибка только переписав функцию dH_dt с разбивкой на вычислительные части.
Оказалось что y_result, то есть результат возврата функции y это кортеж, а не число.
Кортеж, потому что в функции y вас стоит случайная запятая после x22=(-1/c3)*(f(x2-x1)+x3),
Таким образом x22 превращается в кортеж из одного числа.
Кортеж, естественно нельзя умножать или складывать с числами, из-за него менялся тип numpy массива на object и т.д.

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
def dH_dt(H, t=0):
    y_result = y(x1, x2, x3)
    f_result = f(H[1] - H[0])
    f_result2 = f_result + H[2]
    f_result3 = y_result - H[1]
    
    arr = [
        (-c1/c3) * f_result,
        (-1/c3) * f_result2 + (с0/c3) * f_result3,
        c3 * H[1]
    ]
    
    return np.array(arr)

@sqrt322 · 24.04.2020, 21:37 **[ТС]**

Garry Galler, У меня возникла новая проблема с это же программой, не могли бы помочь еще раз, пожалуйста.В это раз сложнее чем с случайной запятой...
Если кратко,то должна быть взаимосвязанная картинка.У меня получилось построить для каждой системы отдельно(как вы прислали, и два графика на одной картинке),но как можно построить взаимосвязанные?Насколько я понял, в конечном счете должно получится что-то такое.То есть по логике, нужно чтобы был массив из функции y(x1,x2,x3),но при этом нужно из этого массива передавать х22(как было до этого) в функцию dH_dt.

@Viktorrus · 24.04.2020, 21:58

Welemir1, После того, как ТС вместо того , что бы поблагодарить, что с ним неучем возятся, он еще в грубой форме заявляет что это Вы плохо понимаете, я бы после этого перестал помогать.

@Welemir1 · 24.04.2020, 22:25

Viktorrus, всегда есть надежда...

Новые блоги и статьи Все статьи Все блоги /
Валидация и контроль данных табличной части документа перед записью Maks 22.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в КА2. Задача: контроль и валидация данных табличной части документа перед записью с учетом регламента компании. . .	Отчёт о затраченных материалах за определенный период с макетом печатной формы Maks 21.04.2026 Отчёт из решения ниже размещён в конфигурации КА2. Задача: разработка отчёта по затраченным материалам за определённый период, с возможностью вывода печатной формы отчёта с шапкой и подвалом. В. . .	Отчёт о спецтехнике находящейся в ремонте Maks 20.04.2026 Отчёт из решения ниже размещен в конфигурации КА2. Задача: отобразить спецтехнику, которая на данный момент находится в ремонте. Есть нетиповой документ "Заявка на ремонт спецтехники" который. . .	Памятка для бота и "визитка" для читателей "Semantic Universe Layer (Слой семантической вселенной)" Hrethgir 19.04.2026 Сгенерировано для краткого описания по случаю сборки и компиляции скелета серверного приложения. И пусть после этого скажут, что статьи сгенерированные AI - туфта и не интересно. И это не реклама -. . .
Запрет удаления строк ТЧ документа при определённом условии Maks 19.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа "Аккумуляторы", разработанного в конфигурации КА2. У данного документа есть ТЧ, в которой в зависимости от прав доступа. . .	Модель заражения группы наркоманов alhaos 17.04.2026 Условия задачи сформулированы тут Суть: - Группа наркоманов из 10 человек. - Только один инфицирован ВИЧ. - Колются одной иглой. - Колются раз в день. - Колются последовательно через. . .	Мысли в слух. Про "навсегда". kumehtar 16.04.2026 Подумалось тут, что наверное очень глупо использовать во всяких своих установках понятие "навсегда". Это очень сильное понятие, и я только начинаю понимать край его смысла, не смотря на то что давно. . .	My Business CRM MaGz GoLd 16.04.2026 Всем привет, недавно возникла потребность создать CRM, для личных нужд. Собственно программа предоставляет из себя базу данных клиентов, в которой можно фиксировать звонки, стадии сделки, а также. . .

@sqrt322 0 / 0 / 0 Регистрация: 03.10.2018 Сообщений: 57
	24.04.2020, 21:37 [ТС]
	Garry Galler, У меня возникла новая проблема с это же программой, не могли бы помочь еще раз, пожалуйста.В это раз сложнее чем с случайной запятой... Если кратко,то должна быть взаимосвязанная картинка.У меня получилось построить для каждой системы отдельно(как вы прислали, и два графика на одной картинке),но как можно построить взаимосвязанные?Насколько я понял, в конечном счете должно получится что-то такое.То есть по логике, нужно чтобы был массив из функции y(x1,x2,x3),но при этом нужно из этого массива передавать х22(как было до этого) в функцию dH_dt. Миниатюры 0

@Viktorrus 1732 / 970 / 199 Регистрация: 22.02.2018 Сообщений: 2,693 Записей в блоге: 6
	24.04.2020, 21:58
	Welemir1, После того, как ТС вместо того , что бы поблагодарить, что с ним неучем возятся, он еще в грубой форме заявляет что это Вы плохо понимаете, я бы после этого перестал помогать. 0

@Welemir1 Автоматизируй это! 7392 / 4819 / 1246 Регистрация: 30.03.2015 Сообщений: 13,695 Записей в блоге: 29
	24.04.2020, 22:25
	Viktorrus, всегда есть надежда... 0

Опции темы

Функция не видит переменные

Решение