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

Функция solve

22.05.2019, 13:42. Показов 9428. Ответов 3
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Здравствуйте.
У меня есть два уравнения из которых мне нужно вывести значения А, ω, x_dot, alpha_dot.
Вот код:
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
%typeset_mode True
import numpy as np
import scipy
import matplotlib.pyplot as plt
import numpy.fft as fft
import matplotlib
import matplotlib.pyplot as plt1
import matplotlib.pyplot as plt2
from scipy import integrate, optimize, interpolate
from scipy.special import gamma
# y = [0..3]
# x=y[0]
# x_dot=y[1]
# alpha=y[2]
# alpha_dot=y[3]
%var angle1, x, x_dot,x1dot, x2dot, epsilon, A, Aa, alpha, alpha_dot, m1, m2, k1,g, l1, Ek, Ep, Lagran1, Lagran2
Ek = m2*((x_dot^2)+(2*x_dot*l1*cos(alpha)*alpha_dot)+(l1^2)*(alpha_dot^2)*((cos(alpha)**2)-(sin(alpha)**2)))/2
Ep = m2*g*(l1-l1*cos(alpha))+(k1*(x^2)/2)
L = Ek - Ep
Ldiff = L.diff(x_dot )
Ldiff2 = L.diff(alpha_dot)
t=var('t')
xdot_tmp = var ('X_Dot')
alphadot_tmp = var ('Alpha_Dot')
alpha_tmp = function ('Alpha')(t)# For substitution alpha
x_tmp = function ('X')(t) #For substitution x
angle1 = var ('theta')
Ldiff3 = Ldiff.subs([x_dot==x_tmp.diff(t), x==x_tmp ,alpha==alpha_tmp,alpha_dot==alpha_tmp.diff(t)]).diff(t)
Ldiff4 = Ldiff2.subs([x_dot==x_tmp.diff(t), x==x_tmp ,alpha==alpha_tmp,alpha_dot==alpha_tmp.diff(t)]).diff(t)
Lagran1 = Ldiff3-L.diff(x)-A*sin(angle1*t)
Lagran2= Ldiff4-L.diff(alpha)
print ('Lagran1')
Lagran1
print ('Lagran2')
Lagran2
omega =var ('omega') #For subs second derivative alpha. It is pendulum acceleration
Aa =var('A') #For subs second derivative x. It is trolley acceleration
eq1=Lagran1.subs([x_tmp.diff(t)==x_dot,x_tmp==x,alpha_tmp.diff(t)==alpha_dot,alpha_tmp==alpha, alpha_tmp.diff(t).diff(t)==omega, x_tmp.diff(t).diff(t)==Aa])
print ('eq1')
eq1
eq2=Lagran2.subs([x_tmp.diff(t)==x_dot,x_tmp==x,alpha_tmp.diff(t)==alpha_dot,alpha_tmp==alpha, alpha_tmp.diff(t).diff(t)==omega, x_tmp.diff(t).diff(t)==Aa])
print ('eq2')
eq2
solve([ eq1,eq2  ],[Aa,omega ])
#x_dot is trolley velocity
#alpha_dot is pendulum velocity
solve([ eq1,eq2  ],[x_dot,alpha_dot ])
Полазив в интернете я решил, что проще всего будет сделать это при помощи функции "solve". Значения А, ω вывелись без проблем, а вот при выводе x_dot, alpha_dot появляются только пустые скобки.
Что я сделал неверно и как исправить код так, чтобы значения выводились?

P.S. На картинке во вложении скриншот работы программы.
Миниатюры
Функция solve  
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
22.05.2019, 13:42
Ответы с готовыми решениями:

solve() takes 0 positional arguments but 3 were given
Есть код: from sympy import * from tkinter import * from tkinter.ttk import Frame, Button, Style import numpy as np from...

Функция solve
Добрый день, господа. У меня возникла проблема с функцией solve при решении системы уравнений. Подскажите в чем проблема. { 2*x1 +...

Функция solve
никак не могу найти функцию solve в 14 маткаде, приходится редактировать файл с примером, можете продиктовать всю последовательность...

3
Эксперт Python
 Аватар для dondublon
4652 / 2072 / 366
Регистрация: 17.03.2012
Сообщений: 10,182
Записей в блоге: 6
22.05.2019, 16:52
Цитата Сообщение от Miclosh Посмотреть сообщение
import matplotlib.pyplot as plt1 import matplotlib.pyplot as plt2
Не надо так.

Вопрос интересный, конечно, ещё бы знать, что у вас за solve.
0
0 / 0 / 0
Регистрация: 21.11.2015
Сообщений: 21
22.05.2019, 23:07  [ТС]
dondublon, а что с import matplotlib.pyplot?

Solve - это функция из Sympy, вот пример: https://docs.sympy.org/latest/... lvers.html.

Добавлено через 2 часа 10 минут
Я разобрался в чем проблема, она не в коде. Вопрос снят.
0
Эксперт Python
 Аватар для dondublon
4652 / 2072 / 366
Регистрация: 17.03.2012
Сообщений: 10,182
Записей в блоге: 6
23.05.2019, 10:13
Miclosh, вы импортируете класс дважды. Так не делается. Импортируйте один раз, создавайте дважды.
matplotlib, конечно, тоже хороши, показывают неправильный шаблон в каждом своём example.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
23.05.2019, 10:13
Помогаю со студенческими работами здесь

Функция solve()
Почему результат расчета функции solve имеет такой вид,будто не до конца решено? >> syms x, X1=solve(0.2*x^3-3*x^2+5*x+25) X1 = ...

Функция solve
нужно найти нуль функции при помощи solve. выдается ошибка такого рода(КАК ИСПРАВИТЬ?): ??? Error using ==> solve>getEqns at 202 ...

Функция solve. Проблема
-3х^5+х^4-2х^2+х+1 solve,x-> Только степени выглядят обычно как в математике без этих значков.Я в Маткаде с помощью функции solve хочу...

Не работает функция solve
Доброй ночи! Набираю готовую программу, но не получается сделать так, как решено. В данном фрагменте программы не работает функция solve....

Функция solve в MathCad не работает
Помогите разобраться с проблемой! В MathCad не находит корни уравнения функция solve. Ниже приложил скриншот и сам файл с расчетом


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

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Новые блоги и статьи
Изучаю 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% до. . .
Модель микоризы: классовый агентный подход 2
anaschu 06.01.2026
репозиторий https:/ / github. com/ shumilovas/ fungi ветка по-частям. коммит Create переделка под биомассу. txt вход sc, но sm считается внутри мицелия. кстати, обьем тоже должен там считаться. . . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru