50946

Екстраполяційний метод Адамса розвязання задачі Коші

Лабораторная работа

Информатика, кибернетика и программирование

Мета. Навчитися знаходити розвязок диференційного рівняння екстраполяційним методом Адамса. Устаткування: папір формату А4, ручка, калькулятор, ПЗ С ++. Хід роботи Правила техніки безпеки Теоретичні дані Індивідуальне завдання. Використовуючи метод Адамса з трьома кінцевими різницями, скласти таблицю наближених значень інтеграла диференційного рівняння, з початковими умовами на відрізку з точністю 0,001. Початковий відрізок встановити методом Рунге-Кутта.

Украинкский

2014-02-02

41 KB

6 чел.

Лабораторна робота №28

Тема. Екстраполяційний метод Адамса розвязання задачі Коші.

Мета. Навчитися знаходити розв’язок диференційного рівняння екстраполяційним методом Адамса.

Устаткування: папір формату А4, ручка, калькулятор, ПЗ  С ++.

Хід роботи

  1.  Правила техніки безпеки
  2.  Теоретичні дані

      Індивідуальне завдання.

Використовуючи метод Адамса з трьома кінцевими різницями, скласти таблицю наближених  значень інтеграла диференційного рівняння , з початковими умовами на відрізку  з точністю 0,001. Початковий відрізок встановити методом Рунге-Кутта.

12)

#include<iostream.h>

#include<math.h>

double f (double x, double y)

{return cos(y)/(1.25+x)-0.1*y*y;}

int i,j;

double h;

double x[20]; 

double y[20];

double dy[20];

double k[20][20];

double q[20][20];

void q_f()

{

int fi,fj;

for(fi=0;fi<20;fi++)

{q[0][fi]=h*f(x[fi],y[fi]);}

for(fj=1;fj<20;fj++)

{for(fi=0;fi<20;fi++){q[fj][fi]=q[fj-1][fi+1]-q[fj-1][fi];}}}

void main()

{cout<<"Input x[0] = ";

cin>>x[0];

cout<<"Input y[0] = ";

cin>>y[0];

cout<<"Input step h = ";

cin>>h;

for(i=0;x[i]<=1;i++){ x[i+1] =x[i]+h;}

int i_max=i;

for(i=0;i<3;i++)

{  k[0][i]=h*f(x[i],y[i]);

 k[1][i]=h*f((x[i]+h/2),(y[i]+k[0][i]/2));

 k[2][i]=h*f((x[i]+h/2),(y[i]+k[1][i]/2));

 k[3][i]=h*f((x[i]+h),(y[i]+k[2][i]));

 dy[i]=(k[0][i]+2*k[1][i]+2*k[2][i]+k[3][i])/6;

 y[i+1] =y[i]+dy[i];

 x[i+1] =x[i]+h;}

for(i=0;i<i_max;i++)

{q_f();

dy[i]=q[0][i]+.5*q[1][i-1]+(5/12)*q[2][i-2]+(3/8)*q[3][i-3];

y[i+1] =y[i]+dy[i];

x[i+1] =x[i]+h;

cout<<i<<"  "<<x[i];}

for(i=0;i<i_max;i++)

{ cout<<"\nx["<<i<<"] = "<<x[i]<<"        y["<<i<<"] = "<<y[i];}}

  1.  Контрольні питання
  •  Для якої задачі використовується метод Адамса?
  •  Які формули методу Адамса?
  •  Як знайти початковий відрізок?
  •  Що являється результатом рішення задачі методом Адамса?
  •  Як оцінити похибку метода на практиці?


 

А также другие работы, которые могут Вас заинтересовать

3259. Техническая эксплуатация систем горячего водоснабжения 64.5 KB
  Техническая эксплуатация систем горячего водоснабжения Назначение, классификация и устройство систем горячего водоснабжения зданий. Техническое обслуживание и ремонт систем горячего водоснабжения зданий. Назначение, классификация и устройство сист...
3260. Обжиг сырьевой смеси 57.5 KB
  Обжиг сырьевой смеси Для обжига сырьевой смеси как при мокром, так и при сухом способе производства почти исключительно применяют вращающиеся печи. При сухом способе производства иногда используют шахтные печи. Длина современных вращающихся печей пр...
3261. Демократия как политический режим 39.66 KB
  Сущность демократического политического режима В переводе с греческого демократия означает власть народа (demos - народ, cratos - власть). Более развернутое определение демократии, ставшее классическим, было дано американским президентом...
3262. Сущность права 59.76 KB
  Одним из важнейших признаков права является формальная определенность, характеризующая его ясность, недвусмысленность, точность, структурированность, устойчивость во времени и по кругу лиц. Если эти качества не обеспечиваются, то возникает ...
3263. Повышение мотивации труда 121.54 KB
  Результаты изучения моделей мотивации не позволяют с психологической точки зрения четко определить, что же побуждает человека к труду. Изучение человека и его поведения в процессе труда дает только некоторые общие объяснения мотивации, но даже они позволяют разрабатывать прагматические модели мотивации работника на конкретном рабочем месте.
3264. Европейский Союз: особенности формирования, этапы и перспективы развития 74.19 KB
  Закономерным результатом развития международных экономических отношений, а именной международной торговли и международного движения факторов производства, стала экономическая интеграция, являющаяся особым этапом интернационализации хозяйст...
3265. Возникновение планетных систем и земли 96.38 KB
  Согласитесь, сегодня человек, в какой бы самой отдаленной области науки или народного хозяйства он ни работал, должен иметь представления, хотя бы общее, о нашей Солнечной системе, звездах и современных достижениях астрономии. Сравнительное...
3266. Агрессия как социально-психологический феномен 60.36 KB
  Ни одно общество не свободно от таких явлений, как убийство и разрушение. Наше непосредственное окружение полно красноречивых сцен: омерзительная пьяная драка на улице, убийство из ревности, ограбление. Насилие устрашает, повергает в смятен...
3267. Природно-ресурсный потенциал мирового хозяйства 53.74 KB
  Современный этап развития мирового хозяйства отличается всевозрастающими масштабами потребления природных ресурсов, резким усложнением процесса взаимодействия природы и общества, интенсификацией и расширением сферы проявления специфических...