Главная // Актуальные документы // Методические указания
СПРАВКА
Источник публикации
М.: Издательство стандартов, 1986
Примечание к документу
Название документа
"МИ 641-84. Методические указания. Расчет значений критериев качества поверки средств измерений методами программного моделирования"
(утв. протоколом НТС ВНИИМС от 30.11.1984 N 16)

"МИ 641-84. Методические указания. Расчет значений критериев качества поверки средств измерений методами программного моделирования"
(утв. протоколом НТС ВНИИМС от 30.11.1984 N 16)


Содержание


Утверждены
Научно-техническим советом
Всесоюзного научно-исследовательского
института метрологической службы
30 ноября 1984 г.
(протокол N 16)
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
РАСЧЕТ ЗНАЧЕНИЙ КРИТЕРИЕВ КАЧЕСТВА ПОВЕРКИ СРЕДСТВ
ИЗМЕРЕНИЙ МЕТОДАМИ ПРОГРАММНОГО МОДЕЛИРОВАНИЯ
МИ 641-84
РАЗРАБОТАНЫ Всесоюзным научно-исследовательским институтом метрологической службы (ВНИИМС)
ИСПОЛНИТЕЛИ: Н.Н. Вострокнутов, канд. техн. наук (руководитель темы), Б.А. Френкель (отв. исполнитель).
УТВЕРЖДЕНЫ Научно-техническим советом Всесоюзного научно-исследовательского института метрологической службы (ВНИИМС) 30 ноября 1984 г. (протокол N 16)
Настоящие методические указания предназначены для использования при разработке и экспертизе нормативно-технических документов (НТД), удовлетворяющих требованиям ГОСТ 8.042-83 "ГСИ.НТД на методики поверки средств измерений" и разработке стандартов на поверочные схемы (ГОСТ 8.061-80)".
ИС МЕГАНОРМ: примечание.
МИ 187-79 утратили силу в связи с введением в действие МИ 187-86 "Методические указания. Государственная система обеспечения единства измерений. Достоверность и требования к методикам поверки средств измерений", утв. ВНИИМС 25.07.1986.
ИС МЕГАНОРМ: примечание.
МИ 188-79 утратили силу в связи с введением в действие МИ 188-86 "Методические указания. Государственная система обеспечения единства измерений. ГСИ. Средства измерений. Установление значений параметров методик поверки", утв. ВНИИМС 25.07.1986.
Методические указания устанавливают методы и программы расчета на универсальных цифровых ЭВМ значений критериев качества поверки, регламентированных МИ 187-79, МИ 188-79 (приложение 1) в заданных точках диапазона измерений без предварительного установления аналитических зависимостей между параметрами методики поверки и значениями критериев качества.
1. ОБЩИЕ УКАЗАНИЯ
1.1. В соответствии с настоящими методическими указаниями производят расчет значений критериев качества проверки следующих характеристик поверяемого СИ: основной погрешности; систематической составляющей основной погрешности; среднего квадратического отклонения (с.к.о.) случайной составляющей основной погрешности; вариации.
1.2. Расчет значений критериев качества поверки СИ производят на ЭВМ серии ЕС под управлением операционной системы ОС версии 4 0 или более совершенной.
1.3. Текст комплекса программ расчета значений критериев качества поверки СИ (ниже - комплекса программ) на алгоритмическом языке Фортран-IV приведен в приложении 2.
Объем оперативной памяти, занимаемой комплексом программ, не превышает 80 кБайт.
Для решения задачи требуется один накопитель на магнитных дисках типа ЕС 5050 или ЕС 5061 для хранения загрузочного модуля комплекса программ.
1.4. Комплекс программ производит расчет: значений критериев качества поверки; максимальной вероятности необнаруженного брака Pнм; максимального выхода за допуск ; средней вероятности фиктивного брака Pф; максимальной вероятности фиктивного брака Pфм; значений оперативной характеристики (приложение 1) исследуемой методики поверки.
1.5. Комплекс программ обеспечивает вывод на печать следующих сообщений: общей информации пользователя о задании (текст); комментариев пользователя к методике поверки (текст); таблицы исходных данных и значений критериев качества поверки; значений оперативной характеристики; графика оперативной характеристики; контрольно-диагностических сообщений и комментариев комплекса программ (приложение 3).
Общая информация пользователя о задании, таблица значений критериев качества поверки и контрольно-диагностические сообщения выводятся на печать всегда. Необходимость печати остальных сообщений должна быть указана особо.
1.6. Работа пользователей с комплексом программ производится в три этапа: этап 1 - подготовка исходных данных; этап 2 - программирование задания; этап 3 - решение задачи на ЭВМ.
Работы по этапу 1 выполняет специалист-метролог, принимающий участие в разработке или экспертизе НТД, указанных во вводной части методических указаний. Указания по выполнению работ, относящихся к этапу 1, даны в разд. 2.
Работы по этапам 2 и 3 выполняет специалист по программированию на языке Фортран (программист). Указания по выполнению работ, относящихся к этапам 2 и 3, даны соответственно в разд. 3 и 4.
Возможно выполнение всех этапов работы специалистом-метрологом, обладающим достаточной подготовкой в программировании на языке Фортран.
В приложении 4 приведены развернутые примеры выполнения работ по всем этапам и результаты расчета значений критериев качества поверки для рассмотренных примеров.
2. ПОДГОТОВКА ИСХОДНЫХ ДАННЫХ
2.1. Для проведения работы с комплексом программ следует подготовить исходные данные, перечисленные в пп. 2.1.1 - 2.1.9.
2.1.1. Формула или алгоритм, описывающие связь показаний (выходного сигнала) Y поверяемого СИ со значением измеряемой величины X и составляющими основной погрешности СИ:
, (2.1)
где - систематическая составляющая погрешности СИ; - случайная составляющая погрешности СИ; - составляющая погрешности, обусловленная гистерезисом; - составляющие погрешности СИ, обусловленные отклонениями влияющих величин и неинформативных параметров входного сигнала от значений, соответствующих нормальным условиям.
Вывод формулы или алгоритма (2.1) производится на основании анализа принципа действия СИ конкретного типа. Выражение (2.1) записывают в общем виде (без подстановки значений параметров).
Составляющие , , погрешности поверяемого СИ представляют в (2.1) указанными ниже способами в зависимости от того, для какой из характеристик погрешности рассчитываются значения критериев качества поверки.
2.1.1.1. При расчете значений критериев качества проверки погрешности поверяемого СИ в формуле (2.1): представляют в виде простой переменной, обозначаемой символом ; - в виде выражения:
, (2.2)
где q - номинальная ступень квантования поверяемого цифрового прибора или аналого-цифрового преобразователя (для аналоговых приборов или преобразователей принимают q = 0); - с.к.о. случайной составляющей погрешности поверяемого СИ; H - вариация поверяемого СИ; - представляют в виде выражения:
, (2.3)
где - безразмерная случайная величина, распределенная по трапецеидальному закону с.к.о., равным 1 (см. рисунок, а).
Примечания:
1. Случайная величина задается в комплексе программ с помощью подпрограмм - генераторов последовательностей случайных чисел, распределенных по различным законам распределения.
2. Предусмотрена возможность наращивания подпрограмм - генераторов случайных чисел.
Законы распределения последовательностей случайных чисел
а
б
Значение в (2.2) и (2.3) принимают равным: пределу , если он нормирован для поверяемого СИ; максимальному значению с.к.о. , которое может иметь место у исправных экземпляров СИ данного типа, если предел не нормирован.
представляют в виде выражения
, (2.4)
где RH = 1, если входной сигнал при поверке регулируется в сторону его увеличения; RH = -1, если входной сигнал регулируется в сторону его уменьшения; D = 1, если петля гистерезиса характеристики преобразования поверяемого СИ имеет направление по часовой стрелке; D = -1 - в противном случае.
Значения H в (2.2) и (2.4) принимают равными: пределу Hop вариации погрешности поверяемого СИ, если он нормирован; максимальному значению Homax составляющей , которое может иметь место у исправных экземпляров СИ данного типа, если Hop не нормирован, а составляющая может иметь место по принципу действия поверяемого СИ.
2.1.1.2. При расчете значений критериев качества проверки систематической составляющей погрешности в формуле (2.1): исключают из формулы; представляют в виде простой переменной, обозначаемой символом ; представляют в виде выражения (2.3); представляют в виде выражения (2.4).
Значения и H выбирают в соответствии с указаниями п. 2.1.1.1.
2.1.1.3. При расчете значений критериев качества проверки с.к.о. случайной составляющей погрешности СИ в формуле (2.1): исключают из формулы; принимают равной нулю, если имеется уверенность в том, что не влияет на погрешность оценивания (обычно это условие выполняется, если поверяемое СИ - аналоговый прибор или преобразователь), в противном случае принимают , где - такое значение , при котором следует ожидать наибольшие погрешности оценивания (зависимость этой оценки от часто имеет место у цифровых СИ); представляют в виде выражения (2.3), где - простая переменная; представляют в виде выражения (2.4). Значение H выбирают в соответствии с указаниями п. 2.1.1.1.
2.1.1.4. При расчете критериев качества контроля вариации поверяемого СИ в (2.1): исключают из формулы; представляют в соответствии с указаниями п. 2.1.1.3; представляют в виде выражения (2.3); значение выбирают в соответствии с указаниями п. 2.1.1.1; представляют в виде выражения (2.4), где H - простая переменная.
2.1.1.5. Составляющие включают в (2.1) только в тех случаях, когда возникает необходимость исследования влияния соответствующих величин на значения критериев качества поверки или имеются основания ожидать заметное влияние какой-либо конкретной влияющей величины на погрешность поверки.
2.1.2. Формула или алгоритм, описывающие связь fe показаний (выходного сигнала) Ye образцового средства измерений с номинальным значением Nx воспроизводимой величины для мер или значением Xe входной величины для измерительных преобразователей и приборов и составляющими погрешности образцового СИ:
для мер
; (2.5)
для приборов и преобразователей
, (2.5а)
где - систематическая составляющая погрешности образцового СИ; - случайная составляющая погрешности образцового СИ; и представляют в (2.5) указанными ниже способами в зависимости от того, какие характеристики погрешности нормированы для образцового СИ конкретного типа.
2.1.2.1. Нормирован предел допускаемых значений погрешности образцового СИ. Принимают ; , представляют в (2.5) или (2.5а) в виде
, (2.6)
где - безразмерная случайная величина, имеющая закон распределения равной плотности (см. рисунок б) в пределах от -1 до 1 для ансамбля образцовых СИ данного типа и принимающая неизменное значение для конкретного экземпляра образцового СИ.
2.1.2.2. Нормированы предел систематической составляющей и предел с.к.о. случайной составляющей погрешности образцового СИ. представляют в (2.5) в виде выражения (2.6); представляют в (2.5) в виде:
, (2.7)
где - безразмерная случайная величина, распределенная по закону (см. рисунок, а) с с.к.о., равным 1.
2.1.2.3. Нормированы математическое ожидание (м.о.) и с.к.о. систематической составляющей погрешности для образцовых СИ данного типа, и предел допускаемых значений с.к.о. случайной составляющей погрешности образцового СИ. представляют в (2.6) в виде выражения (2.7); представляют в (2.6) в виде:
, (2.8)
где - безразмерная случайная величина, распределенная по закону равной плотности с с.к.о., равным 1, но принимающая неизменное значение для конкретного экземпляра образцового СИ.
2.1.3. Описание алгоритма проведения поверки. Алгоритм выполнения поверки описывают путем перечисления каждого действия, выполняемого при поверке, в необходимой последовательности.
Все эти действия, выполняемые при проведении измерений и обработке их результатов, должны быть записаны последовательно даже в том случае, когда в действительности они выполняются одновременно.
Формулы для вычислений (программы вычислений) и логические операции принимаются при моделировании без измерений по сравнению с проводимыми в исследуемой методике поверки. Рекомендуется указывать число знаков, с которым производятся вычисления и правила округления результатов вычислений.
По описанию операций и перечню расчетных формул разрабатывают блок-схему алгоритма поверки по ГОСТ 19.002-80.
Примеры разработки описания алгоритма поверки приведены в приложении 4.
Примечание. В одном задании для комплекса программ допускается одновременно рассчитывать значения критериев качества поверки, выполняемой по разным (до девяти) алгоритмам. В этом случае разрабатывают модели каждого алгоритма и присваивают им условный номер Nм от 1 до 9.
2.1.4. Числовые значения перечисленных ниже параметров:
а) предел Aд допускаемых значений контролируемой характеристики погрешности (абсолютное значение); б) отношение контрольного допуска Aк к Aд; в) граница фиктивного бракования (см. приложение 1); г) объем выборки Nр для оценки вероятностей признания годным (для расчета значений оперативной характеристики по методу Монте-Карло); д) условные номера Nм алгоритмов поверки (см. примечание к п. 2.1.3).
Рекомендуемое значение .
Значение Nр выбирают таким, чтобы обеспечить выполнение неравенства
, (2.9)
где - половина интервала допускаемого отклонения оценки Pнм от ее истинного значения (в % к ожидаемому значению Pнм); - квантиль интеграла вероятности для доверительной вероятности .
Рекомендуется принимать Nр > 100.
2.1.5. Число вариантов исходных данных.
Вариантом исходных данных называют совокупность исходных данных по пп. 2.1.1 - 2.1.4. В одном задании для комплекса программ возможен расчет значений критериев качества поверки СИ для различных (не более девяти) вариантов исходных данных.
2.1.6. Текст общей информации о задании - текст (не более 1000 знаков) произвольного содержания, отражающий специфику решаемой задачи. Указанный текст предназначен для распечатки в заголовке листа результатов расчета и служит для облегчения последующей классификации и хранения распечаток. Пример текста общей информации о задании дан в приложении 4.
2.1.7. Комментарии к методике поверки (приводить не обязательно) - текст (не более 1000 знаков) произвольного содержания, отражающий специфику исследуемой методики поверки. Назначение комментариев - то же, что у текста общей информации о задании (п. 2.1.6).
Программисту должно быть указано, что помимо общей информации по п. 2.1.6 должны быть распечатаны комментарии к методике поверки.
2.1.8. Требуемый объем выходных данных: требуется ли распечатка таблицы исходных данных (для проверки правильности их кодирования); требуется ли распечатка значений оперативных характеристик методики поверки; требуется ли построение графиков оперативных характеристик методики поверки.
Примечание. Рекомендуется указывать программисту требование распечатки значений и графиков оперативных характеристик методики поверки в следующих случаях:
а) при расчете критериев качества поверки, который проводится с целью сопоставления разрабатываемой и известных методик. В этом случае графики оперативных характеристик являются полезными для принятия решения о введении (изменении) контрольного допуска, числа отсчетов при контроле или других характеристик методики поверки;
б) если при расчете критериев качества поверки (без распечатки значений и графиков оперативных характеристик) получено сообщение комплекса программ (приложение 3), свидетельствующее о возникновении ошибки расчета. В этом случае для отыскания причин возникновения ошибки следует повторить расчет, указав требование распечатать значения и графики оперативных характеристик.
2.1.9. Условный номер задания (для последующей классификации и облегчения хранения распечаток).
2.2. Оформление исходных данных по п. 2.1.
2.2.1. Обозначения и наименования характеристик, входящих в (2.1), (2.6), и алгоритма поверки, принимающих числовые значения, записывают в таблицу вида (табл. 1).
Таблица 1
Обозначение
Наименование характеристики
Классификационный признак
Имя
1
2
3
4
В графу 1 табл. 1 записывают условное обозначение характеристики, входящей в (2.1) и (2.6); в графу 2 - наименование характеристики; в графу 3 - классификационный признак характеристики. Каждой характеристике присваивают один из следующих признаков: "контролируемый параметр" - контролируемой погрешности или характеристике ее составляющей поверяемого СИ; "контрольный допуск" - величине, с которой сравнивается при контроле оценка контролируемой характеристики погрешности поверяемого СИ; "параметр целого типа" - детерминированным характеристикам, принимающим неслучайные и только целые значения; "параметр вещественного типа" - характеристикам, способным принимать не только целые, но и дробные или только дробные неслучайные значения; "случайная величина" - характеристикам, принимающим случайное, причем при каждом измерении новое значение; "случайная величина по ансамблю СИ" - характеристикам, принимающим случайное, однако для одного экземпляра СИ неизменяемое значение; "промежуточная переменная" - характеристикам, являющимся функциями других.
Графа 4 табл. 1 заполняется программистом (см. п. 3.3.3).
2.2.2. Наименования и числовые значения параметров вещественного и целого типов записывают в раздельные таблицы параметров вещественного и целого типов, составленные по форме таблицы (табл. 2).
Таблица 2
NN п/п
Наименование
Номер варианта
1
2
3
1
2
3
4
5
В графу 2 таблицы параметров вещественного типа сначала записывают сокращенные (до 40 знаков) наименования параметров по п. 2.1.4 а, б, и в в порядке их перечисления в п. 2.1.4, а в графы 3, 4, 5 ... - числовые значения этих параметров для каждого варианта исходных данных. Затем в графу 2 таблицы параметров вещественного типа записывают сокращенные наименования характеристик из табл. 1, имеющие признак "параметр вещественного типа", а в графы 3, 4, 5 ... - числовые значения этих характеристик для каждого варианта исходных данных.
В графу 2 таблицы параметров целого типа сначала записывают сокращенные наименования параметров по п. 2.1.4 г и д в порядке их перечисления в п. 2.1.4, а в графы 3, 4, 5 ... - числовые значения этих параметров для каждого варианта исходных данных. Затем в графу 2 таблицы целого типа записывают сокращенные наименования характеристик из табл. 1, имеющие признак "параметр целого типа", а в графы 3, 4, 5 ... - числовые значения этих характеристик для каждого варианта исходных данных.
2.3. Формулы, указанные в пп. 2.1.1 - 2.1.3, блок-схемы алгоритмов по п. 2.1.4, таблицы по п. 2.2, тексты и указания по п. 2.1.5 - 2.1.9 передают программисту для программирования задания и решения задачи на ЭВМ.
3. ПРОГРАММИРОВАНИЕ ЗАДАНИЯ
3.1. Программирование задания включает: формирование пакета исходных данных; программирование исследуемой методики поверки СИ (составление программной модели методики поверки); подготовку пакета инструкций операционной системе ЭВМ для управления выполнением задания (п. 3.4).
3.1.1. Формирование пакета исходных данных осуществляется на основе материалов, разработанных в соответствии с разд. 2 и включает: кодирование управляющей карты; подготовку блоков текстовой информации и числовых исходных данных.
Правила формирования пакета исходных данных изложены в п. 3.2.
3.1.2. Программирование методики поверки СИ заключается в написании программы методики поверки (ПМП) СИ по алгоритму, разработанному по указаниям разд. 2, и сопряжении ПМП с комплексом программ расчета значений критериев качества поверки СИ (ниже - комплексом программ). Правила программирования методики поверки СИ изложены в п. 3.3.
3.2. Формирование пакета исходных данных.
3.2.1. В пакет исходных данных к заданию включают: управляющую карту (УК); блоки текстовой информации; блоки числовых исходных данных.
3.2.2. УК предназначена для сообщения комплексу программ указаний пользователя по пп. 2.1.6 и 2.1.9. В УК задают следующие коды:
а) условный пятизначный номер задания. На УК он кодируется в позициях с 1 по 5;
б) код типа решаемой задачи (цифры 02) в позициях 6 и 7;
в) число вариантов исходных данных. Допустимое число вариантов исходных данных в одном задании от 1 до 9. Код занимает позицию 9 УК;
г) код объема распечатываемой текстовой информации указывается в позиции 13 УК. Применяется следующая кодировка: 0 - печать текстовой информации в полном объеме; распечатываются: текст общей информации о задании по п. 2.1.7, комментарии к методике поверки по п. 2.1.8 и комментарии комплекса программ; 1 - печать только текстовой информации по пп. 2.1.7 и 2.1.8; 2 - печать только общей информации о задании по п. 2.1.7;
д) код объема распечатываемых таблиц указывается в позиции 15 УК. Применяется следующая кодировка: 0 - печать таблиц в полном объеме: 1 - не печатаются таблицы значений оперативных характеристик; 2 - не печатаются таблицы значений исходных данных и оперативных характеристик;
е) код распечатки графиков оперативной характеристики указывается в позиции 17 УК. Применяется следующая кодировка: 0 - графики распечатываются; 1 - графики не распечатываются.
3.2.3. Состав, назначение и правила кодирования блоков текстовой информации и числовых исходных данных (БТЧИ). Состав БТЧИ: блок 1 - общая информация о задании п. 2.1.6; блок 2 - числовые значения параметров вещественного типа и их наименования по п. 2.2.2; блок 3 - числовые значения параметров целого типа и их наименования по п. 2.2.2; блок 4 - комментарии к методике поверки по п. 2.1.7 (необязательный блок).
БТЧИ кодируют на стандартных перфокартах; каждый из них может состоять из произвольного числа перфокарт.
Началом блока считается отдельная перфокарта, в позиции 1 которой закодирован номер блока. Последняя карта каждого блока должна содержать символ "9" в позиции 1.
Текст блоков 1 и 4 должен занимать на каждой карте позиции с 2 по 69. Блок 4 составляется только при задании кодов 0 или 1 объема распечатываемой текстовой информации (п. 3.2.2).
В блок 2 включают наименования и числовые значения параметров, записанные в таблице параметров вещественного типа, по п. 2.2.2. На каждый параметр отводится две карты: карта наименования и карта числовых значений.
На карте наименования в первых 40 позициях в виде алфавитных символов кодируется наименование параметра.
На карте числовых значений кодируют значения указанного параметра, число которых соответствует числу вариантов исходных данных. На каждое значение отводится поле в пять позиций, начиная со второй. Присутствие десятичной точки в числовых значениях переменных обязательно. Нулевое значение переменной может кодироваться в виде пробелов в соответствующем поле. Последняя карта блока должна содержать символ "9" в позиции 1.
В блок 3 включают наименования и числовые значения параметров целого типа, записанные в таблице по п. 2.2.2. Правила формирования и кодирования блока 3 те же, что для блока 2 за исключением требования о наличии десятичной точки в кодах чисел. Формат значений переменных - I5.
3.3. Программирование методики поверки СИ.
3.3.1. Программу методики поверки (ПМП) СИ составляют на стандартном алгоритмическом языке Фортран-IV.
3.3.2. Написание ПМП производят в три этапа:
I - присвоение имен характеристикам, записанным в табл. 1 (п. 3.3.3);
II - описание наименования ПМП (п. 3.3.4) и имен характеристик табл. 1 (пп. 3.3.5 - 3.3.8), выходной переменной результата контроля (п. 3.3.9) и аргумента программы генератора случайных чисел (п. 3.3.10);
III - программирование алгоритма поверки (пп. 3.3.12 и 3.3.13).
3.3.3. В графу 4 табл. 1 записывают имена характеристик СИ и алгоритма поверки. Имена характеристик должны состоять из букв латинского алфавита (число букв в имени - от 1 до 6). Имена характеристик целого типа должны начинаться с одной из букв: I, J, K, L, M, N, а имена характеристик вещественного типа - с любой из остальных букв алфавита.
3.3.4. Первой инструкцией ПМП должен быть оператор наименования вида: SUBROUTINE MET "I", где I = 1, 2, ..., 9 - номер методики поверки СИ. Пример. SUBROUTINE MET1.
Если в задании требуется рассчитать значения критериев качества для нескольких методик поверки СИ, соответствующие подпрограммы должны иметь те же номера, что и номер схемы алгоритма поверки по п. 2.1.4.
3.3.5. Имена характеристик из табл. 1, отмеченных признаком "параметр вещественного типа" должны быть в ПМП описаны с помощью оператора вида: COMMON/REAL/ "список имен". Пример. COMMON/REAL/ETAL, SIGMA, POINT.
3.3.6. Имена характеристик из табл. 1, отмеченных признаком "параметр целого типа", должны быть описаны в ПМП оператором COMMON/INTEG/ "список имен".
3.3.7. Порядок следования имен в списках должен соответствовать порядку следования наименований характеристик в табл. 1 и таблицах по п. 2.2.2.
3.3.8. Имена характеристик из табл. 1 с признаками "контролируемый параметр" и "контрольный допуск" в ПМП должны быть описаны с помощью оператора: COMMON/CONTRL/ "имя контролируемого параметра", "имя контрольного допуска".
Пример. Контролируется с.к.о. случайной составляющей погрешности СИ (имя SIGMA1) на соответствие контрольному допуску (имя SCONTR). Оператор описания этих параметров должен иметь вид: COMMON/CONTRL/SIGMA1, SCONTR.
Примечание. Текущие значения контролируемого параметра задаются комплексом программ в соответствии с алгоритмом решаемой задачи.
3.3.9. В тексте ПМП должны присутствовать операторы принятия решения по результатам контроля вида:
GODEN = TRUE. - СИ признано годным;
GODEN = FALSE. - СИ забраковано.
Выходная переменная результата контроля должна иметь имя GODEN и подлежит описанию в ПМП с помощью операторов:
LOGICAL GODEN
COMMON/LOGIC/GODEN
3.3.10. Аргумент программы ГСЧ, предназначенной для моделирования случайных составляющих погрешности СИ, подлежит описанию в ПМП вида: COMMON/RND/IX. Начальное значение аргумента IX, равное 137, задается комплексом программ.
3.3.11. Характеристики в табл. 1 с признаками "случайная величина" и "промежуточная переменная" описания в ПМП не требуют, за исключением промежуточных переменных, представляющих собой массивы данных, описание которых следует производить по правилам алгоритмического языка Фортран-IV.
3.3.12. Программирование алгоритма поверки производят по блок-схеме алгоритма, разработанной по указаниям п. 2.1.3.
Операции установления выходного сигнала образцовой меры (установление проверяемой точки), определения показания образцового (поверяемого) СИ программируют путем записи на языке Фортран выражения (2.6) и (2.1) с использованием в качестве характеристик имен из табл. 1. Все текущие значения характеристик, входящих в выражения (2.1) и (2.6), должны быть определены. Для этого предварительно должны быть присвоены значения всем характеристикам, обозначенным в табл. 1 признаками "случайная величина", "случайная величина по ансамблю СИ" и "промежуточная переменная". Правила задания значений случайным величинам даны в п. 3.3.13.
Операции, включаемые в цикл измерений, рекомендуется оформлять в виде цикла DO.
До начала участка ПМП, соответствующего циклу измерений, должны быть присвоены значения всем характеристикам типа "случайная величина по ансамблю СИ". Характеристикам типа "случайная величина" значения присваивают внутри указанного выше участка.
Допускается использование дополнительных переменных для хранения результатов промежуточных вычислений. Выбор и описание имен этих переменных производят по правилам языка Фортран-IV.
Пример программирования алгоритма поверки дан в приложении 4.
3.3.13. Программирование характеристик, приведенных в табл. 1 с признаками "случайная величина", производят с применением программ генераторов случайных чисел (ГСЧ) (приложение 2).
Обращение к программе ГСЧ производят по выражению вида: R - "имя ГСЧ" (IX), где R - имя величины, получающей случайное значение с законом распределения, соответствующим имени программы ГСЧ; IX - аргумент программы ГСЧ (п. 3.3.10).
В табл. 3 указаны принятые имена программ ГСЧ и их основные характеристики. При пользовании данными табл. 3 следует иметь в виду, что в приведенной в приложении 2 версии комплекса программ использована расширенная библиотека подпрограмм генераторов случайных чисел.
Таблица 3
Условный номер программы
Имя ГСЧ
Закон распределения
Пределы распределения
С.к.о.
1
RAND
Равномерный
от -1 до 1
0,58
2
SIMP
Симпсона
от -1 до 1
0,41
3
TRAP
Трапецеидальный
от -1 до 1
0,45
4
GAUSS
Нормальный
-
1
5
AMOD1
Антимодальный 1
от -1 до 1
0,71
6
AMOD2
Антимодальный 2
от -1 до 1
0,79
7
RANDS
Равномерный
от -1,73 до 1,73
1
8
SIMPS
Симпсона
от -2,44 до 2,44
1
9
TRAPS
Трапецеидальный
от -2,19 до 2,19
1
10
AMODS1
Антимодальный 1
от -1,41 до 1,41
1
11
AMODS2
Антимодальный 2
от -1,27 до 1,27
1
Для реализации изложенного метода расчета достаточно иметь ГСЧ NN 1, 3, 7, 9.
Косвенное обращение к программе ГСЧ производят по выражению вида: R = RASP (Nд), где Nд - условный номер программы ГСЧ из табл. 3; R - имя величины, получающей случайное значение с законом распределения с номером Nд.
Примечание. Условный номер программы ГСЧ может быть записан в таблицу по п. 2.2.2 и введен в перфокарт в числе других параметров целого типа.
3.4. Для запуска комплекса программ под управлением ОС ЕС следует подготовить пакет задания, включающий операторы на языке управления заданиями.
Пакет задания должен иметь вид:
//OPTIVER JOB
// EXEC FORTGCLG,
// PARM. LKED = LET,
// COND.GO = (8, LT, LKED)
//FORT. SYSIN DD*
Пакет подпрограмм методик поверки СИ
//LKED. ADDLIB DD DSN = OPTLIB, DISP-SHR
//LKED. SYSIN DD*
INCLUDE ADDLIB (OPTIVER)
ENTRY MAIN
//GO. SYSIN DD*
УПРАВЛЯЮЩАЯ КАРТА
БЛОКИ ТЕКСТОВОЙ ИНФОРМАЦИИ И ЧИСЛОВЫХ ИСХОДНЫХ ДАННЫХ
ПУСТАЯ КАРТА
//
Примечания:
1. Указаны только основные операторы и их операнды.
2. OPTLIB - имя библиотеки загрузочных модулей, разделом которых является комплекс программ;
3. OPTIVER - имя раздела OPTLIB, соответствующего комплексу программ.
4. РЕШЕНИЕ ЗАДАЧИ НА ЭВМ
4.1. Решение задачи на ЭВМ производится в пакетном режиме. Вмешательство оператора ЭВМ в ход выполнения задания не требуется.
4.2. Результаты работы комплекса программ распечатываются на стандартном АЦПУ. Описание сообщений комплекса программ дано в приложении 3. Там же приведено описание возможных причин появления контрольно-диагностических сообщений.
Приложение 1
Справочное
КРИТЕРИИ
КАЧЕСТВА ПОВЕРКИ СИ И ОПИСАНИЕ АЛГОРИТМА ИХ РАСЧЕТА
МЕТОДОМ ПРОГРАММНОГО МОДЕЛИРОВАНИЯ
ИС МЕГАНОРМ: примечание.
МИ 187-79 утратили силу в связи с введением в действие МИ 187-86 "Методические указания. Государственная система обеспечения единства измерений. Достоверность и требования к методикам поверки средств измерений", утв. ВНИИМС 25.07.1986.
1.1. Применяют следующие критерии для оценки качества поверки СИ (по МИ 187-79): наибольшую вероятность принятия любого негодного экземпляра СИ в качестве годного (максимальная вероятность необнаруженного брака Pнм); отношение наибольшего возможного значения Aм характеристики погрешности СИ, признанного по результатам поверки годным, но в действительности негодного, к пределу Aд ее допускаемых значений (максимальный выход за допуск ); отношение числа годных, но забракованных СИ, к числу всех в действительности годных СИ (средняя вероятность фиктивного брака Pф); максимальная вероятность фиктивного брака Pфм.
Применяемые критерии качества поверки СИ иллюстрирует рис. 1, где x = |A/Aд| - отношение действительного значения контролируемой характеристики A погрешности СИ к Aд; P(x) - вероятность признания поверяемого СИ годным для заданного x (оперативная характеристика); - вспомогательный параметр - верхняя граница зоны значений x, в которой бракование СИ считается фиктивным (граница фиктивного бракования); .
Рис. 1
1.2. Для расчета значений критериев качества поверки СИ применяется комплекс программ в составе: основных программ - MAIN, QUALMT, PROB; вспомогательных программ - PFB, PRINOP, PROBIN, ERPROC, TXOUT, RPROUT, IPROUT, PRIM1, QUER, PRRES, OPCHR, GRAF; программ генераторов случайных чисел.
Алгоритмы основных программ комплекса даны в п. 1.4.
1.3. Назначение основных и вспомогательных программ комплекса.
1.3.1. Задачами программы MAIN (главная программа) являются: организация ввода и размещения в памяти исходных данных к заданию; построение последовательности вызова других программ; организация диагностики и вывода результатов расчета.
1.3.2. Подпрограмма QUALMT предназначена для вычисления значений критериев качества поверки СИ путем расчета оперативной характеристики исследуемой методики поверки.
1.3.3. Подпрограмма PROB предназначена для оценки условной вероятности признания СИ годным при поверке при условии, что контролируемая характеристика погрешности поверяемого СИ имеет конкретное значение.
1.3.4. Назначение вспомогательных программ отражено в комментариях к операциям вызова этих программ из основных, помеченным соответствующими именами в алгоритмах п. 1.4.
1.4. Алгоритмы основных программ комплекса.
На рис. 2 - 4 представлены схемы алгоритмов основных программ комплекса - MAIN, QUALMT, PROB
Программа MAIN
Лист 1
Лист 2
Лист 3
Лист 4
Рис. 2
Подпрограмма QUALMT
Рис. 3
Подпрограмма PROB
Рис. 4
1.5. Показатели точности расчета критериев качества поверки по программе QUALMT.
Абсолютная погрешность расчета критериев качества поверки Pнм и Pфм определяется по формуле , где - оценка Pнм или Pфм.
Абсолютная погрешность расчета Pф определяется по формуле , где - оценка Pф.
За оценку принимается минимальное из множества значений контролируемой характеристики погрешности поверяемого СИ, при которых вероятность признания последнего годным не превышает 0,01.
Приложение 2
C
КОМПЛЕКС ПРОГРАММ РАСЧЕТА КРИТЕРИЕВ КАЧЕСТВА ПОВЕРКИ
C
СРЕДСТВ ИЗМЕРЕНИЙ
C
C
РАЗРАБОТЧИКИ
Н.Н. ВОСТРОКНУТОВ, Б.А. ФРЕНКЕЛЬ (ВНИИМС)
C
КОМПЛЕКСА
C
C
C
ГЛАВНАЯ ПРОГРАММА
C
C
ВЫЗЫВАЕМЫЕ ПОДПРОГРАММЫ
C
C
ОСНОВНЫЕ
ВСПОМОГАТЕЛЬНЫЕ
МЕТОДИКИ ПОВЕРКИ
C
C
QUALMT
PRINOP
MET1
C
PROBIN
ERPROC
MET2
C
PRRES
TXOUT
MET3
C
RPROUT
MET4
C
IPROUT
MET5
C
DISPET
MET6
C
QUER
MET7
C
OPCHR
MET8
C
GRAF
MET9
EXTERNAL MET1, MET2, MET3, MET4, MET5, MET6, MET7, MET8,
*MET9
DIMENSION TEXT (3000), RES (10, 9), OPX (101, 9), NBEG (4),
*LTX(4)
DIMENSION RPARM (20, 9), IPARM (20, 9)
DIMENSION CROUT (3)
COMMON/REAL/RP (50)/REALP/RPP (50)
COMMON/INTEG/IP (20)/INTP/IPP (20)
COMMON/RND/IX
COMMON/ARG/X (101)
ITER = 1
IX = 137
C
C
ЧТЕНИЕ УПРАВЛЯЮЩЕЙ КАРТЫ
C
READ 11, NJOB, NPR, NEXP, NOUS, MTEXT, MTABL, MGRAF,
*MOD
PRINT 7
PRINT 9
PRINT 10
PRINT 9
PRINT 7
PRINT 12, NJOB
5
PRINT 6, NEXP
IND = 4
8
PRINT 7
CALL PRINOP (MTEXT, MTABL MGRAF, NOUS, NPR)
PRINT 7
C
C
ВВОД ИСХОДНЫХ ДАННЫХ И ОБРАБОТКА ОШИБОК ВВОДА
C
CALL PROBIN (RPARM, LR, IPARM, LI, TEXT, NBEG, LTX, NEXP,
*IER)
CALL ERPROC (IER, MTEXT)
PRINT 7
C
C
ВЫВОД ТАБЛИЦЫ ИСХОДНЫХ ДАННЫХ
C
PRINT 15
CALL TXOUT (TEXT, NBEG (1), LTX (1), MOD)
IF (MTEXT.NE.O) GO TO 40
PRINT 7
PRINT 16
CALL TXOUT (TEXT, NBEG (4), LTX (4), MOD)
40
IF (MTABL.EQ.2) GO TO 41
PRINT 7
PRINT 17
CALL RPROUT (TEXT, NBEG (2), LTX(2), RPARM, LR, NEXP,
*MOD)
CALL IPROUT (TEXT, NBEG (3), LTX (3), IPARM, LI, NEXP, MOD)
41
IF (MTEXT.NE.2) CALL PRIM1
PRINT 7
DO 50 I = 1, 101
50
X (1) = 0.02* (I - 1)
C
C
НАЧАЛО ЦИКЛА РЕШЕНИЯ ЗАДАЧ
C
DO 60 N = 1, NEXP
C
C
ПЕРЕДАЧА ПАРАМЕТРОВ МЕТОДИКЕ ПОВЕРКИ
C
DO 55 L = 1, LR
55
RPP (L) = RPARM (L, N)
DO 56 L = 1, LI
56
IPP (L) = IPARM (L, N)
DO 57 L = IND, LR
57
RP (L - IND + 1) = RPARM (L, N)
DO 58 L = 3, LI
58
IP (L - 2) = IPARM (L, N)
C
C
ОБРАЩЕНИЕ, ЕСЛИ НЕОБХОДИМО, К ДОПОЛНИТЕЛЬНОЙ ПРОГРАММЕ ПОЛЬЗОВАТЕЛЯ
C
И ПЕРЕХОД К РЕШЕНИЮ ЗАДАЧИ
C
IF (N.LE.NOUS) CALL DISPET (N)
31
METOD = IPP (2)
PRINT 19
GO TO (61, 62, 63, 64, 65, 66, 67, 68, 69) METOD
61
CALL QUALMT (MET1, X, OPX (1, N), RES (1, N), IER)
GO TO 70
62
CALL QUALMT (MET 2, X, OPX, (1, N), RES (1, N), IER)
GO TO 70
63
CALL QUALMT (MET3, X, OPX (1, N), RES (1, N), IER)
GO TO 70
64
CALL QUALMT (MET4, X, OPX, (1, N), RES (1, N), IER)
GO TO 70
65
CALL QUALMT (MET5, X, OPX (1, N), RES (1, N), IER)
GO TO 70
66
CALL QUALMT (MET6, X, OPX (1, N), RES (1, N), IER)
GO TO 70
67
CALL QUALMT (MET7, X, OPX (1, N), RES (1, N), IER)
GO TO 70
68
CALL QUALMT (MET8, X, OPX (1, N), RES (1, N), IER)
GO TO 70
69
CALL QUALMT (MET9, X, OPX (1, N), RES (1, N), IER)
70
PRINT 20, N, IER, (RES (I, N), I = 1, 4)
IF (IER.EQ.J) GO TO 60
CALL QUER (IER)
60
CONTINUE
PRINT 7
IF (NPR.LT.10.AND.NPR.NE.2) GO TO 34
PRINT 9
PRINT 7
C
C
ВЫВОД РЕЗУЛЬТАТОВ ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ СИ
C
PRINT 26
CALL PRRES (RES, NEXP, MOD)
PRINT 7
PRINT 9
PRINT 7
IF (MTABL.NE.0) GO TO 81
PRINT 27
CALL OPCHR (X, OPX, NEXP, 1)
PRINT 7
CALL OPCHR (X, OPX, NEXP, 51)
PRINT 7
81
IF (MGRAF. EQ. 0) CALL GRAF (NEXP, OPX, X, RPARM (3, 1),
RPP (2))
PRINT 7
34
PRINT 28
PRINT 9
STOP
C
C
СПИСОК ФОРМАТОВ
C
6
FORMAT (T25, 'ОЦЕНКА КАЧЕСТВА ПОВЕРКИ СИ'/Т25, 'ЧИСЛО ОЦЕНОК', 12)
7
FORMAT (1HO)
9
FORMAT (T10, 100 ('*'))
10
FORMAT (T10, '*** КОМПЛЕКС ПРОГРАММ РАСЧЕТА КРИТЕРИЕВ',
*
'КАЧЕСТВА ПОВЕРКИ СРЕДСТВ ИЗМЕРЕНИЙ *ВЕРСИЯ 4/1982/***')
11
FORMAT (15, 712)
12
FORMAT (T10, 'ЗАДАНИЕ', 16)
15
FORMAT (T20, 'ОБЩАЯ ИНФОРМАЦИЯ О ЗАДАНИИ')
16
FORMAT (T20, 'КОММЕНТАРИИ К МЕТОДИКАМ')
17
FORMAT (T20, 'ТАБЛИЦА ИСХОДНЫХ ДАННЫХ')
19
FORMAT (T10, 'ДИАГНОСТИКА ОБРАЩЕНИЯ К П/П ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ')
20
FORMAT (T10, 15, 'КОД ОШИБКИ - ',12, 4X, 'PHM =',F5.3, 3X, 'ДМ =', F5.3,
*
3X, 'РФ =', F5.3, 3X, 'РФМ =', F5.3)
26
FORMAT (T20, 'ТАБЛИЦА РЕЗУЛЬТАТОВ')
27
FORMAT (T20, 'ТАБЛИЦА ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ')
28
FORMAT (T20, 'ЗАДАНИЕ ВЫПОЛНЕНО')
END
SUBROUTINE QUALMT (MET, X, PX, RES, IER)
C
C
C
ПОДПРОГРАММА - QUALMT
C
C
C
C
C
НАЗНАЧЕНИЕ
C
C
ПРОГРАММА ПРОИЗВОДИТ ОЦЕНКУ КАЧЕСТВА МЕТОДИКИ ПОВЕРКИ СИ,
C
ЗАДАВАЕМОЙ ПОЛЬЗОВАТЕЛЕМ, ПО КРИТЕРИЯМ МИ 187-79.
C
ИСПОЛЬЗУЕТСЯ МЕТОД СТАТИСТИЧЕСКИХ ИСПЫТАНИЙ
C
ИМИТАЦИОННОЙ МОДЕЛИ ПОВЕРКИ.
C
C
C
ВЫЗЫВАЕМЫЕ ПОДПРОГРАММЫ
C
C
PROB - П/П ОЦЕНКИ ВЕРОЯТНОСТИ ПРИЗНАНИЯ СИ ГОДНЫМ
C
PFB - ИНТЕГРИРОВАНИЕ МЕТОДОМ ТРАПЕЦИЙ
C
C
C
C
ОПИСАНИЕ ФОРМАЛЬНЫХ ПАРАМЕТРОВ
C
MET
- ИМЯ ПРОГРАММЫ МЕТОДИКИ ПОВЕРКИ
C
X
- МАССИВ АРГУМЕНТОВ ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
C
PX
- МАССИВ ЗНАЧЕНИЙ ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
C
RES
- НАБОР ЗНАЧЕНИЙ КРИТЕРИЕВ КАЧЕСТВА ПОВЕРКИ СИ
C
IER
- КОД ОШИБКИ РЕШЕНИЯ
C
C
C
C
ОПИСАНИЕ ПЕРЕМЕННЫХ В БЛОКАХ COMMON
C
C
EDOP
- ПРЕДЕЛ ДОП. ЗНАЧЕНИЙ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
C
GAMMA
- КОНТРОЛЬНЫЙ ДОПУСК
C
BETA
- ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
C
NPOV
- ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТИ СОБЫТИЙ
C
C
DIMENSION X (101), PX (101), RES (10)
COMMON/REALP/EDOP, GAMMA, BETA
COMMON/INTP/NPOV
COMMON/CONTRL/ERSP, ERC
DO 35 I = 1, 101
35
PX (I) = 0.
ERC = GAMMA*EDOP
IER = 0
DO 8 I = 1, 101
ERSP = X (I)*EDOP
PX (I) = PROB (NPOV, MET)
IF (PX(I).LE..01) GO TO 7
NST = 0
IF (X(I).GE.1.01.0R.X(I).LE..99) GO TO 6
RES (1) = PX (I)
IF (RES(1).LT..99) GO TO 8
IER = 4
GO TO 8
6
IF (PX(I).GE..99) GO TO 8
IF (I.EQ.1) IER = 1
GO TO 8
7
NST= NST + 1
IF (NST.EQ.3) GO TO 4
8
CONTINUE
IER = 3
GO TO 5
4
RES (2) = X (I - 2)
IF (RES(2).GT.1) GO TO 5
RES (1) = 0.
IER = 2
5
DELTA = X(2) - X(1)
CALL PFB (BETA, DELTA, PX, RES(3), RES (4))
IF (RES(4).GE..99) IER = 5
RETURN
END
SUBROUTINE QUER (IER)
C
C
C
ПОДПРОГРАММЫ - QUER, PRRES, OPCHR, GRAF
C
C
C
C
C
НАЗНАЧЕНИЕ
C
QUER:
ОБРАБОТКА ОШИБОК ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ
C
PRRES:
РАСПЕЧАТКА ТАБЛИЦЫ РЕЗУЛЬТАТОВ ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ
C
GRAF:
ПОСТРОЕНИЕ СВОДНОГО ГРАФИКА ОПЕРАТИВНЫХ
ХАРАКТЕРИСТИК
C
OPCHR:
РАСПЕЧАТКА СВОДНОЙ ТАБЛИЦЫ ЗНАЧЕНИЙ ОХ
C
C
C
C
C
ОПИСАНИЕ ФОРМАЛЬНЫХ ПАРАМЕТРОВ
C
C
IER
- КОД ОШИБКИ РЕШЕНИЯ
C
RES
- НАБОР ЗНАЧЕНИЙ КРИТЕРИЕВ КАЧЕСТВА ПОВЕРКИ СИ
C
NEXP
- ЧИСЛО СТОЛБЦОВ ТАБЛИЦЫ
C
MOD
- КОД ИНТЕРВАЛА РАСПЕЧАТКИ
C
X
- МАССИВ АРГУМЕНТОВ ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
C
OPX, F
- МАССИВ ЗНАЧЕНИЙ ОХ
C
NBEG
- НОМЕР НАЧАЛЬНОЙ СТРОКИ РАСПЕЧАТКИ
C
NGRAF
- ЧИСЛО ВЫВОДИМЫХ ГРАФИКОВ
C
B, G
- ХАРАКТЕРНЫЕ ТОЧКИ (ПОМЕЧАЮТСЯ НА ГРАФИКЕ)
C
C
GO TO (71, 72, 73, 74, 75), IER
71
PRINT 21
RETURN
72
PRINT 22
RETURN
73
PRINT 23
RETURN
74
PRINT 24
RETURN
75
PRINT 25
RETURN
21
FORMAT (T12,'* БРАК ПРИ МИНИМАЛЬНОМ ЗНАЧЕНИИ КОНТР. ПАРАМЕТРА')
22
FORMAT (T12,'* НУЛЕВОЙ НЕОБНАРУЖЕННЫЙ БРАК')
23
FORMAT (T12,'* МАКС. ВЫХОД ЗА ДОПУСК > 2')
24
FORMAT (T12,'* ВЕРОЯТНОСТЬ Н. Б. РАВНА 1')
25
FORMAT (T12,'* ВЕРОЯТНОСТЬ Ф. Б. РАВНА 1')
END
SUBROUTINE PRRES (RES, NEXP, MOD)
DIMENSION RES (10, 9)
1
FORMAT (1H)
2
FORMAT (T10, 'ПОКАЗАТЕЛЬ КАЧЕСТВА', T40, 9I8/)
3
FORMAT (T2, 'МАКС. ВЕРОЯТНОСТЬ НЕОБНАРУЖ. БРАКА', T40, 9F8.3)
4
FORMAT (T2, 'МАКС. ВЫХОД ЗА ДОПУСК', T40, 9F8.3)
5
FORMAT (T2, 'СРЕДН. ВЕРОЯТНОСТЬ ФИКТ. БРАКА', T40, 9F8.3)
6
FORMAT (T2, 'МАКС. ВЕРОЯТНОСТЬ ФИКТ. БРАКА', T40, 9F8.3)
PRINT 1
PRINT 2, (I, I = 1, NEXP)
PRINT 1
IF (MOD.EQ.l) PRINT 1
PRINT 3, (RES (1, I), I = 1, NEXP)
IF (MOD.EQ.l) PRINT 1
PRINT 4, (RES (2, I), I = 1, NEXP)
IF (MOD.EQ.l) PRINT 1
PRINT 5, (RES (3, I), I = 1, NEXP)
IF (MOD.EQ.l) PRINT 1
PRINT 6, (RES (4, I), I = 1, NEXP)
RETURN
END
SUBROUTINE OPCHR (X, OPX, NEXP, NBEG)
DIMENSION X (101), OPX (101, 9)
1
FORMAT (1H)
2
FORMAT (T10, 'КОНТРОЛИРУЕМЫЙ ПАРАМЕТР', T40, 9I8)
3
FORMAT (T2, I10, F10, 2, T40, 9F8.3)
PRINT 1
PRINT 2, (I, I = 1, NEXP)
PRINT 1
DO 10 I = 1, 50
K = I + NBEG - 1
PRINT 3, K, X (K), (OPX (K, N), N = 1, NEXP)
10
CONTINUE
RETURN
END
SUBROUTINE GRAF (NGRAF, F, X, B, G)
DIMENSION F (101, NGRAF), X (101), FOR (104), MS (101)
DATA PRB/4HF5.1/, BL1/3H, 1X/, BL5/3H, 5X/, TRK/4H, '-'/,
*
CLM/4H, ':'/, DT/3H, I1/, FOR (I)/4H(T8,/, FOR (104)/1H)/,AXOR/4H, 'I'/
10
FORMAT (T4, 11 (9X, 'I'))
11
FORMAT (T6, 11 (7X, F3.1))
12
FORMAT (/T13, 'КОНТРОЛИРУЕМЫЙ ПАРАМЕТР (В ЕД. ПРЕДЕЛА ДОП. ЗНАЧЕНИЙ)')
13
FORMAT (///T25, 'ГРАФИК ОПЕРАТИВНЫХ ХАРАКТЕРИСТИК'//)
C
ВЫЧИСЛЕНИЕ КООРДИНАТ ХАРАКТЕРНЫХ ТОЧЕК
KB = INT (B*50. + .1) + 3
KG = INT (G*50. + .1) + 3
C
ПЕЧАТЬ ЗАГОЛОВКА
PRINT 13
C
ЦИКЛ ФОРМИРОВАНИЯ И ПЕЧАТИ ФОРМАТА СТРОКИ
DO 1 NROW = 1, 101
NPX = 101 - NROW
C
ПОСТРОЕНИЕ ОСИ ОРДИНАТ
FOR (3) = AXOR
C
ПОИСК СТРОК, КРАТНЫХ 10
DO 2 N = 1, 11
IF (NROW.EQ.10*(N - 1) + 1) GO TO 3
2
CONTINUE
C
ПОСТРОЕНИЕ НЕКРАТНЫХ СТРОК
FOR (2) = BL5
ASSIGN 14 TO LABEL
ASSIGN 16 TO METKA
C
ОЧИСТКА СТРОКИ
DO 6 I = 4, 103
6
FOR (I) = BL1
GO TO 5
C
ПОСТРОЕНИЕ КРАТНЫХ СТРОК
3
P = .01*NPX
FOR (2) = PRB
ASSIGN 15 TO LABEL
ASSIGN 17 TO METKA
DO 20 I = 4, 103
20
FOR (I) = TRK
C
ПОСТРОЕНИЕ ХАРАКТЕРНЫХ ТОЧЕК
5
FOR (KB) = CLM
FOR (KG) = CLM
FOR (53) = CLM
C
ЦИКЛ ЗАПОЛНЕНИЯ ПОЗИЦИЙ СТРОКИ
NS = 1
DO 7 I = 1, 101
K = I + 2
C
ЦИКЛ ПОИСКА ТОЧКИ ГРАФИКА
DO 8 N = 1, NGRAF
IF (NPX.NE.INT (F (I, N)*100. + .5)) GO TO 8
IF (FOR (K).EQ.DT) GO TO 9
FOR (K) = DT
MS (NS) = N
GO TO 8
C
ТОЧКА ПЕРЕСЕЧЕНИЯ ГРАФИКОВ
9
MS (NS) = 0
GO TO 18
8
CONTINUE
IF (FOR (K).NE.DT) GO TO 7
18
NS = NS + 1
7
CONTINUE
C
ПЕЧАТЬ, ЕСЛИ ТОЧКА ПУСТАЯ
IF (NS.EQ.l) GO TO METKA, (16, 17)
NS = NS - 1
GO TO LABEL, (14, 15)
C
ПЕЧАТЬ НЕПУСТОЙ СТРОКИ
14
PRINT FOR, (MS(I), I = 1, NS)
GO TO 1
C
ПЕЧАТЬ НЕПУСТОЙ КРАТНОЙ СТРОКИ
15
PRINT FOR, P, (MS(I), I = 1, NS)
GO TO 1
C
ПЕЧАТЬ ПУСТОЙ СТРОКИ
16
PRINT FOR
GO TO 1
C
ПЕЧАТЬ ПУСТОЙ КРАТНОЙ СТРОКИ
17
PRINT FOR, P
1
CONTINUE
C
ПЕЧАТЬ ОСИ АБСЦИСС
PRINT 10
PRINT 11, (X(I), I = 1, 101, 10)
PRINT 12
RETURN
END
FUNCTION RASP (N)
C
C
C
КОМПЛЕКС ПРОГРАММ - ГЕНЕРАТОРОВ СЛУЧАЙНЫХ ЧИСЕЛ
C
C
C
C
C
НАЗНАЧЕНИЕ ПРОГРАММ
C
C
ГЕНЕРАТОРЫ СЛУЧАЙНЫХ ЧИСЕЛ (ГСЧ) ПРЕДНАЗНАЧЕНЫ ДЛЯ
C
МОДЕЛИРОВАНИЯ СЛУЧАЙНЫХ СОСТАВЛЯЮЩИХ ПОГРЕШНОСТЕЙ СРЕДСТВ ИЗМЕРЕНИЙ
C
RASP (N) - КОСВЕННОЕ ОБРАЩЕНИЕ К ГСЧ (ПО УСЛОВНОМУ НОМЕРУ)
C
C
НОМЕР
ОБОЗН.
РАСПРЕДЕЛЕНИЕ
ХАРАКТЕРИСТИКИ
C
C
1.
RAND
РАВНОМЕРНОЕ
РАЗМАХ ОТ -1 до 1
C
2.
SIMP
СИМПСОНА
ТО ЖЕ
C
3.
TRAP
ТРАПЕЦЕИДАЛЬНОЕ
ТО ЖЕ
C
4.
GAUSS
ГАУССА
С.К.О. РАВНО 1
C
5.
AMOD1
АНТИМОДАЛЬНЫЙ 1
РАЗМАХ ОТ -1 ДО 1
C
6.
AMOD2
АНТИМОДАЛЬНЫЙ 2
ТО ЖЕ
C
7.
RANDS
РАВНОМЕРНОЕ
С.К.О. РАВНО 1
C
8.
SIMPS
СИМПСОНА
ТО ЖЕ
C
9.
TRAPS
ТРАПЕЦЕИДАЛЬНОЕ
ТО ЖЕ
C
10.
AMODS1
АНТИМОДАЛЬНЫЙ 1
ТО ЖЕ
C
11.
AMODS2
АНТИМОДАЛЬНЫЙ 2
ТО ЖЕ
C
C
C
COMMON/RND/IX
IF (N) 20, 20, 30
20
RASP = 0.
RETURN
30
GO TO (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11), N
1
RASP = RAND (IX)
RETURN
2
RASP = SIMP (IX)
RETURN
3
RASP = TRAP (IX)
RETURN
4
RASP = GAUSS (IX)
RETURN
5
RASP = AMOD1 (IX)
RETURN
6
RASP = AMOD2 (IX)
RETURN
7
RASP = RANDS (IX)
RETURN
8
RASP = SIMPS (IX)
RETURN
9
RASP = TRAPS (IX)
RETURN
10
RASP = AMODS1 (IX)
RETURN
11
RASP = AMODS2 (IX)
RETURN
END
FUNCTION RAND (IX)
IY = IX*65539
IF (IY) 1, 1, 2
1
IY = IY + 2147483647 + 1
2
Y = IY
Y = Y*.4656613E - 9
IX = IY
RAND = 2.*Y - 1.
RETURN
END
FUNCTION SIMP (IX)
T = RAND (IX)/2.
SIMP = T + RAND (IX)/2.
RETURN
END
FUNCTION TRAP (IX)
T = 0,75*RAND (IX)
TRAP = T + 0,25*RAND (IX)
RETURN
END
FUNCTION GAUSS (IX)
B = SQRT (-2.*ALOG ((RAND (IX) + 1.)/2.))
A = 6.283185*RAND (IX)
GAUSS = B*COS (A)
RETURN
END
FUNCTION AMOD1 (IX)
T = SIMP (IX)
IF (T) 1, 2, 2
1
AMOD1 = - 1. - T
RETURN
2
AMOD1 = 1. - T
RETURN
END
FUNCTION AMOD2 (IX)
T = SIMP (IX)
IF (T) 1, 2, 2
1
AMOD2 = - 1. - 2.*T/3.
RETURN
2
AMOD2 = 1. - 2.*T/3.
RETURN
END
FUNCTION RANDS (IX)
IY = IX*65539
IF (IY) 1, 1, 2
1
IY = IY + 2147483647 + 1
2
Y = IY
Y = Y*.4656613E - 9
IX = IY
RANDS = (2.*Y - 1.)*1.73
RETURN
END
FUNCTION SIMPS (IX)
T = RAND (IX)/2.
SIMPS = (T + RAND (IX)/2.)*2.44
RETURN
END
FUNCTION TRAPS (IX)
T = 0.75*RAND (IX)
TRAPS = (T + 0.25*RAND (IX)) *2.19
RETURN
END
FUNCTION AMODS1 (IX)
T = SIMP (IX)
IF (T) 1, 2, 2
1
AMODS1 = (- 1. - T)*1.41
RETURN
2
AMODS1 = (1. - T)*1.41
RETURN
END
FUNCTION AMODS2 (IX)
T = SIMP (IX)
IF (T) 1, 2, 2
1
AMODS2 = (- 1. - 2.*T/3.)*1.27
RETURN
2
AMODS2 = (1. - 2.*T/3.)*1.27
RETURN
END
SUBROUTINE PRINOP (MTEXT, MTABL, MGRAF, NOUS, NPP)
C
C
C
СЛУЖЕБНЫЕ ПОДПРОГРАММЫ -
PRINOP, PROBIN, ERPROC, TXOUT, RPROUT, IPROUT,
C
PRIM1, DISPET
C
C
C
C
C
НАЗНАЧЕНИЕ
C
C
PRINOP:
ОБРАБОТКА ДОПОЛНИТЕЛЬНЫХ ЗАДАНИЙ ПОЛЬЗОВАТЕЛЯ
C
PROBIN
ВВОД ЧИСЛОВЫХ ДАННЫХ И ТЕКСТОВОЙ ИНФОРМАЦИИ
C
ERPROC:
ОБРАБОТКА ОШИБОК ВВОДА
C
TXOUT:
РАСПЕЧАТКА ТАБЛИЦ ИСХОДНЫХ ДАННЫХ И КОММЕНТАРИЕВ К НИМ
C
RPROUT:
C
IPROUT
C
PRIM1:
C
DISPET:
УПРАВЛЕНИЕ ДОПОЛНИТЕЛЬНЫМИ ПРОГРАММАМИ
C
ПОЛЬЗОВАТЕЛЯ
C
10
FORMAT (/T10, '* ДОПОЛНИТЕЛЬНЫЕ ВОЗМОЖНОСТИ *')
11
FORMAT (T10, '* ПЕЧАТЬ ТЕКСТА ПОЛНАЯ *')
12
FORMAT (T10, ' * ОБЩАЯ ИНФОРМАЦИЯ И КОММЕНТАРИИ *')
13
FORMAT (T10, ' * ОБЩАЯ ИНФОРМАЦИЯ *')
14
FORMAT (T10, '* ПЕЧАТЬ ТАБЛИЦ ПОЛНАЯ*')
15
FORMAT (T10, '* ТАБЛИЦЫ РЕЗУЛЬТАТОВ И ИСХ. ДАННЫХ*')
16
FORMAT (T10, '* ТАБЛИЦА РЕЗУЛЬТАТОВ*')
17
FORMAT (T10, '* ПЕЧАТЬ ГРАФИКОВ ПОЛНАЯ *')
18
FORMAT (T10, '* ГРАФИКОВ НЕ ПЕЧАТАТЬ *')
19
FORMAT (T10, '* ',I1,' ДОПОЛН. ПРОГРАММ ПОЛЬЗОВАТЕЛЯ*')
PRINT 10
PRINT 19, NOUS
IF (NPR.NE.0.AND.NPR.NE 3) GO TO 30
MTABL = 0
MGRAF = 1
30
IF (MTEXT - 1) 20, 21, 22
20
PRINT 11
GO TO 23
21
PRINT 12
GO TO 23
22
PRINT 13
23
IF (MTABL - 1) 24, 25, 26
24
PRINT 14
GO TO 27
25
PRINT 15
GO TO 27
26
PRINT 16
27
IF (MGRAF) 28, 28, 29
28
PRINT 17
RETURN
29
PRINT 18
RETURN
END
SUBROUTINE PROBIN (RPARM, LR, IPARM, LI, TEXT, NBEG, LTX, NEXP, IER)
DIMENSION RPARM (20, 9), IPARM (20, 9)
DIMENSION TEXT (1), LTX (1), NBEG (1)
10
FORMAT (I1)
11
FORMAT (I1, 17A4)
33
FORMAT (10A4)
34
FORMAT (I1, 9F5.0)
35
FORMAT (I1, 915)
IER = 2
NE = 0
5
READ 10, NBL
IF (NBL) 30, 31, 30
31
RETURN
30
NBEG (NBL) = NE + 1
GO TO (1, 2, 2, 4), NBL
2
IF (IER - 1) 21, 22, 23
21
IER = 3
23
RETURN
22
DO 7 I = 1, 20
N = NBEG (NBL) + 10* (I - 1)
L = N + 9
READ 33, (TEXT (J), J = N, L)
IF (NBL.EQ.3) GO TO 24
READ 34, NEND, (RPARM (I, J), J = 1, NEXP)
IF (NEND - 9) 7, 26, 7
26
LR = I
GO TO 32
24
READ 35, NEND, (IPARM (I, J), J = 1, NEXP)
IF (NEND - 9) 7, 36, 7
36
LI = I
GO TO 32
7
CONTINUE
IER = 4
RETURN
32
LTX (NBL) = 10*I
GO TO 16
1
NE = 0
4
DO 6 I = 1, 100
N = NBEG (NBL) + 17*(I - 1)
L = N + 16
READ 11, NEND, (TEXT (J), J = N, L)
IF (NEND - 9) 6, 15, 6
6
CONTINUE
IER = 5
RETURN
15
IER = IER - 1
LTX (NBL) = 17*I
16
NE = L
GO TO 5
END
SUBROUTINE ERPROC (IER, MTEXT)
50
FORMAT (T20, '***ВВОД ДАННЫХ - КОД ОШИБКИ = ',I1,'***')
51
FORMAT (T20, '***НЕОБЯЗАТЕЛЬНЫЙ БЛОК ОТСУТСТВУЕТ***')
52
FORMAT (T20, '***БЛОК ДАННЫХ ОТСУТСТВУЕТ***')
53
FORMAT (T20, '***НАРУШЕН ПОРЯДОК ВВОДА***')
54
FORMAT (T20, '***ПРЕВЫШЕНА ДЛИНА БЛОКА ДАННЫХ***')
55
FORMAT (T20, '***ПРЕВЫШЕНА ДЛИНА БЛОКА ТЕКСТА***')
56
FORMAT (T20, '***ВЫВОД БЛОКА НЕ ОБЕСПЕЧЕН***')
PRINT 50, IER
IF (IER - 1) 60, 61, 62
61
PRINT 51
IF (MTEXT) 63, 63, 60
63
PRINT 56
MTEXT = 1
GO TO 60
62
IF (IER - 3) 64, 65, 66
64
PRINT 52
STOP
65
PRINT 53
STOP
66
IF (IER - 5) 67, 68, 68
67
PRINT 54
STOP
68
PRINT 55
STOP
60
RETURN
END
SUBROUTINE TXOUT (TEXT, NBEG, LTX, MOD)
DIMENSION TEXT (1)
1
FORMAT (1H0)
2
FORMAT (T20, 17A4)
3
FORMAT (T20, 17A4/)
L = LTX + NBEG - 1
PRINT 1
IF (MOD) 20, 20, 30
20
PRINT 2, (TEXT (I), I = NBEG, L)
RETURN
30
PRINT 3, (TEXT (I), I = NBEG, L)
RETURN
END
SUBROUTINE RPROUT (TEXT, NBEG, LTX, RPARM, LR, NEXP, MOD)
DIMENSION TEXT (1)
DIMENSION RPARM (20, 9)
1
FORMAT (1H0)
2
FORMAT (T2, 10A4, 9F8.3)
3
FORMAT (T20, 'ПАРАМЕТР', I1X, 9I8)
4
FORMAT (1H)
PRINT 1
PRINT 3, (I, I = 1, NEXP)
PRINT 1
N = NBEG
L = NBEG + 9
DO 10 K = 1, LR
IF (MOD.EQ.l) PRINT 4
PRINT 2, (TEXT (I), I = N, L), RPARM (K, J), J = 1, NEXP)
N = L + 1
10
L = N + 9
RETURN
END
SUBROUTINE IPROUT (TEXT, NBEG, LTX, IPARM, LI, NEXP, MOD)
DIMENSION TEXT (1)
DIMENSION IPARM (20, 9)
2
FORMAT (T2, 10A4, 918)
4
FORMAT (1H)
N = NBEG
L = NBEG + 9
DO 10 K = 1, LI
IF (MOD.EQ.l) PRINT 4
PRINT 2, (TEXT (I), I = N, L), IPARM (K, J), J = 1, NEXP)
N = L + 1
10
L = N + 9
RETURN
END
SUBROUTINE PRIM1
18
FORMAT (T20, 'ПРИМЕЧАНИЕ: ПРИНЯТАЯ КОДИРОВКА ЗАКОНОВ РАСПРЕД.:',
1/T33,' 1, 7 - ЗАКОН РАВНОЙ ПЛОТНОСТИ',
2/T33,' 2, 8 - ЗАКОН СИМПСОНА (ТРЕУГОЛЬНЫЙ)',
3/T33,' 3, 9 - ТРАПЕЦЕИДАЛЬНЫЙ ЗАКОН',
4/T33,' 4, НОРМАЛЬНЫЙ ЗАКОН',
5/T33,' 5, 10 - АНТИМОДАЛЬНЫЙ 1'
6/T33,' 6, 11 - АНТИМОДАЛЬНЫЙ 2')
PRINT 18
RETURN
END
SUBROUTINE DISPET (N)
GO TO (1, 2, 3, 4, 5, 6, 7, 8, 9), N
1
CALL USPGM1
RETURN
2
CALL USPGM2
RETURN
3
CALL USPGM3
RETURN
4
CALL USPGM4
RETURN
5
CALL USPGM5
RETURN
6
CALL USPGM6
RETURN
7
CALL USPGM7
RETURN
8
CALL USPGM8
RETURN
9
CALL USPGM9
RETURN
END
FUNCTION PROB (N, MET)
C
C
C
ПОДПРОГРАММЫ - PROB, PFB
C
C
C
C
C
НАЗНАЧЕНИЕ
C
C
PROB:
ПРОИЗВОДИТ ОЦЕНКУ ВЕРОЯТНОСТИ ПРИЗНАНИЯ СИ ГОДНЫМ
C
PFB:
ВЫЧИСЛЯЕТ СРЕДНЮЮ И МАКСИМАЛЬНУЮ
C
ВЕРОЯТНОСТИ ФИКТИВНОГО БРАКОВАНИЯ
C
C
C
C
C
ОПИСАНИЕ ФОРМАЛЬНЫХ ПАРАМЕТРОВ
C
C
N
- ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ СОБЫТИЙ
C
MET
- ИМЯ ПРОГРАММЫ МЕТОДИКИ ПОВЕРКИ
C
BETA
- ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
C
STEP
- ШАГ АРГУМЕНТА ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
C
P
- ОТСЧЕТЫ ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
C
C
LOGICAL GOD
COMMON/LOGIC/GOD
P = 0.
DO 1 I = 1, N
CALL MET
IF (GOD) P = P + 1.
1
CONTINUE
PROB = P/N
RETURN
END
SUBROUTINE PFB (BETA, STEP, P, PF, PFM)
DIMENSION P (1)
NSTEP = (BETA + .01)/STEP + 1
PF = (P(1) + P(NSTEP))/2.
K = NSTEP - 1
DO 1 I = 2, K
1
PF = PF + P (I)
PF = BETA - STEP*PF
PFM = 1. - P (NSTEP)
RETURN
END
Приложение 3
Справочное
СООБЩЕНИЯ
КОМПЛЕКСА ПРОГРАММ
3.1. Сообщения комплекса программ включают в себя: контрольно-диагностические сообщения; сообщения о результатах расчета; информационные сообщения.
Контрольно-диагностические сообщения предназначены для контроля правильности кодирования исходных данных в пакете задания и обработки основных программ комплекса в процессе расчета.
Результаты расчета включают таблицу значений критериев качества поверки СИ, таблицы и графики оперативных характеристик.
Информационные сообщения служат для облегчения расшифровки результатов расчета, классификации и хранения распечаток.
3.2. Контрольно-диагностические сообщения.
3.2.1. Сообщения о дополнительных указаниях пользователя "ДОПОЛНИТЕЛЬНЫЕ ВОЗМОЖНОСТИ". Вслед за данным сообщением распечатываются возможности, предоставляемые пользователю в текущем задании. Если представляемые возможности не соответствуют указаниям пользователя по п. 2.1.9, следует проверить правильность их кодирования в УК.
3.2.2. Сообщения об ошибках ввода данных:
а) "ВВОД ДАННЫХ - КОД ОШИБКИ n". Сообщение распечатывается после завершения ввода исходных данных. n = 0 указывает на отсутствие ошибок ввода. При других кодах вслед за последним сообщением распечатывается диагностическое сообщение о характере ошибки;
б) "НЕОБЯЗАТЕЛЬНЫЙ БЛОК ОТСУТСТВУЕТ". Причина. Отсутствует необязательный блок комментариев к методике поверки (блок 4). Если в УК пользователем указана распечатка текстовой информации в полном объеме, вслед за данным сообщением появляется сообщение "ВЫВОД БЛОКА НЕ ОБЕСПЕЧЕН". Обработка задания продолжается;
в) "БЛОК ДАННЫХ ОТСУТСТВУЕТ". Причина. Отсутствует один или несколько блоков исходных данных или символы, идентифицирующие начало или конец блока. Обработка задания прекращается;
г) "НАРУШЕН ПОРЯДОК ВВОДА". Причина. Блоки исходных данных следуют в пакете задания не по порядку возрастания их номеров. Обработка задания прекращается;
д) "ПРЕВЫШЕНА ДЛИНА БЛОКА ДАННЫХ". Причина. Число параметров хотя бы в одном из блоков числовых исходных данных превышает 20. Обработка задания прекращается.;
е) "ПРЕВЫШЕНА ДЛИНА БЛОКА ТЕКСТА". Причина. Общий объем текстовой информации превышает 3000 знаков. Обработка задания прекращается.
3.2.3. Сообщения об исходных данных. "ТАБЛИЦА ИСХОДНЫХ ДАННЫХ". Вслед за данным сообщением распечатывается содержание блоков 2 и 3 исходных данных к заданию. Таблица включает в себя колонку наименований параметров и колонки значений, соответствующие вариантам решения задачи. Таблица служит для контроля правильности кодирования наименований и числовых значений параметров в пакете задания. Последовательно распечатываются параметры вещественного, затем целого типов.
3.2.4. По завершении расчета значений критериев качества поверки СИ по каждому из вариантов исходных данных распечатывается сообщение "ДИАГНОСТИКА ОБРАЩЕНИЯ К П/П ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ", "КОД ОШИБКИ - n", n = 0 указывает на нормальное завершение работы комплекса программ. Другие коды сопровождаются следующими сообщениями:
а) "БРАК ПРИ МИНИМАЛЬНОМ ЗНАЧЕНИИ КОНТР. ПАРАМЕТРА". Причина. Возможно бракование поверяемых СИ, у которых контролируемая характеристика погрешности равна нулю. При отсутствии ошибок в задании исходных данных указанное событие возможно, если методика поверки предусматривает чрезмерно жесткие условия бракования или недостаточную точность образцового СИ;
б) "НУЛЕВОЙ НЕОБНАРУЖЕННЫЙ БРАК". Причина. Не наблюдается случаев признания годным фактически негодных СИ. Возможно наличие чрезмерно жестких условий бракования поверяемых СИ;
в) "МАКС. ВЫХОД ЗА ДОПУСК > 2". Причина. Возможно признание годным СИ, контролируемая характеристика погрешности которого превышает удвоенное значение предела допускаемых значений. Фактическое значение максимального выхода за допуск не может быть вычислено. Возможна недостаточная точность образцового СИ, применяемого при поверке, или чрезмерно большое значение с.к.о. случайной составляющей погрешности поверки, или недостаточное число отсчетов погрешности при поверке;
г) "ВЕР-ТЬ Н. Б. РАВНА 1". Причина. Практически невероятно бракование СИ, значение контролируемой характеристики погрешности которого равно пределу допускаемых значений. Возможно наличие слишком мягких условий бракования СИ;
д) "ВЕРОЯТНОСТЬ Ф. Б. РАВНА 1". Причина. Максимальная вероятность фиктивного брака близка к 1.
Возможно наличие чрезмерно жестких условий бракования, в результате чего возможно бракование большого числа фактически годных СИ.
3.3. Сообщения о результатах расчета содержат: таблицы значений критериев качества поверки СИ; таблицы оперативных характеристик (ОХ); график ОХ.
3.3.1. В таблице значений критериев качества поверки СИ каждая колонка значений соответствует заданному варианту исходных данных.
3.3.2. Таблицы ОХ включают в себя действительные значения контролируемой характеристики погрешности СИ, отнесенные к пределу допускаемых значений, и оценки вероятности признания годным СИ, контролируемая характеристика погрешности которого имеет соответствующее значение. Номер колонки значений вероятности соответствует номеру варианта исходных данных.
3.3.3. График ОХ наглядно интерпретирует таблицу ОХ. Шаг дискретизации оси абсцисс составляет 0,02 предела допускаемых значений контролируемой характеристики погрешности СИ, по оси ординат - 0,01.
Кривые на графике кодируются цифрами, соответствующими варианту исходных данных. Совпадение двух или нескольких кривых в одной точке кодируется символом "О".
3.4. Информационные сообщения.
3.4.1. "ЗАДАНИЕ n". В сообщении печатается условный номер задания, кодируемый пользователем в УК. Вслед за этим сообщением распечатывается заголовок: "оценка качества поверки СИ".
3.4.2. "ОБЩАЯ ИНФОРМАЦИЯ О ЗАДАНИИ". Под этим заголовком распечатывается блок 1 исходных данных в том виде, в котором он был закодирован пользователем.
3.4.3. "КОММЕНТАРИИ К МЕТОДИКАМ". Под этим заголовком распечатывается блок 4 исходных данных в том виде, в котором он был закодирован пользователем.
Помимо указанных сообщений в данной версии комплекса программ в качестве комментария распечатывается принятая кодировка программ генераторов случайных чисел.
Приложение 4
Справочное
ПРИМЕРЫ
ВЫПОЛНЕНИЯ РАСЧЕТОВ КРИТЕРИЕВ КАЧЕСТВА ПОВЕРКИ СИ МЕТОДОМ
ПРОГРАММНОГО МОДЕЛИРОВАНИЯ
ПРИМЕР 1 (цифровое средство измерений).
ИС МЕГАНОРМ: примечание.
МИ 187-79 утратили силу в связи с введением в действие МИ 187-86 "Методические указания. Государственная система обеспечения единства измерений. Достоверность и требования к методикам поверки средств измерений", утв. ВНИИМС 25.07.1986.
Постановка задачи. Необходимо рассчитать значения установленных МИ 187-79 критериев качества проверки ( - максимальный выход за допуск; Pнм - максимальная вероятность необнаруженного брака; Pф - средняя вероятность фиктивного брака; Pфм - максимальная вероятность фиктивного брака) систематической составляющей основной погрешности цифрового вольтметра, который обладает перечисленными ниже характеристиками, по указанной ниже методике.
1. Принцип действия - время-импульсный, погрешность из-за гистерезиса отсутствует по принципу действия прибора.
2. Номинальная ступень квантования qsf = 5 мВ.
3. Пределы допускаемых значений систематической составляющей основной погрешности .
4. Предел допускаемых значений с.к.о. случайной составляющей основной погрешности в зависимости от модификации прибора; 2,5 или 10 мВ.
Поверка такого вольтметра производится в соответствии с указаниями п. 4.4.11.2 методики МИ 118-77, так как с.к.о. случайной составляющей основной погрешности превышает 0,25 qsf.
Принимаем, что в качестве образцового средства измерений используется калибратор (многозначная мера), для которого нормированы пределы допускаемых значений погрешности.
ЭТАП 1. Подготовка исходных данных (см. п. 1.6)
1. Формула, описывающая связь показаний Y поверяемого вольтметра (см. п. 2.1.1 основного текста), обладающего указанными выше свойствами, со значением измеряемой величины имеет вид <*>
--------------------------------
<*> Н.Н. Вострокнутов. Испытания и поверка цифровых измерительных устройств, - М., Изд-во стандартов, 1977, 148 с.
,
где Xвх - входное напряжение поверяемого СИ; Int - операция выделения целой части числа; - систематическая составляющая погрешности поверяемого СИ, представляемая в соответствии с указаниями п. 2.1.1.1 в виде простой переменной; - случайная составляющая погрешности поверяемого СИ, представляемая в соответствии с указаниями п. 2.1.1.2 формулой (2.3): , где - безразмерная случайная величина, распределенная по трапецеидальному закону с с.к.о., равным 1 (см. п. 2.1.1.1).
2. Формула, описывающая связь выходного сигнала Ye образцового средства измерений (калибратора) выбирается в соответствии с указаниями п. 2.1.2: , где Nx - номинальное значение выходного сигнала калибратора; - погрешность калибратора (см. п. 2.1.2.1);
- безразмерная случайная величина, распределенная по закону равной плотности в пределах от -1 до +1 для ансамбля образцовых СИ и принимающая неизменное значение для конкретного экземпляра образцового СИ.
3. Описание алгоритма выполнения поверки (см. 2.1.3) производят на основании п. 4.4.11.2 МИ 118-77.
3.1. Устанавливают выходной сигнал калибратора, соответствующим поверяемой точке Nx.
3.2. Регистрируют показание Yi поверяемого СИ.
3.3. Повторяют операцию п. 3.2 N = 10 раз (при i = 1, 2, ... 10).
3.4. Вычисляют оценку систематической составляющей погрешности поверяемого СИ по формуле: .
3.5. Если выполняется неравенство , где - контрольный допуск, поверяемый прибор бракуют, в противном случае - признают годным.
По приведенному описанию алгоритма разрабатывают его блок-схему (см. п. 2.1.3), приведенную на рис. 1.
Рис. 1
ИС МЕГАНОРМ: примечание.
МИ 187-79 утратили силу в связи с введением в действие МИ 187-86 "Методические указания. Государственная система обеспечения единства измерений. Достоверность и требования к методикам поверки средств измерений", утв. ВНИИМС 25.07.1986.
4. На основании приведенных исходных данных, полученных из нормативной документации на поверяемое СИ, задаем значения параметров, указанных в 2.1.4: ; ; (граница фиктивного бракования - см. МИ 187-79); Nр = 500 (см. рекомендации п. 2.1.4); Nм = 1 (по заданию исследованию подлежит только один алгоритм поверки).
5. Задаем число вариантов исходных данных (см. п. 2.1.5). В нашем случае это число равно двум, так как в соответствии с заданием значения критериев качества нужно вычислить при двух различных значениях .
6. Составляем и записываем текст общей информации о задании (см. п. 2.1.6): "Расчет критериев качества проверки систематической составляющей погрешности цифрового вольтметра типа ... при наличии существенной случайной составляющей погрешности поверяемого СИ. Расчет производится при двух различных пределах допускаемых значений с.к.о. случайной составляющей погрешности поверяемого СИ".
Дополнительные комментарии к методике поверки не приводятся (см. п. 2.1.7), на печать следует выводить только общую информацию по п. 2.1.7.
7. Формулируем указания программисту по требуемому объему печати выходных данных (информации, выводимой на печать по окончании расчета) - см. п. 2.1.8. Требуется напечатать: таблицы исходных данных (для проверки правильности кодирования задания); значения оперативных характеристик (для оценки достаточности выбранного числа Nр); графики оперативных характеристик (для наглядного сравнения результатов расчета и диагностики возможных ошибок).
8. Назначаем условный номер задания (см. п. 2.1.9) - 1.
ИС МЕГАНОРМ: примечание.
В официальном тексте документа, видимо, допущена опечатка: имеется в виду п. 2.2.1.
9. Оформляем подготовленные исходные данные в виде табл. 1 (см. п. 2.2.9).
Таблица 1
Обозначения
Наименование параметра
Признак
Имя
Систематическая составляющая погрешности поверяемого СИ
Контролируемый параметр
DSYST
Контрольный допуск
Контрольный допуск
DCONTR
Предел с.к.о. случайной составляющей погрешности поверяемого СИ
Параметр вещественного типа
SIGMA
Предел допускаемых значений погрешности калибратора
То же
DETAL
Nx
Номинальное значение выходного сигнала калибратора (поверяемая точка)
"
PNOM
qsf
Ступень квантования поверяемого СИ (номинальная)
"
QUANT
N
Число регистрируемых показаний поверяемого СИ
Параметр целого типа
NMES
Xвх
Входной сигнал поверяемого СИ
Промежуточная переменная
POINT
Y
Показание поверяемого СИ
То же
SHOW
Случайная составляющая погрешности поверяемого СИ (трапецеидальный закон распределения с с.к.о. равным )
Случайная величина
DRND
Погрешность калибратора (равномерное распределение по ансамблю СИ в пределах )
Случайная величина по ансамблю СИ
ERROR
Примечание. Графу 4 заполняет программист при выполнении этапа 3 (см. ниже).
10. В соответствии с указаниями п. 2.2.2 составляем таблицы параметров вещественного типа (табл. 2) и целого типа (табл. 3).
Таблица 2
N п/п
Наименование параметра
Номер варианта
1
2
1
ПДЗ контролируемого параметра
10,0
10,0
2
Отношение гамма
0,9
0,9
3
Граница фиктивного бракования
0,8
0,8
4
С.к.о. случайной составляющей погрешности поверяемого СИ
2,5
10,0
5
ПДЗ погрешности образцового СИ
3,33
3,33
6
Поверяемая точка
995,0
995,0
7
Ступень квантования поверяемого СИ
5,0
5,0
Примечание. Здесь и ниже ПДЗ - предел допускаемых значений.
Таблица 3
N п/п
Наименование параметра
Номер варианта
1
2
1
Объем выборки для оценивания вероятности
500
500
2
Номер методики поверки
1
1
3
Число измерений в цикле
10
10
11. Все материалы, разработанные на этапе 1, передаются программисту для программирования задания и выполнения работ по этапам 2 и 3 (см. п. 1.6).
ЭТАП 2. Программирование задания (см. п. 1.6).
1. Кодируем УК по указаниям пользователя (см. п. 3.2.2):
2. Кодируем блок 1 (см. п. 3.2.3) пакета исходных данных:
1
РАСЧЕТ ЗНАЧЕНИЙ КРИТЕРИЕВ КАЧЕСТВА КОНТРОЛЯ СИСТЕМАТИЧЕСКОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ЦИФРОВОГО ВОЛЬТМЕТРА ТИПА ... ПРИ НАЛИЧИИ СУЩЕСТВЕННОЙ СЛУЧАЙНОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ПОВЕРЯЕМОГО СИ. РАСЧЕТ ПРОИЗВОДИТСЯ ДЛЯ ДВУХ РАЗЛИЧНЫХ ЗНАЧЕНИЙ
9
С.К.О. СЛУЧАЙНОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ПОВЕРЯЕМОГО СИ.
3. Кодируем блоки 2 и 3 исходных данных:
2
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
10.0.10.0
ОТНОШЕНИЕ ГАММА
0.9 0.9
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0.8 0.8
С.К.О. СЛ. СОСТ. ПОГРЕШНОСТИ ПОВЕР. СИ
2.5 10.0
ПДЗ ПОГРЕШНОСТИ ОБРАЗЦОВОГО СИ
3.3333.333
ПОВЕРЯЕМАЯ ТОЧКА
995. 995.
СТУПЕНЬ КВАНТОВАНИЯ ПОВЕРЯЕМОГО СИ
9
5.0 5.0
3
ЧИСЛО ПОВЕРОК ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500 500
НОМЕР МЕТОДИКИ ПОВЕРКИ
1
1
ЧИСЛО ИЗМЕРЕНИЙ ПРИ ПОВЕРКЕ
9
10 10
4. Записываем имена параметров в табл. 1 и на основании блок-схемы алгоритма поверки составляем программу (описание) методики поверки (см. пп. 3.1.2 и 3.3):
C
SUBROUTINE MET1
C
ОПИСАНИЕ ПАРАМЕТРОВ
C
LOGICAL GODEN
COMMON/LOGIC/GODEN
COMMON/CONTRL/DSYST, DCONTR
COMMON/INTEG/NMES
COMMON/REAL/SIGMA, DETAL, PNOM, QUANT
COMMON/RND/IX
DIMENSION SHOW (10)
C
C
ОПИСАНИЕ ПРОЦЕДУРЫ КОНТРОЛЯ
C
C
УСТАНАВЛИВАЕТСЯ ВЫХОДНОЙ СИГНАЛ КАЛИБРАТОРА
C
ERROR = DETAL*RAND (IX)
POINT = PNOM + ERROR
C
C
ЦИКЛ РЕГИСТРАЦИИ ПОКАЗАНИЙ ПОВЕРЯЕМОГО СИ
C
DO 1 I = 1, NMES
DRND = SIGMA*TRAPS (IX)
SHOW (I) = QUANT*
1
AINT ((POINT + DRND + DSYST)/QUANT + 0.5)
1
CONTINUE
C
C
ВЫЧИСЛЕНИЕ СРЕДНЕГО ЗНАЧЕНИЯ ПОКАЗАНИЙ
C
SUMMA = 0.0
DO 2 I = 1, NMES
SUMMA = SUMMA + SHOW (1)
2
CONTINUE
YMED = SUMMA/NMES
C
C
СРАВНЕНИЕ С КОНТРОЛЬНЫМ ДОПУСКОМ И
C
ПРИНЯТИЕ РЕШЕНИЯ О ГОДНОСТИ СИ
C
IF (ABS (YMED - PNOM) - DCONTR) 3, 3, 4
3
GODEN = .TRUE.
RETURN
4
GODEN = .FALSE.
RETURN
END
5. Подготавливаем пакет инструкций по операционной системе ЭВМ (см. 3.4).
ЭТАП 3. Решение задачи на ЭВМ (см. п. 1.6)
ИС МЕГАНОРМ: примечание.
Здесь и далее в официальном тексте документа, видимо, допущены опечатки: имеется в виду раздел 4, а не пункт 4.
Задача решается на ЭВМ в соответствии с указаниями п. 4. В результате работы ЭВМ по составленному заданию и в соответствии с материалами, подготовленными на этапах 1 и 2, получим приводимую ниже распечатку.
ЗАДАНИЕ 1
ОЦЕНКА КАЧЕСТВА ПОВЕРКИ СИ
ЧИСЛО ОЦЕНОК 2
* ДОПОЛНИТЕЛЬНЫЕ ВОЗМОЖНОСТИ*
* 0 ДОПОЛН. ПРОГРАММ ПОЛЬЗОВАТЕЛЯ *
* ОБЩАЯ ИНФОРМАЦИЯ *
* ПЕЧАТЬ ТАБЛИЦ ПОЛНАЯ *
* ПЕЧАТЬ ГРАФИКОВ ПОЛНАЯ *
ВВОД ДАННЫХ - КОД ОШИБКИ = 1
НЕОБЯЗАТЕЛЬНЫЙ БЛОК ОТСУТСТВУЕТ
ОБЩАЯ ИНФОРМАЦИЯ О ЗАДАНИИ
РАСЧЕТ ЗНАЧЕНИЙ КРИТЕРИЕВ КАЧЕСТВА КОНТРОЛЯ СИСТЕМАТИЧЕСКОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ЦИФРОВОГО ВОЛЬТМЕТРА ТИПА ... ПРИ НАЛИЧИИ СУЩЕСТВЕННОЙ СЛУЧАЙНОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ПОВЕРЯЕМОГО СИ РАСЧЕТ ПРОИЗВОДИТСЯ ДЛЯ ДВУХ РАЗЛИЧНЫХ ЗНАЧЕНИЙ С.К.О. СЛУЧАЙНОЙ СОСТАВЛЯЮЩЕЙ ПОГРЕШНОСТИ ПОВЕРЯЕМОГО СИ.
ТАБЛИЦА ИСХОДНЫХ ДАННЫХ
ПАРАМЕТР
1
2
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
10,000
10,000
ОТНОШЕНИЕ ГАММА
0,900
0,900
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0,800
0,800
С.К.О. СЛ. СОСТ. ПОГРЕШНОСТИ ПОВЕР. СИ
2,500
10,000
ПДЗ ПОГРЕШНОСТИ ОБРАЗЦОВОГО СИ
3,333
3,333
ПОВЕРЯЕМАЯ ТОЧКА
995 000
995 000
СТУПЕНЬ КВАНТОВАНИЯ ПОВЕРЯЕМОГО СИ
5,000
5,000
ЧИСЛО ПОВЕРОК ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500
500
НОМЕР МЕТОДИКИ ПОВЕРКИ
1
1
ЧИСЛО ИЗМЕРЕНИЙ ПРИ ПОВЕРКЕ
10
10
ДИАГНОСТИКА ОБРАЩЕНИЯ К П/П ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ
1 КОД ОШИБКИ - 0 РНМ = 0,366 ДМ = 1,320 РФ = 0,050 РФМ = ,374
ДИАГНОСТИКА ОБРАЩЕНИЯ К П/П ОЦЕНКИ КАЧЕСТВА ПОВЕРКИ
2 КОД ОШИБКИ - 0 РНМ = 0,392 ДМ = 1,720 РФ = 0,116 РФМ = 0,430
ТАБЛИЦА РЕЗУЛЬТАТОВ
ПОКАЗАТЕЛЬ КАЧЕСТВА
1
2
МАКС. ВЕРОЯТНОСТЬ НЕОБНАРУЖ. БРАКА
0,366
0,392
МАКС. ВЫХОД ЗА ДОПУСК
1,320
1,720
СРЕДН. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,050
0,116
МАКС. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,374
0,430
ТАБЛИЦА ОПЕРАТИВНОЙ ХАРАКТЕРИСТИКИ
КОНТРОЛИРУЕМЫЙ ПАРАМЕТР
1
2
1
0,0
1,000
0,990
2
0,02
1,000
0,974
3
0,04
1,000
0,994
4
0,06
1,000
0,974
5
0,08
1,000
0,980
6
0,10
1,000
0,976
7
0,12
1,000
0,970
8
0,14
1,000
0,982
9
0,16
1,000
0,982
10
0,18
1,000
0,960
11
0,20
1,000
0,968
12
0,22
1,000
0,958
13
0,24
1,000
0,958
14
0,26
1,000
0,940
15
0,28
1,000
0,950
16
0,30
1,000
0,922
17
0,32
1,000
0,936
18
0,34
1,000
0,918
19
0,36
0,994
0,922
20
0,38
0,996
0,918
21
0,40
1,000
0,898
22
0,42
0,998
0,894
23
0,44
0,982
0,854
24
0,46
0,986
0,860
25
0,48
0,998
0,822
26
0,50
0,984
0,852
27
0,52
0,970
0,816
28
0,54
0,944
0,812
29
0,56
0,940
0,822
30
0,58
0,934
0,770
31
0,60
0,918
0,774
32
0,62
0,884
0,738
33
0,64
0,874
0,712
34
0,66
0,860
0,712
35
0,68
0,812
0,698
36
0,70
0,762
0,684
37
0,72
0,784
0,680
38
0,74
0,714
0,664
39
0,76
0,700
0,614
40
0,78
0,672
0,578
41
0,80
0,626
0,570
42
0,82
0,604
0,548
43
0,84
0,576
0,530
44
0,86
0,558
0,476
45
0,88
0,498
0,494
46
0,90
0,478
0,494
47
0,92
0,452
0,430
48
0,94
0,420
0,414
49
0,96
0,432
0,414
50
0,98
0,384
0,376
КОНТРОЛИРУЕМЫЙ ПАРАМЕТР
1
2
51
1,00
0,366
0,392
52
1,02
0,290
0,344
53
1,04
0,274
0,310
54
1,06
0,248
0,290
55
1,08
0,232
0,308
56
1,10
0,228
0,300
57
1,12
0,140
0,224
58
1,14
0,154
0,240
59
1,16
0,100
0,216
60
1,18
0,104
0,208
61
1,20
0,058
0,210
62
1,22
0,044
0,196
63
1,24
0,044
0,168
64
1,26
0,022
0,156
65
1,28
0,016
0,172
66
1,30
0,014
0,102
67
1,32
0,002
0,118
68
1,34
0,004
0,108
69
1,36
0,010
0,094
70
1,38
0,0
0,086
71
1,40
0,0
0,070
72
1,42
0,0
0,080
73
1,44
0,0
0,080
74
1,46
0,0
0,048
75
1,48
0,0
0,052
76
1,50
0,0
0,034
77
1,52
0,0
0,046
78
1,54
0,0
0,030
79
1,56
0,0
0,036
80
1,58
0,0
0,026
81
1,60
0,0
0,026
82
1,62
0,0
0,032
83
1,64
0,0
0,018
84
1,66
0,0
0,018
85
1,68
0,0
0,012
86
1,70
0,0
0,016
87
1,72
0,0
0,006
88
1,74
0,0
0,010
89
1,76
0,0
0,008
90
1,78
0,0
0,0
91
1,80
0,0
0,0
92
1,82
0,0
0,0
ПРИМЕР 2 (аналоговое средство измерений)
Постановка задачи. Рассчитать значения критериев качества проверки погрешности и вариации стрелочного манометра МТИ-25, обладающего указанными ниже характеристиками и свойствами:
а) пределы измерения 0 - 2,4 МПа (0 - 25 кгс/см2);
б) пределы допускаемых значений основной погрешности Xk = +/- 147 МПа (Xk - верхний предел измерения);
в) предел допускаемых значений вариации Hop = 0,006 Xk;
г) цена деления шкалы (0,2 кгс/см2);
д) случайная составляющая погрешности практически отсутствует.
В качестве образцового средства измерений используется задатчик давления, пределы допускаемых значений погрешности которого Nx, где Nx - номинальное значение установленного давления.
Погрешность и вариацию в каждой проверяемой точке шкалы определяют по изложенной ниже методике (ГОСТ 15614-70):
а) поднимая давление, устанавливают его равным номинальным значениям давления Nx для каждой проверяемой отметки шкалы и записывают показания Y' поверяемого прибора, округляя результат до 0,1 деления шкалы;
б) снижая давление, устанавливают его равным номинальным значениям давления Nx для каждой проверяемой отметки шкалы и записывают показания Y" поверяемого прибора, округляя до 0,1 деления шкалы;
в) подсчитывают оценки погрешности и вариации по формулам: 
где Max - операция выделения максимального числа из перечисленных в скобках: ;
г) если хотя бы в одной точке шкалы выполняется одно или оба неравенства , , где - контрольный допуск для погрешности; - контрольный допуск для вариации, поверяемый манометр бракуют, в противном случае признают годным.
Значение в манометрии принимают равным единице.
ЭТАП 1. Подготовка исходных данных (см. п. 1.6)
1. В соответствии с указаниями п. 1.1 с помощью комплекса программ можно рассчитывать значения критериев качества поверки только для одной из проверяемых характеристик. Поэтому для решения поставленной выше задачи надо произвести два расчета и, соответственно, подготовить два задания на расчет для ЭВМ; первое - для расчета критериев качества проверки погрешности; второе - для расчета критериев качества проверки вариации.
Исходные данные для этих заданий приводятся в табл. 4.
Таблица 4
Задание N 1. Расчет критериев качества проверки погрешности
Задание N 2. Расчет критериев качества проверки вариации
1. Формула, описывающая связь показаний поверяемого манометра (см. п. 2.1.1) со значением измеряемой величины имеет вид:
,
где вычисляется по формуле (2.2) (п. 2.1.1.1), в которой принимают:
(см. постановку задачи); q = 0 (поверяемый манометр - аналоговое СИ), т.е. , в которой - простая переменная; H = Hop; - см. п. 2.1.1.1; H = Hop; D = 1; RH = 1 - при изменении давления в сторону увеличения и RH = -1 - при изменении давления в сторону уменьшения;
где - см. п. 2.1.1.3; - см. п. 2.1.1.1; H - простая переменная - см. п. 2.1.1.4; D = 1 - см. п. 2.1.1.1; RH = l - при изменении давления в сторону увеличения и RH = -1 - при изменении давления в сторону уменьшения;
- погрешность считывания показаний. Показания манометра считываются со шкалы манометра оператором. Из литературы по эргономике известно, что при считывании показаний наблюдатель допускает ошибку в оценке дробной части деления, которая может быть представлена случайной величиной, распределенной по закону Симпсона (треугольному) в пределах +/- 0,1 деления шкалы. Таким образом, может быть представлена в виде: , где - безразмерная случайная величина, распределенная по закону Симпсона в пределах от -1 до +1; R0 - операция округления при считывании показаний. Округление показаний производится наблюдателем до 0,1 деления шкалы.
В общем случае операция округления показаний до 10 деления шкалы может быть представлена в виде: , где Y0 - неокругленное показание поверяемого СИ; Int - операция выделения целой части числа, записанного в скобках.
Таким образом, формула, описывающая связь показаний поверяемого манометра со значением измеряемой величины приобретает вид:
2. Формула, описывающая связь выходного сигнала образцового задатчика давления (калибратора) выбирается в соответствии с указаниями п. 2.1.2: , где Nx - номинальное значение выходного сигнала задатчика; - погрешность задатчика, которая в соответствии с указаниями п. 2.1.2.1 и с учетом свойств погрешности задатчика (см. постановку задачи) должна быть равной: , где - безразмерная случайная величина, распределенная по закону равной плотности в пределах от -1 до +1 для ансамбля задатчиков и принимающая неизменное значение для конкретного экземпляра задатчика давления.
3. Описание алгоритма выполнения поверки приведено выше (см. постановку задачи). Установим, что значения критериев качества следует рассчитывать для поверяемой точки, соответствующей пределу диапазона измерений, т.е. Nx = Xk = 2,4 МПа (25 кгс/см2), так как именно в этой точке получается наибольшая допускаемая погрешность образцового задатчика.
Блок-схемы алгоритмов проверки погрешности и вариации приведены соответственно на:
4. Задаем значения параметров, указанных в п. 2.1.4.
,
,
Nм = 1, Nр = 500
Hop = 0,15 кгс/см2
Nм = 1, Nр = 500
5. Число вариантов исходных данных в каждом задании - 1.
6. Составляем и записываем текст общей информации о задании (см. п. 2.1.6)
"Расчет критериев качества проверки погрешности манометра МТИ-25"
"Расчет критериев качества проверки вариации манометра МТИ-25"
Дополнительные комментарии к методике поверки не приводятся, следует печатать только общую информацию о задании по п. 2.1.7.
7. Даем указания по требуемому объему выходных данных (дополнительную информацию, выводимую на печать после проведения расчета на ЭВМ - см. п. 2.1.8):
Требуется напечатать таблицы исходных данных (для проверки правильности кодирования задания).
8. Назначаем условный номер задания
1
2
9. Оформляем подготовленные числовые исходные данные в виде таблиц (см. п. 2.2.1)
10. Составляем таблицы параметров (см. п. 2.2.1):
вещественного типа - табл. 7
целого типа - табл. 8
вещественного типа - табл. 9
целого типа - табл. 8
Рис. 2
Рис. 3
Таблица 5
Обозначение
Наименование характеристики
Классификационный признак
Имя
Y'
Показание манометра при регулировании давления в сторону увеличения
Промежуточная переменная
SHOW1
Y"
То же - в сторону уменьшения
То же
SHOW2
Цена деления
Параметр вещественного типа
CDEL
X
Входное давление
Промежуточная переменная
PINP
Погрешность манометра
Контролируемый параметр
DELTA
Nx
Номинальное значение выходного сигнала
Параметр вещественного типа
PNOM
Ye
Выходной сигнал задатчика
Промежуточная переменная
POUT
Погрешность задатчика
То же
DELET
Предел допускаемых значений погрешности манометра
Контрольный допуск
DCONTR
Hop
Предел допускаемых значений вариации
Параметр вещественного типа
HOP
Оценка погрешности манометра
Промежуточная переменная
D
RH
Направление регулирования давления
То же
RH
Безразмерная случайная величина, распределенная по закону Симпсона в пределах от -1 до +1
Случайная величина
SIMP (IX)
Безразмерная случайная величина, распределенная по закону равной плотности в пределах от -1 до +1
Случайная величина по ансамблю СИ
RAND (IX)
Таблица 6
Обозначения
Наименования характеристик
Классификационный признак
Имя
Y'
Показание манометра при регулировании давления в сторону увеличения
Промежуточная переменная
SHOW1
Y"
То же, при уменьшении давления
То же
SHOW2
Цена деления
Параметр вещественного типа
CDEL
X
Входное давление
Промежуточная переменная
PINP
H
Вариация
Контролируемый параметр
HVAR
Безразмерная случайная величина, распределенная по закону Симпсона в пределах +/- 1
Случайная величина
SIMP (IX)
Nx
Номинальное значение выходного сигнала задатчика
Параметр вещественного типа
PNOM
Ye
Выходной сигнал задатчика
Промежуточная переменная
POUT
Погрешность задатчика
То же
DELET
Безразмерная случайная величина, распределенная по закону равной плотности в пределах +/- 1
Случайная величина по ансамблю СИ
RAND (IX)
Hop
Предел допускаемых значений вариации
Контрольный допуск
HCONTR
Оценка вариации
Промежуточная переменная
H
RH
Направление регулирования
То же
RH
Таблица 7
N п/п
Параметр
Номер варианта
1
1
ПДЗ контролируемого параметра
0,15
2
Отношение гамма
1,0
3
Граница фиктивного бракования
0,8
4
Цена деления шкалы
0,2
5
Номинальное значение выходного сигнала задатчика
25,0
6
ПДЗ вариации
0,15
Таблица 8
N п/п
Параметр
Номер варианта
1
1
Объем выборки для оценки вероятностей
500
2
Условный номер алгоритма
1
Таблица 9
N п/п
Параметр
Номер варианта
1
1
ПДЗ контролируемого параметра
0,15
2
Отношение гамма
1,0
3
Граница фиктивного бракования
0,8
4
Цена деления шкалы
0,2
5
Номинальное значение выходного сигнала задатчика
25,0
11. Все разработанные на этапе 1 материалы передают программисту для выполнения этапов 2 и 3.
ЭТАП 2. Программирование задания (см. п. 1.6)
Задание N 1
1. Кодируем УК по указаниям пользователя (см. п. 3.2.2): _ _ _ _1_2_1_ _ _ 2_1_1
2. Кодируем блок 1 (см. п. 3.2.3) пакета исходных данных:
1
РАСЧЕТ КРИТЕРИЕВ КАЧЕСТВА ПРОВЕРКИ
9
ПОГРЕШНОСТИ МАНОМЕТРА МТИ-25
3. Кодируем блоки 2 и 3 исходных данных:
2
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
0,15
ОТНОШЕНИЕ ГАММА
1,0
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0,8
ЦЕНА ДЕЛЕНИЯ ШКАЛЫ
0,2
НОМ. ЗНАЧЕНИЕ ВЫХ. СИГНАЛА ЗАДАТЧИКА
25,0
ПДЗ ВАРИАЦИИ
9
0,15
3
ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500
УСЛОВНЫЙ НОМЕР АЛГОРИТМА
9
1
4. Записываем имена параметров в табл. 5 и на основании блок-схемы алгоритма (рис. 2) составляем программу проверки погрешности манометра (см. пп. 3.1.2 и 3.3):
SUBROUTINE MET1
LOGICAL GODEN
COMMON/LOGIC/GODEN
COMMON/RND/IX
COMMON/CONTRL/DELTA, DCONTR
COMMON/REAL/CDEL, PNOM, HOP
C
C
УСТАНОВИТЬ ВХОДНОЙ СИГНАЛ В СТОРОНУ УВЕЛИЧЕНИЯ И
C
ЗАПИСАТЬ ПОКАЗАНИЕ ПОВЕРЯЕМОГО СИ
C
DELET = 0.0005*RAND (IX)*PNOM
POUT = PNOM + DELET
PINT = POUT
RH = 1.
SHOW1 = CDEL*AINT (10.*((PINP + DELTA -
*
HOP/2. + HOP*RH/2.)/CDEL + 0.1*SIMP (IX)) + 0.5)/10.
C
C
УСТАНОВИТЬ ВХОДНОЙ СИГНАЛ В СТОРОНУ УМЕНЬШЕНИЯ И
C
ЗАПИСАТЬ ПОКАЗАНИЕ ПОВЕРЯЕМОГО СИ
C
RH = -1.
SHOW2 = CDEL*AINT (10.*((PINP + DELTA -
*
HOP/2. + HOP*RH/2.)/CDEL + 0.1*SIMP (IX)) + 0.5)/10.
C
C
ОЦЕНКА ПОГРЕШНОСТИ И ПРИНЯТИЕ РЕШЕНИЯ О ГОДНОСТИ СИ
C
D = AMAX1 (ABS(SHOW1 - PNOM), ABS (SHOW2 - PNOM))
IF (D - DCONTR) 1, 1, 2
1
GODEN = .TRUE.
RETURN
2
GODEN = .FALSE.
RETURN
END
5. Подготавливаем пакет инструкций операционной системы ЭВМ (см. п. 3.4)
Задание N 2
1. Кодируем УК по указаниям пользователя (см. п. 3.2.2): _ _ _ _2_2_1_ _ _2_1_1
2. Кодируем блок 1 (см. п. 3.2.3) пакета исходных данных:
1
РАСЧЕТ КРИТЕРИЕВ КАЧЕСТВА ПРОВЕРКИ
9
ВАРИАЦИИ МАНОМЕТРА МТИ-25
3. Кодируем блоки 2 и 3 исходных данных:
2
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
0,15
ОТНОШЕНИЕ ГАММА
1,0
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0,8
ЦЕНА ДЕЛЕНИЯ ШКАЛЫ
0,2
НОМ. ЗНАЧЕНИЕ ВЫХ. СИГНАЛА ЗАДАТЧИКА
9
25,0
3
ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500
УСЛОВНЫЙ НОМЕР АЛГОРИТМА
9
1
4. Записываем имена параметров в табл. 6 и на основании блок-схемы алгоритма составляем программу проверки вариации манометра (см. пп. 3.1.2 и 3.3):
SUBROUTINE MET1
LOGICAL GODEN
COMMON/LOGIC/GODEN
COMMON/RND/IX
COMMON/CONTRL/HVAR, HCONTR
COMMON/REAL/CDEL, PNOM
C
C
УСТАНОВИТЬ ВХОДНОЙ СИГНАЛ В СТОРОНУ УВЕЛИЧЕНИЯ И
C
ЗАПИСАТЬ ПОКАЗАНИЕ ПОВЕРЯЕМОГО СИ
C
DELET = 0.0005*RAND (IX)*PNOM
POUT = PNOM + DELET
PINT = POUT
RH = 1.
SHOW1 = CDEL*AINT(10.*((PINP +
*
HVAR*RH/2.)/CDEL + 0.1*SIMP (IX)) + 0.5)/10.
C
C
УСТАНОВИТЬ ВХОДНОЙ СИГНАЛ В СТОРОНУ УМЕНЬШЕНИЯ И
C
ЗАПИСАТЬ ПОКАЗАНИЕ ПОВЕРЯЕМОГО СИ
C
RH = -1.
SHOW2 = CDEL*AINT (10.*((PINP +
*
HVAR*RH/2.)/CDEL + 0.1*SIMP (IX)) + 0.5)/10.
C
C
ОЦЕНКА ВАРИАЦИИ И ПРИНЯТИЕ РЕШЕНИЯ О ГОДНОСТИ СИ
C
H = ABS (SHOW1 - SHOW2)
IF (H - HCONTR) 1, 1, 2
1
GODEN = .TRUE.
RETURN
2
GODEN = .FALSE.
RETURN
END
5. Используем пакет инструкций операционной системы ЭВМ, подготовленный для задания N 1 данного примера.
ЭТАП 3. Решение задачи на ЭВМ (см. п. 1.6).
Задача решается на ЭВМ в соответствии с указаниями п. 4. В результате работы ЭВМ по составленным заданиям и в соответствии с материалами, подготовленными на этапах 1 и 2, получим приводимые ниже распечатки.
ЗАДАНИЕ
1
ОЦЕНКА КАЧЕСТВА ПОВЕРКИ СИ
ЧИСЛО ОЦЕНОК 1
* ДОПОЛНИТЕЛЬНЫЕ ВОЗМОЖНОСТИ *
* О ДОПОЛН. ПРОГРАММ ПОЛЬЗОВАТЕЛЯ
* ОБЩАЯ ИНФОРМАЦИЯ
* ТАБЛИЦЫ РЕЗУЛЬТАТОВ И ИСХ. ДАННЫХ
* ГРАФИКОВ НЕ ПЕЧАТАТЬ
ВВОД ДАННЫХ = КОД ОШИБКИ = 1
НЕОБЯЗАТЕЛЬНЫЙ БЛОК ОТСУТСТВУЕТ
ОБЩАЯ ИНФОРМАЦИЯ О ЗАДАНИИ
РАСЧЕТ КРИТЕРИЕВ КАЧЕСТВА
ПРОВЕРКИ
ПОГРЕШНОСТИ МАНОМЕТРА МТИ-25
ТАБЛИЦА ИСХОДНЫХ ДАННЫХ
ПАРАМЕТР
1
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
0,150
ОТНОШЕНИЕ ГАММА
1,000
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0,800
ЦЕНА ДЕЛЕНИЯ ШКАЛЫ
0,200
НОМ. ЗНАЧЕНИЕ ВЫХ. СИГНАЛА ЗАДАТЧИКА
25,000
ПДЗ ВАРИАЦИИ
0,150
ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500
УСЛОВНЫЙ НОМЕР АЛГОРИТМА
1
ТАБЛИЦА РЕЗУЛЬТАТОВ
ПОКАЗАТЕЛЬ КАЧЕСТВА
1
МАКС. ВЕРОЯТНОСТЬ НЕОБНАРУЖ. БРАКА
0,512
МАКС. ВЫХОД ЗА ДОПУСК
1,160
СРЕДН. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,0
МАКС. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,0
ЗАДАНИЕ
2
ОЦЕНКА КАЧЕСТВА ПОВЕРКИ СИ
ЧИСЛО ОЦЕНОК 1
* ДОПОЛНИТЕЛЬНЫЕ ВОЗМОЖНОСТИ *
* О ДОПОЛН. ПРОГРАММ ПОЛЬЗОВАТЕЛЯ *
* ОБЩАЯ ИНФОРМАЦИЯ *
* ТАБЛИЦЫ РЕЗУЛЬТАТОВ И ИСХ. ДАННЫХ *
* ГРАФИКОВ НЕ ПЕЧАТАТЬ *
ВВОД ДАННЫХ - КОД ОШИБКИ = 1
НЕОБЯЗАТЕЛЬНЫЙ БЛОК ОТСУТСТВУЕТ
ОБЩАЯ ИНФОРМАЦИЯ О ЗАДАНИИ
РАСЧЕТ КРИТЕРИЕВ КАЧЕСТВА ПРОВЕРКИ
ВАРИАЦИИ МАНОМЕТРА МТИ-25
ТАБЛИЦА ИСХОДНЫХ ДАННЫХ
ПАРАМЕТР
1
ПДЗ КОНТРОЛИРУЕМОГО ПАРАМЕТРА
0,150
ОТНОШЕНИЕ ГАММА
1,000
ГРАНИЦА ФИКТИВНОГО БРАКОВАНИЯ
0,800
ЦЕНА ДЕЛЕНИЯ ШКАЛЫ
0,200
НОМ. ЗНАЧЕНИЕ ВЫХ. СИГНАЛА ЗАДАТЧИКА
25,000
ОБЪЕМ ВЫБОРКИ ДЛЯ ОЦЕНКИ ВЕРОЯТНОСТЕЙ
500
УСЛОВНЫЙ НОМЕР АЛГОРИТМА
1
ТАБЛИЦА РЕЗУЛЬТАТОВ
ПОКАЗАТЕЛЬ КАЧЕСТВА
1
МАКС. ВЕРОЯТНОСТЬ НЕОБНАРУЖ. БРАКА
0,552
МАКС. ВЫХОД ЗА ДОПУСК
1,200
СРЕДН. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,000
МАКС. ВЕРОЯТНОСТЬ ФИКТ. БРАКА
0,004