До початку сторінка 10 з 14
Таблица 1
Ординаты аппроксимирующих кривых обеспеченности сумм жидких осадков и глубин промерзания на водотоке Б
IСЕ1
А, В - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) СУМ.
Ж. ОСАД.
Е - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ
А=6.15400Е-03
В=2.56446Е-03
E=5.61809E-04
А, В - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) ГЛУБ. ПРОМ.
Е - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ
А= .29688
Р ЭМП
B= .400021
СУМ. Ж. ОСАД. НАБЛ.
E= .0542944 СУМ. Ж. ОСАД. АППР.
ГОД
ГЛУБ. ПРОМ. PACЧ
ГЛУБ. ПРОМ. АППР.
ГОД
.0384615
270.4
356.996
1982
2.33395
2.43049
1966
.0769231
270.1
329.181
1972
2.32161
2.36486
1964
.115385
263.4
305.387
1977
2.23754
2.30268
1979
.153846
263.3
284.8
1973
2.17201
2.24368
1977
.192308
257.5
266.814
1966
2.12654
2.18763
1962
.230769
256.5
250.965
1985
2.11576
2.13432
1972
.269231
248.7
236.893
1983
2.09443
2.08354
1975
.307692
240.1
224.316
1978
2.06688
2.03512
1981
.346154
236
213.007
1975
2.04223
1.9889
1984
.384615
215.2
202.783
1974
1.87887
1.94474
1985
.423077
210.9
193.496
1984
1.87116
1.90249
1970
.461538
201.1
185.022
1961
1.86676
1.86204
1973
.5
198.2
177.259
1969
1.80822
1.82327
1968
.538462
176
170.121
1979
1.78721
1.78608
1982
.576923
174.4
163.536
1962
1.78648
1.75039
1978
.615385
165
157.442
1967
1.77471
1.71609
1963
.653846
162.5
151.786
1968
1.772
1.6831
1986
.692308
146.1
146.522
1970
1.743
1.65137
1976
.730769
145.5
141.611
1976
1.72919
1.6208
1983
.769231
143.8
137.018
1980
1.69469
1.59135
1980
.807692
142.8
132.714
1981
1.66125
1.56295
1971
.846154
138.2
128.672
1963
1.64034
1.53555
1974
.884615
117.8
124.869
1965
1.5986
1.50909
1965
.923077
107.3
121.284
1971
1.56892
1.48352
1969
.961538
99.2
117.9
1964
1.07274
1.45881
1967
Рис. 1. Кривые обеспеченности сумм жидких осадков и глубин промерзания на водотоке Б (питание грунтовыми водами)
Площадь фильтрационной зоны до постройки моста принимается равной пятикратной длине отверстия моста. При толщине фильтрационного слоя 4,9 м площадь фильтрации после промерзания в естественных условиях
wФ=5[3×17(4,90-2,26)]=672 м2.
Площадь стеснения фильтрационного потока постройки моста (см. рис. 3 и табл. 1):
wС=nОП(b-2b1)(d-df)+(d'f-df)(l-2b1)nПР=4(0,5-2×0,235)(4,90-2,26)+(2,56-2,26)(51-2×0,235)3=
=55,7 м2;
d'f =2,56 - глубина промерзания под мостом.
Коэффициент, учитывающий стеснение,
По формуле (6) определяем прогнозируемый объем наледи
VНР=gnVНИkxkfkC=1,2×12400×1,52×1,17×1,08=28580 м3.
Средняя расчетная толщина наледи на наледной поляне
м.
Максимальная толщина наледи
Zmax p=YZСР Р=1,13×1,32=1,5 м,
где
ПРИЛОЖЕНИЕ 4 Программа I СЕ 2 (Лед-II) расчетов осадков, глубин промерзания и ординат аппроксимирующих кривых обеспеченности при прогнозировании наледей поверхностных вод II. ИНСТРУКЦИЯ ПО ЭКСПЛУАТАЦИИ ПРОГРАММЫ I СЕ 2
Программу I СЕ 2 применяют для расчета жидких осадков, глубин промерзания поверхностного потока и ординат аппроксимирующих кривых при прогнозировании наледей поверхностных вод.
Программа I СЕ 2 составлена на алгоритмическом языке БЕЙСИК и выполняется в режиме интерпретации под управлением операционной системы реального времени с разделением функций (РАФОС). Объем оперативной памяти, занимаемой интерпретатором, не более 14 К слов. Минимальная конфигурация, необходимая для работы интерпретатора, это базовый комплект УВК СМ 4 (или микро-ЭВМ «Электроника-60»).
1. Инструкция по заполнению исходных данных
Исходные данные в программу I СЕ 2 заносят двумя способами:
1) путем предварительной организации блока данных n; g0; (lC)j, j=1, ... , 6; (hC)ij, i=1, ... , n; j=1, ... , 6; (tB)ij, i=1, ..., n; j=1, ..., 6; (О0)i, (О1)i, (О2)i, (О3)i, i=1, ..., n; в тексте программы I СЕ 2 занесением данных в виде операторов. В этом случае операторы READ и DATA совместно попользуются для ввода численной информации;
2) путем предварительного изменения (если в этом есть необходимость) правой части операторов LET, которые определяют некоторые физические константы, характеризующие процесс lЛ, a, Lv.
Для организации блока данных занести n и g0 в строку 26. Например, 26 DATA, 25, 1961
(lC)1 в строку 450
(lC)2 в строку 500
(lC)3 в строку 550
(lC)4 в строку 600
(lC)5 в строку 650
(lC)6 в строку 700
Например, 450 DATA, 2031
(hC)i1 i=1, ..., n в строки 460, 465
(hC)i2 i=1, ..., n в строки 510, 615
(hC)i3 i=1, ..., n в строки 560, 565
(hC)i4 i=1, ..., n в строки 610, 615
(hC)i5 i=1, ..., n в отроки 660, 665
(hC)i6 i=1, ..., n в строки 710, 715
Например,
460 DATA 2, 7.67, 4, 7.33, 4.67, 2, .33, 1.33, 3.33, 8.67
465 DATA 1.67, 5, .67, 5, 5.33, 1.67, 7.33, 2.33, 9, .67, 0, 1, .33, 8.33, 8.33
(tВ)i1 i=1, ..., n в строки 490¸497
(tВ)i2 i=1, ..., n в строки 540¸547
(tВ)i3 i=1, ..., n в строки 590¸597
(tВ)i4 i=1, ..., n в строки 640¸647
(tВ)i5 i=1, ..., n в строки 690¸697
(tВ)i6 i=1, ..., n в строки 740¸747
Например,
490 DATA -11.6, -13,2, -9.8, -14.5, -12.1
495 DATA -10.2, -4.0, -8.8, -9.9, -11.7
496 DATA -6, -5.2, -3.9, -4.4, -3.3, -6, -11.4, -5.4, -8, -8.8
497 DATA -7, -5.7, -7, -2.7, -5.3
2. Расшифровка результатов
Результаты расчета выводятся на печать в виде таблицы:
1 столбец - эмпирические вероятности превышения сумм жидких осадков и глубин промерзания за n лет наблюдений;
2 столбец - суммы жидких осадков по данным гидрометеослужбы за n лет наблюдений;
3 столбец - ординаты аппроксимирующей кривой сумм жидких осадков;
4 столбец - годы наблюдений осадков;
5 столбец - расчетные глубины промерзания поверхностного потока;
6 столбец - ордината аппроксимирующей кривой глубин промерзания поверхностного потока;
7 столбец - годы, для которых вычислялась глубина промерзания.
3. Инструкция оператору
а) Загрузить операционную систему РАФОС.
б) Запустить интерпретатор БЕЙСИК, например, под управлением si - монитора, т. е.
Набрать команду BASIC
Печатается
BASIC/RAFOS V.02
OPTIONAL FUNCTIONS (ALL, NONE, OR INDIVIDUAL)?
Набрать отпет А (ответ А включает все средства расширения функций интерпретатора).
Печатается
READY
в) Предварительно занести в текст программы исходные данные, пользуясь инструкцией по заполнению исходных данных.
г) Запустить программу I СЕ 2:
Если программа I СЕ 2 уже находится в памяти, то набрать команду RUN.
Если программа I СЕ 2 находится на устройстве файловой структуры (диски), то набрать команду RUN STR, где STR - спецификация файла, например набрать команду RUN DX1 : I СЕ 2. BAS.
д) Печатается таблица результатов, описанная в расшифровке результатов.
III. ТЕКСТ ПРОГРАММЫ I СЕ 2 (ЛЕД-II)
4 КЕМ ICE2
5 REM ПРОГРАММА РАСЧЕТА ГЛУБИН ПРОМЕРЗАНИЯ И
6 REM СУММ ЖИДКИХ ОСАДКОВ ЗА N ЛЕТ НАБЛЮДЕНИЙ
7 REM АППРОКСИМАЦИЯ КРИВЫХ ПО МЕТОДУ НАИМЕНЬШИХ КВАДРАТОВ
8 W9=SYS(7,0)
9 OPEN «LP:» FOR OUTPUT AS FILE #1
10 PRINT #1, «ICE2»\PRINT #1,
11 DIM L(6), H(30,6), Т1(30,6)
15 DIM F(30), G1(30), F2(30)
16 DIM F5(30,6)
22 DIM O0(30), O1(30), O2(30), O3(30), O(30), G2(30)
23 DIM X(30), Y(30), U(30), V(30), V1 (30)
24 DIM Y1(30), Y2(30)
25 READ N, GO
26 DATA 25, 1961
27 REM
33 Z1=2.22
34 A=23.33
38 Z2=3.07280E+08
39 REM ЦИКЛ ВВОДА ПО 6 МЕСЯЦАМ ЗИМЫ
40 FOR J=1 TO 6
49 REM ВВОД КОЭФФИЦИЕНТА ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА КАЖДОГО МЕСЯЦА
50 READ L(J)
59 REM ВВОД СРЕДНЕМЕСЯЧНОЙ ТОЛЩИНЫ СНЕЖНОГО ПОКРОВА КАЖДОГО ГОДА
60 FOR I=1 TO N
70 READ H(I, J)
80 NEXT I
82 REM ВВОД СРЕДНЕМЕСЯЧНЫХ ТЕМПЕРАТУР КАЖДОГО ГОДА
84 FOR I=1 TO N
86 READ T1(I, J)
88 NEXT I
90 REM ОКОНЧАНИЕ ЦИКЛА ВВОДА ПО 6 МЕСЯЦАМ ЗИМЫ
92 NEXT J
100 REM WATER
110 FOR I=1 TO N
210 F5(I, 1)=0
219 FOR J=2 TO 6
230 H5=Z1/A+H(I, J)*01*Z1/L(J)
240 F5(I, J)=.8*(SQR((H5+F5(I, J-1))^2-2*Z1*T1(I, J)*720*3600/Z2)-H5)
280 NEXT J
360 F(I)=F5(I,6)
364 REM ФОРМИРОВАНИЕ РЯДА ЛЕТ НАБЛЮДЕНИЙ ТОЛЩИНЫ СНЕЖНОГО ПОКРОВА
365 G1(I)=G0+1
390 NEXT I
399 REM РАНЖИРОВАНИЕ РЯДА ГЛУБИН ПРОМЕРЗАНИЙ
400 К=0
402 FOR I=1 TO N-1
404 IF F(1)>F(I+1) THEN 420
406 R1=F(I+1)
408 R2=G1(I+1)
410 F(I+1)=F(I)
412 G1(I+1)=G1(I)
414 F(I)=R1
416 G1(I)=R2
418 K=1
420 NEXT I
422 IF K=1 THEN 400
449 REM КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА ОКТЯБРЯ
459 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА ОКТЯБРЯ
460 DATA 2, 7.67, 4, 7.33, 4.67, 2, .33, 1.33, 3.33, 8.67
465 DATA 1.67, 5, .67, 5, 5.33, 1.67, 7.33, 2.33, 9, .67, 0, 1, .33, 8.33, 8.33
489 REM СРЕДНЯЯ ТЕМПЕРАТУРА ОКТЯБРЯ
490 DATA -11.6, -13.2, -9.8, -14.5, -12.1
495 DATA -10.2, -4.0, -8.8, -9.9, -11.7
496 DATA -6, -5.2, -3.9, -4.4, -3.3, -6, -11.4, -5.4, -8, -8.8
497 DATA -7 -5.7, -7, -2.7, -5.3
499 REM КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА НОЯБРЯ
509 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА НОЯБРЯ
510 DATA 15.67, 8.67, 8.67, 12, 6.33, 14.67, 2, 3.67, 14.33, 14,33
515 DATA 3.33, 5.33, 7.33, 4.67, 10.67, 6.33, 14.67, 5, 0, 4.33, 9.33, 10.33, .33, 5.67, 11.67
539 REM СРЕДНЯЯ ТЕМПЕРАТУРА НОЯБРЯ
540 DATA -23.6, -22.7, -18.6, -19.6, -26.1
545 DATA -19.2, -12.0, -14.3, -15.6, -21.3
546 DATA -15.1, -21.6, -14.7, -19, -17.9, -15.3, -16.1, -15.8, -19.1, -17.5
547 DATA -14.3, -20.4, -17.4, -15.6, -11.7
549 REM КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА ДЕКАБРЯ
550 DATA .2031
559 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА ДЕКАБРЯ
560 DATA 11.33, 11, 8.67, 19.67, 6.33, 35, 15, 18, 11, 15.33
565 DATA 5, 19.33, 1.2.33, 4.67, 8.67, 6, 11.33, 0, 15, 6.67, 12, 9, 3.33, 8, 10.33
589 REM. СРЕДНЯЯ ТЕМПЕРАТУРА ДЕКАБРЯ
590 DATA -20.6, -23.1, -27.2, -26.7, -28.9
595 DATA -25.6, -27.8, -29.6, -17.3, -16
596 DATA -20.1, -19.2, -16.7, -20.5 -15.9, -20.9, -18.6, -28.1, -21.9, -15.2
597 DATA -25.0, -25.2, -18.2, -25.8, -20.1
599 RЕМ КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА ЯНВАРЯ
600 DATA .349
609 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА ЯНВАРЯ
610 DATA 8, 24, 8, 32, 7, 39.67, 13.33, 20.67, 10.33, 18.33
615 DATA 5, 14, 11, 5, 8.33, 5.67, 13, 8, 18.67, 6.33, 14.67, 19, 13.33, 8, 7.67
639 REM СРЕДНЯЯ ТЕМПЕРАТУРА ЯНВАРЯ
640 DATA -26.6, -19.6, -31, -26.7, -27.1
645 DATA -26.8, -21.2, -19.8, -20.2, -18.6
646 DATA -28.4, -30.4, -17.8, -23.7, -19.5, -19.2, -23.4, -21.9, -14.9, -20.7
647 DATA -26.1, -16.3, -20.6, -21.7, -19.9
649 REM КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА ФЕВРАЛЯ
650 DATA .349
659 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА ФЕВРАЛЯ
660 DATA 8, 24.67, 6.67, 28, 6, 59.33, 10, 26, 10, 13.33
665 DATA 8.33, 13, 14, 6.67, 12.67, 4.67, 8.33, 9.33, 20, 9, 17.33, 18.33, 11, 14.67, 8
689 REM СРЕДНЯЯ ТЕМПЕРАТУРА ФЕВРАЛЯ
690 DATA -27.1, -23.6, -25.8, -22.4, -29.4
695 DATA -18.7, -19.3, -18.6, -23.6, -21.1
696 DATA -14.8, -21.5, -18.7, -16.7, -18.6, -23.4, -21.3, -15, -19.9, -19
697 DATA -16.2, -21.4, -20.2, -15, -17
699 REM КОЭФФИЦИЕНТ ТЕПЛОПРОВОДНОСТИ СНЕЖНОГО ПОКРОВА МАРТА
700 DATA .349
709 REM СРЕДНЯЯ ТОЛЩИНА СНЕЖНОГО ПОКРОВА МАРТА
710 DATA 6, 26, 6.67, 27, 15, 64, 10.67, 25, 8, 12.33
715 DATA 5.33, 12.33, 13.33, 8.33, 12.67, 3.67, 8.33, 4, 18, 7.33, 16, 16, 9.67, 14.67, 8.33
739 REM СРЕДНЯЯ ТЕМПЕРАТУРА МАРТА
740 DATA -19.2, -24.4, -23.2, -20.2, -9.2
745 DATA -13.4, -17.4, -15.4, -17.2, -16.1
746 DATA -12.2, -16.7, -14, -14.4, -17, -18.6, -11.2, -16.4, -17, -16.9
747 DATA -18.9, -14.5, -13.1, -20.5, -16.8
1000 REM РАСЧЕТ СУММ ЖИДКИХ ОСАДКОВ
1043 REM ВВОД КОЛ-BA ЖИДКИХ ОСАДКОВ ЗА ИЮЛЬ КАЖДОГО ГОДА
1044 FOR I=1 ТО N
1046 READ O0(I)
1048 NEXT I
1049 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА АВГУСТ КАЖДОГО ГОДА
1050 FOR I=1 ТО N
1052 READ O1(I)
1054 NEXT I
1055 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА СЕНТЯБРЬ КАЖДОГО ГОДА
1056 FOR I=1 ТО N
1058 READ O2(I)
1060 NEXT I
1061 REM ВВОД КОЛ-ВА ЖИДКИХ ОСАДКОВ ЗА ОКТЯБРЬ КАЖДОГО ГОДА
1062 FOR I=1 ТО N
1064 READ O3(I)
1066 NEXT I
1499 REM РАСЧЕТ СУММ ЖИДКИХ ОСАДКОВ ЗА 3 МЕСЯЦА
1500 FOR I=1 TO N
1510 O(I)=0
1520 NEXT I
1530 FOR I=1 TO N
1540 O(I) =O0(I)+O1(I)+O2(I)+O3(I)
1560 NEXT I
1609 REM ФОРМИРОВАНИЕ РЯДА ЛЕТ НАБЛЮДЕНИЙ ЖИДКИХ ОСАДКОВ
1610 FOR I=1 TO N
1620 G2(I)=G0-1+I
1630 NEXT I
1639 REM РАНЖИРОВАНИЕ РЯДА СУММ ЖИДКИХ ОСАДКОВ
1640 K=0
1650 FOR I=1 TO N-1
1660 IF O(I)>O(I+1) THEN 1740
1670 R1=O(I+1)
1680 R2=G2(I+1)
1690 O(I+1)=O(I)
1700 G2(I+1)=G2(I)
1710 O(I)=R1
1720 G2(I)=R2
1730 K=1
1740 NEXT I
1750 IF K=1 THEN 1640
1760 FOR I=1 TO N
1764 REM ВЫЧИСЛЕНИЕ ЭМПИРИЧЕСКОЙ ВЕРОЯТНОСТИ
1765 F2(I)=I/(N+1)
1780 NEXT I
1949 REM ЖИДКИЕ ОСАДКИ ЗА ИЮЛЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ
1950 DATA 83.1, 78.8, 49.5, 24, 35.8
1955 DATA 74.1, 50.8, 73.3, 98.2, 44.1, 17.4, 79.3, 15.5, 66.2, 67.7
1956 DATA 41.4, 93.9, 149.3, 64, 67.7, 62.7, 34.9, 94.3, 93.8, 163.2
1959 REM ЖИДКИЕ ОСАДКИ ЗА АВГУСТ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ
1960 DATA 37.6, 51.2, 21.8, 73.9
1970 DATA 60, 76.6, 17.9, 30.1, 61.6, 61, 46.9
1975 DATA 160.9, 59.2, 18.2, 117.5, 102.2, 75.7, 59.3
1976 DATA 73.3, 33.7, 24.3, 111.4, 103.0, 84.6, 64.7
1979 REM ЖИДКИЕ ОСАДКИ ЗА СЕНТЯБРЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ
1980 DATA 80.4, 44.4, 8.4
1990 DATA .5, 22, 106.8, 96.3, 59.1, 38.4, 41, 43
1995 DATA 28.1, 98.6, 97.4, 50, 1.9, 93.8, 15.2
1996 DATA 38.7, 42.4, 38.9, 47.7, 51.4, 12.9, 28.1
1998 REM ЖИДКИЕ ОСАДКИ ЗА ОКТЯБРЬ ПО КАЖДОМУ ГОДУ НАБЛЮДЕНИЙ
1999 DATA 0, 0, 58.5, .8, 0, 0, 0, 0, 0, 0, 0
2000 DATA 1.8, 90, 33.4, .8, 0, 0, 16.3, 0, 0, 16.9, 76.4, 0, 19.6, .5
2100 REM АППРОКСИМАЦИЯ КРИВЫХ СУММ ЖИДКИХ ОСАДКОВ И ГЛУБИН ПРОМЕРЗАНИЯ
2110 FOR I=1 TO N
2120 X(I)=F2(I)
2130 Y(I)=O(I)
2140 NEXT I
2215 PRINT #1, «А, В - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) СУМ. Ж. ОСАД.»
2146 PRINT #1. «Е – ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ»
2150 GOSUB 3000
2160 А1=А9
2170 В1=В9
2180 Е1=Е
2190 FOR I=1 TO N
2200 Y(I)=F(I)
2210 NEXT I
2215 PRINT #1, «А, В - КОЭФФИЦИЕНТЫ УРАВНЕНИЯ КРИВОЙ Y=1/(A*X+B) ГЛУБ. ПРОМ.»
2216 PRINT #1, «E - ТОЧНОСТЬ - УСРЕДНЕННОЕ СРЕДНЕКВАДРАТИЧНОЕ ОТКЛОНЕНИЕ»
2220 GOSUB 3000
2230 А2=А9
2240 В2=В9
2250 Е2=Е
2255 REM ПЕЧАТЬ ЗАГОЛОВКА ТАБЛИЦЫ РЕЗУЛЬТАТОВ
2280 PRINT #1, «Р ЭМП», «СУМ. Ж. ОСАД.», «СУМ. Ж. ОСАД.», «ГОД», «ГЛУБ. ПРОМ.», «ГЛУБ. ПРОМ.», «ГОД»
2265 PRINT #1, « », «НАБЛ.», «АППР.», « », «РАСЧ.», «АППР.»
2270 FOR I=1 ТО N
2280 Y1(I)=1/(A1*F2(I)+B1)
2290 Y2(I)=1/(A2*F2(I)+B2)
2295 REM ПЕЧАТЬ ТАБЛИЦЫ РЕЗУЛЬТАТОВ
2300 PRINT #1, F2(I), O(I), Y1(I), G2(I), F(I), Y2(I), G1(I)
2310 NEXT I
2315 CLOSE #1
2350 STOP
3000 REM ПОДПРОГРАММА АППРОКСИМАЦИИ КРИВОЙ ПО МЕТОДУ НАИМЕНЬШИХ КВАДРАТОВ
3020 S=0
3021 S1=0
3622 S2=0
3023 S3=0
3024 S4=0
3070 FOR I=1 TO N
3075 U(I)=X(I)
3080 V(1)=1/Y(I)
3085 NEXT I
3165 FOR I=1 TO N
3170 S1=S1+U(I)*U(I)
3175 S2=S2+U(I)
3180 S3=S3+U(I)*V(I)
3185 S4=S4+V(I)
3190 NEXT I
3195 D=S1*N-S2*S2
3200 D1=S3*N-S4*S2
3205 D2=S1*S4-S2*S3
3210 A9=D1/D
3211 B9=D2/D
3215 FOR I=1 TO N
3220 V1(I)=A9*U(I)+B9
3225 S=S+(V1(I)-V(I))*(V1(I)-V(I))
3230 NEXT I
3235 E=SQR(S/N)
3255 PRINT #1, «A=»; A9, «B=»; B9,«E=»; E
3260 RETURN
3300 END
ПРИЛОЖЕНИЕ 5 Пример прогнозирования расчетных параметров наледи поверхностных вод на реке В Река находится в первой климатической зоне. Долина реки в месте перехода имеет ширину 840 м, ограниченную пологими склонами гор. Русло сложено из гравия и гальки с песком, также сложены поймы. На глубине 3,9 м залегают глинистые сланцы. Площадь водосбора выше перехода 870 км2, расход воды 246 м3/с с вероятностью превышения 2 %. Уклон реки 0,006. Средняя глубина воды в русле при осеннем уровне 2,7 м.
Побачили розбіжність з офіційним текстом?