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

Форум пользователей пакета Maple

Модератор: Admin

Vladislove1976
Сообщения: 6
Зарегистрирован: Вт дек 08, 2015 11:53 pm

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

Сообщение Vladislove1976 » Ср дек 09, 2015 7:46 pm

Добрый день.Требуется консультация по след. проблеме. Имеется система трансцендентных уравнений вида
h=w1(x,y), h=w2(x,y). Необходимо, меняя h с заданным шагом в интервале, скажем, (0,1), выводить решения этой системы x(i)=x(h(i)), y(i)=y(h(i)) в заданном интервале (x1..x2,y1..y1) при каждом h(i), а также значения заданной функции q(i)=q(h(i))=q(x(h(i)),y(h(i))). Желательно вывести на график полученные зависимости x(h), y(h), q(h).
Есть также нюанс - выводить нужно ближайшее по одной переменной к нулю значение, проверяя выполнение условия устойчивости s(x(i),y(i))>0 в противном случае - переходить к следующему корню.
Насколько знаю, для решения можно использовать встроенную в Maple процедуру численного решения fsolve, но вопрос в том, как ее правильно встроить в цикл. Заранее благодарен

алексей_алексей
Сообщения: 1776
Зарегистрирован: Вс май 01, 2005 9:02 pm

Re: Решение системы трансцендентных уравнений в цикле

Сообщение алексей_алексей » Ср дек 09, 2015 9:39 pm

w
Последний раз редактировалось алексей_алексей Чт июн 16, 2016 7:07 pm, всего редактировалось 1 раз.
не педераст

Vladislove1976
Сообщения: 6
Зарегистрирован: Вт дек 08, 2015 11:53 pm

Сообщение Vladislove1976 » Ср дек 09, 2015 9:45 pm

Немного не так. Здесь h - не неизвестная, а параметр, который нужно изменять и на каждом шаге получать решения системы (x,y)

алексей_алексей
Сообщения: 1776
Зарегистрирован: Вс май 01, 2005 9:02 pm

Сообщение алексей_алексей » Ср дек 09, 2015 9:51 pm

f
Последний раз редактировалось алексей_алексей Чт июн 16, 2016 7:08 pm, всего редактировалось 1 раз.
не педераст

Vladislove1976
Сообщения: 6
Зарегистрирован: Вт дек 08, 2015 11:53 pm

Сообщение Vladislove1976 » Ср дек 09, 2015 10:12 pm

Алексей Алексей, спасибо за ссылочку, сразу не заметил.

Vladislove1976
Сообщения: 6
Зарегистрирован: Вт дек 08, 2015 11:53 pm

Сообщение Vladislove1976 » Ср дек 09, 2015 10:17 pm

Честно говоря, у меня подобные идеи решения этой системы возникали:
∆h=w_1θ (θ,φ)∆θ+w_1φ (θ,φ)∆φ
∆h=w_2θ (θ,φ)∆θ+w_2φ (θ,φ)∆φ
Отсюда
∆θ=(w_2φ-w_1φ)/(w_1θ w_2φ-w_2θ w_1φ ) ∆h,∆φ=(w_1θ-w_2θ)/(w_1θ w_2φ-w_2θ w_1φ ) ∆h
Но, может быть, в Maple есть уже готовые процедуры, которые мне нужно просто встроить в цикл?