50947

Метод прогонки розвязання крайової задачі. Складання алгоритму

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

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

Мета. Навчитися використовувати метод прогонки розв’язання крайової задачі звичайного диференційного рівняння. Скласти алгоритм. Устаткування: папір формату А4, ручка, калькулятор, С++.

Украинкский

2014-02-02

40.5 KB

3 чел.

Лабораторна робота №29-30

Тема. Метод прогонки розв’язання крайової задачі. Складання алгоритму. 

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

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

Хід роботи

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

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

12)

#include<iostream.h>

void main()

{

int n=20,i;

double h,b,a,A,B,alpha0,alpha1,beta0,beta1,n1;

double*x=new double[n],*y=new double[n],*q=new double[n];

double*m=new double[n];

double*p=new double[n];

double*N=new double[n];

double*F=new double[n];

double*f=new double[n];

double*c=new double[n];

double*d=new double[n];

cout<<"INPUT a,b,A,B,alpha0,alpha1,beta0,beta1,h : \n";

cin>>a>>b>>A>>B>>alpha0>>alpha1>>beta0>>beta1>>h;

n=(b-a)/h;

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

 {  x[i]=a+i*h;}

c[0]=alpha1/(alpha0*h-alpha1);

d[0]=A*h/alpha1;

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

 {p[i]=1/x[i];

 q[i]=2;

 f[i]=x[i];

 m[i]=(2*h*h*q[i]-4)/(2+h*p[i]);

 N[i]=(2-h*p[i])/(2+h*p[i]);

 F[i]=(2*f[i])/(2+h*f[i]);

 c[i]=1.0/(m[i]-N[i]*c[i-1]);

 d[i]=F[i]*h*h-N[i]*c[i-1]*d[i-1];}

y[n]=(B*h+beta1*c[n-1]*d[n-1])/(beta0*h+beta1*(c[n-1]+1));

for(i=n-1;i>=0;i--)

 {y[i]=c[i]*(d[i]-y[i+1]);}

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

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

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

  •  Сформулюйте крайову задачу, що вирішується методом прогонки
  •  Які перетворення виконують над областю задання диференційного рівняння?
  •  Які перетворення виконуються над крайовими умовами?
  •  Що є кінцевою метою методу сіток?
  •  З чого складається відповідь до задачі?


 

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

3380. Элементы кинематики 359 KB
  Лекция 1. Элементы кинематики. Введение. Основные кинематические понятия и характеристики. Нормальное, тангенциальное и полное ускорения. Угловая скорость, угловое ускорение. 1. Введение. Физику можно назвать наукой о н...
3381. Обработка результатов измерений в физическом практикуме 303.5 KB
  Содержит изложение методики обработки результатов измерений, получаемых во время практических занятий в учебной лаборатории. Предназначен для оказания помощи студентам технических специальностей всех форм обучения при подготовке к лабораторным работ...
3382. Потенциал электростатического поля. Диэлектрики в электростатическом поле 350 KB
  Потенциал электростатического поля. Диэлектрики в электростатическом поле. Работа при перемещении заряда в электростатическом поле. Потенциал. Разность потенциалов. Связь между напряженностью и потенциалом электростатического поля. Эквип...
3383. Изучение свободных колебаний математического и пружинного маятников 398.5 KB
  Изучение свободных колебаний математического и пружинного маятников Цель работы: изучение физических основ свободных незатухающих колебаний, определение ускорения свободного падения с помощью математического маятника и коэффициента упругости пружины...
3384. Книжные социальные сети как канал распространения издательской и книготорговой библиографической информации 104 KB
  Книжные социальные сети как канал распространения издательской и книготорговой библиографической информации На сегодняшний день в мире существует более 130 миллионов компьютеров и более 80 % из них объединены в различные информационные сети от малых...
3385. Проектирование фундаментов производственных зданий 606.5 KB
  ИСХОДНЫЕ ДАННЫЕ - тип объекта – производственное здание с подвалом и гибкой конструктивной схемой    - район строительства – г. Магнитогорск; -величины нагрузок на фундаменты представлены в табл. 1. Таблица 1 Величины нагрузок...
3386. Рассчитать и спроектировать резец для обработки наружной поверхности детали 8.03 MB
  Оправки в шпинделе закрепляют штоком, проходящим через шпиндель станка. Шток имеет на конце захватное устройство. Инструментальные оправки имеют соответствующие этому устройству наружные, внутренние или резьбовые поверхности захвата...
3387. Реконструкция многоквартирного крупноблочного дома серии 1-439А 888.5 KB
  При модернизации и реконструкции жилых зданий массовой застройки предусматривается решение следующих задач: приведение планировочной структуры здания в соответствие с требованиями к потребительским и эксплуатационным качествам современного ...
3388. Проект разработки роторного снегоочистителя 857.5 KB
  В процессе подготовки будущего инженера к самостоятельному решению технических и производственных задач одно из ведущих мест принадлежит курсовому проектированию. Цель данного курсового проекта – закрепить и обобщить теоретический мате...