35257

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

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

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

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

Украинкский

2013-09-09

29.5 KB

1 чел.

еревко О.В.                                                                                                              2ПМС-07

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

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

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

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

Хід роботи

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

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

16)

Программа на С++

#include<iostream.h>

#include<math.h>

double PX(double x)

{

return -3;

}

double QX(double x)

{

return 1/x;

}

double FX(double x)

{

return 1;

}

void main()

{

int n;

double h,a,b,LO,L1,A,BO,B1,B;

cout<<"Vvedite shag h:";

cin>>h;

cout<<"Vvedite a: ";

cin>>a;

cout<<"Vvedite b: ";

cin>>b;

n=(b-a)/h;

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

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

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

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

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

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

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

cout<<"Vvedite L[0]: ";

cin>>L0;

cout<<"Vvedite L[1]: ";

cin>>L1;

cout<<"Vvedite A: ";

cin>>A;

cout<<"Vvedite B[0]: ";

cin>>B0;

cout<<"Vvedite B[1]: ";

cin>>B1;

cout<<"Vvedite B: ";

cin>>B;

x[0]=a;

c[0]=L1/(L0*h-L1);

d[0]=A*h/Ll;

//Prjamoi hod

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

{

x[i]=x[0]+i*h;

m[i]=2*h*h*QX(x[i])-4)/(2+h*PX(x[i]));

nn[i]=(2-h*PX(x[i]))/(2+h*PX(x[i]));

f[i]=2*FX(x[i])/(2+h*PX(x[i]));

if(i!=0 {c[i]=1/(m[i]-nn[i]*c[i-1]);

d[i]=FX(x[i])*h*h-nn[i]*c[i-1]*d[i-1];}

}

//Obratnii hod

y[n]=(B*h+Bl*c[n-1]*d[n-1])/(BO*h+B1*(c[n-1]+1));

for (i=n-1;i>=0;i--) y[i]=c[i]*(d[i]-y[i+1]);

cout<<"Tablichnaja fynkchija:\n";

for(i=0;i<=n;i++) cout<<"x["<<i<<"]="<<x[i]<<"\ty["<<i<<"]="<<y[i]<<"\n";

}

Результат программы


 

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

46338. Обзор деятельности компании Nissan Motor Co 267 KB
  Компания Nissan Motor Co была создана в декабре 1933 года и получила свое современное название в 1934 году. Первый компактный автомобиль марки Datsun сошел с конвейера в апреле 1935 года. Сегодня заводы компании, производящие автомобили, расположены в 16 странах, официальные продажи осуществляются в 190 странах мира. В компании работает более 180 000 человек. В глобальном модельном ряду компании – более 60 различных моделей.
46339. Разработка и реализация экономико-математической модели расчета стоимости автоперевозок на примере OOO «Heineken» 184.85 KB
  Особенно бурное развитие она получила в период второй мировой войны когда была применена для решения стратегических задач и четкого взаимодействия оборонной промышленности типовых и снабженческих баз и транспорта с целью своевременного обеспечения армии вооружением ГСМ и продовольствием. Подразделения логистики созданы на предприятиях промышленности аграрнопромышленного комплекса транспорта в аппарате НАТО они включаются в состав организационных комитетов по проведению крупных международных соревнований и т. К концу 20 века...
46340. Расчёт и конструирование перекрытия здания 1.51 MB
  Расчёт и конструирование сборной плиты перекрытия Расчётная схема: задаёмся размерами сечения ригеля: Сбор нагрузок: нормативные и расчётные нагрузки на 1 м2 перекрытия: Вид нагрузки нормативная Н м2 коэффициент надёжности по нагрузке расчётная нагрузка Н м2 Постоянная: собственный вес ребристой плиты; вес пола: ц. Высота сечения ребристой предварительно напряжённой плиты
46342. Модернизация бульдозерного оборудования 2.98 MB
  Земляные работы являются самым распространенным видом работ строительного производства. Они применяются в транспортном, гидротехническом, промышленном и городском строительстве, а также в сельском хозяйстве. Вместе с тем они принадлежат к весьма важным видам работ, от качества которых во многом зависит устойчивость и продолжительность службы устроенных на грунтовых основаниях таких инженерных сооружений, как автомобильные и железные дороги, плотины, дамбы.
46343. Обзор конструкций гидротрансформаторов 4.18 MB
  Для предохранения от самопроизвольного включения или выключения передачи в механизмах управления применяются стержневые или шариковые фиксаторы. Специальные клапаны позволяют не только плавно переключать ступени но и производить эти переключения с определенным сдвигом по фазе. В...
46345. Разработка привода цепного конвейера 986.97 KB
  Требуемая частота вращения вала электродвигателя где общее передаточное число кинематических пар изделия: где предварительное значение передаточного числа планетарного механизма u56 предварительное значение передаточного числа цепной передачи. Частоту вращения приводного вала вычисляем по формуле: об мин об мин Подбираем электродвигатель мощнстью 15 кВт и синхронной частотой вращения 750 об мин по табл. Определим чатоты вращения валов:...
46346. История автомобилестроения 207 KB
  В конце XIX столетия развитие автомобилестроения двигалось медленно и неповоротливо: с современным его роднил разве что двигатель внутреннего сгорания. Однако в 1914 году, с началом первой мировой войны, — спустя всего 29 лет с момента появления первой незатейливой модели Бенца