Решение систем обыкновенных дифференциальных уравнений методом Рунге_Кутта 4-го порядка

@Turin464 · Регистрация: 18.01.2022

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

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

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
import math
import numpy as np
import matplotlib.pyplot as plt
import sympy as sp
import scipy
 
delta=0.12
n=3000
w1=2*math.pi*n
E=210000
mu=0.3
ro=8000
y1=np.array([[1],[0]]) #для однородной системы
r_int=np.linspace(80,640,1200)
p2=40000000
y0=np.array([[0],[0]]) #для неоднородной системы
r1=80
r2=640
 
def h(r):
    if r <= 160 and r >= 80:
        return 140
    if r > 160 and r < 600:
        return 80 - 13 * (r - 160) / 88
    if r >= 600 and r <= 640:
        return 30
def Fy0(r,Y,h): #для неоднородного ДУ
    Fy1=-((1+mu)/r)*Y[0]+((1-mu**2)/(r*E*h(r)))*Y[1]
    Fy2=(E*h(r)/r)*Y[0]+((mu-1)/r)*Y[1]-ro*h(r)*(w1**2)*r
    Fy0=np.array([[Fy1],[Fy2]])
    return Fy0
 
def Fy_n(r,Y,h): #для однородного ДУ
    Fy1=-((1+mu)/r)*Y[0]+((1-mu**2)/(r*E*h(r)))*Y[1]
    Fy2=(E*h(r)/r)*Y[0]+((mu-1)/r)*Y[1]
    Fy_ne0=np.array([[Fy1],[Fy2]])
    return Fy_ne0
 
y10=scipy.integrate.RK45(Fy_0,r1,y1,r2,0.5,1e-3,1e-6,True)

@Royal_X · 05.12.2023, 19:20

Turin464, не ответили в теме Решение систем обыкновенных дифференциальных уравнений методом Рунге_Кутта 4-го порядка, решили еще одну создать? А может вам стоит задуматься, почему никто не ответил в той теме? Для начала напишите все ваши дифференциальные уравнения в общепринятой математической нотации, используя редактор формул.

@Turin464 · 05.12.2023, 19:33 **[ТС]**

Конкретно по этому вопросу не понимаю, зачем вам нужны эти дифференциальные уравнения? Я не разобрался, как работает эта встроенная функция, почему она выдаёт ошибку и как эту ошибку решить.Вопрос, скорее,общего характера. Но, если с красиво вбитыми графическими формулами вам будет проще мне помочь, то без проблем.

@Royal_X · 05.12.2023, 20:05

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

Я не разобрался, как работает эта встроенная функция

тогда нужно читать документацию https://docs.scipy.org/doc/sci... .RK45.html

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

@Royal_X 6222 / 2922 / 1046 Регистрация: 01.06.2021 Сообщений: 10,819
	05.12.2023, 19:20
	Turin464, не ответили в теме Решение систем обыкновенных дифференциальных уравнений методом Рунге_Кутта 4-го порядка, решили еще одну создать? А может вам стоит задуматься, почему никто не ответил в той теме? Для начала напишите все ваши дифференциальные уравнения в общепринятой математической нотации, используя редактор формул. 0

@Turin464 0 / 0 / 0 Регистрация: 18.01.2022 Сообщений: 40
	05.12.2023, 19:33 [ТС]
	Конкретно по этому вопросу не понимаю, зачем вам нужны эти дифференциальные уравнения? Я не разобрался, как работает эта встроенная функция, почему она выдаёт ошибку и как эту ошибку решить.Вопрос, скорее,общего характера. Но, если с красиво вбитыми графическими формулами вам будет проще мне помочь, то без проблем. 0