дого из времен и изделий производится расчет значений выходного параметра. На печать выдаются таблицы предсказанных значений и сообщения о том, находится ли данное значение внутри поля допуска. Кроме того, рассчитывается вероятность безотказной работы для данного времени.

Во втором режиме для всех времен контроля рассчитываются таблицы, содержащие истинные и предсказанные значения выходного параметра, а также погрешности предсказания. Для каждого времени рассчитываются средняя и максимальная погрешности предсказания и среднее квадратическое отклонение предсказанного значения от истинного, строится гистограмма распределения абсолютных значений погрешностей предсказания, отнесенных к полю допуска.

Дополнительно в пакет включены еще две программы, одна из которых осуществляет построение вариационных рядов, а вторая - предсказание и расчет временных зависимостей основных количественных характеристик надежности (время безотказной работы, частоту, интенсивность отказов, гамма-процентный ресурс) по группе выходных параметров.

Пакет программ STAB реализован на ЕС ЭВМ, написан на языке ПЛ/1. Объем оперативной памяти не превышает 130 К байт для оверлейного и 200 К байт для быстродействующего варианта пакета. Время построения математической модели оптимальной сложности пропорционально квадрату числа параметров и количеству изделий в эмпирической таблице и составляет около 40 мин процессорного времени ЭВМ ЕС-1033 для 50 изделий, семи параметров и трех уровней селекции: время работы программ ввода и предсказания составляют 45 и 30 с соответственно.

На рис. 2.41 представлен ряд временных зависимостей отклонения сопротивления резистора, нормированного к номинальному значению от приведенного времени, штриховой линией показаны предсказанные значения.

При решении практических задач предсказания было замечено, что точность предсказания существенно зависит от вида зависимостей выходного параметра от конструктивных параметров. Так, при работе в импульсном режиме и длительности импульса, прикладываемого к резистору, менее критического значения тр, соответствующего максимальной температуре Гдоп в РЭ, резистор работает в недогруженном режиме. При превьппении этого значения деградация резистивного материала резко ускоряется, что позволяет вьщелить две области для множества допустимых значений длительности импульса. Поэтому разобьем ноштель пространства входных параметров на ряд нечетких областей и для каждой области будем строить отдельную модель. Пример такого разбиения приведен на рис. 2.42, а и 6.

Пусть носитель каждого из входных параметров д., к = 1, К разбит нечеткими множествами (g.) на «д. областей, Ид; > 1; каж-

Vuc 2.41. Временные зависимости отклонения сопротивления от номинального для четырех пепроволочных силовых резисторов

0,2 0,4 0,Б 0,6 t

J/f/Л,

Рис. 2.42. Разбиение входного пространства на области для повышения точности математической модели изменения conpoTHsneHMflipesHCTopa:

а - временные зависимости отклонения сопротивления силового резистора от номинального значения, Т - время (7 - изменение сопротивления в режиме недогрузки; 2 - то же в режиме перегрузки); б - разбиение носителя параметра длительности импульса (Ti) с помощью двух нечетких множеств (i - соответствует нечеткому множеству силовой резистор недогружен; 2 - то же силовой резистор перегружен); в - зависимость средней скорости изменения сопротивления от длительности импульса (/ - режим недогрузки; 2 - режим перегрузки)

дой области соответствует модель

Model (а, q, т, к, п), значение выходного параметра R(q, т) определяется по формуле 2 7fc,nModel(a,q,T,A:,n)

R(q. г) = ----, (2.98)

ft,я

Tft я = i""» maxju, „ (q), к п

Следует отметить, что использование разбиения пространства входных параметров обеспечивает "мягкую" аппроксимацию выходных параметров, так, на рис. 2.42, в приведена зависимость средней скорости изменения сопротивления от длительности tj .

Построение моделей вида (2.98) можно осуществить с помощью МНК. При увеличении числа разбиений число частных моделей, равное к

П щ, быстро растет. Разбиение пространства входных параметров ft = l

рекомендуется производить либо на основе экспертных оценок, либо с помощью МГУА, вводя внещний цикл, осуществляющий подбор

„ (q) в процедуру построения модели оптималыюй сложности.

Значительного сокращения перебора частных моделей можно достигнуть вводом в программу эвристической информации о недопустимых сочетаниях значений параметров и областей их значений.

Оперативно, хотя и с меньшей точностью, информацию о деградации параметров силовых резисторов во время эксплуатации можно получить с помощью рассмотренного ранее алгоритма индуктивной интерполяции (ИИ). Алгоритм ИИ реализован в виде диалогового пакета прикладных программ (ДППП) Proect*, который функционирует на ЕС ЭВМ в среде диалоговой системы ПРИМУС. Ддя облегчения работы неподготовленного пользователя ДППП вьтолнен на механизме меню, входные данные и результаты вычислений представляются в процессе работы в квазиграфической форме. Зависимости между параметрами определяются в диалоге; выходные параметры оцениваются как с помощью ИИ, так и с помощью дшамически загружаемых программ пользователя.

*В разработке ДППП участвовал А.А. Березкин.