Подякувати Студентському архіву довільною сумою
Admin
26.02.2023 12:38
Дякуємо, що користуєтесь нашим архівом!
Задача апроксимації.
Інформація про навчальний заклад
ВУЗ:
Національний університет Львівська політехніка
Інститут:
Не вказано
Факультет:
Не вказано
Кафедра:
Програмного забезпечення (ПЗ)
Інформація про роботу
Рік:
2008
Тип роботи:
Лабораторна робота
Предмет:
Чисельні методи в інформатиці
Група:
КН-3
Частина тексту файла (без зображень, графіків і формул):
МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
Національний університет Львівська політехніка
Кафедра ПЗ
До лабораторної роботи № 3
З курсу “Чисельні методи”
На тему:
“ Задача апроксимації”.
Виконав:
Студент групи КН-3
Львів – 2008
Тема:
Задача апроксимації..
Мета:
Оволодіти засобами розв’язування задач апрксимації.
Побудувати поліном Лагранжа для функції, заданої таблично та побудувати графік отриманого многочлена.
Варіанти завдань:
Таблиця 1
Номер точки
Значення функції згідно номера варіанта
1
2
3
4
5
6
7
2
5.06e-10
4.51e-10
5.87e-10
4.75e-10
4.57e-10
4.52e-10
3.51e-10
Таблиця 2
Номер точки
Значення функції згідно номера варіанта
8
9
10
11
12
13
14
2
2.10e-06
2.59e-06
2.11e-06
2.13e-06
2.19e-06
2.14e-06
1.13e-06
В процесі розроблення лабораторної роботи було розроблено наступний програмний продукт:
//---------------------------------------------------------------------------
#include <vcl.h>
#pragma hdrstop
#include "Main.h"
#include <math.h>
#include <string.h>
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TMainForm *MainForm;
//---------------------------------------------------------------------------
__fastcall TMainForm::TMainForm(TComponent* Owner)
: TForm(Owner)
{
}
//---------------------------------------------------------------------------
double *TMainForm::LagrangeGetPolinom()
{
int i, j, k, h;
double *Result = new double[COUNT];
double *Temp = new double[COUNT - 1];
double ResMul, ResAdd, Denom;
int *Counter;
for(i = 0; i < COUNT; i++)
{
Result[i] = 0.0;
}
for(h = 0; h < COUNT; h++)
{
Denom = 1.0;
for(i = 0, j = 0; i < COUNT; i++)
{
if(i != h)
{
Temp[j] = -DATA_X[i];
Denom *= DATA_X[h] - DATA_X[i];
j++;
}
}
for(i = 0; i < COUNT; i++)
{
ResAdd = 0.0;
Counter = new int[i];
for(j = 0; j < i; j++)
{
Counter[j] = j;
}
for(j = 0; j < Factorial(COUNT - 1) / (Factorial(i) * Factorial(COUNT - 1 - i)); j++)
{
ResMul = 1.0;
for(k = 0; k < i; k++)
{
ResMul *= Temp[Counter[k]];
}
if(i != 0)
{
IncCounter(Counter, i, COUNT - 1, i - 1);
}
ResAdd += ResMul;
}
if(Denom != 0)
{
Result[i] += ResAdd / Denom * DATA_Y[h];
}
delete[]Counter;
}
}
return Result;
}
void TMainForm::LagrangeDrawGraph(const double *Data)
{
int i;
double min, max, d;
min = max = DATA_X[0];
for(i = 1; i < COUNT; i++)
{
if(DATA_X[i] < min)
{
min = DATA_X[i];
}
if(DATA_X[i] > max)
{
max = DATA_X[i];
}
}
for(d = min; d <= max; d+=(max-min) / 100)
{
Chart->Series[0]->AddXY(d, Polynom(d, Data), "", TColor(5));
}
}
__int64 Factorial(__int64 x)
{
if(x < 2)
{
return 1;
}
return x * Factorial(x - 1);
}
int IncCounter(int *Counter, int Size, int Capacity, int Number)
{
Counter[Number]++;
if(Counter[Number] > Capacity - Size + Number)
{
if(Number != 0)
{
Counter[Number] = IncCounter(Counter, Size, Capacity, Number - 1);
Counter[Number]++;
}
else
{
Counter[Number] = Capacity - Size + Number;
}
}
return Counter[Number];
}
double Polynom(double x, const double *Data)
{
int i;
double result = 0.0;
for(i = 0; i < COUNT; i++)
{
result += pow(x, COUNT - 1 - i) * Data[i];
}
return result;
}
void __fastcall TMainForm::FormCreate(TObject *Sender)
{
double *Data = LagrangeGetPolinom();
AnsiString str;
LagrangeDrawGraph(Data);
for(int i = 0; i < COUNT; i++)
{
str += " ";
if(Data[i] > 0 && i != 0)
{
str += "+";
}
str += FloatToStrF(Data[i], ffFixed, 3, 3) + "x^" + (COUNT - 1 - i);
}
str += " =y";
RichEdit->Clear();
RichEdit->Text = str;
delete[]Data;
}
//---------------------------------------------------------------------------
Висновок:
На цих лабораторній роботі я оволодів засобами розв’язування задач апроксимації.
Ділись своїми роботами та отримуй миттєві бонуси!
0%
Оголошення від адміністратора