ГОСТ 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 |
- расчет () и () по формулам (3) и (4). |
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 |
********** ПРИ РАСЧЕТЕ ЧИСЛО ТОЧЕК СТАЛО МЕНЬШЕ ДОПУСТИМОГО N0=ХХХХ |
Сообщение 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