Всем доброго времени суток!
Необходим универсальный алгоритм, который позволит разбивать ломаную, состоящую из некоего числа отрезков, на части не кратные числу отрезков.
Как делаю сейчас (приведу код на Delphi):
| Код |
for Interval:=1 to CatenaryArray[CatenaryNumber].VertNumber-1 do //цикл от второй до последней точки ломаной
begin
XApp1:=CatenaryArray[CatenaryNumber].X[Interval-1]; //координата X первой точки отрезка "1-2"
YApp1:=CatenaryArray[CatenaryNumber].Y[Interval-1]; //координата Y первой точки отрезка "1-2"
ZApp1:=CatenaryArray[CatenaryNumber].Z[Interval-1]; //координата Z первой точки отрезка "1-2"
XApp2:=CatenaryArray[CatenaryNumber].X[Interval]; //координата X второй точки отрезка "1-2"
YApp2:=CatenaryArray[CatenaryNumber].Y[Interval]; //координата Y второй точки отрезка "1-2"
ZApp2:=CatenaryArray[CatenaryNumber].Z[Interval]; //координата Z второй точки отрезка "1-2"
LOtr:=LOtr+sqrt(sqr(XApp2-XApp1)+sqr(YApp2-YApp1)+sqr(ZApp2-ZApp1)); //Считаем длину новой ломаной
if LOtr<=dLA then //если длина новой ломаной меньше требуемой (переменная dLA), то проводим операции для данного отрезка новой ломаной
begin
//проводим операции для этого отрезка
end
else //если же сумма длин предыдущих отрезков и текущего превышает требуемую, то разбиваем отрезок "1-2" на отрезок "1-3" и "3-2", где 3-это точка, лежащая на отрезке "1-2"
begin
prirost:=dLA-(LOtr-sqrt(sqr(XApp2-XApp1)+sqr(YApp2-YApp1)+sqr(ZApp2-ZApp1))); // высчитываем какую часть нам надо взять из текущего отрезка для достижения требуемой длины ломаной
XApp3:= XApp1+prirost*(XApp2-XApp1)/sqrt(sqr(XApp2-XApp1)+sqr(YApp2-YApp1)+sqr(ZApp2-ZApp1)); //координата X для нужной части отрезка
YApp3:= YApp1+prirost*(YApp2-YApp1)/sqrt(sqr(XApp2-XApp1)+sqr(YApp2-YApp1)+sqr(ZApp2-ZApp1)); //координата Y для нужной части отрезка
ZApp3:= ZApp1+prirost*(ZApp2-ZApp1)/sqrt(sqr(XApp2-XApp1)+sqr(YApp2-YApp1)+sqr(ZApp2-ZApp1)); //координата Z для нужной части отрезка
//проводим нужные операции с новым отрезком "1-3", суммируем все, что получили в предыдущих, считаем что новая ломаная обработана и обнуляем все ненужные переменные
//отрезок "3-2" зачисляем в следующую ломаную
LOtr:=sqrt(sqr(XApp3-XApp2)+sqr(YApp3-YApp2)+sqr(ZApp3-ZApp2)); //считаем длину второй части нового отрезка, которая не вошла в предыдущую часть ломаной
//проводим операции для этого отрезка
end;
end;
|
В общем-то все работает, есть конечно небольшие доработки, которые не стал тут описывать, т.к. это все детали, связанные с точностью.
Но вот в чем вопрос - все это будет работать, если длина требуемой ломаной больше чем длина любого из отрезков исходной ломаной (например, длина всей ломаной 170 метров, разбить надо на 10 частей по 17 метров, длина каждого отрезка в ломаной примерно 6-6.5 метров) . А как написать универсальный алгоритм, который будет работать независимо от требуемой длины новой ломаной? Например, если мне надо будет разбить ломаную на 1000 отрезков, хотя сама она состоит из 50, хотя для моей задачи это не совсем нужно, но в тоже время в моей задаче может возникнуть необходимость разбить ломаную из 10 отрезков на 20 частей, а это тоже самое...
Очень прошу помощи!!! |
