Лабораторный практикум по современным компьютерным технологиям. Часть 3. MathCAD. Ершова Е.Е - 29 стр.

UptoLike

ВУЗ: 

НГАСУ «Сибстрин» | Новосибирск

Составители: 

Рубрика: 

g
i 0,
f0 x
i
()
:= g
i 1,
f1 x
i
(
)
:= g
i 2,
f2 x
i
(
)
:=
g
i 3,
f3 x
i
()
:= g
i 4,
f4 x
i
()
:=
Сформируем матрицу системы уравнений для нахождения коэф-
фициентов пробной функции. Для этого зададим параметры:
i1 0 n..:= j1 0 n..:=
Формируем матрицу по формуле:
M
i1 j1,
0
N
k
g
ki1,
g
kj1,
⋅
()
∑
=
:=
Далее строим вектор правых частей для данной системы уравне-
ний:
f
i1
0
N
k2
y
k2
g
k2 i1,
⋅
()
∑
=
:=
Решаем нашу систему уравнений
M*a=f (см. лабораторную рабо-
ту № 3).
После чего получаем требуемую аппроксимацию:
PX() a
0
f0 X()⋅ a
1
f1 X()⋅+ a
2
f2 X()⋅+ a
3
f3 X()⋅+ a
4
f4 X()⋅+:=
График аппроксимации
P(X) и экспериментальные точки
представлены на рисунке: 
             gi, 0 := f0( xi) gi, 1 := f1( xi)            gi, 2 := f2( xi)

             gi, 3 := f3( xi)          gi, 4 := f4( xi)
     Сформируем матрицу системы уравнений для нахождения коэф-
 фициентов пробной функции. Для этого зададим параметры:

                           i1 := 0 .. n j1 := 0 .. n
        Формируем матрицу по формуле:
                                        N
                      Mi1, j1 :=       ∑       ( gk, i1 ⋅ gk, j1)
                                   k     =0
        Далее строим вектор правых частей для данной системы уравне-
 ний:
                                       N
                        f i1 :=    ∑          ( yk2 ⋅ gk2, i1)
                                  k2   =0
     Решаем нашу систему уравнений M*a=f (см. лабораторную рабо-
 ту № 3).
     После чего получаем требуемую аппроксимацию:

P ( X) := a0 ⋅ f0( X) + a1 ⋅ f1( X) + a2 ⋅ f2( X) + a3 ⋅ f3( X) + a4 ⋅ f4( X)
        График аппроксимации P(X) и экспериментальные точки
 представлены на рисунке:

Страницы