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.  Контрольні питання

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


 

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

25912. Автоматические выключатели. Основные типы. Принцип действия. Основные типы расцепителей. Механизм свободного расцепления. Условия выбора 42 KB
  Конструктивная схема автомата На рисунке дана условная конструктивная схема универсального автомата в упрощенном изображении. Создается усилие которое перемещая рычаги 4 и 5 вправо будет поворачивать основную несущую деталь 6 автомата вокруг неподвижной оси О по часовой стрелке. Замыкаются и включают цепь тока вначале дугогасительные 8 и 10 а затем главные 7 и 11 контакты автомата. Отключающая пружина 2 взводится при включении автомата.
25913. Выключатели переменного тока высокого напряжения. Назначение, основные требования. Номинальный ток отключения. Классификация высоковольтных выключателей по виду дугогасительной среды и изоляции межконтактного промежутка. Использование АПВ. Условия выбора 45.5 KB
  Выключатели переменного тока высокого напряжения. Параметры В соответствии с ГОСТ Р 525652006 выключатели характеризуются следующими параметрами: номинальное напряжение Uном напряжение сети в которой работает выключатель; номинальный ток Iном ток через включённый выключатель при котором он может работать длительное время; номинальный ток отключения Iо.ном наибольший ток короткого замыкания действующее значение который выключатель способен отключить при напряжении равном наибольшему рабочему напряжению при заданных условиях...
25914. Конструкция и принцип действия малообъемных масленых выключателей. Достоинства и недостатки. Условия выбора. Сравнение с другими высоковольтными выключателями 329 KB
  Масляные выключатели В дугогасительных устройствах масляных выключателей гашение дуги осуществляется путем эффективного ее охлаждения в потоке газопаровой смеси вырабатываемой дугой в результате разложения и испарения масла. В зависимости от назначения масла можно выделить 2 основные группы масляных выключателей: 1. В состав газопаровой смеси возникающей в результате разложения масла под действем дуги входит до 70 водорода обладающего по сравнению с воздухом в 8 раз более высокой теплопроводностью но меньшей предельной электрической...
25915. Конструкция и принцип действия вакуумных выключателей. Достоинства и недостатки. Условия выбора. Сравнение с другими высоковольтными выключателями 22.5 KB
  Таким образом дуга в вакууме существует изза ионизации паров контактного материала вначале за счет материала контактного мостика а затем в результате испарения материала электродов под воздействием энергии дуги. Поэтому если поступление паров контактного материала будет недостаточно вакуумная дуга должна погаснуть.
25916. Тепловые процессы в электрических аппаратах. Источники теплоты. Поверхностный эффект и эффект близости. Способы распространения теплоты в пространстве. Термическая стойкость электрических аппаратов 292 KB
  ТЕРМИЧЕСКАЯ СТОЙКОСТЬ ЭЛЕКТРИЧЕСКИХ АППАРАТОВ НАГРЕВ ПРИ КЗ Под термической стойкостью аппарата понимают его способность выдерживать без повреждений и перегрева свыше норм термическое действие токов короткого замыкания определенной длительности. Термическая стойкость характеризуется током термической стойкости т. Для выключающих аппаратов термическая стойкость задается обычно 10секундным током и...
25917. Контактные явления в электрических аппаратах. Классификация контактов их конструкция и материал исполнения. Понятия переходного сопротивления контакта 49 KB
  Классификация контактов их конструкция и материал исполнения. Чем больше контактов в цепи тем сильнее сопротивление. При точечном контакте контактные нажатия небольшие и для уменьшения сопротивления контактов применяют драгоценные металлы не образующие окиси. Для этих контактов применяют медь.
25919. Электромеханические реле. Принцип действия. Виды электромеханических реле, их назначение. Основные характеристики, требования 25.5 KB
  Электромеханические реле. Виды электромеханических реле их назначение. Электромеханическое релекоммутационное устройство предназначенное производить скачкообразные изменения в управляющих цепях. реле подразделяются на 2 класса: электромеханические статические Эл.
25920. Электромеханические реле времени. Тепловые реле. Принцип работы. Область применения 24 KB
  Электромеханические реле времени. Тепловые реле. реле времени. Схема защиты реле автоматикичасто требуется выдержка времени когда выдержка устанавливается для предотвращения срабатывания защиты от пусковых токов.