17081

Формули Н’ютона через кінцеві різниці

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

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

Лабороторна робота №6 Тема. Формули Нютона через кінцеві різниці Мета. Навчитися обчислити значення функції при даному значенні аргумента використовуючи формули Нютона через кінцеві різниці. Обладнання. Лист формату А4 ручка олівець програмне забезпечення С...

Украинкский

2013-06-29

50 KB

5 чел.

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

Тема. Формули Н’ютона через кінцеві різниці

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

Обладнання. Лист формату А4, ручка, олівець, програмне забезпечення С++.

Хід роботи

  1.  Правила ТБ
  2.  Теоретичні відомості

Друга інтерполяційна формула Ньютона(використовується, коли х знаходиться в кінці таблиці):

     3. Індивідуальна робота

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

x

y

0,101

1,26183

0,106

1,27644

0,111

1,29122

0,116

1,30617

0,121

1,32130

0,126

1,33660

0,131

1,35207

0,136

1,36773

0,141

1,38357

0,146

1,39959

0,151

1,41579

№ варіанта

х1

х2

2

0,1035

0,1492

#include<iostream.h>

#include<math.h>

double f1(double xn[],double yn[],double Xn,int nt)

{int i,j;

double q,h,s,qn,raz,fakt;

s=yn[0];

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

{h=xn[i+1]-xn[i];

q=(Xn-xn[0])/h;

qn=1.0;

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

{qn=qn*(q-j+1);

}

fakt=1.0;

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

{fakt=fakt*j;

}

raz=pow(yn[1]-yn[0],i);

s=s+raz*qn/fakt;

}

return s;}

double f2(double xn[],double yn[],double Xn,int nt)

{int i,j;

double q,h,s,qn,raz,fakt;

s=yn[nt];

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

{h=xn[i+1]-xn[i];

q=(Xn-xn[0])/h;

qn=1.0;

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

{qn=qn*(q-j+1);

}

fakt=1.0;

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

{fakt=fakt*j;

}

raz=pow(yn[nt-i+1]-yn[nt-i],i);

s=s+raz*qn/fakt;

}

return s;}

void main()

{double X,a,b,Pnx;

int i,n;

cout<<"Введите X\n";

cin>>X;

cout<<"Введите n\n";

cin>>n;

double*x=new double[n];

double*y=new double[n];

cout<<"Введите массив x\n";

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

cin>>x[i];

cout<<"Введите массив y\n";

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

cin>>y[i];

a=fabs(X-x[0]);

b=fabs(X-x[n]);

if(a<b)Pnx=f1(x,y,X,n);

else Pnx=f2(x,y,X,n);

cout<<"Pn(x)="<<Pnx<<"\n";}

  1.  Контрольні питання:
  2.  Дати визначення кінцевої різниці 1-го, к-го порядку
  3.  Поставте задачу інтерполяції функції
  4.  Запишіть 1-у інтерполяційну формулу Н’ютона. В якому випадку вона застосовується?
  5.  Запишіть 2-у інтерполяційну формулу Н’ютона. В якому випадку вона застосовується?

5. Звіт.        

6. Захист роботи  


 

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

6792. Изучение методов вывода иллюстраций, создания битовых образов, создания движущихся и изменяющихся рисунков 70 KB
  Цель работы: Изучение методов вывода иллюстраций, создания битовых образов, создания движущихся и изменяющихся рисунков. Приобретение навыков создания программ  мультипликации. Постановка задачи:Изобразить на экране прямоугольник, враща...
6793. ВЫПРЯМИТЕЛЬ СО СГЛАЖИВАЮЩИМИ ФИЛЬТРАМИ 160.5 KB
  ВЫПРЯМИТЕЛЬ СО СГЛАЖИВАЮЩИМИ ФИЛЬТРАМИ Цель работы - изучение принципа действия и экспериментальные исследования однофазного маломощного выпрямителя с простейшими сглаживающими фильтрами. 1.1 Основные положения Для питания постоянным током разл...
6794. УПРАВЛЯЕМЫЙ ВЫПРЯМИТЕЛЬ 139 KB
  УПРАВЛЯЕМЫЙ ВЫПРЯМИТЕЛЬ Цель работы - ознакомление с принципом работы и экспериментальное исследование характеристик однофазного управляемого выпрямителя на тиристорах при различных видах нагрузки. 2.1 Основные положения Для регулирования выпря...
6795. ПОЛУПРОВОДНИКОВЫЕ СТАБИЛИЗАТОРЫ 126 KB
  ПОЛУПРОВОДНИКОВЫЕ СТАБИЛИЗАТОРЫ Цель работы - исследование рабочего участка вольт-амперной характеристики полупроводникового стабилитрона, изучение основных характеристик параметрического и компенсационного стабилизаторов напряжения и определен...
6796. Усилители постоянного тока 121.5 KB
  Усилители постоянного тока Цель работы - изучение принципа работы простейшей схемы дифференциального усилителя и экспериментальное определение его параметров ознакомление с работой операционного усилителя (ОУ) и экспериментальное определение е...
6797. Усилительные каскады на транзисторах 135.5 KB
  Усилительные каскады на транзисторах Цель работы - исследование одиночных усилительных каскадов RC-типа на биполярном и полевом транзисторах исследование характеристик и определение параметров изучаемых усилителей. 5.1 Основные положения Усилители ...
6798. Аналоговые схемы на операционных усилителях 146.5 KB
  Аналоговые схемы на операционных усилителях Цель работы - изучение некоторых схем включения операционного усилителя для обработки аналоговых сигналов определение характеристик и параметров инвертирующего и неинвертирующего усилителей, сумматор...
6799. Импульсные схемы на операционных усилителях 123.5 KB
  Импульсные схемы на операционных усилителях Цель работы - изучение принципа работы компаратора и триггера Шмитта на операционном усилителе исследование и определение параметров схем мультивибратора, одновибратора и генератора треугольных импульсов...
6800. Генераторы импульсов на транзисторах и интегральных микросхемах 98 KB
  Генераторы импульсов на транзисторах и интегральных микросхемах Цель работы - ознакомление с принципом работы мультивибратора на биполярных транзисторах и определение его параметров изучение принципа работы мультивибратора, одновибратора и генерато...