35234

Метод Рунге-Кутта вирішення задачі Коші. Складання алгоритму

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

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

Навчитися вирішувати задачу Коші методом Рунге-Кутта; скласти алгоритм.

Украинкский

2013-09-09

37.5 KB

2 чел.

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

 Вариант 11

Лабораторна робота №26-27

Тема: Метод Рунге-Кутта вирішення задачі Коші. Складання алгоритму.

Мета: Навчитися вирішувати задачу Коші методом Рунге-Кутта; скласти алгоритм.

Устаткування: папір формату А4, ПК, програмне забезпечення Borland С++.

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

Вирішити задачу, використовуючи метод Рунге-Кутта для диференціального рівняння  с початковою умовою , обчислити з точністю до  рішення цього рівняння при , . Обчислення провести із двома запасними знаками.

11)

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

#include<iostream.h>

#include<math.h>

void main()

{int i,n;

double h,a,b,dely;

double*k = new double [4];

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(i=0;i<=n;i++)

{

k[0]=h*(0.3*x[i]+y[i]*y[i]);

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

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

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

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

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

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

}

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

cout<<endl;

}

  1.  Контрольні питання
  •  Поставте задачу, розв’язувану методом Рунге-Кутта.
  •  Сформулюйте алгоритм метода Рунге-Кутта.
  •  Який порядок точності методу Рунге-Кутта?
  •  Як практично досягти заданого ступеня точності?
  1.  Звіт про роботу.
  2.  Захист роботи.


 

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

49495. Определение в планируемом периоде количества ТО и КР 287.5 KB
  Каждому типу машин присуще свое определенное распределение трудоемкости по видам работ. Удельный вес видов работ в общем, объеме трудоемкость остается без существенных изменений, несмотря на совершенствование технологии ремонта и снижение общих трудозатрат на ремонт машин данного типа.
49496. Разработка стенда для диагностирования системы охлаждения 1.15 MB
  Ремонтно-механические мастерские, как правило, работают в одну смену, и только при большой загрузке и в целях лучшего использования дорогостоящего оборудования механические отделения и некоторые другие участки иногда работают в две смены.
49497. Проект ОКС 7 на ГТС с УВС и УИС 717 KB
  ОКС7 предоставляет универсальную структуру для организации сигнализации сообщений сетевого взаимодействия и технического обслуживания телефонной сети. SS5 и более ранние версии использовали принцип сигнализации в линии где информация необходимая для соединения передавалась специальными тонами DTMF в телефонной линии известной как Bканал. Такой тип сигнализации создавал уязвимость в безопасности протокола поскольку злоумышленник мог эмулировать набор служебных...
49498. Система автоматического регулирования частоты вращения ДПТ 857 KB
  Область применения системы. Принцип работы системы. Передаточные функции системы. Анализ структурной устойчивости САР 20 Коэффициент усиления системы в разомкнутом состоянии.
49499. РАЗРАБОТКА ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА ИЗГОТОВЛЕНИЯ ДЕТАЛИ 671 KB
  Физическая сущность процесса сварки заключается в образовании прочных связей между атомами или молекулами на соединяемых поверхностях заготовок. Экономическая эффективность применения сварки по сравнению с механическими способами соединения деталей и литьем заключается в экономии металла снижении трудоемкости работ и технологической гибкости процесса.д Все способы сварки условно делятся на две группы. К первой относятся способы сварки при которых соединение получается за счет расплавления металла.
49501. РАЗРАБОТКА ТЕХНОЛОГИЧЕСКОГО ПРОЦЕССА ИЗГОТОВЛЕНИЯ ДЕТАЛИ 416 KB
  Технологический процесс изготовления заготовок. Технологический процесс изготовления детали. технологический процесс изготовления детали.
49502. Микропроцессорная система управления объектом с переменным запаздыванием 1.36 MB
  В тех случаях, когда к качеству переходных процессов в системе регулирования предъявляют высокие требования, целесообразно использовать специальные регуляторы или алгоритмы, обеспечивающие компенсацию чистого запаздывания. Примерами такого регулятора является регулятор Смита. Он представляет собой замкнутую систему, состоящую из обычного регулятора и модели объекта, включенной в линию обратной связи к регулятору.