Минимизация погрешности решения системы ОДУ

@Sawarez · Регистрация: 20.09.2021

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

Здравствуйте! Операция минимизации функции A (сумма погрешностей решения системы ОДУ) выдаёт ошибку и не находит решения. С начальными приближениями констант k и показателей степеней a,b,c и тд решалась система диффуров, результаты сравнивались с экспериментальными данными и складывались в погрешность (ф-ия А). Необходимо минимизировать погрешность, оптимизируя все эти параметры k и показатели степени. Может ли кто-то помочь? Может быть я просто допустил какую-то банальную синтаксическую ошибку или маткад вообще не умеет минимизировать таким образом? Менять начальные приближения я пытался, это не помогло. Если проблема сводится к этому, то каким способом лучше всего их приблизительно найти?

Минимизация погрешности решения системы ОДУ

@Krasme · 21.09.2021, 11:23

1. приложите книгу
2. (наскоро глянула) вы сначала находите z как массив, потом пытаетесь оптимизировать расхождение двух готовых таблиц...

@Sawarez · 21.09.2021, 17:50 **[ТС]**

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

1. приложите книгу
2. (наскоро глянула) вы сначала находите z как массив, потом пытаетесь оптимизировать расхождение двух готовых таблиц...

Книга - Вы имеете ввиду задание? Это не книжное задание, это реальные экспериментальные данные кинетики реакции образование комплекса CeNi6Ile12 3+. Был предложен механизм, в соответствии с ним была предложена система диффуров. V1, V2, V3 - концентрации определенных промежуточных продуктов. Константы k1 k2 и показатели a b c и тд - это то что нужно в конечном счете найти. То есть у нас есть реальные экспериментальные данные - таблица 1, в начале файлы, в столбцах концентрации разнесены по времени и есть массив который мы получаем, решая систему диффуров при задании констант какими-то числами. Погрешность рассчитывается разностью соответствующих чисел из двух массивов. Планировалось минимизировать эту самую погрешность для всех данных, тем самым получив правильные значения параметров k и других. Но собсна функция A оптимизироваться не хочет, а почему - не знаю...

@Krasme · 21.09.2021, 17:54

книга - Книга1.xlsx
я ведь не могу экспериментировать с вашей задачей, не имея внешние данные

@Sawarez · 21.09.2021, 17:56 **[ТС]**

А данные исчезают из маткадовского файла? Извиняюсь, не знал. Прикрепляю данные экселевским файлом)

@Krasme · 22.09.2021, 15:47

т.к. у вас 6-ой прайм, у меня 7-й, могу дать только pdf-результат

kinetika.pdf

считает долго, надо набраться терпения

@Sawarez · 22.09.2021, 16:15 **[ТС]**

Спасибо огромное! У меня на самом деле тоже 7й, просто тут в форуме в списке я видел только 6й максимум. Если это все переписать, он по идее будет работать? По графикам конечно видно что они сильно расходятся... Не знаете, можно ли как-то улучшить оптимизацию или это просто невозможно? Вероятно, что могут происходить некоторые иные процессы в системе, соответственно, систему диффуров надо усложнять или как-то еще менять + эксперимент. данные не самые лучшие) Но если сама "программа" будет работать, то не составит труда заменить данные или уравнения. В любом случае - спасибо огромное еще раз, очень выручаете!

@Krasme · 22.09.2021, 16:23

раз седьмой, держите файл
kinetika.zip

по части оптимизации, маткад выдал результат по заданным уравнениям. вероятно, в системе диффуров не те уравнения должны быть.
я ограничения по k выключила, маткад дал те же результаты.
с начальными приближениями попробовала поиграть, терпения не хватило дожидаться))

@Sawarez · 14.10.2021, 23:32 **[ТС]**

Здравствуйте! Периодически вожусь всё с той же задачей и так и не пришел к решению. Нашел в уравнении ошибки, исправил их, но только маткад не может минимизировать мою функцию A) Даже если задать один параметр, а мне надо 8. Не может - в смысле процесс не заканчивается в принципе, можно ждать сутки... и не дождаться. Может всё-таки что-то принципиально не так? Можно ли в маткаде шаги решения посмотреть?

@Krasme · 16.10.2021, 16:22

Не по теме:

для меня кинетика - темный лес..
я не могу оценить, насколько достоверны ваши формулы кинетики, помимо них еще надо знать, в каких диапазонах должны быть искомые параметры, не только k1,k2,k3.
есть вероятность, что встроенные функции оптимизации не смогут найти вам оптимальное решение.
надо искать другие способы, возможно, даже перестроить саму постановку задачи. это возможно сделать, только будучи сведущим в теме кинетики. здесь я вам помочь ничем не могу.

Новые блоги и статьи Все статьи Все блоги /
Как дизайн сайта влияет на конверсию: 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-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .	Реалии Hrethgir 01.03.2026 Нет, я не закончил до сих пор симулятор. Эта задача сложнее. Не получилось уйти в плавсостав, но оно и к лучшему, возможно. Точнее получалось - но сварщиком в палубную команду, а это значит, в моём. . .	Ритм жизни kumehtar 27.02.2026 Иногда приходится жить в ритме, где дел становится всё больше, а вовлечения в происходящее — всё меньше. Плотный график не даёт вниманию закрепиться ни на одном событии. Утро начинается с быстрых,. . .

@Krasme 7247 / 5100 / 2151 Регистрация: 02.02.2014 Сообщений: 13,461
	21.09.2021, 11:23
	1. приложите книгу 2. (наскоро глянула) вы сначала находите z как массив, потом пытаетесь оптимизировать расхождение двух готовых таблиц... 0

@Krasme 7247 / 5100 / 2151 Регистрация: 02.02.2014 Сообщений: 13,461
	21.09.2021, 17:54
	книга - Книга1.xlsx я ведь не могу экспериментировать с вашей задачей, не имея внешние данные 0

@Krasme 7247 / 5100 / 2151 Регистрация: 02.02.2014 Сообщений: 13,461
	22.09.2021, 15:47
	т.к. у вас 6-ой прайм, у меня 7-й, могу дать только pdf-результат kinetika.pdf считает долго, надо набраться терпения 0

@Sawarez 0 / 0 / 0 Регистрация: 20.09.2021 Сообщений: 11
	22.09.2021, 16:15 [ТС]
	Спасибо огромное! У меня на самом деле тоже 7й, просто тут в форуме в списке я видел только 6й максимум. Если это все переписать, он по идее будет работать? По графикам конечно видно что они сильно расходятся... Не знаете, можно ли как-то улучшить оптимизацию или это просто невозможно? Вероятно, что могут происходить некоторые иные процессы в системе, соответственно, систему диффуров надо усложнять или как-то еще менять + эксперимент. данные не самые лучшие) Но если сама "программа" будет работать, то не составит труда заменить данные или уравнения. В любом случае - спасибо огромное еще раз, очень выручаете! 0

@Krasme 7247 / 5100 / 2151 Регистрация: 02.02.2014 Сообщений: 13,461
	22.09.2021, 16:23
	раз седьмой, держите файл kinetika.zip по части оптимизации, маткад выдал результат по заданным уравнениям. вероятно, в системе диффуров не те уравнения должны быть. я ограничения по k выключила, маткад дал те же результаты. с начальными приближениями попробовала поиграть, терпения не хватило дожидаться)) 0

@Krasme
	16.10.2021, 16:22
	Не по теме: для меня кинетика - темный лес.. я не могу оценить, насколько достоверны ваши формулы кинетики, помимо них еще надо знать, в каких диапазонах должны быть искомые параметры, не только k1,k2,k3. есть вероятность, что встроенные функции оптимизации не смогут найти вам оптимальное решение. надо искать другие способы, возможно, даже перестроить саму постановку задачи. это возможно сделать, только будучи сведущим в теме кинетики. здесь я вам помочь ничем не могу. 0