2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки


Правила форума


В этом разделе нельзя создавать новые темы.



Начать новую тему Ответить на тему
 
 Как обойти неустойчивость
Сообщение13.02.2012, 15:45 


29/12/09
366
Попытаюсь кратко, описать суть вопроса.
Если откинуть все детали в уравнении(дополнительные коэффициенты и т.п.), то получается нужно решить такое уравнение с краевыми условиями:
$y''(x)=F(y)$
$ \phi(y(a),y'(a))=0;$
$\psi(y(b),y'(b))=0;$
Задача, состоит в том чтобы решить это уравнение численно.
Записываю для него разносную схему, причем линеаризую, правую часть$F(y^{i+1}_n)=F(y^i_n)+F_y(y^i_n)(y^{i+1}_n-y^{i}_n)$:
$y_{n-1}-2y_n-y_{n+1}=h^2(F(y^i_n)+F_y(y^i_n)(y^{i+1}_n-y^{i}_n))$
Далее преобразую к диагональному виду, чтобы решить методом прогонки:
$a_ny_{n-1}+b_ny_n+c_ny_{n+1}=d_n$ ,где коэффициенты получаются такими:
$a_n=1;$
$b_n=-(2+h^2F_y(y^i_n));$
$c_n=1;$
$d_n=h^2(F(y^i_n)-F_y(y^i_n)y^{i}_n);$
условие корректности прогонки $|b_n|\ge|a_n|+|c_n|$
Если $F_y(y^i_n)<0$ и $-4<h^2F_y(y^i_n)<0$, то получается, что метод прогонки будет неустойчив, как обойти эту неприятность? Помогите пожалуйста!

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ 1 сообщение ] 

Модераторы: Модераторы Математики, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group