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

Решение системы дифференциальных уравнений

20.02.2021, 17:29. Показов 3509. Ответов 1
Метки нет (Все метки)

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

Построить фазовый портрет системы дифференциальных уравнений

https://www.cyberforum.ru/cgi-bin/latex.cgi?x'' +y'' - 3*x' - y' -2*x +2*y = 0

https://www.cyberforum.ru/cgi-bin/latex.cgi?2*x'' +3*y'' - 9*x' - 2*y' -4*x +6*y = 0

В конце концов получил общее решение:
https://www.cyberforum.ru/cgi-bin/latex.cgi?x(t) = \frac{1}{10}*{e}^{-2t} * ({C}_{1}*(3* {e}^{4*t} + 4*{e}^{2*t} * cos(t) + 3) + {C}_{2}*(3* {e}^{4*t} - 2*{e}^{2*t} * sin(t) -3) - {C}_{3}*({e}^{4*t} - 4*{e}^{2*t} * sin(t) -1) + {C}_{4}*({e}^{4*t} - 2*{e}^{2*t} * cos(t) + 1))

https://www.cyberforum.ru/cgi-bin/latex.cgi?y(t) = \frac{1}{10}*{e}^{-2t} * (3*{C}_{1}*({e}^{4*t} - 4*{e}^{2*t} * sin(t) - 1) + 3*{C}_{2}*({e}^{4*t} - 2*{e}^{2*t} * cos(t) + 1) - {C}_{3}*({e}^{4*t} - 12*{e}^{2*t} * cos(t) + 1) + {C}_{4}*({e}^{4*t} + 6*{e}^{2*t} * sin(t) - 1))

Проверил при помощи онлайн калькулятора, вроде как, сошлось. *скрин прилагаю*

Однако не могу понять, как считать далее.

Наверное можно взять случайные начальные условия, откуда найти Cn и затем подставить их и идти циклом, изменяя t.

В диссонанс вводит такой момент, что все заданий, что выше моего по методичке, можно спокойно решить при помощи метода Эйлера для задачи Коши.
Миниатюры
Решение системы дифференциальных уравнений  
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
20.02.2021, 17:29
Ответы с готовыми решениями:

Численное решение дифференциальных уравнений первого порядка
1.Метод Эйлера. 2.Метод Рунге-Кутта четвертого порядка уравнение - у'=sec(x)- y tg(x) # * - coding: utf-8 - *- # подключение...

Решение системы уравнений: LinAlgError: Singular matrix
Коллеги, выполняю такой код: import matplotlib.pyplot as plt %matplotlib inline import math import numpy as np def f(x):...

Методы решения дифференциальных уравнений
Здравствуйте. По практике необходимо найти решение дифференциального уравнения: \dot{\vartheta }={(1+e*\cos \vartheta...

1
 Аватар для regio1961
602 / 294 / 178
Регистрация: 06.06.2016
Сообщений: 552
24.02.2021, 22:31
Цитата Сообщение от Jagon Посмотреть сообщение
Наверное можно взять случайные начальные условия, откуда найти Cn и затем подставить их и идти циклом, изменяя t.
Да, если важны именно начальные условия. А можно использовать Cn как параметры и строить фазовый портрет при их различных значениях

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
from sympy import *    
from sympy.plotting import plot_parametric
import math
 
C1 = 1
C2 = 2
C3 = 3
C4 = 4
 
t = symbols( 't' )
# сюда надо подставить свои выражения, мне лень набирать
x = C1*sin(t) + C2*sin(2*t) + C3*sin(3*t) + C4*sin(4*t) 
y = C1*cos(t) - C2*cos(-2*t) + C3*cos(-3*t) + C4*cos(-4*t)
 
init_printing()
plot_parametric( x, y, ( t, -10, 10 ) )
При других значениях Cn кривые будут ещё красивее
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
24.02.2021, 22:31
Помогаю со студенческими работами здесь

Система нелинейных дифференциальных уравнений с гран. усл. solver_bvp
У меня есть система нелинейных ДУ. 1) Ошибка с f1f2_model(y, u1, u2, x1, x2, x22) в spi.solve_bvp- если я подаю функцию как lambda y,...

Численное решение системы дифференциальных уравнений методом Хэмминга
Доброе время суток, помогите, пожалуйста, написать программу для решения системы дифференциальных методом Хэмминга на python: dx/dt=y ...

Решение систем дифференциальных уравнений
Доброго времени суток!)) Разбирался с scipy и numpy, понял как решать системы дифференциальных уравнений , если у нас есть фиксированное...

Численное решение дифференциальных уравнений
Дано: y' = y * cos(x), y(0) = 1 Точное решение: y = {e}^{sin(x)} Найти решения с помощью методов Эйлера и Рунге - Кутты 2-го...

Численное решение дифференциальных уравнений первого порядка
Здравствуйте. Собственно, само задание: С помощью численных методов решить дифференциальное уравнение первого порядка у'=f(x,y). ...


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

Или воспользуйтесь поиском по форуму:
2
Ответ Создать тему
Новые блоги и статьи
Знаешь почему 90% людей редко бывают счастливыми?
kumehtar 14.04.2026
Потому что они ждут. Ждут выходных, ждут отпуска, ждут удачного момента. . . а удачный момент так и не приходит.
Фиксация колонок в отчете СКД
Maks 14.04.2026
Фиксация колонок в СКД отчета типа Таблица. Задача: зафиксировать три левых колонки в отчете. Процедура ПриКомпоновкеРезультата(ДокументРезультат, ДанныеРасшифровки, СтандартнаяОбработка) / / . . .
Настройки VS Code
Loafer 13.04.2026
{ "cmake. configureOnOpen": false, "diffEditor. ignoreTrimWhitespace": true, "editor. guides. bracketPairs": "active", "extensions. ignoreRecommendations": true, . . .
Оптимизация кода на разграничение прав доступа к элементам формы
Maks 13.04.2026
Алгоритм из решения ниже реализован на нетиповом документе, разработанного в конфигурации КА2. Задачи, как таковой, поставлено не было, проделанное ниже исключительно моя инициатива. Было так:. . .
Контроль заполнения и очистка дат в зависимости от значения перечислений
Maks 12.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: реализовать контроль корректности заполнения дат назначения. . .
Архитектура слоя интернета для сервера-слоя.
Hrethgir 11.04.2026
В продолжение https:/ / www. cyberforum. ru/ blogs/ 223907/ 10860. html Знаешь что я подумал? Раз мы все источники пишем в голове ветки, то ничего не мешает добавить в голову такой источник, который сам. . .
Подстановка значения реквизита справочника в табличную часть документа
Maks 10.04.2026
Алгоритм из решения ниже реализован на примере нетипового документа "ПланированиеПерсонала", разработанного в конфигурации КА2. Задача: при выборе сотрудника (справочник Сотрудники) в ТЧ документа. . .
Очистка реквизитов документа при копировании
Maks 09.04.2026
Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru