Поиск:

Ответ в темуСоздание новой темы Создание опроса
> Программа для решения систем уравнений численно 
:(
    Опции темы
MathFortran
  Дата 22.4.2014, 13:14 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 4
Регистрация: 22.4.2014

Репутация: нет
Всего: нет



Добрый день!
Есть система алгебраических уравнений, которую надо решать численно (методом Ньютона, например).
Я нашел программу здесь: http://www.zib.de/en/numerik/software/ant/nleq1.html
Но никак не могу разобраться как внутри ввести нужные мне уравнения.
Буду рад любой помощи. В частности, если кто-нибудь подскажет, где найти более прозрачную (user-friendly) программку.

PM MAIL   Вверх
Фантом
Дата 22.4.2014, 16:09 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 5
Всего: 49



Вообще говоря, метод Ньютона проще написать самому, это исходник на пару килобайт. 

Сколько уравнений решать надо?
PM   Вверх
MathFortran
Дата 22.4.2014, 16:15 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 4
Регистрация: 22.4.2014

Репутация: нет
Всего: нет



Там система уравнений химической кинетики для нахождения равновесия.
Такие системы получаются жесткие, поэтому в лоб обычно метода не проходит.
В принципе, судя по статьям солвер в ссылке достаточно хороший просто бестолково документация написана   smile 
PM MAIL   Вверх
Фантом
Дата 22.4.2014, 16:46 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 5
Всего: 49



Цитата(MathFortran @  22.4.2014,  17:15 Найти цитируемый пост)
Там система уравнений химической кинетики для нахождения равновесия.
Такие системы получаются жесткие, поэтому в лоб обычно метода не проходит.

Уравнения химической кинетики - дифференциальные, а не алгебраические. Жесткие системы, естественно, тоже. Поэтому применение к ним метода Ньютона - это что-то, мягко говоря, странное.
PM   Вверх
MathFortran
Дата 22.4.2014, 17:18 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 4
Регистрация: 22.4.2014

Репутация: нет
Всего: нет



Чтобы найти равновесные концентрации, производные от концентраций по времени приравнивают к нулю.
Отсюда получается система алгебраических уравнений, которая в целом тоже достаточно неприятная в вычислительном плане.
PM MAIL   Вверх
Фантом
Дата 22.4.2014, 22:26 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 5
Всего: 49



Цитата(MathFortran @  22.4.2014,  18:18 Найти цитируемый пост)
Чтобы найти равновесные концентрации, производные от концентраций по времени приравнивают к нулю.


Тогда это не совсем кинетика.  smile Ну ладно, с этим понятно. А все-таки - сколько уравнений в системе?
PM   Вверх
MathFortran
Дата 23.4.2014, 08:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 4
Регистрация: 22.4.2014

Репутация: нет
Всего: нет



немного: 5-6 уравнений
PM MAIL   Вверх
Фантом
Дата 23.4.2014, 15:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 5
Всего: 49



Тогда это стрельба из пушки по воробьям, проще написать свою. Если есть возможность, дайте ссылку на уравнения и величины параметров.
PM   Вверх
tzirechnoy
Дата 24.4.2014, 23:08 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
***


Профиль
Группа: Завсегдатай
Сообщений: 1173
Регистрация: 30.1.2009

Репутация: нет
Всего: 16



Честно говоря, по-моему проще как раз наоборот не писать свою, а взять какой-нибудь octave. Или scilab. В котором это делается в число уравнений плюс 3 строчек, и описано в примерах для начинающих.
PM MAIL   Вверх
Фантом
Дата 25.4.2014, 00:50 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Вы это прекратите!
***


Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 5
Всего: 49



Цитата(tzirechnoy @  25.4.2014,  00:08 Найти цитируемый пост)
Честно говоря, по-моему проще как раз наоборот не писать свою, а взять какой-нибудь octave.

Своя будет не намного длиннее.  smile 
PM   Вверх
kovalkos
Дата 30.4.2014, 07:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 2
Регистрация: 29.4.2014
Где: Луганск

Репутация: нет
Всего: нет



Решай методом Гаусса:)) Если метод Ньютона не подходит:))
PM MAIL WWW ICQ Skype   Вверх
Santik
Дата 1.3.2016, 07:53 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Новичок



Профиль
Группа: Участник
Сообщений: 45
Регистрация: 13.3.2012
Где: Мирный (Якутия)

Репутация: нет
Всего: нет



А библиотеку IMSL почему не используем???
PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Fortran | Следующая тема »


 




[ Время генерации скрипта: 0.1285 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.