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

@ProBat9 · Регистрация: 20.06.2022

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

Здравствуйте пишу вам с проблемой. Мне дано было задание решение системы дифференциальных уравнений с помощью метода Рунге-Кутты 4 порядка с помощью функции rkfixed. Но постоянно выдаёт ошибку у D: значение должно быть действительным и должно быть функцией. Прикрепляю файл. До этого задания никогда не использовал маткад. Прошу понимания

@nick55782012 · 21.06.2022, 07:33

ProBat9, мало начальных условий,еще корень квадратный извлекаете из отрицательного числа,да и в самом уравнении,возможно,есть ошибки.

@mathmichel · 21.06.2022, 08:57

Что-то с системой дифференциальных уравнений не так - четыре уравнения и только три интегрируемые функции! Приведите эту систему в оригинале! Да и под квадратными корнями не могут быть отрицательные выражения!

@ProBat9 · 21.06.2022, 14:00 **[ТС]**

Уважаемые @nick55782012, @mathmichel мне надо было решить вот систему состоящую из этих уравнений
$https://www.cyberforum.ru/cgi-bin/latex.cgi?(d[I])/dt=-k_d [I] (2.36) <br /> (d[R^*])/dt=2fk_d [I] (2.37) <br /> (d[S])/dt=-K(r_s [S]^2+[S][A]) (2.38) <br /> (d[A])/dt=-K(r_A [A]^2+[S][A])$

где K = $https://www.cyberforum.ru/cgi-bin/latex.cgi?\inline [(d[R^* ]/dt)/(k_tss/(k_psa^2 ) [S]^2+2 k_tss/(k_psa k_pas ) [S][A]+k_taa/(k_pas^2 ) [A]^2 )]^(1/2)$
В файле что я кидал x1,x2,x3 это [I][S][A]
Преподаватель мне выдал придал задание ввиде скриншота, прикрепляю его. Я его воспроизвёл. Но что то пошло не так. По вашему совету я вынес знак минуса за скобки, и выражение стало действительным. Но жаловалось на малое количество исходных данных. В примере учителя были значения x1,x3,x4. Но не было x2. И я вот сделал вектор из 4 значений где х2 нулевой но видимо что то не то

@mathmichel · 21.06.2022, 14:52

У Вас уже в начальной точке t=0 производные (которые в матрице D) с огромными порядками чисел при обычных начальных значениях x₁, х₂, х₃ и х₄, поэтому система сразу идёт в разнос. Проверьте выражения, которые стоят в матрице D!

@nick55782012 · 21.06.2022, 15:12

ProBat9, kd слишком велико,если уменьшить и добавить нач.условие,то вот.

Новые блоги и статьи Все статьи Все блоги /
Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .	SDL3 для Web (WebAssembly): Сборка SDL3 из исходников с помощью CMake и Emscripten 8Observer8 27.02.2026 Недавно вышла версия 3. 4. 2 библиотеки SDL3. На странице официальной релиза доступны исходники, готовые DLL (для x86, x64, arm64), а также библиотеки для разработки под Android, MinGW и Visual Studio. . . .	SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами 8Observer8 20.02.2026 Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .	Конвертировать закладки radiotray-ng в m3u-плейлист damix 19.02.2026 Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI Eddy_Em 18.02.2026 Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»

@nick55782012 2744 / 1920 / 954 Регистрация: 25.12.2016 Сообщений: 5,463
	21.06.2022, 07:33
	ProBat9, мало начальных условий,еще корень квадратный извлекаете из отрицательного числа,да и в самом уравнении,возможно,есть ошибки. 0

@mathmichel $Эксперт по математике/физике$ 11073 / 7374 / 3990 Регистрация: 14.01.2014 Сообщений: 16,810
	21.06.2022, 08:57
	Что-то с системой дифференциальных уравнений не так - четыре уравнения и только три интегрируемые функции! Приведите эту систему в оригинале! Да и под квадратными корнями не могут быть отрицательные выражения! 0

@mathmichel $Эксперт по математике/физике$ 11073 / 7374 / 3990 Регистрация: 14.01.2014 Сообщений: 16,810
	21.06.2022, 14:52
	У Вас уже в начальной точке t=0 производные (которые в матрице D) с огромными порядками чисел при обычных начальных значениях x₁, х₂, х₃ и х₄, поэтому система сразу идёт в разнос. Проверьте выражения, которые стоят в матрице D! 1

@nick55782012 2744 / 1920 / 954 Регистрация: 25.12.2016 Сообщений: 5,463
	21.06.2022, 15:12
	ProBat9, kd слишком велико,если уменьшить и добавить нач.условие,то вот. Миниатюры 1