35235

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

Практическая работа

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

h double Fdouble x double y { return cos2xy1.5xy; } void min {int n; double hb; doublek=new double [4]; doubleq=new double[n1]; doubledq1=new double[n1]; doubledq2=new double[n1]; doubledq3=new double[n1]; doublex=new double[n1]; doubley=new double[n1]; cout Vvedite bh endl; cin ; cin b; cin h; cout Vvedite y[0] endl; cin y[0]; n=b h; x[0]=; cout x y ; cout endl; cout ; cout endl; for int i=0; i =2; i { k[0]=hFx[i]y[i]; k[1]=hFx[i]h 2y[i]k[0] 2;...

Украинкский

2013-09-09

42 KB

0 чел.

7.04.09  Мишура Денис, гр. 1ПМ-06

 Вариант 11

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

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

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

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

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

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

11)

Текст програми:

#include <iostream.h>

#include <math.h>

double F(double x, double y)

{

return (cos(2*x+y)+1.5*(x-y));

}

void main()

{int n;

double h,a,b;

double*k=new double [4];

double*q=new double[n+1];

double*dq1=new double[n+1];

double*dq2=new double[n+1];

double*dq3=new double[n+1];

double*x=new double[n+1];

double*y=new double[n+1];

cout<<"Vvedite a,b,h "<<endl;

cin>>a;

cin>>b;

cin>>h;

cout<<"Vvedite y[0] "<<endl;

cin>>y[0];

n=(b-a)/h;

x[0]=a;

cout<<"    x    |   y  ";

cout<<endl;

cout<<"-------------------";

cout<<endl;

for (int i=0; i<=2; i++)

{

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

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

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

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

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

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

 q[i]=h*F(x[i],y[i]);

 cout<<"  "<<x[i]<<"    |  "<<y[i]<<"\n";

}

cout<<endl;

for (i=3; i<=n+1;i++)

{

q[i]=h*F(x[i],y[i]);

for (int k=0;k<i-1;k++) dq1[k]=q[k+1]-q[k];

for (k=0;k<i-2;k++) dq2[k]=dq1[k+1]-dq1[k];

for (k=0;k<i-3;k++) dq3[k]=dq2[k+1]-dq2[k];

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

y[i+1]=y[i]+q[i]+dq1[i-1]/2+5*dq2[i-2]/12+3*dq3[i-3]/8;

cout<<"  "<<x[i]<<"    |  "<<y[i]<<"\n";

}

cout<<endl;

}

  1.  Контрольні питання

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

PAGE  1


 

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

39363. Тарифная политика страховщика на региональном рынке страхования (по различным видам) 44.17 KB
  Только правильно сформированная и реализованная тарифная политика сможет обеспечить выполнение целей и задач страховой компании: рентабельность деятельности и общедоступность страхования для населения, организаций и государства.
39364. Расчет и основные характеристики привода ленточного конвейера 941 KB
  Мощность двигателя зависит от требуемой мощности рабочей машины а его частота вращения от частоты вращения приводного вала рабочей машины.2 Определение передаточного числа привода и его ступеней Передаточное число привода определяется отношением номинальной частоты вращения двигателя к частоте вращения приводного вала рабочей машины при номинальной нагрузке и равно произведению передаточных чисел закрытой и открытой передач.1 Частота вращения приводного вала рабочей машины 6 где: скорость тягового органа м с; диаметр...
39365. Проект детского кафе на 50 мест в Торгово- офисном центре г. Пушкино 1.34 MB
  Разработка проекта детского кафе проведена в здании реально существующего Торгово-офисного центра «ВИТ», расположенного в г. Пушкино, ул. Чехова 12. Привлекательность реализации данного проекта обосновывается положительными прогнозами экспертов относительно роста численности целевой аудитории детских кафе в России в последующие годы.
39366. Социальная организация понятие, признаки и функции 84.5 KB
  Социальная организация — это социальная система, которая характеризуется определенной коллективной тождественностью (идентичностью), имеет точный список членов, программу деятельности и процедуру перемещения (или замещения) членов.
39367. Основные данные и расчет привода ленточного конвейера 905 KB
  2 Определяем общий коэффициент полезного действия КПД привода по формуле 3 где: коэффициент полезного действия цилиндрической передачи 096; коэффициент полезного действия червячной передачи 08; коэффициент полезного действия открытой муфты 098; коэффициент полезного действия пары подшипников 099 Определяем общий КПД 2. Примем стандартное передаточное число червячной передачи тогда 9 где: передаточное число червячной передачи 20 2.1 Определяем мощности а двигателя б быстроходного вала цилиндрической...
39368. Привод электрической лебедки 852.5 KB
  Мощность двигателя зависит от требуемой мощности рабочей машины а его частота вращения от частоты вращения приводного вала рабочей машины.2 Определение передаточного числа привода и его ступеней Передаточное число привода определяется отношением номинальной частоты вращения двигателя к частоте вращения приводного вала рабочей машины при номинальной нагрузке и равно произведению передаточных чисел закрытой и открытой передач.1 Определяем мощности а двигателя б быстроходного вала редуктора 7 в тихоходного вала редуктора 8 г рабочей...
39369. Россия в начале XX века: революция или реформы 78 KB
  Россия не являлась одним из промышленных или финансовых гигантов, противоречия между которыми привели к войне Россия была в целом заинтересована в сохранении территориально-политического раздела мира Россия обладала огромным военным потенциалом, значительными сырьевыми ресурсами и пользовалась большим авторитетом на международной арене
39370. Проектирование и расчет цилиндрического редуктора 564 KB
  Расчет электродвигателя и кинематический расчет привода. Расчет зубчатой передачи редуктора. Предварительный расчет валов редуктора. Конструктивные размеры шестерни и колеса. Конструктивные размеры корпуса редуктора. Первый этап компоновки редуктора
39371. Основы генетики и селекций 29.85 KB
  Ген (от греч. genos — род, происхождение) – участок молекулы геномной нуклеиновой кислоты, характеризуемый специфической для него последовательностью нуклеотидов, представляющий единицу функции, отличной от функций других генов, и способный изменяться путем мутирования.