(Sin (b2 * pi/180)/Sin ((180 - d2 - b2) * pi/180))=AC3 * (Sin (b3 * pi/180)/Sin ((180 - d3 - b3) * pi/180))=AC4 * (Sin (b4 * pi/180)/Sin ((180 - d4 - b4) * pi/180))=(AB1 + AB2 + AB3 + AB4)/4 AB (середнє) = raquo ;; ABSub
Рис. 5 Робоче вікно Visual Basic при обчисленні неприступного відстані
Результат обчислення неприступного відстані в Basic 6.0 приблизно АВ=223.408 м
2. Рішення прямого кутової засічки за формулою Юнга
У цьому завданню необхідно було розрахувати координати, віддаленої точки P, за формулами Юнга (рис. 5). Дано точки 1 і 2 з координатами і відповідно. Крім цього за допомогою кутомірного інструменту виміряні кути і. Варто звернути увагу, що при вирішенні геодезичних задач на відміну від загальноприйнятих обчислень осі координат 0x і 0y переставлені місцями.
Рис. 6. Схема розрахунку координат за формулами Юнга
Для обчислення координат віддаленої точки P використовувалися такі розрахункові формули:
,
.
Спочатку задача була вирішена засобами МS Excel. Результат представлений на рис. 7.
Рис. 7. Робоче вікно MS Excel при розрахунку координати за формулами Юнга
Результат обчислення координат точки Р в МS Excel м, м.
Потім рішення здійснювалося засобами MathCAD. В результаті чого було отримано наступний вигляд обчислення координат точки Р:
Рис. 8. Робоче вікно MathCad при розрахунку координати за формулами Юнга
Результат обчислення координат точки Р в MathCAD м, м.
Далі рішення виконувалося засобами MatLab. Результат представлений на рис. 7.
Програмний код:=input(laquo;x1-gt;raquo;);=input(laquo;y1-gt;raquo;);=input(laquo;x2-gt;raquo;);=input(laquo;y2-gt;raquo;);g=input(laquo;b1(grad)-gt;raquo;);m=input(laquo;b1(min)-gt;raquo;);s=input(laquo;b1(sec)-gt;raquo;);g=input(laquo;b2(grad)-gt;raquo;);m=input(laquo;b2(min)-gt;raquo;);s=input(laquo;b2(sec)-gt;raquo;);=b1g+b1m/60+b1s/3600;=b2g+b2m/60+b2s/3600;=b1*pi/180;=b2*pi/180;= (x1 * cot (b2) + x2 * cot (b1) -y1 + y2)/(cot (b1) + cot (b2))=(y1 * cot (b2) + y2 * cot (b1) + x1-x2 )/(cot (b1) + cot (b2))
Рис. 9. Робоче вікно MatLab при розрахунку координати за формулами Юнга
Результат обчислення координат точки Р в MatLab м, м.
Рішення даного завдання завершувалося обчисленням координат точки Р засобами Visual Basic 6.0. Було отримано наступний вид рішення, результат якого представлений на рис. 8.
Програмний код:
ctg (x)=Cos (x)/Sin (x) FunctionSub Command1_Click ()=Val (Text1 (0) .Text): Y1=Val (Text2.Text)=Val (Text6.Text): Y2= Val (Text7.Text)=Val (Text3.Text) + Val (Text4.Text)/60 + Val (Text5.Text)/3600=Val (Text8.Text) + Val (Text9.Text)/60 + Val ( Text10.Text)/3 600=3.1415926=b1 * pi/180=b2 * pi/180=(X1 * ctg (b2) + X2 * ctg (b1) - Y1 + Y2)/(ctg (b1) + ctg (b2 ))=(Y1 * ctg (b2) + Y2 * ctg (b1) + X1 - X2)/(ctg (b1) + ctg (b2)) xp= raquo ;; xp, yp= raquo ;; ypSub
Рис. 10. Робоче вікно Visual Basic при розрахунку координати за формулами Юнга
Результат обчислення координат точки Р в Basic 6.0 приблизно м, м.
Таким чином, середнє значення координат точки Р, розрахованих за формулами Юнга, виявилися рівні xp=52.225 м yp=45.559 м.
. Рішення прямого кутової засічки за формулами Гаусса
У цьому завданню необхідно було вирахувати координати віддаленої точки P, за координатами двох заданих точок 1 і 2 і дирекційний кутах і (рис. 11).
Рис. 11. Схема для обчислення координат за формулами Гаусса
Для обчислення координат віддаленої точки P використовувалися формули Гауса:
,
.
Спочатку задача була вирішена засобами МS Excel. Результат представлений на рис. 12
Рис. 12. Робоче вікно MS Excel при розрахунку прямої кутової засічки за формулами Гаусса
Результат обчислення координат точки Р в МS Excel м, м.
Потім рішення здійснювалося засобами MathCAD. В результаті чого було отримано наступний вигляд обчислення координат точки Р:
Рис. 12. Робоче вікно MathCad при розрахунку прямої кутової засічки за формулами Гаусса
Результат обчислення координат точки Р в MathCAD м, м.
Далі рішення виконувалося засобами MatLab. Результат представлений на рис. 11.
Програмний код:
=input(laquo;x1-gt;raquo;);=input(laquo;y1-gt;raquo;);=input(laquo;x2-gt;raquo;);=input(laquo;y2-gt;raquo;);=input(laquo;a(grad)-gt;raquo;);=input(laquo;a(min)-gt;raquo;);=input(laquo;b(grad)-gt;raquo;);=input(laquo;b(min)-gt;raquo;);=ag+am/60;=bg+bm/60;=b*pi/180;=a*pi/180;= (x1 * tan (a) -x2 * tan (b) -y1 + y2)/(tan (a) -tan (b))=y1 + (xp-x1) * tan (a)
Рис. 12. Робоче вікно MatLab при розраху...
