close

Вход

Забыли?

вход по аккаунту

?

Моделирование рельефа с использованием ортонормированных полиномов Чебышева.

код для вставкиСкачать
УДК 528.72
В.С. Коркин
СГГА, Новосибирск
МОДЕЛИРОВАНИЕ РЕЛЬЕФА С ИСПОЛЬЗОВАНИЕМ ОРТОНОРМИРОВАННЫХ
ПОЛИНОМОВ ЧЕБЫШЕВА
Предлагается алгоритм построения цифровой модели рельефа (ЦМР) по
нерегулярной сетке исходных точек, с соблюдением непрерывности
моделируемой поверхности.
Пусть в некоторой области  R в N точках с координатами
R R
XK
, YK   R измерены значения неизвестной функции рельефа f(x,у),
определенной на всей области моделирования  
n
 i ,
i 1
где i   j  0
для i, j .
Будем аппроксимировать данную функцию бикубическим полиномом:
3
3
  aijRTˆi ( x)Tˆ j ( y) ,
R
P ( x, y ) 
(1)
i 0 j 0
где aijR - неизвестные коэффициенты; Tˆi ( x),Tˆ j ( y) - ортонормированные
полиномы Чебышева, получены путем нормирования полиномов Чебышева lго рода: T0 ( x)  1, T1 ( x)  x , T2 ( x)  2 x 2  1 , T3 ( x)  4 x 3  3x , номера i и j
указывают на степень полинома.
Согласно методу наименьших квадратов запишем функционал

N

K 1
3 3
( f ( xk yk )    aijRTˆi ( xk )Tˆ j ( yk ))2  min .
R
(2)
i 1 j 1
Дифференцируя функционал (2) по неизвестным, получим систему
нормальных уравнений, после решения которой найдем коэффициенты aij .
Для решения поставленной задачи разобьем область  на ячейки R ,
имеющие форму квадратов (прямоугольников).
В каждой ячейке сетки аппроксимация выполняется независимо друг от
друга. После нахождения неизвестных коэффициентов aij переходим ко
второму этапу моделирования. На этом этапе выполняется аппроксимация по
алгоритму, изложенному в [1], полиномом вида
P ( x, y ) 
3
3
  aij Tˆi ( x)Tˆ j ( y) .
(3)
i 0 j 0
Отличие состоит только в том, что аппроксимировать будем не функцию
рельефа, представленную дискретно, а уже полученную функцию (1), которая
непрерывна на каждом элементарном участке  R .
Запишем минимизируемый функционал и условие непрерывности на
границах элементарных участков  R
1 1 3
2

Rˆ
Rˆ

ˆ
ˆ
      aij Ti ( x)T j ( y )    aij Ti ( x)T j ( y )  dxdy  min


i 0 j 0
1 1 i  0 j  0

3
3
3
(4)



i 0

3
Rˆ
l 1, r
 aij Ti (1)  b j  0 ( j  0,...,3)
i 0

.
(5)
3
Rˆ
l,r

 aij T j (1)  ai  0

j 0

3

Rˆ
l , r 1
 0 (i  0,...3) 
 aij T j (1)  ai

j 0
В выражениях (4) и (5):
aij - коэффициенты полинома, полученные из первой аппроксимации;
3
 aijRTˆi (1)  b,jr
0
aij - искомые коэффициенты полинома; Tˆi - 1,Tˆ j 1 - значения полиномов
Чебышева на границах элементарного участка; b j , ai - коэффициенты
полиномов профилей, полученных по линиям сетки.
Перед тем, как решать задачу вторичной аппроксимации,
предварительно строиться сетка профилей. Полиномы, аппроксимирующие
функцию рельефа по линиям сетки, получают из выражения (1) путем
фиксирования одной из координат x  d e или y  cr .
По линиям сетки, параллельным осям x и y имеем
3
P( x)   aiTˆi ( x) , P( y ) 
i 0
3
 b jTˆ j ( y) ,
(6)
j 0


(i  0,...3; j  0,...,3) .
где ai  aij T j (cr ) , b j  aij Ti (d e ) ,
Теперь задача состоит в том, чтобы обеспечить условие непрерывности
по линиям сетки в ее узлах. Для этого вычисляются значения функции (1) в
узлах.
Далее приступаем к нахождению новых коэффициентов ai и bj полиномов
(6) при условиях непрерывности в узлах сетки.
Задача решается следующим образом. По линиям сетки параллельным
оси у, необходимо найти
1 3
2

  b  l , r Tˆ ( y )   bl , r Tˆ ( y )  dy  min
j
j

j
j

i 0
1 j  0

при граничных условиях
3
(7)

l,r ˆ
l,r
b
T
(

1
)

Z

0
 j j


j 0
(8)
.
3

 blj, r Tˆ j (1)  Z l , r 1  0
j 0

Если проинтегрировать левую часть выражения (7) с учетом свойств
ортонормированных полиномов Чебышева, получим функционал
3
3
3
(  b jl ,r   b lj,r ) 2  min .
j 0
(9)
i 0
Для решения (9) и (8) составляем функцию Лагранжа


3
3
3
j 0
i 0
i 0
L b j , 1, 2  (  b j l , r   b lj, r ) 2  1l , r (  b lj, r T (1)  Z l , r ) 
3
 1l , r (  b lj, r T (1)  Z l , r 1 ),
i 0
(10)
где 1 и  2 - множители Лагранжа, которые являются дополнительными
неизвестными; Z l , r , Z l , r 1 - усредненные значения функции (1) в узлах сетки
R .
Находим частные производные функции L по неизвестным b lj,r , 1 ,  2 ,
приравниваем их к нулю и после преобразований получим систему из двух
уравнений
d1l ,r  el2,r  l l ,r 
(11)
,
l ,r
l ,r
l ,r 1
e1  d 2  l

3
3
j 0
3
j 0
где d   Tˆ j (1)Tˆ j (1)   Tˆ j (1)Tˆ j (1), e 
3
Tˆ j (1)Tˆ j (1),
j 0
3
l l , r  2(  b jl , r Tˆ j (1)  Z l , r ), l l , r 1  2(  b jl , r Tˆ j (1)  Z l , r 1 ).
j 0
j 0
Решая (11) находим 1 и 2 и вычисляем значения b j по формуле
b jl , r  (1l , r Tˆ j (1)  l2, r Tˆ j (1))
l,r
bj 
.
2
Построение профилей, параллельных оси х, и нахождение
коэффициентов ai выполняется аналогично. Теперь решается задача (4) и (5)
для определения коэффициентов aij каждого элементарного участка с
соблюдением условий непрерывности на границах.
Запишем функцию Лагранжа
2
3 3
 3 3

L(aij , v j , w j , i , i )     aijl , r    aijl , r  


i 0 j 0
 i 0 j 0

 3

 3

 v lj, r   aijl , r Tˆi (1)  b lj, r   wlj1, r   aijl , r Tˆi (1)  b lj1, r  




 j 0

 j 0

 3

 3

il , r   aijl , r Tˆi (1)  a lj, r      aijl , r Tˆi (1)  a lj, r 1 , ( j  0,...,3), (i  0,...,3).




 j 0

 j 0

(12)
Найдем частные производные (12) по переменным (aij , v j , w j ,  i ,  i ) и
прировняем их к нулю

L
 2aijl , r  2aijl , r  v lj, r Tˆi (1)  wlj1, r Tˆi (1)    il.r Tˆ j (1)   il , r 1Tˆ j (1)  0
b j


3

L
l,r ˆ
l,r
  aij Ti (1)  b j  0,

l,r
v j
j 0

.
3
L

l,r ˆ
l 1, r
  aij Ti (1)  b j
 0, ( j  0,...,3)

l 1, r
w j

j 0

3
L

l,r ˆ
l,r
  a j Ti (1)  a  0,

 il , r j  0


3
L
l,r ˆ
l , r 1

  a j Ti (1)  a
 0, (i  0,...,3)
l , r 1

 i
j 0

(13)
Выразим из первого уравнения (13) aij
aijl , r  (v lj, r Tˆi (1)  wlj1, r Tˆi (1)  il.r Tˆ j (1)   il , r 1Tˆ j (1)
l,r
aij 
2
(14)
и подставим в следующие уравнения системы (13).
В результате получим систему нормальных уравнений с шестнадцатью
неизвестными множителями Лагранжа. Коэффициенты d, l и c0 ,...,c9
системы нормальных уравнений имеют постоянные значения и зависят
только от значений ортонормированных полиномов Чебышева на концах
интервала [-1,1].
3
d   Tˆi (1)Tˆi (1) ;
i 0
l
3
Tˆi (1)Tˆi (1) ,
i 0
c0  Tˆ0 (1)Tˆ0 (1),
c1  Tˆ0 (1)Tˆ1 (1),
c2  Tˆ0 (1)Tˆ2 (1),
c5  Tˆ1 (1)Tˆ2 (1),
c6  Tˆ1 (1)Tˆ3 (1),
c7  Tˆ2 (1)Tˆ2 (1),
c3  Tˆ0 (1)Tˆ3 (1),
c4  Tˆ1 (1)Tˆ1 (1),
c8  Tˆ2 (1)Tˆ3 (1),
c9  Tˆ3 (1)Tˆ3 (1),
 3

u1lj, r  2  aijl , r Tˆi (1)  b lj, r ,


 i 0

 3

u lj1, r  2  aijl , r Tˆi (1)  b lj1, r ,


 i 0

 3
 1li , r  2



aijl , r Tˆ j (1)  ail , r ,

 j 0

 3

 il , r  2  aijl , r Tˆ j (1)  ail , r 1 ,


 j 0

где u1lj, r , u lj1, r ,  1li , r ,  il , r - свободные члены уравнений.
Решая полученную систему, находим множители Лагранжа и вычисляем
для каждого элементарного участка aijl ,r коэффициенты аппроксимирующего
полинома (3). Необходимо заметить, что обратная матрица коэффициентов
нормальных уравнений постоянна и все решения сводятся к вычислению
свободных членов и умножению обратной матрицы на вектор свободных
членов преобразованной системы.
Таким образом, получаем аппроксимирующую поверхность рельефа,
описанную полиномом (3) в каждом элементарном участке при соблюдении
условий непрерывности по всему полю моделирования.
Точность аппроксимации будет зависеть от размера элементарного
участка, плотности расположения исходных точек на поверхности рельефа и
сложности самого рельефа.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Журкин И. Г., Коркин В.С. Об одном подходе к построению цифровой модели
рельефа // Изв. Вузов. Геодезия и аэрофотосъемка. – 1983. – № 4. – с. 66-71.
 В.С. Коркин, 2005
Документ
Категория
Без категории
Просмотров
5
Размер файла
411 Кб
Теги
ортонормированным, моделирование, чебышева, использование, рельеф, полиномов
1/--страниц
Пожаловаться на содержимое документа