3855

Інтерполяційна схема Ейткена

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

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

Інтерполяційна схема Ейткена Мета роботи:Засвоїти теоретичний матеріалі методи апроксимації функцій, набути практичні навики знаходження наближених значень функцій. Короткі теоретичні відомості На практиці зустрічаються випадки, коли потрібно ...

Украинкский

2012-11-09

324 KB

15 чел.

Інтерполяційна схема Ейткена

Мета роботи: Засвоїти теоретичний матеріалі методи апроксимації функцій, набути практичні навики знаходження наближених значень функцій.

Короткі теоретичні відомості

На практиці зустрічаються випадки, коли потрібно мати значення інтерполяційного багаточлена Лагранжа в деякій точці х, а не загальний його вигляд. Тоді зручно користуватись інтерполяційною схемою Ейткена.

Обчислювальний алгоритм має такий вигляд:

Таблиця 1.

причому

  

– інтерполяційний багаточлен Лагранжа за вузлами хо,х,...,хп. Кожен із   отримують з  та   шляхом перехресного множення та ділення. Із застосуванням схеми Ейткена поступово можна залучати щораз нові значення вузлів хк доти, поки обчислення не засвідчать, що точність уже не зростає.

Варіант 17

Користуючись таблицею значень cos(x), xє[0.75k, 0.8k] з кроком h=0.01. Знайти наближене значення cos(x) при x=0.775, де k-списковий номер студента.

Результати обчислень подати таблично.

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

//---------------------------------------------------------------------------

#include <vcl.h>

#include <math.h>

#pragma hdrstop

#include "EytkinUn.h"

//---------------------------------------------------------------------------

#pragma package(smart_init)

#pragma resource "*.dfm"

TForm1 *Form1;

AnsiString str="L", str1;

float x0,xn,h,x,k,xk;

float masxkx[1000]={0}, masL[1000][1000]={0};

//---------------------------------------------------------------------------

__fastcall TForm1::TForm1(TComponent* Owner)

 : TForm(Owner)

{}

//---------------------------------------------------------------------------

void __fastcall TForm1::FormCreate(TObject *Sender)

{

 StringGrid1->Cells[0][0]="xk";

 StringGrid1->Cells[1][0]="fk";

 StringGrid1->Cells[2][0]="xk-x";

}

//---------------------------------------------------------------------------

void __fastcall TForm1::Button1Click(TObject *Sender)

{

 x0=StrToFloat(LabeledEdit1->Text);

 xn=StrToFloat(LabeledEdit2->Text);

 h=StrToFloat(LabeledEdit3->Text);

 x=StrToFloat(LabeledEdit4->Text);

 k=StrToFloat(LabeledEdit5->Text);

 StringGrid1->RowCount=2;

 x0*=k;   xn*=k;   x*=k;  xk=x0;   int i=0;    float min=fabs(xk-x);

 int t=0;

 while(xk<xn)

 {

   StringGrid1->Cells[0][i+1]=FloatToStrF(xk,ffFixed,8,2);

   masL[0][i]=cos((float)xk);

   StringGrid1->Cells[1][i+1]=FloatToStrF(masL[0][i],ffFixed,16,8);

   masxkx[i]=xk-x;

   if(min>=fabs(masxkx[i])){t=i;min=fabs(masxkx[i]);}

   StringGrid1->Cells[2][i+1]=FloatToStrF(masxkx[i],ffFixed,8,4);

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

   {

     masL[j][i]=1/(xk-x0)*(masL[j-1][i-1]*(xk-x)-masL[j-1][i]*(x0-x));

     StringGrid1->ColCount=2+i;

     StringGrid1->Cells[2+j][i+1]=FloatToStrF(masL[j][i],ffFixed,16,8);

   }

   xk+=h;

   StringGrid1->RowCount++;   i++;

 }

 Label1->Caption=(AnsiString)"sin("+

                 FloatToStrF(x,ffFixed,8,4)+

                 (AnsiString)") = "+

                 FloatToStrF(masL[t][t],ffFixed,16,8);

}

//---------------------------------------------------------------------------

Результат виконання програми:

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


 

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

33398. Структура РПП и ВПП КР1816ВЕ51 28.5 KB
  Организация памяти в микроконтроллере иллюстрируется рисунке Память программ имеет 16битовую адресную шину ее элементы адресуются с использованием счетчика команд PC или инструкций которые вырабатывают 16разрядные адреса. Память программ доступна только по чтению. ОМЭВМ не имеют команд и управляющих сигналов предназначенных для записи в память программ.
33399. Структура РПД и ВПД КР1816ВЕ51 27.5 KB
  Организация памяти в микроконтроллере иллюстрируется рисунке Память данных делится на внешнюю и внутреннюю каждая из них имеет свое пространство адресов. В архитектуре MК51 пространство адресов внутренней памяти данных объединяет все внутренние программно доступные ресурсы. Это пространство размером 256 байт в свою очередь делится на пространство адресов внутреннего ОЗУ резидентная память данных РПД размером 128 байт и пространство адресов регистров специальных функций.
33400. Порты ввода-вывода КР1816ВЕ51 34.5 KB
  Каждый порт содержит управляемые регистрзащелку входной буфер и выходной драйвер. Выходные драйверы портов 0 Р0 и 2 Р2 а также входной буфер Р0 используются при обращении к внешней памяти ВПП и ВПД.5 Вход таймера счетчика 1 или тест вход.4 Вход таймера счетчика 0 или тест вход.
33401. Таймер/Счетчики КР1816ВЕ51 30 KB
  Для управления режимами работы Т С и для организации взаимодействия таймеров с системой прерывания используются два регистра специальных функций РРТС и РУСТ. РРТС определяет включение и выключение T C а также их режимы работы. Используются четыре режима работы Т С. Режим 0.
33402. Образ организации-работодателя 15.43 KB
  Профиль восприятия или имидж предприятия как работодателя в персоналмаркетингеmix называют персоналимиджем personlimge. Персоналимидж организации управляет поведением кандидата на наем ищущего работу. Сначала персоналимидж влияет на решение кандидата: можно ли вообще рассматривать организацию как работодателя. Если организация в итоге потенциальной для найма определяется то персоналимидж свое решающее влияние может оказать на претендента в процессе собеседования в результате которого у кандидата на должность сформируется...
33403. КОНСТРУКТИВНЫМ КОНФЛИКТ 45 KB
  Чтобы конфликт начал разрастаться нужен инцидент т. По значению конфликты делятся на конструктивные созидательные и деструктивные неконструктивные разрушительные. КОНСТРУКТИВНЫМ КОНФЛИКТ бывает тогда когда оппоненты имея собственную позицию не выходят за рамки этических норм деловых отношений и разумных аргументов.
33404. Етикет подарка 35.5 KB
  Ваш подарок может и должен зависеть от характера Ваших отношений с тем кому Вы собираетесь его преподнести. Такой подарок разрушает имидж слабого пола . В портмоне или кошелек надо обязательно положить монетку иначе подарок сулит безденежье. Этот подарок может смутить и вызвать неприятные ощущения если будет подарен в присутствии смешанной компании.
33405. ПРИЁМ ПОСЕТИТЕЛЕЙ И ОБЩЕНИЕ С НИМИ 35.5 KB
  Хозяин офиса принимает посетителя сидя ЗА РАБОЧИМ СТОЛОМ. Он отвечает на приветствие и приглашает сесть перед своим столом. Расположение участников беседы в офисе за стандартным прямоугольным столом при четырех положениях вашего собеседника: 1. Люди не желающие взаимодействовать за столом друг с другом занимают независимую позицию.
33406. Психологическая культура делового разговора 12.94 KB
  Создание благоприятного психологического климата. Задавайте вопросы, на которые вашему партнеру приятно будет отвечать. Поощрите его к тому, чтобы он рассказывал и о достижениях своей фирмы.