Ну, как-то доказалось, наконец...
Утверждение начального поста можно изложить так.
Плоскую кривую длины
с кривизной, монотонно изменяющейся от значения
до
, можно заменить кривой, составленной из трёх гладко сoпряжённых круговых дуг с кривизнами
, и имеющей ту же длину, те же граничные точки и касательные в них.Рис. 1:
![Изображение](https://dl.dropboxusercontent.com/u/14367396/dxdy/Close/Fig1.jpg)
На Рис. 1 в левой колонке вверху показан график функции
![$k=k(s)$ $k=k(s)$](https://dxdy-02.korotkov.co.uk/f/1/9/f/19fd1a0a0dc15bb1ce69007a9ac2dee082.png)
,
![$0\le s\le L$ $0\le s\le L$](https://dxdy-03.korotkov.co.uk/f/a/6/d/a6d0eb316e2a719b088d984814f89d6b82.png)
. Это натуральное уравнение некой плоской кривой длины
![$L$ $L$](https://dxdy-02.korotkov.co.uk/f/d/d/c/ddcb483302ed36a59286424aa5e0be1782.png)
с кривизной, монотонно возрастающей от
![$k(0)=k_1$ $k(0)=k_1$](https://dxdy-04.korotkov.co.uk/f/b/b/6/bb631a278a2be3f26fcdce23d6f88c1382.png)
до
![$k(L)=k_2$ $k(L)=k_2$](https://dxdy-03.korotkov.co.uk/f/e/e/b/eebfe85e107bc5d15ca480c9f247a9b082.png)
. Чисто для понта добавлена лёгкая волнистость.
Под ним прорисован график функции
![$\tau(s)$ $\tau(s)$](https://dxdy-01.korotkov.co.uk/f/4/f/5/4f5fef9a9edf93045bf307e49d136cd182.png)
, описывающей угол наклона касательной к кривой:
![$$ \tau(s)=\tau_0+\int_{0}^{s}{k(t)}\,{dt}. \eqno(1) $$ $$ \tau(s)=\tau_0+\int_{0}^{s}{k(t)}\,{dt}. \eqno(1) $$](https://dxdy-01.korotkov.co.uk/f/4/2/9/4292f6378f1fba74e7b2aee43c8e4e1482.png)
Монотонность
![$k(s)$ $k(s)$](https://dxdy-04.korotkov.co.uk/f/f/f/d/ffd004fe212e0bda8a367b0fc4febf2882.png)
влечёт выпуклость графика
![$\tau(s)$ $\tau(s)$](https://dxdy-01.korotkov.co.uk/f/4/f/5/4f5fef9a9edf93045bf307e49d136cd182.png)
. Наклоны граничных касательных к графику равны
![$k_1$ $k_1$](https://dxdy-03.korotkov.co.uk/f/a/a/9/aa90653a26bc63b138fb304972d8158982.png)
и
![$k_2$ $k_2$](https://dxdy-03.korotkov.co.uk/f/a/8/e/a8ebf8c468236800b8ed78d42ddbfa5782.png)
.
Eщё ниже прорисована сама кривая
![$[x(s),y(s)]$ $[x(s),y(s)]$](https://dxdy-01.korotkov.co.uk/f/0/1/5/015673eee64bb636d65221db9e7fe39a82.png)
:
![$$
x(s)=x_0+\int_{0}^{s}{\cos\tau(t)}\,{dt},\quad
y(s)=y_0+\int_{0}^{s}{\sin\tau(t)}\,{dt}. \eqno(2) $$ $$
x(s)=x_0+\int_{0}^{s}{\cos\tau(t)}\,{dt},\quad
y(s)=y_0+\int_{0}^{s}{\sin\tau(t)}\,{dt}. \eqno(2) $$](https://dxdy-04.korotkov.co.uk/f/f/4/3/f4385b38690ea82c819046b0732be51382.png)
Кривая выпущена из начала координат
![$(x_0,y_0)=(0,0)$ $(x_0,y_0)=(0,0)$](https://dxdy-03.korotkov.co.uk/f/a/3/7/a37d0ac43143163edcd6837a43ff515082.png)
, касательная в начальной точке направлена вдоль оси абсцисс (
![$\tau_0=0$ $\tau_0=0$](https://dxdy-02.korotkov.co.uk/f/d/d/c/ddc9743bc038eee94d6f3c36644a8ad082.png)
).
В правой колонке исходная функция
![$k(s)$ $k(s)$](https://dxdy-04.korotkov.co.uk/f/f/f/d/ffd004fe212e0bda8a367b0fc4febf2882.png)
подменена кусочно-постоянной функцией
![$\bar{k}(s)$ $\bar{k}(s)$](https://dxdy-04.korotkov.co.uk/f/f/6/1/f61c258f8f883f85faff530774a5129982.png)
, также монотонной. Положение и кривизна
![$k_0$ $k_0$](https://dxdy-02.korotkov.co.uk/f/1/b/6/1b63f35b880a6307974273a3ff6063c982.png)
промежуточного участка вычислены так, чтобы обеспечить те же граничные условия (координаты и наклоны касательных), что и у исходной кривой. График
![$\bar{\tau}(s)$ $\bar{\tau}(s)$](https://dxdy-01.korotkov.co.uk/f/8/8/7/8871cb032aa57514ea028b3c109fca4282.png)
(бывшая
![$\lambda(s)$ $\lambda(s)$](https://dxdy-02.korotkov.co.uk/f/d/3/e/d3efbfa10592569f18e492377c72bfcc82.png)
) становится при этом кусочно-линейным (и по-прежнему выпуклым), а сама кривая --- "кусочно-окружностной". Различия между оригиналом и его окружностной аппроксимацией слабо заметны: точечки, которыми представлена оригинальная кривая, лишь слегка отклонены от круговых дуг.
Рис. 2:
![Изображение](https://dl.dropboxusercontent.com/u/14367396/dxdy/Close/Fig2.jpg)
На Рис. 2 приведены ещё пять примеров фунций
![$k(s)$ $k(s)$](https://dxdy-04.korotkov.co.uk/f/f/f/d/ffd004fe212e0bda8a367b0fc4febf2882.png)
, монотонно возрастающих от
![$k_1$ $k_1$](https://dxdy-03.korotkov.co.uk/f/a/a/9/aa90653a26bc63b138fb304972d8158982.png)
до
![$k_2$ $k_2$](https://dxdy-03.korotkov.co.uk/f/a/8/e/a8ebf8c468236800b8ed78d42ddbfa5782.png)
, и их 3-ступенчатые аналоги. Им соответствуют 5 кривых во второй колонке.
У любого скучающего человека естественно возникает вопрос: а как локализованы конечные точки этих кривых?
Возможность подмены всего
![$\infty$ $\infty$](https://dxdy-04.korotkov.co.uk/f/f/7/a/f7a0f24dc1f54ce82fecccbbf48fca9382.png)
-мерного пространства таких кривых трёхпараметрическим семейством тридуг позволяет скучающему человеку отмоделировать множество концевых точек тремя вложенными циклами: если
![$k_1,l_1$ $k_1,l_1$](https://dxdy-01.korotkov.co.uk/f/0/5/f/05f2ee87016ce5f9c5c5d3111d9df67482.png)
--- кривизна и длина первой дуги,
![$k_0,l_0$ $k_0,l_0$](https://dxdy-01.korotkov.co.uk/f/0/d/8/0d8e06ab95c7782dd5dc2a92aed0bf2982.png)
--- кривизна и длина промежуточной дуги,
![$k_2,l_2$ $k_2,l_2$](https://dxdy-03.korotkov.co.uk/f/6/4/4/644cebc2869ec74035fb1a1d6a30bcc582.png)
--- третьей, то параметрами семейства (и переменными циклов) служат
![$l_1$ $l_1$](https://dxdy-01.korotkov.co.uk/f/4/6/9/469f525d671e1e96713a0a17a13f246882.png)
,
![$k_0$ $k_0$](https://dxdy-02.korotkov.co.uk/f/1/b/6/1b63f35b880a6307974273a3ff6063c982.png)
и
![$l_2$ $l_2$](https://dxdy-04.korotkov.co.uk/f/7/2/5/7252ad06a4944da2b6628a58281cb88782.png)
с ограничениями
![$$
l_1>0,\quad l_2>0,\quad l_1+l_2\le L,\quad(l_0=L-l_1-l_2\ge0),
\quad k_1<k_0<k_2.\eqno(3)
$$ $$
l_1>0,\quad l_2>0,\quad l_1+l_2\le L,\quad(l_0=L-l_1-l_2\ge0),
\quad k_1<k_0<k_2.\eqno(3)
$$](https://dxdy-01.korotkov.co.uk/f/c/2/6/c26f75daa3586c5454989f3ebafdd69182.png)
Любая современная ЭВМ способна быстренько выполнить три вложенных цикла с разумным шагом по каждой
из переменных. Справа эти пять кривых нарисованы вместе, а серые точки показывают результат моделирования,
т.е. концы тысяч других кривых этого семейства. Границы этого множества, кривые
![$\Gamma_1$ $\Gamma_1$](https://dxdy-04.korotkov.co.uk/f/f/a/b/fabaf4ddb27af5131de71a9579a32e7982.png)
и
![$\Gamma_2$ $\Gamma_2$](https://dxdy-01.korotkov.co.uk/f/4/d/9/4d91c09273763aee35f2224b6aa8f6b982.png)
, угадываются лёгким напряжением мозга.
Если кривые
![$\Gamma_1$ $\Gamma_1$](https://dxdy-04.korotkov.co.uk/f/f/a/b/fabaf4ddb27af5131de71a9579a32e7982.png)
и
![$\Gamma_2$ $\Gamma_2$](https://dxdy-01.korotkov.co.uk/f/4/d/9/4d91c09273763aee35f2224b6aa8f6b982.png)
продлить за естественные границы изменения параметра, то получится Вау!, показанное на Рис. 3:
Рис. 3:
![Изображение](https://dl.dropboxusercontent.com/u/14367396/dxdy/Close/Fig3.jpg)