ГОСТ 9.803-88
Группа Т97
ГОСУДАРСТВЕННЫЙ СТАНДАРТ СОЮЗА ССР
Единая система защиты от коррозии и старения
ФУНГИЦИДЫ
Метод определения эффективности
Unified system of corrosion and ageing protection.
Fungicides. Efficiency determining method
ОКСТУ 0009
Дата введения 1990-01-01
ИНФОРМАЦИОННЫЕ ДАННЫЕ
1. РАЗРАБОТЧИКИ
Х.Н.Фидлер, канд. техн. наук; Ю.В.Моисеев, д-р хим. наук; Э.Г.Африкян, д-р биол. наук; Б.А.Чепенко, канд. техн. наук; З.С.Боголюбова; А.А.Малама, канд. биол. наук (руководители темы); О.А.Хачатурова; В.А.Габечава; Р.Э.Хведелидзе; С.А.Семенов, канд. техн. наук; А.А.Рыжков, канд. хим. наук; К.З.Гумаргалиева, канд. хим. наук; Л.С.Хачатрян, канд. биол. наук; С.Н.Миронова, канд. биол. наук; Т.В.Филимонова
2. УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Постановлением Государственного комитета СССР по стандартам от 28.11.88 N 3857
3. Срок проверки - 1996 г., периодичность проверки - 5 лет
4. ВВЕДЕН ВПЕРВЫЕ
5. ССЫЛОЧНЫЕ НОРМАТИВНО-ТЕХНИЧЕСКИЕ ДОКУМЕНТЫ
|
Обозначение НТД, на который дана ссылка |
Номер раздела, пункта, приложения |
|
1.1; 2.1; Разд.3; 4.2; 4.3; 4.4; 4.5; 4.9; 4.10, разд.7 | |
|
2.1 | |
|
2.1 | |
|
1.5 | |
|
2.1 | |
|
Приложение 3 | |
|
Приложение 1 | |
|
Разд.3 | |
|
Приложение 1 | |
|
Разд.3 | |
|
Приложение 1 | |
|
ГОСТ 24104-88 |
Разд.3 |
|
Разд.3, приложение 1 | |
|
Разд.3 |
6. ПЕРЕИЗДАНИЕ
Настоящий стандарт распространяется на вещества, обладающие фунгицидной и (или) фунгистатической активностью (далее - фунгициды), предназначенные для защиты материалов и изделий от повреждения плесневыми грибами (далее - грибами) и устанавливает метод испытаний эффективности фунгицидов.
Стандарт не распространяется на газообразные фунгициды, а также водонерастворимые фунгициды, разлагающиеся при температуре менее 90 °С.
Сущность метода заключается в культивировании грибов на питательной среде, содержащей фунгициды, и оценке их эффективности по кинетическим параметрам развития грибов.
1. ОТБОР ОБРАЗЦОВ
1.1. В качестве образцов для испытаний водорастворимых фунгицидов применяют их растворы в среде Чапека-Докса. Среду Чапека-Докса готовят по ГОСТ 9.048. Концентрации водорастворимых фунгицидов выбирают из ряда: 2х10
; 4х10
; 8х10
; 2х10
; 4х10
; 8х10
; 2х10
; 8х10
; 2х10
; 4х10
; 2х10
; 4х10
; 8х10
моль/дм
.
1.2. В качестве образцов для испытаний водонерастворимых фунгицидов применяют пленки этих веществ, осажденных на синтезированных гидрогелевых подложках.
Концентрация водонерастворимых фунгицидов на поверхности подложки должна составлять: 1х10
; 2х10
; 5х10
; 8х10
; 1х10
; 2х10
; 5х10
; 8х10
; 1х10
; 2х10
; 5х10
моль/см
.
Концентрация водонерастворимых фунгицидов в растворе органических растворителей, применяемых для осаждения на поверхность подложки, соответственно выбирают из ряда: 1,57х10
; 3,14х10
; 7,85х10
; 10,56х10
; 1,57х10
; 3,14х10
; 7,85х10
; 12,50х10
; 1,57х10
; 3,14х10
; 7,85х10
моль/дм
.
Метод изготовления гидрогелевой подложки с осажденной на ней пленкой водонерастворимого фунгицида приведен в приложении 1.
1.3. Количество исследуемых концентраций каждого фунгицида должно быть не менее 10.
1.4. Контрольными образцами для водорастворимых фунгицидов является среда Чапека-Докса, не содержащая фунгицидов, а для водонерастворимых - гидрогелевая подложка без пленки фунгицидов.
1.5. Количество проб (N) для каждого из отборов каждой из концентраций фунгицидов и для контрольных образцов рассчитывают по ГОСТ 9.707, приложение 3. Если относительная ошибка и вероятность попадания среднего арифметического значения в доверительный интервал не заданы, количество проб на один отбор должно быть не менее семи.
2. ВИДЫ ГРИБОВ
2.1. Виды грибов для испытаний фунгицидов, применяемых для защиты:
полимерных материалов - по ГОСТ 9.049;
лакокрасочных покрытий - по ГОСТ 9.050;
тканей - по ГОСТ 9.802;
других типов материалов и изделий - по ГОСТ 9.048.
3. АППАРАТУРА, МАТЕРИАЛЫ И РЕАКТИВЫ
Аппаратура, материалы и реактивы по ГОСТ 9.048.
Подложки гидрогелевые (см. приложение 1).
Колба с тубусом исполнения 1 номинальной вместимостью 1000 см
по ГОСТ 25336.
Фильтр обеззоленный - синяя лента диаметром 5,5 см.
Воронка Бюхнера N 1 по ГОСТ 9147.
Насос водоструйный по ГОСТ 25336.
Эксикатор исполнения 1 диаметром корпуса 250 мм по ГОСТ 25336.
Лупа четырехкратная по ГОСТ 25706.
Пипетка вместимостью 2 см
.
Пипетка вместимостью 5 см
.
Весы лабораторные общего назначения 2-го класса точности с наибольшим пределом взвешивания 200 г по ГОСТ 24104.
Чашка исполнения 2 типа ЧБН по ГОСТ 25336.
Стаканчик типа СВ с взаимозаменяемым конусом 34/12 по ГОСТ 25336.
Силикагель технический прокаленный по ГОСТ 3956.
Колба коническая по ГОСТ 25336 вместимостью 100 см
из стекла группы ТС.
4. ПОДГОТОВКА К ИСПЫТАНИЯМ
4.1. Для проведения испытаний составляют программу, в которой указывают:
объект испытаний (наименование, марка фунгицида и НТД), завод-изготовитель;
растворимость в воде и других растворителях, режим испытаний;
перечень используемых при испытаниях НТД;
метрологическое обеспечение испытаний;
обозначение настоящего стандарта.
4.2. Посуду, применяемую при испытании, готовят по ГОСТ 9.048.
4.3. Питательные среды для выращивания и хранения культур плесневых грибов готовят по ГОСТ 9.048.
4.4. Пересев, выращивание, хранение культур плесневых грибов и оценка их жизнеспособности - по ГОСТ 9.048.
4.5. Суспензию спор каждого вида грибов в воде готовят по ГОСТ 9.048, при этом концентрация их должна быть 1-2 млн/см
.
4.6. Готовят растворы водорастворимых фунгицидов в среде Чапека-Докса по п.1.1.
4.7. Готовят растворы водонерастворимых фунгицидов в соответствующих органических растворителях квалификации ч. и ч.д.а. по п.1.2. Подготовленные растворы разливают в чашки ЧБН по (5,00±0,05) см
и выпаривают органический растворитель, помещая ЧБН в термостат, полки которого устанавливают горизонтально по уровню. Термостат устанавливают в вытяжном шкафу. Устанавливают температуру в термостате на (5±2) °С ниже температуры кипения соответствующего растворителя и выпаривают до постоянной массы.
4.8. В чашках ЧБН синтезируют гидрогелевые подложки (см. приложение 1).
4.9. Гидрогелевые подложки стерилизуют по ГОСТ 9.048 в чашках ЧБН, устанавливая их в боксы под ртутными или ртутно-кварцевыми лампами на расстоянии (50±10) см и выдерживают (240±6) мин.
4.10. Для испытаний водонерастворимых фунгицидов подготавливают чашки ЧБН, дно которых полностью покрывают отрезками стеклянных трубок и стерилизуют по ГОСТ 9.048.
5. ПРОВЕДЕНИЕ ИСПЫТАНИЙ
5.1. Проведение испытаний водорастворимых фунгицидов
5.1.1. В конические колбы вместимостью 100 см
разливают по 40 см
растворов, приготовленных по п.4.6, и контрольную пробу (см. п.1.4). Количество параллельных проб в соответствии с п.1.5.
5.1.2. Из приготовленных суспензий спор грибов (п.4.5) берут по 2 см
каждого из видов и смешивают между собой.
5.1.3. В конические колбы, подготовленные по п.5.1.1, наливают по 1 см
суспензии спор грибов, подготовленной по п.5.1.2.
5.1.4. Конические колбы, подготовленные по пп.5.1.1-5.1.3, закрывают ватными пробками в стерильных условиях и помещают в термостат при температуре (29±2) °С.
5.1.5. Промежуток времени, в течение которого не наблюдается развитие грибов (далее - лаг-фаза
), при всех концентрациях водорастворимого фунгицида и в контрольных пробах устанавливают по появлению хотя бы на одной пробе биомассы грибов, видимой невооруженным глазом. Периодичность осмотров для выявления лаг-фазы 24 ч.
Если разность минимального и максимального значений лаг-фазы в пробах в любом из отборов превышает 75% ее минимального значения, испытания повторяют на новой серии образцов.
5.1.6. Для количественного определения биомассы при всех концентрациях отборы проводят один раз в 3 сут, а для контрольных проб - один раз в сутки. Общее количество отборов должно быть не менее семи. Продолжительность испытаний - не более 42 сут.
5.1.7. Содержимое колбы сливают на предварительно высушенный до постоянной массы при температуре (105±2) °С и взвешенный обеззоленный фильтр, установленный в приспособлении для вакуумного фильтрования, состоящем из воронки Бюхнера, колбы с тубусом и водоструйного насоса. Биомассу отмывают дистиллированной водой в количестве 100 см
.
5.1.8. После вакуумного фильтрования фильтр с биомассой помещают в предварительно высушенный при температуре (105±2) °С и взвешенный до постоянной массы стаканчик СВ и сушат в термостате с открытой крышкой при температуре (105±2) °С с периодичностью взвешивания 2 ч. Стаканчик с фильтром после каждого из периодов сушки помещают в эксикатор с силикагелем на 1 ч до полного остывания, а затем взвешивают. Сушку прекращают, когда масса стаканчика с фильтром при последующем взвешивании отличается от его массы в предыдущем не более чем на ±0,4 мг.
5.1.9. Если в контрольном опыте по истечении 5 сут не наблюдается развития грибов, споры, использованные для заражения, считают нежизнеспособными.
5.1.10. Вычисляют биомассу каждой пробы каждого отбора каждой из концентраций фунгицида.
5.1.11. Определяют среднее арифметическое значение биомассы (
) в каждом из отборов.
5.1.12. Испытания при каждой из концентраций фунгицидов прекращают, если количество биомассы в последующем отборе не увеличивается по сравнению с предыдущим.
5.1.13. Данные испытаний заносят в протокол, форма которого приведена в приложении 2.
5.2. Проведение испытаний водонерастворимых фунгицидов
5.2.1. Испытания проводят на гидрогелевых подложках с осажденным на них водонерастворимым фунгицидом различных поверхностных концентраций (пп.4.7-4.9).
5.2.2. В чашки ЧБН, подготовленные по п.4.10, наливают питательную среду Чапека-Докса, в них укладывают гидрогелевые подложки (п.5.2.1) так, чтобы поверхность подложки, свободная от фунгицида, контактировала с питательной средой, а на поверхность, содержащую фунгицид, равномерно наносят 1 см
суспензии спор грибов, полученной по п.5.1.2.
5.2.3. Чашки ЧБН накрывают крышками, помещают в камеру или эксикатор, на дно которого налита вода. Эксикатор закрывают крышкой. Испытания проводят при температуре (29±2) °С и относительной влажности более 90%. В камере (эксикаторе) не допускают конденсации влаги, принудительной вентиляции и воздействия прямого естественного или искусственного освещения.
Не реже одного раза в 7 сут эксикаторы или камеры открывают для притока воздуха.
5.2.4. Лаг-фазу при всех концентрациях водонерастворимого фунгицида и в контрольных пробах определяют по времени, в течение которого хотя бы на одной из проб появляется биомасса грибов, видимая при четырехкратном увеличении. Периодичность осмотров - 24 ч.
5.2.5. Для количественного определения биомассы периодичность отборов и их количество - по п.5.1.6. При каждом отборе подложку вынимают из чашки ЧБН, биомассу полностью снимают с поверхности подложки шпателем или лопаткой на предварительно высушенный до постоянной массы при температуре (105±2) °С беззольный фильтр с синей полосой, установленный в приспособлении для вакуумного фильтрования. Проводят вакуумное фильтрование. Биомассу отмывают дистиллированной водой в количестве 100 см
.
5.2.6. В дальнейшем испытания проводят в соответствии с пп.5.1.8- 5.1.13.
6. ОБРАБОТКА РЕЗУЛЬТАТОВ
6.1. Эффективность фунгицидов определяют по зависимости кинетических параметров уравнения, описывающего рост биомассы, от концентрации этих веществ.
6.2. Кинетические параметры (коэффициенты
и
) уравнения, связывающего текущее среднее арифметическое значение удельной биомассы (
) и среднее максимальное (
) значение удельной биомассы при различных концентрациях фунгицида вычисляют по формуле
, (1)
где
- средняя максимальная удельная биомасса, достигнутая при развитии грибов в присутствии фунгицида каждой из концентраций и в контрольной пробе, мг/см
или мг/дм
;
- минимальная, экспериментально полученная лаг-фаза для данной концентрации фунгицида и в контрольной пробе, ч;
- продолжительность выдержки, после которой произведен отбор, ч;
- коэффициент, характеризующий способность споры образовывать биомассу при данной концентрации фунгицида и в контрольной пробе;
- коэффициент, характеризующий удельную скорость развития микроорганизмов на данной питательной среде (при
-й концентрации фунгицида), ч
.
Преобразуем уравнение (1)
(2)
и методом наименьших квадратов вычисляют коэффициенты
и
по формулам:
; (3)
, (4)
где
- порядковый номер отбора (
=1, 2, ...,
).
6.3. Определенные по п.6.2 коэффициенты
и
подставляют в формулу (1) и находят расчетные значения
.
6.4. Рассчитывают средние квадратические отклонения (
) экспериментальных значений
от расчетных по формуле
.
Если в эксперименте имеются значения
, которые не попадают в интервал 
(
), их исключают из данных эксперимента и коэффициенты
и
вновь определяют по п.6.2 для оставшихся точек до тех пор, пока все оставшиеся экспериментальные точки будут попадать в интервал
. Количество их должно быть не менее пяти.
6.5. Оставшиеся для расчета экспериментальные данные биомассы
, а также
и полученные по п.6.2 коэффициенты
и
подставляют в формулу (2) и получают продолжительность лаг-фазы (
) для каждой из концентраций фунгицида и в контрольной пробе (
).
6.6. Определяют константу
, характеризующую зависимость лаг-фазы от концентрации фунгицида, по формуле
, (5)
где
- концентрация фунгицида от
=0 до
;
- минимальная продолжительность лаг-фазы при каждой концентрации фунгицида, ч;
- минимальная продолжительность лаг-фазы в контрольной пробе, ч.
Преобразуют уравнение (5)
(6)
и методом наименьших квадратов рассчитывают константу
по формуле
, (7)
где
- порядковый номер концентрации (
= 1, 2,...,
).
6.7. Используя полученную по п.6.6 константу
, рассчитывают
по формуле (6).
6.8. Рассчитывают среднее квадратическое отклонение (
) полученных по п.6.5 значений
от их расчетных значений.

Если есть вычисленные по п.6.5 значения
, не попадающие в интервал 
(
), их исключают и константу
рассчитывают по оставшимся значениям
до тех пор, пока все оставшиеся значения
будут попадать в интервал
. Количество их должно быть не менее пяти.
6.9. Если экспериментальные данные не удовлетворяют требованиям п.6.8, испытания повторяют.
6.10. По полученным (пп.6.2-6.4) коэффициентам для различных концентраций фунгицида находят константу
, характеризующую зависимость удельной скорости развития грибов от концентрации фунгицида по формуле
, (8)
где
- коэффициент уравнения (1) для концентрации фунгицида (
);
- коэффициент уравнения (1) в контрольной пробе;
- концентрация фунгицида, включая нулевое значение, моль/см
или моль/дм
;
- константа, численно равная концентрации фунгицида, при которой коэффициент
уменьшается в 2 раза по сравнению с развитием гриба на среде, не содержащей фунгицида.
Преобразуют уравнение (8) в линейное
(9)
и методом наименьших квадратов вычисляют значение параметра
по формуле
, (10)
где
=1, 2,...,
- порядковый номер концентрации.
6.11. Чем ниже значение
и выше значение
, тем более эффективным является фунгицид.
6.12. Результаты испытаний могут быть обработаны на ЭВМ в соответствии с программой, приведенной в приложении 3.
6.13. Для предварительной оценки эффективности фунгицида определяют лаг-фазу, используя вместо семи отборов, предусмотренных в пп.5.1.6 и 5.2.5, только три отбора, и рассчитывают
с использованием уравнения (1) и п.6.6.
7. ТРЕБОВАНИЯ БЕЗОПАСНОСТИ
Требования безопасности - по ГОСТ 9.048.
ПРИЛОЖЕНИЕ 1
Обязательное
МЕТОД ИЗГОТОВЛЕНИЯ ГИДРОГЕЛЕВОЙ ПОДЛОЖКИ
Гидрогелевую подложку изготовляют полимеризацией монометакрилового эфира этиленгликоля в присутствии инициатора и сшивающего агента.
1. АППАРАТУРА, МАТЕРИАЛЫ И РЕАКТИВЫ
Термостат, обеспечивающий температуру до 200 °С с погрешностью не более ±2 °С.
Мензурка вместимостью 100 см
по ГОСТ 1770.
Пипетка исполнения 4 1-го класса точности вместимостью 2 см
.
Пипетка исполнения 2 1-го класса точности вместимостью 5, 20, 25 и 50 см
.
Чашка типа ЧКЦ исполнения 1 номинальной вместимостью 2500 см
по ГОСТ 25336.
Часовое стекло диаметром 20 мм.
Спектрофотометр СФ-16.
Чашки типа ЧБН исполнения 2 по ГОСТ 25336.
Стакан типа В исполнения 1 номинальной вместимостью 150 см
из термически стойкого стекла группы ТС по ГОСТ 25336.
Этиленгликольмонометакрилата, ч.
Этиленгликольдиметакрилат, ч.
Бензоила перекись техническая 1-го сорта по ГОСТ 14888.
Вода дистиллированная по ГОСТ 6709.
2. ПОДГОТОВКА К ИЗГОТОВЛЕНИЮ
2.1. Взвешивают часовое стекло с погрешностью не более ±0,5 мг.
2.2. Устанавливают в термошкафу температуру (90±2) °С, а полки его выставляют горизонтально по уровню.
3. ИЗГОТОВЛЕНИЕ ПОДЛОЖЕК
3.1. На часовом стекле взвешивают навеску перекиси бензоила массой (70,0±0,5) мг.
3.2. В стеклянный стакан вместимостью 150 см
наливают 70 см
дистиллированной воды, опускают в нее часовое стекло с перекисью бензоила и перемешивают до полного ее растворения.
3.3. Затем в стакан наливают 30 см
этиленгликольмонометакрилата, содержимое стакана тщательно перемешивают в течение (3,0±0,1) мин, из пипетки доливают (1,50±0,05) см
этиленгликольдиметакрилата и тщательное перемешивание продолжают еще в течение (1,0±0,1) мин.
3.4. В чашки ЧБН (п.4.7 настоящего стандарта) наливают по 40 см
смеси компонентов (пп.3.2, 3.3), устанавливают их на полки термостата при температуре (90±2) °С и выдерживают (120±5) мин. При изготовлении большого количества подложек смесь компонентов из расчета 40 см
на одну подложку берут в соотношениях, указанных в пп.3.1-3.3.
3.5. Вынимают чашки ЧБН из термостата и охлаждают до комнатной температуры. Заполимеризованные подложки опускают в чашки ЧКЦ с дистиллированной водой, взятой в количестве (1500±50) см
.
3.6. Отмывают подложки десятикратно со сменой воды два раза в сутки. Полноту отмывки от мономера, димера и перекиси бензоила контролируют спектрофотометрически при длине волны 220 нм.
Оптическая плотность поглощения последней из промывных вод из чашек ЧКЦ относительно дистиллированной воды не должна превышать (0,12±0,05).
3.7. Отмытые подложки хранят в запаянных полиэтиленовых пакетах в холодильнике при температуре не выше 5 °С.
ПРИЛОЖЕНИЕ 2
Рекомендуемое
ПРОТОКОЛ
испытаний эффективности фунгицида
1. Наименование, марка фунгицида
2. Завод-изготовитель
3. НТД, по которому выпускается
4. Вид растворителя
5. Условия испытаний
|
Концентрация фунгицида |
Номер отбора |
Продол- житель- ность выдер- жки до отбора, ч, сут |
Продолжительность лаг-фазы, ч, сут, для пробы |
Среднее ариф- мети- ческое зна- чение лаг-фазы, ч, сут |
Количество биомассы для пробы |
Среднее арифме- тическое значение удельной биомассы, мг/см | ||||||||||||||||||||
|
водо- раство- римого, моль/дм |
водо- нераство- римого, моль/см |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
||||
|
|
|
|
|
|
|
|
мг |
мг/см |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
ПРИЛОЖЕНИЕ 3
Рекомендуемое
ОПИСАНИЕ ПРОГРАММЫ ОБРАБОТКИ РЕЗУЛЬТАТОВ ИСПЫТАНИЙ ПО ОПРЕДЕЛЕНИЮ ЭФФЕКТИВНОСТИ ФУНГИЦИДОВ
1. Функциональное назначение
Настоящая программа предназначена для реализации алгоритма определения эффективности фунгицидов, изложенного в настоящем стандарте, и составлена в соответствии с требованиями ГОСТ 19.402.
Программа написана на языке Фортран и функционирует под управлением операционной системы ОС ЕС версии 6.1 и выше на ЭВМ ЕС.
Программа состоит из двух подпрограмм.
Первая подпрограмма предназначена для ввода контроля и формирования исходных данных экспериментальных наблюдений.
Вторая подпрограмма, которая вызывается из первой, осуществляет действия согласно алгоритму, описанному в настоящем стандарте.
При одном пуске программы может быть определена эффективность неограниченного числа фунгицидов не более 20 различных концентраций при количестве отборов для каждой концентрации не более 20 и количества проб в одном отборе не более 16.
Программа осуществляет ввод экспериментальных данных и контроль наличия ошибок, допущенных при подготовке их к вводу, а также формирует величины максимального значения биомассы для каждой концентрации фунгицида.
Программа отсеивает отбор, если число проб в этом отборе меньше минимально допустимого (F) количества; кроме того, осуществляется отсев концентрации фунгицида, если количество отборов в ней окажется меньше минимально допустимого (N0) и прекращается обработка результатов по данному фунгициду, если количество концентраций меньше минимально допустимого (N0) .
Все описанные ошибки выдаются на печать в виде диагностических сообщений с указанием номера отбора, номера концентраций и номера фунгицида, что облегчает поиск и исправление ошибки.
По введенным без ошибок данным программа выполняет расчет средней биомассы для каждого отбора, а затем поиск средней максимальной биомассы (MB), достигнутой при развитии грибов в присутствии фунгицида и соответствующего ему интервала времени (Т), считая с начала эксперимента.
Эти данные подаются на вход второй подпрограммы, которая реализует алгоритм определения эффективности фунгицида, приведенный в настоящем стандарте.
Процесс счета сопровождается подробной диагностической печатью, способствующей его анализу.
Данная программа выполняет обработку результатов эксперимента для всех концентраций каждого фунгицида.
Вначале методом наименьших квадратов по формулам (3) и (4) настоящего стандарта для каждой концентрации фунгицида вычисляют значения (
) и (
), затем отсеивают те экспериментальные точки, которые оказались за пределами интервала
(
- среднее квадратическое отклонение экспериментальных данных от расчетных).
Вычисления (
) и (
) повторяются после каждого отсева точек каждой концентрации.
Если в результате отсева значений биомассы какой-либо концентрации фунгицида количество оставшихся отборов станет меньше заданного, то данные этой концентрации фунгицида в дальнейшем не обрабатывают и не учитывают.
Если количество концентраций, для которых удалось рассчитать (
) и (
), равно или превышает заданное, то рассчитывают значение (
) по формуле (7), затем аналогично описанному выше отсеивают концентрации фунгицида, оказавшиеся вне интервала (
) от расчетной.
Если в результате отсева число оставшихся концентраций меньше чем заданное, то прекращается обработка данных по этому фунгициду с выдачей соответствующего сообщения на печать и программа возвращает управление в первую подпрограмму.
Если число оставшихся концентраций фунгицида обеспечивает вычисление (
), значение (
) вычисляют по формуле (10) и печатают значения концентраций фунгицида, лаг-фаз, (
) и (
).
Затем программа переходит к обработке данных следующего фунгицида.
2. Логическая структура вызывающей подпрограммы

3. Описание логической структуры вызывающей программы (первая подпрограмма)
|
А0 |
- пуск программы. |
|
В0 |
- ввод данных, определяющих количество обрабатываемых фунгицидов, минимально допустимое количество отборов (Р0) и минимально допустимое количество концентраций фунгицида (N0). |
|
C0 |
- формирование цикла по количеству фунгицидов (N). |
|
D0 |
- ввод лаг-фазы и количества концентраций для конкретного фунгицида. |
|
Е0 |
- проверка числа введенных концентраций. Если оно больше или равно минимально допустимому количеству, то переход на (Н0), если меньше, то на (F0). |
|
F0 |
- печать сообщения 1. Список выдаваемых сообщений приводится в п.6 настоящего приложения. |
|
G0 |
- прекращение работы программы. |
|
Н0 |
- формирование цикла по количеству концентраций. |
|
I0 |
- ввод параметров конкретной концентрации. |
|
J0 |
- печать параметров конкретной концентрации. |
|
A1 |
- проверка достижения минимально допустимого значения числа отборов. |
|
B1 |
- формирование цикла по количеству отборов данной концентрации. |
|
C1 |
- ввод времени взятия отбора и количества проб в данном отборе. |
|
D1 |
- контроль значения количества проб в данном отборе (не равно нулю). |
|
E1 |
- контроль значения количества проб в данном отборе (не более 16). |
|
G1 |
- печать сообщения 4. |
|
H1 |
- печать сообщения 3. |
|
A2 |
- печать сообщения 2, содержащего значения концентрации фунгицида, лаг-фазы, количества отборов и количества проб в каждом отборе. |
|
B2 |
- формирование цикла по количеству отборов |
|
C2 |
- ввод и распечатка значений биомасс в пробах данного отбора. |
|
D2 |
- проверка окончания цикла по количеству отборов. |
|
E2 |
- наращивание параметра цикла отборов. |
|
F2 |
- ввод и печать значений биомасс в пробах данного отбора и уменьшение количества отборов на единицу. |
|
G2 |
- проверка наличия достижения максимально допустимого количества отборов. Если количество отборов меньше минимально допустимого, то переход на C1, если количество отборов равно или больше минимально допустимого, то переход на 12. |
|
I2 |
- уменьшение количества концентраций на единицу. |
|
J2 |
- проверка достижения минимально допустимого числа концентрации фунгицида. Если число концентраций меньше минимально допустимого количества, то переход к печати сообщений 1 (F0), если число концентраций больше или равно минимально допустимому количеству, то переход на I0. |
|
B3 |
- ввод и печать значений биомассы в пробах данного отбора. |
|
C3 |
- расчет среднего значения биомассы в данном отборе. |
|
D3 |
- проверка окончания цикла по количеству отборов, если число отборов исчерпано, то переход на G3, если не исчерпано, то переход на C1. |
|
E3 |
- увеличение параметра цикла на единицу для перехода к обработке данных следующей серии проб следующего отбора. |
|
G3 |
- поиск и определение предельного значения биомассы, при котором уже не наблюдается ее рост. |
|
H3 |
- проверка достижения минимально допустимого количества отборов данной концентрации. Если число отборов меньше минимально допустимого, то переход на I3, если оно больше или равно минимально допустимому, то переход на А4. |
|
I3 |
- печать сообщения 5. |
|
А4 |
- проверка окончания цикла по числу концентраций данного фунгицида. |
|
B4 |
- увеличение значения параметра цикла по концентрациям на единицу. |
|
D4 |
- вызов второй подпрограммы счета по алгоритму, приведенному в настоящем стандарте. |
|
E4 |
- проверка окончания цикла по числу фунгицидов. |
|
F4 |
- увеличение значения параметра цикла по фунгицидам. |
|
H4 |
- окончание работы программы. |
4. Логическая структура второй подпрограммы

5. Описание логической структуры второй подпрограммы
|
А0 |
- печать сообщения 8. |
|
В0 |
- осуществляется "сжатие" строки массива данных, содержащего средние значения биомасс по отборам данной концентрации фунгицида, а также связанные с этим формирования. |
|
С0 |
- проверка достижения предельно допустимого значения количества отборов для каждой концентрации. |
|
D0 |
- печать сообщения 9. |
|
Е0 |
- уменьшение количества концентраций на единицу. |
|
F0 |
- проверка достижения минимально допустимого количества концентраций. |
|
Если количество концентраций фунгицида меньше, чем минимально допустимое, то переход на G0, если число концентраций больше или равно минимально допустимому, то переход на I0. | |
|
I0 |
- осуществляется сжатие всех столбцов массива, содержащего значения биомасс по концентрациям фунгицида. |
|
А1 |
- формирования, связанные с получением управления из вызывающей программы. |
|
В1 |
- формирование цикла по количеству концентраций исследуемого фунгицида. |
|
С1 |
- печать поступивших на вход программы исходных данных. |
|
D1 |
- формирования для расчета коэффициентов ( |
|
Е1 |
- печать значений исходных данных, используемых для расчета ( |
|
F1 |
- расчет ( |
|
G1 |
- печать вычисленных значений ( |
|
H1 |
- выявление тех отборов данной концентрации фунгицида, значения биомасс которых оказались за пределами интервала |
|
I1 |
- проверка признака, сформированного на предыдущем шаге H1, о необходимости сжатия строки массива, содержащего данные биомасс каждой концентрации. |
|
А2 |
- проверка окончания цикла по количеству концентраций. |
|
В2 |
- увеличение параметра цикла концентраций на единицу. |
|
С2 |
- печать значений концентраций и соответствующих им значений лаг-фаз после окончания цикла по концентрациям. |
|
D2 |
- вычисление методом наименьших квадратов значения |
|
Е2 |
- поиск тех значений концентраций, для которых экспериментальные значения лаг-фаз лежат за пределами интервала |
|
F2 |
- проверка наличия отсеиваемых точек. |
|
G2 |
- вычисление методом наименьших квадратов величины. |
|
Н2 |
- печать результатов работы программы. |
|
I2 |
- возврат управления в первую подпрограмму. |
|
F3 |
- печать сообщения 10. |
|
G3 |
- "Сжатие" массива значений лаг-фаз и уменьшение числа концентраций на единицу. |
|
H3 |
- проверка количества концентраций фунгицида. Если оно больше или равно минимально допустимому числу концентраций, то переход на С3, если оно меньше, то переход на I3. |
|
I3 |
- печать сообщения 6. |
6. СПИСОК ВЫДАВАЕМЫХ СООБЩЕНИЙ
|
Сообщение 1 |
|
***** ДЛЯ ФУНГИЦИДА ХХ РАСЧЕТ НЕВОЗМОЖЕН (N=ХХ МЕНЬШЕ N0=ХХ) |
|
Сообщение 2 |
|
******** ДЛЯ ДАННОЙ КОНЦЕНТРАЦИИ С=ХХХХХХ ОБРАБОТКА НЕ ПРОИЗВОДИТСЯ (Р=ХХХХХХ МЕНЬШЕ Р0=ХХХХХХ) |
|
Сообщение 3 |
|
** N3 ** ДАННЫЕ ДЛЯ Т=ХХХХХХ НЕ ОБРАБАТЫВАЮТСЯ, Т.К. ЧИСЛО ПАРАЛЛЕЛЬНЫХ ПРОБ РАВНО НУЛЮ |
|
Сообщение 4 |
|
** N4 ** ДАННЫЕ ДЛЯ Т=ХХХХХХ НЕ ОБРАБАТЫВАЮТСЯ, Т.К. ЧИСЛО ПАРАЛЛЕЛЬНЫХ ПРОБ БОЛЬШЕ 16 |
|
Сообщение 5 |
|
********************************************************************************************* |
|
Сообщение 6 |
|
********** ПРИ РАСЧЕТЕ |
|
Сообщение 7 |
|
*************************************** |
|
Сообщение 8 |
|
ОТСЕИВАЕТСЯ ТОЧКА ХХ М=ХХ.ХХ Т=ХХХХ |
|
Сообщение 9 |
|
<<<<<<<<<<<<<<<< ОТСЕИВАЕТСЯ КОНЦЕНТРАЦИЯ ХХ >>>>>>>>>>>>>>>>>>>> |
|
Сообщение 10 |
|
********** ПРИ РАСЧЕТЕ |
7. Входные данные
Входные данные, представляющие собой цифровую информацию, должны быть помещены в файл с последовательной организацией.
Первой записью этого файла, определяющего режим работы программы, является запись, содержащая:
(KF) - количество фунгицидов;
(Р0) - минимально допустимое количество отборов при каждой концентрации фунгицида;
(N0) - минимально допустимое количество концентраций каждого из фунгицидов.
В указанной последовательности эти параметры помещаются в формате 3I3 языка Фортран.
Вслед за первой записью следует KF групп записей, каждая из которых содержит данные о действии того или иного фунгицида, включая и контрольную - не содержащую фунгицид.
Первая запись в каждой группе содержит в формате I3 языка Фортран количество концентраций (N) данного фунгицида, включая данные о контрольных пробах, не содержащих фунгицид (
). Далее следует N подгрупп, каждая из которых содержит данные о каждой концентрации.
Первая запись в подгруппе - заголовок подгруппы, содержит сведения о концентрации фунгицида (С), лаг-фазе (LE), количестве отборов (Р) и количестве проб в каждом отборе (FN). Все перечисленные данные в указанной последовательности содержатся в записи в формате 416 языка Фортран (
).
Вслед за заголовком подгруппы располагаются
пар записей, содержащих данные о каждом из отборов. Первая запись в упомянутой паре записей содержит данные о времени (TR), прошедшем с начала эксперимента до взятия отбора и количества проб (FIR) в данном отборе*. Значения TR и FIR представляют в записи в формате 2I6 языка Фортран. Вторая запись содержит данные о биомассах не более чем в 16 пробах каждого из отборов. Значения биомасс в каждом отборе представляются в записи в формате I6F 5.2.
____________________________
* FN распространяются на всю подгруппу, FIR - на конкретный отбор. Если указаны оба, то предпочтение отдается FIR. Нельзя, чтобы FN и FIR отсутствовали. FN и FIR не должны превышать 16.
8. Выходные данные
Выходом программ является распечатка поступающих на обработку данных, а также сообщения, выдаваемые программой в процессе счета. Все выдаваемые программой сообщения и промежуточные результаты счета не нуждаются в дополнительных комментариях.
9. Текст вызывающей (первой) подпрограммы
|
С |
**** |
*********************************************************************************************** | |
|
INTEGER P(20), C(20), LE(20), T(0, 20), N, KF, N0, P0, TR, FIR, FI, FN, F, P1 | |||
|
С С |
|
| |
|
REAL М(20, 20), МВ(20), R(20) | |||
|
С |
ВВОД: KF - КОЛИЧЕСТВО ФУНГИЦИДОВ; | ||
|
С |
Р0 - ПРЕДЕЛЬНО ДОПУСТИМОЕ ЧИСЛО ОТБОРОВ; (МИНИМАЛЬНОЕ) | ||
|
С |
N0 - ПРЕДЕЛЬНО ДОПУСТИМОЕ ЧИСЛО КОНЦЕНТРАЦИЙ; (МИНИМАЛЬНОЕ) | ||
|
С |
|||
|
С |
**** |
************************************************************************************************ | |
|
31 |
FORMAT (1X, 100('*')/'* ДАННЫЕ ДЛЯ ', 13, ' - ОЙ КОНЦЕНТРАЦИИ', | ||
|
*' ОБРАБАТЫВАТЬСЯ НЕ БУДУТ *', /1X, 100('*')) | |||
|
C |
**** |
************************************************************************************************ | |
|
С |
|||
|
READ (5, 1) KF, P0, N0 | |||
|
1 |
FORMAT (313, F7, 5) | ||
|
PRINT 101 | |||
|
101 |
FORMAT ('KF Р0 N0') | ||
|
PRINT 1, KF, P0, N0 | |||
|
С |
|
| |
|
С С |
|
DO 999 LL=1, KF | |
|
C |
ВВОД: N - ЧИСЛО КОНЦЕНТРАЦИЙ ДЛЯ ДАННОГО ФУНГИЦИДА (N<=20). | ||
|
READ (5, 2) N | |||
|
2 |
FORMAT (13) | ||
|
IF (N.GE. N0) GO TО 3 | |||
|
9 |
PRINT 71, LL, N, N0 | ||
|
GO TO 999 | |||
|
710 |
FORMAT ('КОЛИЧЕСТВО КОНЦЕНТРАЦИЙ N=', 12) | ||
|
|
71 |
FORMAT ('**** ДЛЯ ФУНГИЦИДА ', 12, ' РАСЧЕТ НЕ ВОЗМОЖЕН' | |
|
3 |
PRINT 710, N | ||
|
С |
------------------------------------------------------------------------------------------------------------------------ | ||
|
C |
ЦИКЛ ПО ЧИСЛУ КОНЦЕНТРАЦИЙ | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
D0 99 I=7, N | |||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
С |
ВВОД: С - ЗНАЧЕНИЕ КОНЦЕНТРАЦИИ; | ||
|
C |
LE - ЛАГ-ФАЗА ЭКСПЕРИМЕНТАЛЬНАЯ; | ||
|
C |
Р - ЧИСЛО ОТБОРОВ ДАННОЙ КОНЦЕНТРАЦИИ; | ||
|
C |
FN - ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ (НЕ БОЛЕЕ 16); | ||
|
C |
------------------------------------------------------------------------------------------------------------------------- | ||
|
10 |
READ (5, 4) C(I), LE(I), P(I), FN | ||
|
4 |
FORMAT (416) | ||
|
P1=P(I) | |||
|
|
41 |
FORMAT ('********** КОНЦЕНТРАЦИЯ', 14, 10X, | |
|
PRINT 41, I | |||
|
PRINT 42 | |||
|
PRINT 4, C(I), LE(I), P(I), FN | |||
|
С |
|
| |
|
42 |
FORMAT (4X, 'С LE P FN') | ||
|
IF (P(I).GE.P0) GO TO 5 | |||
|
K1=1 | |||
|
21 |
PRINT 72, C(I), P(I), P0 | ||
|
72 |
FORMAT ('******** ДЛЯ ДАННОЙ КОНЦЕНТРАЦИИ С=', 16, | ||
|
*'ОБРАБОТКА НЕ ПРОИЗВОДИТСЯ', | |||
|
*/ 15X, ' (Р=', 16, ' МЕНЬШЕ Р0=', 16, ') ') | |||
|
C |
-----------------------------------------------------------------------------------------------------------------------+ | ||
|
С |
ПРОГОН ДАННЫХ ЭТОЙ КОНЦЕНТРАЦИИ |
! | |
|
C |
-----------------------------------------------------------------------------------------------------------------------+ | ||
|
DO 6 K=K1, P1 | |||
|
READ (5, 7) TR, FIR | |||
|
7 |
FORMAT (216) | ||
|
PRINT 7, TR, FIR | |||
|
READ (5, 8) (R(J), J=1, 16) | |||
|
8 |
FORMAT (16F5.2) | ||
|
PRINT 8, (R(J), J=1, 16) | |||
|
6 |
CONTINUE | ||
|
C |
------------------------------------------------------------------------------------------------------------------------- | ||
|
C |
УМЕНЬШИТЬ ЧИСЛО КОНЦЕНТРАЦИЙ ДАННОГО ФУНГИЦИДА | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
11 |
N=N-1 | ||
|
IF (N.LT.N0) GO TO 9 | |||
|
GO TO 10 | |||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
С |
ЦИКЛ ПО ЧИСЛУ ОТБОРОВ ДАННОЙ КОНЦЕНТРАЦИИ | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
5 |
DO 12 К=1, Р1 | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
С |
ВВОД: Т - ВРЕМЯ С НАЧАЛА РОСТА БИОМАССЫ | ||
|
C |
FI - ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
19 |
READ (5, 7) T(I, K), FI | ||
|
IF(FI.EQ.0) GO TO 13 | |||
|
F=FI | |||
|
GO TO 14 | |||
|
13 |
IF (FN, EQ.0) GO TO 15 | ||
|
F=FN | |||
|
14 |
IF(F.GT.16) GO TO 16 | ||
|
GO TO 17 | |||
|
15 |
PRINT 18, Т (I, К) | ||
|
|
18 |
FORMAT ('**N3** ДАННЫЕ ДЛЯ Т=', 16, ' НЕ ОБРАБАТЫВАЮТСЯ Т.К. ЧИСЛО ПРОБ В ДАННОМ ОТБОРЕ РАВНО НУЛЮ') | |
|
23 |
READ (5, 8) (R(J), J=1, 16 | ||
|
PRINT 8, (R(J), J=1, 16) | |||
|
P(I)=P(I)-1 | |||
|
IF(P(I), GE.P0) GO TO 19 | |||
|
PRINT 20 | |||
|
20 |
FORMAT (1X, 70('-')) | ||
|
K1=K+1 | |||
|
GO TO 21 | |||
|
16 |
PRINT 22, Т (I, К) | ||
|
|
22 |
FORMAT ('**N4** ДАННЫЕ ДЛЯ Т=', 16, ' НЕ ОБРАБАТЫВАЮТСЯ Т.К.' | |
|
GO TO 23 | |||
|
17 |
PRINT 7, Т (I, К), F | ||
|
READ (5, 8) (R(J), J=1, F) | |||
|
PRINT 24, (R(J), J=1, F) | |||
|
24 |
FORMAT ('**', 16 (1X, F5.2)) | ||
|
M (I, К)=0 | |||
|
DO 25 KK=1, F | |||
|
M(I, K)=M(I, K)+R(KK) | |||
|
25 |
CONTINUE | ||
|
M(I, K)=M(I, K)/F | |||
|
12 |
CONTINUE | ||
|
C |
------------------------------------------------------------------------------------------------------------------------ | ||
|
С |
ЗАПОЛНЕНИЕ МАССИВА MB | ||
|
K=2 | |||
|
28 |
IF ((M(I, K)-M(I, К-1)).LT.0) GO TO 26 | ||
|
IF(K.EQ.P(I) GO TO 27 | |||
|
K=K+1 | |||
|
GO TO 28 | |||
|
27 |
MB(I)=M(I, K) | ||
|
P(I)=P(I)-1 | |||
|
GO TO 99 | |||
|
26 |
MB(I)=M(I, K-1) | ||
|
P(I)=К-2 | |||
|
IF (P(I).LE.P0) GO TO 30 | |||
|
99 |
CONTINUE | ||
|
CALL DATVLA (C, LE, MB, M, T, N, P, P0, N0) | |||
|
999 |
CONTINUE | ||
|
STOP | |||
|
30 |
PRINT 31, I | ||
|
GO TO 99 | |||
|
END | |||
10. Текст второй подпрограммы
|
C |
----------------------------------------------------------------------------------------------------------------- | ||||||
|
C |
С - КОНЦЕНТРАЦИЯ |
| |||||
|
C -------------------------------------------------------------------------------------------------------------------- | |||||||
|
SUBROUTINE DATVLA (С, LE, MB, М, Т, N, Р, Р0, N0) |
|||||||
|
INTEGER MKO(20), C(20), LE(20), T(20, 20), Р(20), ММ(20),ТТ(20) |
|||||||
|
INTEGER N.PI.N0, Р0 |
|||||||
|
REAL F(20), M(20, 20), МВ(20), KL, КС, А(20), В(20), LR(20), SAM(20) |
|||||||
|
PRINT 1013 |
|||||||
|
1012 |
FORMAT (2(/1X, 100('*'))) |
||||||
|
1013 |
FORMAT (3(/1X, 100('*'))) |
||||||
|
102 |
FORMAT (6Х, '1', 5Х, '2', 5Х, '3', 5Х, '4', 5Х, '5' |
||||||
|
103 |
FORMAT ('C', 14, 19 (1Х, 15)) |
||||||
|
104 |
FORMAT ('LE', 14, 19 (1Х, 15)) |
||||||
|
105 |
FORMAT ('MB', F4.2, 19 (1Х, F5.2)) |
||||||
|
С |
C ------------------------------------------- * |
||||||
|
PRINT 103, (C(J), J=1, N) |
|||||||
|
PRINT 104, (LE(J), J=1, N) |
|||||||
|
PRINT 105, (MB(J), J=1,N) |
|||||||
|
C |
C -------------------------------------------------------------------------------------------------------------* |
||||||
|
C |
C |
* |
|||||
|
1020 |
FORMAT (2(/1Х, 100('*')), /50Х, 'РАСЧЕТ KL') |
||||||
|
C |
|
||||||
|
C ---------------------------------------------------------------------------* |
|||||||
|
NKO=0 |
* |
||||||
|
NK=N |
|||||||
|
DO 100 I=1, NK |
|||||||
|
PI=P(I) |
|||||||
|
1000 |
PRINT 1012 |
||||||
|
110 |
FORMAT (/IХ, 20('='), 'ДАННЫЕ ДЛЯ КОНЦЕНТРАЦИИ', 12, 20('=')) |
||||||
|
111 |
FORMAT ('M', F4.2, 19(1Х, F5.2)) |
||||||
|
112 |
FORMAT ('T', 14, 19(1Х, 15)) |
||||||
|
113 |
FORMAT ('М - ЭКСП.МАССА, MR - РАСЧ.МАССА, Т - ВРЕМЯ, CKО - ОТКЛОНЕНИЕ') |
||||||
|
C ------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
С |
* |
||||||
|
С |
ПЕЧАТЬ ИСХОДНЫХ ДАННЫХ, УЧАСТВУЮЩИХ В РАСЧЕТАХ А И В * |
||||||
|
С |
* |
||||||
|
C ------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
151 PRINT 110,I |
|||||||
|
PRINT 102 |
|||||||
|
PRINT 112, (T(I, J), J=1, P I) |
|||||||
|
PRINT 111, (M(I, J), J=1, P I) |
|||||||
|
C -------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
С |
* |
||||||
|
C |
РАСЧЕТ КОЭФФИЦИЕНТОВ А И В |
* |
|||||
|
С |
* |
||||||
|
C -------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
SE=0 |
|||||||
|
SC=0 |
|||||||
|
SR=0 |
|||||||
|
SD=0 |
|||||||
|
SL=LE(I) |
|||||||
|
DO=120 J=1, PI |
|||||||
|
ST=T(I, J) |
|||||||
|
F(J)=ALOG ((MB(I)-M(I, J))/M(I, J)) |
|||||||
|
H=ST-SL |
|||||||
|
SE=SE+F(J) |
|||||||
|
SC=SC+H |
|||||||
|
SD=SD+F(J) * H |
|||||||
|
SR=SR+H * H |
|||||||
|
120 |
CONTINUE |
||||||
|
PT=PI |
|||||||
|
B(I)=(SC*SE-PT*SD)/(PT*SR-SC*SC) |
|||||||
|
A(I)=(SE*SR-SC*SD)/(PT*SR-SC*SC) |
|||||||
|
106 |
FORMAT ('+', 50Х, ' LN(A)=', F12.7, ' B=', F12.7, 'СКО=', Е16.10) |
||||||
|
GO TO 3000 |
|||||||
|
3003 |
PRINT106, A(I), B(I), SIGMA |
||||||
|
PRINT 7002 |
|||||||
|
C -------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
C |
ОТСЕВ ТОЧЕК С ОШИБКОЙ, ПРЕВЫШАЮЩЕЙ 2 * SIGMA |
* |
|||||
|
C -------------------------------------------------------------------------------------------------------------------------* |
|||||||
|
П=0 |
|||||||
|
СКО=2*SIGMA |
|||||||
|
DO 130 J=1, PI |
|||||||
|
ST=ABS(M(I, J)-SAM(J)) |
|||||||
|
IF(ST, GT, CKO) GO TO 131 |
|||||||
|
130 |
CONTINUE |
||||||
|
IF (II, EQ, 0) GO TO 200 |
|||||||
|
GO TO 1400 |
|||||||
|
C |
|||||||
|
C |
|||||||
|
C |
|||||||
|
131 |
PRINT 132, J, M(I, J), T(I, J) |
||||||
|
132 |
FORMAT (1Х, 'ОТСЕИВАЕТСЯ ТОЧКА', 12, ' М=', F5, 2, 'Т', 14) |
||||||
|
II=II+1 |
|||||||
|
MM(II)=J |
|||||||
|
GO TO 130 |
|||||||
|
1400 |
DO 2003 I1=1, II |
||||||
|
KR=MM(I1) |
|||||||
|
DO 2001 I2=KR, PI |
|||||||
|
T(I, I2)=T(I, I2+1) |
|||||||
|
M(I, I2)=M(I, I2+1) |
|||||||
|
2001 |
CONTINUE |
||||||
|
PI=PI-1 |
|||||||
|
DO2002 I3=1, II |
|||||||
|
152 |
FORMAT (1Х, 20('<'), 'ОТСЕИВАЕТСЯ КОНЦЕНТРАЦИЯ', 12, 20('>')) |
||||||
|
MM(I3)=MM(I3)-1 |
|||||||
|
2002 |
CONTINUE |
||||||
|
2003 |
CONTINUE |
||||||
|
IF(PI, GE, P0) GO TO 151 |
|||||||
|
PRINT 152, I |
|||||||
|
NKO=NKO+1 |
|||||||
|
MKO (NKO)=I |
|||||||
|
С ------------------------------------------------ |
|||||||
|
С ------------------------------------ |
|||||||
|
С ------------------------------------ |
|||||||
|
100 |
CONTINUE |
||||||
|
GO TO 6000 |
|||||||
|
160 |
PRINT 163 |
||||||
|
163 |
FORMAT (1Х, 50('*'), /1Х, 'ЭКСПЕРИМЕНТ НЕУДАЧНЫЙ', /1Х, 50 ('*')) |
||||||
|
999 |
RETURN |
||||||
|
501 |
PRlNT1020 |
||||||
|
5010 |
PRINT 102 |
||||||
|
PRINT 103, (C(J), J=1, NK) |
|||||||
|
PRINT 105, (MB(J), J=1, NK) |
|||||||
|
PRINT 105, (LE(J), J=1, NK) |
|||||||
|
DO7000 J=1, NK |
|||||||
|
7000 |
PRINT 7001, J, A(J), B(J) |
||||||
|
7001 |
FORMAT ('I=', 13, 'A=', E16, 10, ' В=', E16.10) |
||||||
|
PRINT 7002 |
|||||||
|
7002 |
FORMAT (1Х, 100('-')) |
||||||
|
302 |
FORMAT ('LR', F4.1, 19(1Х, F5.1)) |
||||||
|
С ----------------------------------------------------------------------------------------------------------------------- |
|||||||
|
C |
РАСЧЕТ KL МЕТОДОМ НАИМЕНЬШИХ |
||||||
|
С ----------------------------------------------------------------------------------------------------------------------- |
|||||||
|
КВАДРАТОВ |
|||||||
|
410 |
FORMAT (1Х, 10('='), 'KL F(J)=', F8.5) |
||||||
|
SD=0 |
|||||||
|
SC=0 |
|||||||
|
SL=0 |
|||||||
|
SK=0 |
|||||||
|
DO 400 J=2, NK |
|||||||
|
DJ = LR(J) |
|||||||
|
D1=LR(1) |
|||||||
|
Q=C(J) |
|||||||
|
F(J)=ALOG(DJ/DI) |
|||||||
|
SD=SD+F(J)*Q |
|||||||
|
SC=SC+Q |
|||||||
|
SL=SL+F(J) |
|||||||
|
SK=SK+Q*Q |
|||||||
|
400 |
CONTINUE |
||||||
|
SN=NK |
|||||||
|
KL=(SN*SD-SC*SL)/(SN*SK-SC*SC) |
|||||||
|
C |
|||||||
|
C |
|||||||
|
С ---------------------------------------------------------------------------------------------------------------- |
|||||||
|
C |
ОТСЕВ ТОЧЕК LE ПРИ РАСЧЕТЕ KL |
||||||
|
С ---------------------ПО ------------2*SIGMA -------------------------------------------------------------- |
|||||||
|
SUM=0 | ||
|
DO 401 J=1, NK | ||
|
F(J)=ABS(LR(J)-LR(1)*EXP(KL*C(J))) | ||
|
401 |
SUM=SUM+F(J)*F(J) | |
|
SIGMA=SQRT(SUM/(PI-1)) | ||
|
PRINT 6402, KL, SIGMA | ||
|
6402 |
FORMAT (1X, 20('<'), 'KL=', E16.10,' SIGMA=', E16.10) | |
|
CKO=2*SIGMA | ||
|
II=0 | ||
|
DO 6400 J=1, NK | ||
|
IF (F(J), LE.CKO) GO TO 6400 | ||
|
II=II+1 | ||
|
MM(II)=J | ||
|
SUM=LR(1)*EXP(KL*C(J)) | ||
|
PRINT 404, C(J), SUM | ||
|
6400 |
CONTINUE | |
|
IF (II, EQ, 0) GO TO 403 | ||
|
DO 6405 I1=1, II | ||
|
KR=MM(I1) | ||
|
DO 6406 I2=KR, NK | ||
|
C(I2)=C(I2+1) | ||
|
MB(I2)=MB(I2+1) | ||
|
LE(I2)=LE(I2+1) | ||
|
LR(I2)=LR(I2+1) | ||
|
B(I2)=B(I2+1) | ||
|
6406 |
A(I2)=A(I2+1) | |
|
NK=NK-1 | ||
|
DO 6407 I2=1, II | ||
|
MM(I2)=MM(I2)-1 | ||
|
6407 |
CONTINUE | |
|
6405 |
CONTINUE | |
|
IF (NK, LT, N0) GO TO 500 | ||
|
GO TO 5010 | ||
|
500 |
PRINT 502, N0 | |
|
GO TO 160 | ||
|
404 |
FORMAT (1X, 10('*'), 'ПРИ РАСЧЕТЕ KL ОТСЕЯЛАСЬ КОНЦЕНТРАЦИЯ' | |
|
C |
||
|
502 |
FORMAT (1X, 10('*'), 'ПРИ РАСЧЕТЕ KL ЧИСЛО ТОЧЕК СТАЛО' | |
|
C |
||
|
C |
||
|
C |
||
|
C |
||
|
С ------------------------------------------------------------------------------------------------------------------------ | ||
|
C |
РАСЧЕТ KC МЕТОДОМ НАИМЕНЬШИХ | |
|
С ------------------------------------------------------------------------------------------------------------------------ | ||
|
КВАДРАТОВ | ||
|
403 |
SD=0 | |
|
SC=0 | ||
|
SK=0 | ||
|
SL=0 | ||
|
DO 700 J=2, NK | ||
|
Q=C(J) | ||
|
H=1/Q | ||
|
W=B(J)/B(1)-B(J) | ||
|
SD=SD+H*W | ||
|
SC=SC+H | ||
|
SL=SL+W | ||
|
SK=SK+H*H | ||
|
700 |
CONTINUE | |
|
SN=NK-1 | ||
|
KC=(SN*SD-SC*SL)/(SN*SK-SC*SC) | ||
|
701 |
FORMAT (16X, 25('*'), /16X, '* РАСЧЕТ ПРОИЗВЕДЕН *' | |
|
С ----------------------------------------------------------------------------------------------------------------------- | ||
|
PRINT 701 | ||
|
PRINT 102 | ||
|
PRINT 103, (C(J), J=1, NK) | ||
|
PRINT 302, (LR(J), J=1, NK) | ||
|
AH=ALOG(KL) | ||
|
AB=ALOG(1/KC) | ||
|
PRINT 702, AH, AB | ||
|
702 |
FORMAT (1X, 10('-'), 5X, 'LN(KL)=', F10.5, 5X, 'LN(1/KC)=', F10.5) | |
|
GO TO 999 | ||
|
C ------------------------------------------------------------------------------+ | ||
|
C |
|
! |
|
C |
РАСЧЕТ ЛАГ-ФАЗ |
! |
|
C |
|
! |
|
C ------------------------------------------------------------------------------+ | ||
|
200 |
SL=0 | |
|
DO 201 J=1, PI | ||
|
ST=T(I, J) | ||
|
SD=ST-(A(I)-ALOG((MB(I)-M(I, J))/M(I, J)))/B(I) | ||
|
C |
PRINT 202, SD | |
|
SL=SL+SD | ||
|
201 |
CONTINUE | |
|
202 |
FORMAT (10X, 10('*'), F8.3) | |
|
LR(I)=SL/PI | ||
|
GO TO 100 | ||
|
С --------------------------------------------------------------------------------------------------------------------- | ||
|
C |
||
|
C |
||
|
С --------------------------------------------------------------------------------------------------------------------- | ||
|
3002 |
FORMAT ('MR', F4.2, 19(1X, F5.2)) | |
|
C |
||
|
C |
||
|
С --------------------------------------------------------------------------------------------------------------------- | ||
|
С --------------------------------------------------------------------------------------------------------------------- | ||
|
3000 |
SUM=0 | |
|
DO 3001 J=1, PI | ||
|
SAM(I)=MB(I)/(1+(EXP(A(I)))*(EXP(B(I)*(LE(I)-T(I, J)))) | ||
|
F(J)=ABS(M(I, J)-SAM(J)) | ||
|
3001 |
SUM=SUM+F(J)*F(J) | |
|
PRINT 3002, (SAM(J), J=1, PI) | ||
|
PRINT 113 | ||
|
SIGMA = SQRT(SUM/PI-1)) | ||
|
GO TO 3003 | ||
|
6000 |
IF (NKO.EQ, 0) GO TO 501 | |
|
DO 6001 I1=1, NKO | ||
|
KR=MKO(I1) | ||
|
DO 161 J=KR, NK | ||
|
C(J)=C(J+1) | ||
|
P(J)=P(J+1) | ||
|
LE(J)=LE(J+1) | ||
|
MB(J)=MB(J+1) | ||
|
LR(J)=LR(J+1) | ||
|
DO 161 K=1, 20 | ||
|
M(J, K)=M(J+1, K) | ||
|
161 |
T(J, K)=T(J+1, K) | |
|
NK=NK-1 | ||
|
DO 6001 I3=1, NKO | ||
|
6001 |
MKO(I3)=MKO(I3)-1 | |
|
IF(NK.GE.N0) GO TO 501 | ||
|
GO TO 160 | ||
|
END | ||
11. Вызов и загрузка программы
Вызов и загрузка программы осуществляется с помощью следующего задания, составленного на языке управления заданиями
//ВВВ JOB MSGLEYEL=(1, 1)
// ЕXЕС FORTGOLG, REGION, FORT=240K, PARM, FORT=(MAP, ID),
// REGION, GO=240K
//FORT, SYSIN DO VOL=SER=XXXXXX, UNIT=5061, DISP=SHR, DSN=ARL(FUNGICID)
//GO, SYSIN DD VOL=SER=MUKMUK, UNIT=5061, DISP=SHR, DSN=DANN
//.
Перед запуском этого задания текст программы должен быть помещен в раздел FUNGICID библиотеки ARL, находящейся на томе прямого доступа с именем XXXXXX, а данные, составленные по правилам, описанным в разд.7 настоящего приложения, должны быть помещены в последовательный набор данных под именем DANN, находящийся на томе прямого доступа с именем MUKMUK.
По этому заданию операционная система осуществит трансляцию, редактирование и выполнение описанной программы, выдав на печать результаты счета.
Если приведенное задание будет помещено в раздел ААА библиотеки ССС, находящейся на томе прямого доступа EEEEEE, то для запуска этого задания нужно набрать на консоли следующий оператор: S_RDRD, 5061, ЕЕЕЕЕЕ, DSN=ССС(ААА).
Текст документа сверен по:
официальное издание
Защита от коррозии.
Часть 6. Защита от биоповреждений. -
М.: Издательство стандартов, 1994