![]() |
Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | |
|
Составить алгоритм, по которому написать и отладить программу аппроксимации функции f (x) = x 3 – 5 x 2 на интервале [-2, 5] многочленом Лагранжа, m – количество точек, в которых известна функция, n – количество рассчитываемых значений.
Вид формы и полученные результаты представлены на рис. 7.2. Тексты функций-обработчиков и функции пользователя будут иметь следующий вид:
double fun(double);
double Mn_Lagr(double*, double, int);
//------------------ Текст функции-обработчика кнопки Вычислить -------------------
double x,h,h1, a, b, *mas_x, *mas_y_t;
int i,n,m;
a = StrToFloat(Edit1->Text); b = StrToFloat(Edit2->Text);
m = StrToInt(Edit3->Text); n = StrToInt(Edit4->Text);
h = (b-a)/(m-1); h1 = (b-a)/(n-1);
mas_x = new double[m+1]; mas_y_t = new double[n+1];
for(x=a, i=0; i<m; i++){
mas_x[i] = x;
x+=h;
}
Memo1->Lines->Add("---- Многочлен Лагранжа ---");
Memo1->Lines->Add("Получили " + IntToStr(n) + " значений:");
for(x=a, i=0; i<n; i++, x+=h1) {
mas_y_t[i] = Mn_Lagr(mas_x,x,m);
Memo1->Lines->Add(" x = "+FloatToStrF(x,ffFixed,8,2)
+" f*(x) = "+FloatToStrF(mas_y_t[i],ffFixed,8,4));
}
//----------------------------- Очистка Графиков -----------------------------------------------
Chart1->Series[0]->Clear(); Chart1->Series[1]->Clear();
Рис. 7.2
//----------------------------- Вывод Графиков -------------------------------------------------
for(x=a-0.1; x<b+0.1; x+=0.01)
Chart1->Series[0]->AddXY(x,fun(x));
for(x=a,i=0; i<n; i++,x+=h1)
Chart1->Series[1]->AddXY(x,mas_y_t[i]);
delete []mas_x;
delete []mas_y_t;
}
//------------------------------ Исходная функция f (x) -----------------------------------------
double fun(double x) {
return pow(x,3) - 5 * x*x;
}
//----------------------------- Многочлен Лагранжа -------------------------------------------
double Mn_Lagr(double *x, double xt, int kol) {
int i, k;
double e, p=0;
for(k=0; k<kol; k++) {
e=1.;
for (i=0;i<kol;i++)
if (i!=k) e *= ((xt-x[i])/(x[k]-x[i]));
p += e*fun(x[k]);
}
return p;
}
Дата публикования: 2015-02-22; Прочитано: 247 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!