2014 dxdy logo

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

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


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


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



Начать новую тему Ответить на тему
 
 Не сходятся некоторые значения в м. Адамса с точным решением
Сообщение02.11.2008, 11:42 


27/09/08
18
Решаю методом Адамса второго порядка по формуле: $y_{i+1}=y_{i}+\frac{h}{2}(3f(x_{i},y_{i})-f(x_{i-1},y_{i-1}))$
$f(x,y)=-2xy+xe^{-x^{2}}$
h = 0.01, x_{i+1} = x_{i} + h
Точное решение:
$u=\frac{e^{-x^{2}}x^{2}}{2}$

Результат
1 - Точное решение
2 - Метод Адамса

Код:
      (1)                (2)
  0             0         
  0.00124688    0.00124736
  0.00495025    0.00495242
*0.0109997     0.0110047 
  0.0192158     0.0192246 
  0.0293567     0.02937   
  0.0411269     0.0411454 
*0.0541882     0.0542122 
*0.0681715     0.0682009 
*0.0826895     0.0827241 
  0.0973501     0.0973894 
*0.111769      0.111812   
*0.125582      0.125628   
*0.138455      0.138503   
*0.150093      0.150143   
*0.160251      0.1603     
  0.168734      0.168781   
  0.1754        0.175446   
*0.180168      0.18021   
  0.183006      0.183045   

* - не совпадают 4 знака после запятой

Может быть в формуле $y_{i+1}$ ошибся?

 Профиль  
                  
 
 
Сообщение04.11.2008, 14:27 
Заслуженный участник


11/05/08
32166
а почему они в точности должны совпадать? 2-й порядок метода означает лишь, что при шаге в одну сотую и решении порядка единицы погрешность должна быть примерно одна десятитысячная, там плюс-минус порядок (и плюс к тому смотря ещё как первый шаг аппроксимировался)

 Профиль  
                  
 
 
Сообщение04.11.2008, 16:20 


27/09/08
18
ewert, первый шаг брался с метода Рунге-Кутта.

 Профиль  
                  
 
 
Сообщение04.11.2008, 16:30 
Заслуженный участник


11/05/08
32166
тогда этот шаг можно условно считать точным, а всё остальное остаётся в силе

 Профиль  
                  
 
 
Сообщение04.11.2008, 18:44 


27/09/08
18
ewert, Спасибо, теперь понятно.
Вопрос решён!

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

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



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

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


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

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