3. Визначаємо значення швидкості за геофізичними моделями за формулою (1).

ω(х)	-9,81E-10
ω(y)	-2,395E-09
ω(z)	3,153E-09

	GUIP	CREI	ENTZ	TSLE	VFCH
Vxм	-0,010348	-0,012096	-0,013168	-0,010851	-0,011581
Vyм	0,017985	0,017855	0,017887	0,018882	0,018234
Vzм	0,010442	0,009799	0,009490	0,010966	0,010247

	LIL2	AUTN	CHIZ	PUYV
Vxм	-0,01245	-0,01213776	-0,01085928	-0,01171955
Vyм	0,017581	0,018264256	0,018448918	0,018608071
Vzм	0,009481	0,010096972	0,010635017	0,010488249

4.Виконуємо інтерполяцію виходячи з формули (2).

Попередньо обчислюємо значення середньої широти B0 та різниці широт Bi-B0 і довготLi-L0.

Для станцій класу А

В0	47,43115	0,827829723
L0	3,545485	0,061880379
Bi-B0
1,0134535	0,017688
1,83179239	0,031971
1,11820906	0,019516
-3,82649733	-0,06678
-0,13695761	-0,00239

За формулою (2) будуємо матриці :

		а1
Матриця Х	=	а2
		а3
		а4

		Lx	Ly	Lz
		-0,001252	-0,000185	0,000158
		-0,000904	-0,000355	0,000801
Матриця L	=	-0,001532	-0,000387	0,000410
		-0,000449	0,000418	0,001234
		-0,000719	-0,000134	0,000653

		1,000000	0,017688	0,010030	0,000177
		1,000000	0,031971	-0,011773	-0,000376
Матриця А	=	1,000000	0,019516	0,047305	0,000923
		1,000000	-0,066785	-0,026556	0,001774
		1,000000	-0,002390	-0,021613	0,000052

Транспонуємо матрицю А :

		1	1	1	1	1
Матриця AT	=	0,0176881	0,031970808	0,01951643	-0,066784977	-0,00239
		0,010030022	-0,011773196	0,047305387	-0,026556355	-0,02161
		0,000177412	-0,000376399	0,000923232	0,001773566	5,17E-05

Обчислюємо матрицю N=ATA :

		5	0	-0,002607081	0,002549474
Матриця N	=	0	0,00618184	0,002549474	-0,000109449
		-0,002607081	0,002549474	0,003649368	1,6687E-06
		0,002549474	-0,000109449	1,6687E-06	4,17371E-06

Обчислюємо матрицю В, для кожної з координат B=ATL(Lx,Ly,Lz):

		-0,00485721	-0,00485721	-0,00485721	-0,00485721
Матриця Bx	=	-4,92205E-05	-4,92205E-05	-4,92205E-05	-4,92205E-05
		-4,69053E-05	-4,69053E-05	-4,69053E-05	-4,69053E-05
		-2,13056E-06	-2,13056E-06	-2,13056E-06	-2,13056E-06
		-0,0006429	-0,0006429	-0,0006429	-0,0006429
Матриця By	=	-4,9813E-05	-4,9813E-05	-4,9813E-05	-4,9813E-05
		-2,4203E-05	-2,4203E-05	-2,4203E-05	-2,4203E-05
		4,78827E-07	4,78827E-07	4,78827E-07	4,78827E-07
		0,00325526	0,00325526	0,00325526	0,00325526
Матриця Bz	=	-4,75517E-05	-4,75517E-05	-4,75517E-05	-4,75517E-05
		-3,53101E-05	-3,53101E-05	-3,53101E-05	-3,53101E-05
		2,327E-06	2,327E-06	2,327E-06	2,327E-06

Знаходимо невідомі величини матриці X=NˉⁱB:

		-0,001154194
X	=	0,005297297
		-0,017533943
		0,340481166

		-0,00011
Y	=	-0,00815
		-0,00101
		-0,02867

		-0,0003623
Z	=	0,043837574
		-0,0414491
		1,944980187

Підставивши в формулу (2) обчислені коефіцієнти отримаємо ΔV:

∆Vx	∆Vy	∆Vz
-0,00083078	-0,000608305	0,002432646
-0,00148237	-9,58897E-05	-0,001339966
-0,00026058	4,98857E-05	0,001184909
-0,00051218	0,000191043	-0,003340576

Використовуючи вираз VсlassB = Vмод – ΔV швидкості перманентних станцій мережі EPN будуть наступні:

	LIL2	AUTN	CHIZ	PUYV
Vclass B	-0,01328123	-0,013620125	-0,011119855	-0,012231731
	0,01697249	0,018168367	0,018498804	0,018799114
	0,01191318	0,008757006	0,011819927	0,007147673