книги из ГПНТБ / Гельфанд Р.Е. Программа уравнивания и предвычисления точности геодезических сетей на плоскости для ЭВМ типа М-220 (Шифр ПГС-УП) [Текст] 1973. - 134, [2] с
.pdf20.
сторонам не измерялись, в исходную информацию вводятся лож ные элементы.
При расчетах по предвычислению точности сетей вся исход ная информация записывается точно также, как при уравнивании, только вместо измеренных элементов в мантиссы вписываются нули. Исключение составляют линейно-угловые сети с двумя группами измеренных направлений. В этих сетях величины 0','01 должна обязательно записываться неуказанном выше месте инфор мации, а в полигонометрическом ходе вместо £ следует за писать любой угол, больший 0','02.
Вконце основного числового массива записывается
"/CS фикт.".
Втом случае, когда рабочие координаты определяемых пунк тов должны быть насчитаны по линейным измерениям, задается дополнительный числовой массив, содержащий шфорицию о последовательности треугольников, из которых могут быть нас читаны рабочие координаты.
Например:
+ |
+ |
+ |
00 |
101 |
102 |
001 |
+ |
+ |
+ |
00 |
001 |
102 |
002 |
|
• |
• |
|
• |
|
|
Х'Хфикт.
Приведенная запись означает, что рабочие координаты пунк та I насчитываются по координатам пунктов 101 и 102 и длинам сторон треугольника I 0 I - I 0 2 - I ; координаты пункта 2 будут нас читаны по координатам пунктов I и 102 и длинам сторон треу гольника I-I02-2 и т.д.
•В указанной записи должны соблюдаться следующие правила: а) . Номер пункта, для которого будет производиться насчет координат, записывается последним.
б) . Обход треугольника в порядке записи номеров происходит против часовой стрелки.
21.
^.Результаты "счета,
Печатается заголовок, определяющий тип сети: "Уравнива ние триалгуляции", или "Уравнивание трилатерации", или "Урав- • нивание линейно-угловой сети". После этого заголовка при предвычислении точности печатается слово "Предвычисления".
Печатается к и вся исходная информация основного массива. Эта печать производится логическими группами, ука занными в предыдущем разделе: общие сведения, информация о дополнительных условиях и т.д.
При наличии второго числового массива в сетях трилатера ции далее печатается к и числовая информация второго массива. В линейно-угловых сетях эта информация печатается после угловых невязок и рабочих координат, насчитанных по угловым измерениям.
Все линейные величины результатов выдаются на печать в том же масштабе, что и исходные данные, все угловые величины- в градусном измерении.
Линейные величины —координаты, длины, ср.кв. ошибки коор динат и длин-округляются по правилу, указанному в предыдущем
.разделе.
Результаты выдаются на печать группами в следующем поряд
ке:
1 . После заголовка "Угловые невязки треугольников" печа таются номера точек, составляющих замкнутые треугольники,и угловые невязки.
Вконце отдельным числом печатается ошибка по Фереро.
Всетях трилатерации и при отсутствии замкнутых треуголь ников эта печать отсутствует.
Влинейно—угловых сетях и сетях триангуляции с двумя группами угловых измерений возможна проверка невязок не всех замкнутых треугольников; соответственно ошибка по Фереро бу дет не верной. То же самое может быть в том случае, когда пункты стоянки в исходной информации заданы не в порядке но меров.
2.Под заголовком "Рабочие координаты" печатаются рабочие координаты всех определяемых пунктов в порядке номеров. При задании всех рабочих координат в исходной информации эта
22.
печать отсутствует.
Признак + + + 3 вместо координаты пункта означает, что координаты этого пункта не насчитываются формулой i t o a .
После заголовка "Результаты" печатаются:
3.Окончательные координаты всех определяемых пунктов.
4.Средняягазадратическаяошибка единицы веса q .
В сетях трилатерации q - в линейном измерении, в остальных-
вградусном.
Впрограмме принимаются значения весов:
В триангуляции рх = |
I , |
pff |
- заданное значение, л |
|||
в линейно-угловых сетяхJ>nl= |
|
>fs= ^ у * |
||||
в трилатерации |
А = |
mSz±(^ |
; |
Л,J 117 <%г'? |
||
где rn^j- |
( I ) - ср. кв. ошибка на I км хода измерений I группы. |
|||||
При предвычислениях |
точности |
q |
задается программой и |
|||
равна: |
триангуляции - |
I " |
, |
|
|
|
в |
|
|
||||
в л1шейно-угловых сетях -п?п1 |
, |
|
||||
в |
трилатерации |
~mfJ{I). |
|
|
||
5. Средниегсвадратическиеошибки координат определяемых |
||||||
пунктов: rn^i f mSi, |
, ntX£ |
, ... |
|
|
||
Результаты, указанные выше в п/пунктах 3,4 и 5 печатают |
||||||
ся при уравнивании дхадвы: для первого и для второго уравни вания.
6. Для каздого указанного в исходной информации элемента
'печатаются 8 чисел: |
+ |
+ + 0 к -1 |
а) . Номера пунктов |
||
б) . Дирекдаонный угол |
<cKi |
|
в) . Длина стороны |
S K i . |
|
Величины *-Ki и S K L насчитываются по округленным оконча |
||
тельным координатам. Поэтому в пределах точности они могут
отличаться от твердых |
и S ° i t если таковые задаваясь. |
||
г) . Ср.кв. ошибка дирекционного угла |
тлкС |
||
д) . Ср.кв. ошибка длины стороны |
тЛк^. |
|
|
Для стороны сети между |
твердыми пунктами вместо/7?,л-и т6к |
||
равных нулю, печатаются укороченные числа |
+ + |
+ 0 0 0 0 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
23. |
+ + + |
0 |
0 |
0 |
|
0 |
0 . О |
(см.контрольные примеры). |
|||||
|
е) . Знаменатель дроби |
|
|
|
|
|
||||||
|
|
|
у |
|
при заданном направлении, |
|
||||||
|
gr. ^ т |
|
|
|||||||||
|
^ |
— |
— г |
- при заданной длине |
|
|
|
|||||
|
|
в масштабе 10 |
|
. Печатаются три первые |
|
|||||||
|
|
значащие цифры, остальные•отбрасываются |
|
|||||||||
|
|
без округления. |
|
|
|
|
|
|
||||
Если m t K i |
или поправка в длину стороны &SKi |
|
равны нулю, |
|||||||||
то вместо знаменателя соответствующей дроби печатается |
|
|||||||||||
укороченное число + |
+ |
+ |
0 |
0 |
О |
О О . |
|
|||||
В. линейно-угловых сетях знаменатель дроби |
s |
1. |
- пе |
|||||||||
чатается с признаком. |
|
|
|
|
|
|
|
«<-/\uSKi\ |
||||
При величине знаменателя большей I 000 ООО печатается |
||||||||||||
число + |
+ |
+ |
0 |
0 |
0 |
0 |
0 |
999 |
999. |
|
|
|
|
ж) .Поправка: |
|
|
|
|
|
|
|
|
|||
|
|
для измеренного направления - поправка в угол |
||||||||||
|
|
между данным направлением к — - |
L |
и предыду |
||||||||
|
|
щим на данном пункте к — — |
/ |
|
|
|
||||||
для измеренной длины - поправка в длину
С признаком печатаются поправки в элементы, содержавшие признак в исходной информации, и поправки в длину в линейноугловых сетях.
+ |
Вместо Sy |
для 1-ого |
направления на пункте печатается |
||||||
+ |
+ |
0 |
0 |
0 |
0 |
0. |
|
|
|
+ |
При предвычислениях вместо поправок печатаются числа |
||||||||
+ |
+ |
0 |
0 |
0 |
0 |
0. |
|
|
|
|
|
|
з) . Признак - разделитель - - |
+ |
0 0 . |
||||
|
|
|
|
Далее печатается каталог результатов. |
|||||
|
В каталоге все линейные величины печатаются в метрах, |
||||||||
все угловые - в градусном измерении. |
|
|
|||||||
|
Величины <<-kL |
nm^KiB |
|
каталоге округляются до десятых |
|||||
долей секунды; поправки |
печатаются до десятых долей без |
||||||||
округления. Угловую величину в каталоге, не содержащую запя-
24.
той, отделяющей десятые секунда, следует читать, как градусы и минуты. Если в каталоге вместо числа напечатана группа
символов /// |
, то это означает, что |
количество значащих цифр |
в числе превосходит предусмотренное |
в каталоге. |
|
Печать X |
в графе каталога "На пункт л/ " означает, что |
|
эта строка каталога соответствует элементу П группы точности; в линейно-угловых сетях такими элементами являются измерен ные длины.
Печать X вместо :в конце строки означает, что строка соответствует элементу, не включенному в расчет - имеющему признак в исходной информации.
Если на пункте задавался только один элемент, то в ката логе печатается строка с координатой у- и пустыми остальны ми графами.
Кроме указанных величин на печать могут выдаваться неко торые значения для выявления грубых ошибок, допущенных в измерениях или при составлении исходной информации.
Печать слов "Расчет не окончен" может означать: 1) . Число заданных элементов ^ 350.
2). Рабочие координаты некоторых пунктов не насчиты вается формулами Юнга.
3). Неправильно составлена информация о базисах, твердых дирекционных углах и привязках.
Если при насчете рабочих координат пункта координаты, насчитанные из трех разных треугольников, различаются более чем на 0,5x10"^, то на печать будут выданы три восьмеричных числа, составленных из "сквозных" номеров пунктов, составляю щих эти_ треугольники
о о о к; к'л к'3 , ...
Сквозные номера для твердых пунктов К'= 2(к- 101); для определяемых К' = 2 ( К +.$,-1).
|
На печать выдаются также свободные члены уравнений пог |
решностей в том случае, если они превосходят: |
|
|
для измеренного направления - I ' , |
I |
для измеренной длины - 1 м . |
25.
Перед печатью результатов в этом случае для каждого такого элемента печатаются номера пунктов и свободные чле ны. В линейно-угловых сетях линейному свободному члену (соответствующему измеренной длине) приписывается признак.
Если в каком-либо пункте на печать выдаются свободные члены уравнений погрешностей для всех направлений, исключая первое в исходной информации, то это может, как правило, означать:
• а). Координаты пункта заданы или насчитаны грубо. В этом случае будут также встречаться на печати большие свободные члены для направлений (и длин) на данный пункт,
б). Ошибка в первом в исходной информации направлешв! на пункте. Свободные члены в этом случае будут одного порядка.
26.
7. Описание программы.
Программа состоит из пяти частей;
I . Ввод и печать исходной информации, перевод и расшиф ровка информации об измеренных элементах, определение и пе чать угловых невязок треугольников и рабочих координат, вы числение констант и коэффициентов подстановок из условных уравнений (3) и ( 4 ) .
П. Формирование нормальной системы.
Ш. Обращение матрицы и решение нормальной системы. Вччисление, печать и запись на МБ окончательных координат оп ределяемых пунктов, средних квадратических ошибок единицы веса и координат.
IV. Вычисление, печать и запись на МБ величин, описанных '
в пункте 6 раздела "Результаты расчета": k-i , оСкС , |
, |
тл , ms , ..'. . |
|
V. Печать каталога результатов. |
|
Вся программа вводится и переписывается на МБ единым массивом. При работе программы части поочередно переписывают ся в МЭЗУ. Команды записи программы на МБ и вызова в ШЗУ располагаются в ячейках OOOI-r-0025.
Расположение программы в памяти ЭВМ.
й части |
|
МБ |
|
ЮЗУ |
||
I |
I |
0030 - |
I 2573 |
5200 |
- 7742 |
|
п |
I |
2574 - |
I 3437 |
7100 |
- 7742 |
|
ш |
I |
3440 |
- 1.4303 |
7100 |
- 7742 |
|
и |
I |
4304 |
- |
I 5147 |
7100 |
- 7742 |
У |
I |
5150 |
- |
I 7223 |
5200 |
- 7252 |
Свободные ячейки в конце части У могут использоваться для изменений и дополнений .программы. Последнее число каждой части - контрольная сумма части.
Для хранения констант л промежуточных результатов исполь-
27.
зуготся в Ш'ЗУ ячейки 0026 ~ 7742; 7750-4- 7767. Рабочие ячейки всех подпрограмм 7743 — 41-YJ.
Исходная ындхэрмация по окончании работы части I хранится в :.ЮЗУ или на :.Д> в преобразованном заде. Игажгоенис состав ляет шяоормацкя о дополнительных условиях. Она переписывается •на X в ячейки 2 6556 -г 2 GS02 без измеиепий и используется
•при повторном уравшизагши. |
|
|
Постояпнне места в кОШ отведены |
|
|
рабочим координатам пунктов - яч.0030 -f- 0173, |
|
|
шуормации о подстановках - яч.0174 ч-0251, |
|
|
ин(оормацгш о весах (или ошибках) - яч.0262 -f- 0267, |
|
|
константам - 7753 - f 7767. |
|
|
Для сетей с числом определяемых пунктов |
^ I 3 F q |
paci!23> |
рованная информация об измеренных цементах хранится в |
-..DSf |
|
в ячейках 1662 Ч- 7076 при работе частей I - |
I / . Прц числе |
|
определяемых пунктов > 13 ста ши)оршция хданитоя па X |
з |
|
ячейках 2 5260-4- 2 6555. |
|
|
В расшнфрованкой шбормацки кавдоцу измеренному элементу соответствуют две ячейки: в перво.: во 2-сь: :•: 3-егл адресе на ходятся сквозные номера точек к ц «1 , во второй - изме ренная величина, is первом адресе первой ячейки шгюрпацик
о первом измеренном на пункте элементе хреиштся удвоенное
число измерению: на этом пушете |
элементов. "Сквозные" коме; г. |
|
для тверда |
пушетов К' = 2( к. - |
I C I ) , а для ошэеделяомс; |
IC~Z{K+fc - |
I ) . |
|
Матрица нормальной систем: занимает ячейки .'£>,j', начшг.л с 1066. На этом «е месте прк работе частей U и Ы распола гается обратная матрица.
Результаты тасчета, округлешше :: переверните и ,;есч- .
тпчную систему дат печати каталога хранятся |
.1-.: |
|
|
координаты всех пу:и-:тов - яч.2 СОЗи — 2 0x74., |
|
||
cp.:3.ocn!diai едашяя веса г. ^оо^-шит - .т-:.2 Ы7у*а |
03--i. |
||
.»-1Я Ксч.;дого пзмегкзнного элемента велпчша: |
, °CKi |
, |
|
5^, т А , rr>gt |
знаменатель :•: ;ю::^л?.в::а - яч.2 0-134-f-2 0257. |
||
При предварительном расчете (без печати ::.••-г-•мох^-.) atiimcb |
|||
с-езульгатов на X |
не щюиззодится.' |
|
|
28.
Укажем место нахождения некоторых констант при работе программы, которые могут быть изменены с пульта:
Ср.кв. ошибка единицы веса при предвычислешш точности перед началом работы части П находится в ячейке 0027.
Допуски для печати величин больших свободных членов урав нений погрешностей находятся в части П программы:
яч.7466 |
I ' |
яч.7471 |
10' |
яч.7472 |
10' |
яч.7473 |
,-5 |
10' |
Величина I м при двух, трех, и четырех значавдс цифрах в исходнш информации о лгшейншс величинах после запятой.
29.
Блок-схема программы.
Ввод программы и запись на МБ
Останов яч.ОООЗ |
_ |
|
|
|
|
||
Перепись части I в ШЗУ |
|
|
|
Ввод и печать исходной информации. Перевод |
|
||
констант 10 — 2 |
_г _ |
|
|
Расшифровка информации об измеренных элементах |
Щ |
||
и связях. |
. |
|
во |
Масштабирование заданных ср.кв.ошибок |
£1* |
||
измерении. |
|
|
|
Вычисление рабочих координат
Вычисление констант и коэффициентов подстановок по условным уравнениям
Определение числа избыточных наблюдений
|
! t=tm |
Перепись части П в Ш З У |
||
|
|
|
||
Й |
|
| |
Формирование нормальной системы |
|
|
|
Перепись части Ш в Ш З У |
||
|
|
Обращение матрицы нормальной системы |
||
1 |
in |
Решение нормальной системы |
||
Вычисление ср.кв.ошибки ед.веса |
||||
•-ч |
||||
I |
|
lea |
Вычислениенормальной системе.ср.кв.ошибок неизвестных в |
|
(Н |
|
|||
оВычисление неизвестныхисключенных из
онормальной системней их ср.кв.ошибок
|
|
|
X |
|
|
I |
а |
Вычисление окончательных координат |
|||
Округление,перевод.печать и запись на |
|||||
|
|||||
|
oot- |
МЕР координат, д |
и ср.кв. ошибок координат |
||
|
Перевод окончательных координат 10 — 2 |
||||
|
0> |
Перепись части 1У в ЮЗУ |
|
||
|
Вычисление,округление,перевод, печать |
||||
|
|
||||
|
|
и запись на МБ величин к-чs<t, |
... |
||
Перепись части У в Ю З У Печать каталога результатов