VPF::Метод Зейделя для системы НЕлинейных уравнений

norman

Дата 30.4.2006, 15:21 (ссылка)

(нет голосов)

Загрузка ...

Новичок

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

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

Писал программу на TP, работать она естественно отказывалась...думаю, что всё дело в алгоритме, помогите найти правильное решение...

Этот метод основан на процессе итерации...как организовать процесс для одного нелинейного уравнения я знаю:

Код

while abs(f(x))>=e do
begin
if f(x)>0 then x:=x-0.1*f(x)
else x:=x+0.1*f(x);
end;

А как для системы...понятия не имею smile

Читал книгу для дипломированных специалистов по прикладной математике, мало чё понял...там вобщем так:
Дана система:
F1(x,y)=0
F2(x,y)=0
Потом надо выразить одно уравнение через Х, а другое через Y...вот примерно так:
x=G1(x,y)
y=G2(x,y)
Вот, я дальше ничего не ясно...надеюсь на вашу помощь...
Благодарен за любую информацию...

Это сообщение отредактировал(а) norman - 30.4.2006, 15:28

bsa

Дата 1.5.2006, 16:33 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Модератор
Сообщений: 9185
Регистрация: 6.4.2006
Где: Москва, Россия

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

Есть несколько интересных ссылок:
http://www.google.com/search?rls=ru&q=...-8&oe=utf-8

--------------------

Правильно заданный вопрос - половина ответа
Отзывы о работе модератора.

hypocrisy

Дата 5.5.2006, 11:11 (ссылка)

(нет голосов)

Загрузка ...

Новичок

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

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

если задано
x=G1(x,y)
y=G2(x,y)
то:

Код


G1_temp:=1;
G2_temp:=1;
while abs(G1_temp*G2_temp)>=e do    
begin    

G1_temp=G1(x,y);
if f1_temp>0 then x:=x-0.1*G1_temp;
else x:=x+0.1*G1_temp;

G2_temp=G2(x,y);
if G2_temp>0 then y:=y-0.1*G2_temp;
else y:=y+0.1*f2_temp;

end;

Модератор: пользуемся кнопкой "Код"

Это сообщение отредактировал(а) maxim1000 - 5.5.2006, 12:49

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Алгоритмы \| Следующая тема »