86294 (612714), страница 3
Текст из файла (страница 3)
Для роботи програми треба виконати наступні дії:
-
у поля "2" "3" потрібно ввести межі інтервалу на якому функція визначена і диференційовна;
-
у поле "1" треба ввести похибку;
-
вибрати вид сплайна 4;
-
вибрати функцію 5;
-
далі натиснути кнопку "6", яка викликає функцію, що будує балансне наближення із заданою похибкою;
-
для побудови балансного наближення із заданою кількістю ланок треба натиснути кнопку "8";
-
вивід результатів буде у полі "7";
-
щоб переглянути графіки функції і сплайну потрібно натиснути на закладку 9;
-
щоб переглянути графік похибки наближення потрібно натиснути на закладку 10;
Опис основних функцій програми:
-
void ermit_1(double *a, double x0, double x1); - функція пошуку коефіцієнтів ермітового сплайна з ланками виду
![]()
; -
void ermit_2(double *a, double x0, double x1); - функція пошуку коефіцієнтів ермітового сплайна з ланками виду
![]()
; -
void ermit_3(double *a, double x0, double x2); - функція пошуку коефіцієнтів ермітового сплайна з ланками виду
![]()
; -
void ermit_4(double *b, double x0, double x2); - функція пошуку коефіцієнтів ермітового сплайна з ланками виду
![]()
; -
void ermit_5(double *a, double zl, double zp); - функція пошуку коефіцієнтів ермітового сплайна з ланками виду
![]()
; -
void ermit_rp(double *a, double a1, double b1, double nyu, int n); - рівномірне наближення ермітовими сплайнами із заданною похибкою;
-
void ermit_rl (double*a, double zl, double b1, int n, int p, int r); - рівномірне наближення ермітовими сплайнами із заданною кількістю ланок.
;
;
;
;
;Література
-
Пізюр Я.В., Попов Б.О. Рівномірне наближення ермітовими сплайнами з парною кількістю параметрів.// Контрольно-вимірювальна техніка.- 1993. – Вип. 50. – С. 8-13
-
Пізюр Я.В. Наближення функцій ермітовими сплайнами з експоненціальними ланками// Вісник НУ "Львівська політехніка". "Фізико-математичні науки" №566, 2006, – С. 68-75.
-
Зав’ялов Ю.С., Квасов Б.И., Мірошниченко В.Л. Методи сплайн функцій. – М.: Наука, 1980. – 352 с.
Додаток
Код програми.
private: System::Void button2_Click(System::Object^ sender, System::EventArgs^ e) {
double zl, zp, x;
zl=Double::Parse(textBox3->Text);
zp=Double::Parse(textBox4->Text);
x =(zl+zp)/2;
int p=2,n=5;
if((s1->Checked==true)||(s2->Checked==true)||(s5->Checked==true)||(s7->Checked==true)||(s9->Checked==true)){
n=4;p=1;
}
double *a = new double[n];
if (s1->Checked == true) {
ermit_1(a, zl, zp, p);
}
if (s2->Checked == true) {
ermit_2(a, zl, zp, p);
}
if (s3->Checked == true) {
ermit_3(a, zl, zp, p);
}
if (s4->Checked == true) {
ermit_4(a, zl, zp, p);
}
if (s5->Checked == true) {
ermit_5(a, zl, zp, p);
}
if (s6->Checked == true) {
ermit_6(a, zl, zp, p);
}
if (s7->Checked == true) {
ermit_8(a, zl, zp, p);
}
if (s8->Checked == true) {
ermit_9(a, zl, zp, p);
}
if (s9->Checked == true) {
ermit_7(a, zl, zp, p);
}
chart1->Series["Function"]->Points->Clear();
chart1->Series["Spline"]->Points->Clear();
chart2->Series["error"]->Points->Clear();
double t,ch1,ch2,ch3;
for(t = zl; t <= zp; t += 0.01)
{
ch1 = f(t,1);
ch2 = sp(a,t);
ch3 = epx_p(a,t,1);
chart1->Series["Function"]->Points->AddXY(t, ch1);
chart1->Series["Spline"]->Points->AddXY(t, ch2);
chart2->Series["error"]->Points->AddXY(t, ch3);
}
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
double epx_p(double *a, double x, int p) {
return Math::Abs((f(x, p) - sp(a, x)) / w(x));
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
public: double f(double x, int p) {
if(f1->Checked==true){
return Math::Sin(x);
}
if(f2->Checked==true){
return Math::Cos(x);
}
if(f3->Checked==true){
return 1 / (2 + x * x);
}
if(f4->Checked==true){
return Math::Log(x + 1);
}
if(f5->Checked==true){
return Math::Exp(x);
}
if(f6->Checked==true){
return x*x;
}
return -1;
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
public: double fp(double x, int p) {
if(f1->Checked==true){
return Math::Cos(x);
}
if(f2->Checked==true){
return -Math::Sin(x);
}
if(f3->Checked==true){
return (-2* x ) / ((2 + x * x) * (2 + x * x));
}
if(f4->Checked==true){
return 1 / (x + 1);
}
if(f5->Checked==true){
return Math::Exp(x);
}
if(f6->Checked==true){
return 2*x;
}
return -1;
}
//*****************************************
public: void ermit_1(double *a, double x0, double x1, int p) {
a[3] = (fp(x1, p) / f(x1, p) + fp(x0, p) / f(x0, p) + 2 * (Math::Log(f(x0, p)
/ f(x1, p)) / (x1 - x0))) / ((x1 - x0) * (x1 - x0));
a[2] = ((2* x0 * x0 - x0 * x1 - x1 * x1) * a[3] - (fp(x0, p) / f(x0, p))
- (Math::Log(f(x0, p) / f(x1, p)) / (x1 - x0))) / (x1 - x0);
a[1] = (1 / (x1 - x0)) * (Math::Log(f(x1, p) / f(x0, p)) - a[3] * (x1 * x1 * x1
- x0 * x0 * x0) - a[2] * (x1 * x1 - x0 * x0));
a[0] = f(x0, p) * Math::Exp(-(a[1] * x0 + a[2] * x0 * x0 + a[3] * x0 * x0 * x0));
}
/////////////////////////////////////////////////////////////////////////////
void ermit_2(double *b, double x0, double x1, int p) {
double h;
h = x1 - x0;
b[3] = 2 * ((fp(x1, p) + fp(x0, p)) / 2 - (f(x1, p) - f(x0, p)) / h) / (h
* h);
b[2] = 0.5 * ((fp(x1, p) - fp(x0, p)) / h - 3* b [3] * (x0 + x1));
b[1] = fp(x1, p) - 2* b [2] * x1 - 3* b [3] * x1 * x1;
b[0] = 0.5 * (f(x1, p) + f(x0, p) - b[3] * (x0 * x0 * x0 + x1 * x1 * x1)
- b[2] * (x0 * x0 + x1 * x1) - b[1] * (x0 + x1));
}
////////////////////////////////////////////////////////////////
void ermit_3(double *a, double x0, double x2, int p) {
double a1, a2, a3, b1, b2, b3, z1, z2, z3, g1, g2, g3;
double x1;
x1 = (x0 + x2) / 2;
a1 = Math::Log(f(x1, p) / f(x0, p)) / (x1 - x0) - Math::Log(f(x2, p) / f(x0, p)) / (x2 - x0);
b1 = (x1 + x0) * (x1 * x1 + x0 * x0) - (x2 + x0) * (x2 * x2 + x0 * x0);
z1 = x1 * x1 + x1 * x0 - x2 * x2 - x2 * x0;
g1 = x1 - x2;
a2 = Math::Log(f(x1, p) / f(x0, p)) / (x1 - x0) - fp(x0, p) / f(x0, p);
b2 = x1 * x1 * x1 + x0 * x1 * (x0 + x1) - 3* x0 * x0 * x0;
z2 = x1 * x1 + x0 * x1 - 2* x0 * x0;
g2 = x1 - x0;
a3 = Math::Log(f(x1, p) / f(x0, p)) / (x1 - x0) - fp(x2, p) / f(x2, p);
b3 = (x0 + x1) * (x0 * x0 + x1 * x1) - 4* x2 * x2 * x2;
z3 = x1 * x1 + x0 * x1 + x0 * x0 - 3* x2 * x2;
g3 = x1 + x0 - 2* x2 ;
a[2] = ((a2 * b1 - a1 * b2) * (z3 * b1 - z1 * b3) - (z2 * b1 - z1 * b2)
* (a3 * b1 - a1 * b3)) / ((g1 * b2 - g2 * b1) * (z1 * b3 - z3 * b1)
+ (z2 * b1 - z1 * b2) * (g1 * b3 - g3 * b1));
a[3] = (a2 * b1 - a1 * b2 + (g1 * b2 - b1 * g2) * a[2]) / (b1 * z2 - z1 * b2);
a[4] = (a1 - z1 * a[3] - g1 * a[2]) / b1;
a[1] = (1 / (x2 - x0)) * (Math::Log(f(x2, p) / f(x0, p)) - a[4] * (x2 * x2 * x2
* x2 - x0 * x0 * x0 * x0) - a[3] * (x2 * x2 * x2 - x0 * x0 * x0)
- a[2] * (x2 * x2 - x0 * x0));
a[0] = f(x0, p) * Math::Exp(-(a[1] * x0 + a[2] * x0 * x0 + a[3] * x0 * x0 * x0 + a[4] * x0 * x0 * x0 * x0));
}
///////////////////////////////////////////////////////////////////////////
void ermit_4(double *a, double x0, double x2, int p) {
double a1, a2, a3, b1, b2, b3, z1, z2, z3, g1, g2, g3;
double x1, h;
x1 = (x0 + x2) / 2;
h = x1 - x0;
a1 = (f(x2, p) - f(x0, p)) / (x2 - x0) - (f(x1, p) - f(x0, p)) / (x1 - x0);
b1 = (x2 + x0) * (x2 * x2 + x0 * x0) - (x1 + x0) * (x1 * x1 + x0 * x0);
z1 = x2 * x2 + x2 * x0 - x1 * x1 - x1 * x0;
g1 = x2 - x1;
a2 = (f(x1, p) - f(x0, p)) / (x1 - x0) - fp(x0, p);
b2 = x1 * x1 * x1 + x0 * x1 * (x0 + x1) - 3* x0 * x0 * x0;
z2 = x1 * x1 + x0 * x1 - 2* x0 * x0;
g2 = x1 - x0;
a3 = (f(x1, p) - f(x0, p)) / (x1 - x0) - fp(x2, p);
b3 = (x0 + x1) * (x0 * x0 + x1 * x1) - 4* x2 * x2 * x2;
z3 = (x1 * x1 + x0 * x1 + x0 * x0) - 3* x2 * x2;
g3 = (x1 + x0 - 2* x2 );
a[2] = ((a2 * b1 - a1 * b2) * (z3 * b1 - z1 * b3) - (z2 * b1 - z1 * b2)
* (a3 * b1 - a1 * b3)) / ((g1 * b2 - g2 * b1) * (z1 * b3 - z3 * b1)
+ (z2 * b1 - z1 * b2) * (g1 * b3 - g3 * b1));
a[3] = (a2 * b1 - a1 * b2 + (g1 * b2 - b1 * g2) * a[2]) / (b1 * z2 - z1
* b2);
a[4] = (a1 - z1 * a[3] - g1 * a[2]) / b1;
a[1] = (f(x1, p) - f(x0, p)) / h - a[2] * (x1 + x0) - a[3] * (x1 * x1 + x1
* x0 + x0 * x0) - a[4] * (x1 + x0) * (x1 * x1 + x0 * x0);
a[0] = f(x0, p) - a[1] * x0 - a[2] * x0 * x0 - a[3] * x0 * x0 * x0 - a[4]
* x0 * x0 * x0 * x0;
}
////////////////////////////////////////////////////////////////
void ermit_5(double *a, double zl, double zp, int p) {
double h, x1;
h = (zp - zl) / 2;
x1 = (zp + zl) / 2;
a[3] = ((h * h * (fp(zp, p) - fp(zl, p))) / (f(zp, p) - f(zl, p) - h * (fp(
zp, p) + fp(zl, p)))) - x1;
a[2] = ((fp(zp, p) - fp(zl, p)) * (((x1 + a[3]) * (x1 + a[3]) - h * h)
* ((x1 + a[3]) * (x1 + a[3]) - h * h))) / (4* h * (x1 + a[3]));
a[0] = fp(zl, p) + a[2] / ((x1 + a[3] - h) * (x1 + a[3] - h));
a[1] = f(zl, p) - a[0] * (x1 - h) - a[2] / (x1 - h + a[3]);
}
////////////////////////////////////////////////////////////////////
void ermit_8(double *a, double zl, double zp, int p) {
double j0, j1, j2, m0, m1;
j0 = (Math::Exp(f(zp, p)) - Math::Exp(f(zl, p))) / (zp - zl);
j1 = (zl * zl - zp * zp) / (zp - zl);
j2 = (zl * zl * zl - zp * zp * zp) / (zp - zl);
m0 = (f(zl, p) * Math::Exp(f(zl, p)) - j0) / (j1 - 2* zl );
m1 = (2* zl * zl - j2) / (j1 - 2* zl );
a[3] = (-fp(zp, p) * Math::Exp(f(zp, p)) + j0 + m0 * j1 - 2* zp * m0) / (3* zp
* zp + 2* zp * m1 - j2 - m1 * j1);
a[2] = m0 + a[3] * m1;
a[1] = j0 + j1 * a[2] + a[3] * j2;
a[0] = Math::Exp(f(zl, p)) - a[1] * zl - a[2] * zl * zl - a[3] * zl * zl * zl;
}
///////////////////////////////////////////////////////////////////
void ermit_9(double *a, double zl, double zp, int p) {
double j1, j2, j3, j4, m1, m2, m3, k1, k2, z;
z = (zp + zl) / 2;
j1 = (zp * zp - z * z) / (z - zp);
j2 = (zp * zp * zp - z * z * z) / (z - zp);
j3 = (zp * zp * zp * zp - z * z * z * z) / (z - zp);
j4 = (Math::Exp(f(z, p)) - Math::Exp(f(zp, p))) / (z - zp);
m1 = (Math::Exp(f(zp, p)) - Math::Exp(f(zl, p)) + j4 * (zl - zp)) / (zp * zp - zl * zl + j1 * (zp - zl));
m2 = (zl * zl * zl - zp * zp * zp + j2 * (zl - zp)) / (zp * zp - zl * zl
+ j1 * (zp - zl));
m3 = (zl * zl * zl * zl - zp * zp * zp * zp + j3 * (zl - zp)) / (zp * zp
- zl * zl + j1 * (zp - zl));
k1 = (Math::Exp(f(zl, p)) * fp(zl, p) - j4 - m1 * j1 - 2* zl * m1) / (j2 + 3* zl
* zl + j1 * m2 + 2* m2 * zl);
k2 = (j3 + 4* zl * zl * zl + j1 * m3 + 2* m3 * zl) / (j2 + 3* zl * zl + j1
* m2 + 2* m2 * zl);
a[4] = (Math::Exp(f(zp, p)) * fp(zp, p) - j4 - m1 - j1 - 2* zp * m1 - k1 * j2
- 3* zp * zp * k1 - k1 * j1 * m2 - 2* m2 * zp * k1) / (j3 + 4* zp
* zp * zp + j1 * m3 + 2* m3 * zp - k2 * j2 - 3* zp * zp * k2 - k2
* j1 * m2 - 2* m2 * zp * k2);
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
