Московский государственный открытый университет им. В.С. Черномырдина
Опубликован: 20.01.2011 | Доступ: свободный | Студентов: 2310 / 498 | Оценка: 4.14 / 3.69 | Длительность: 27:06:00
Лекция 4:

Разработка "сквозных" математических моделей технологического процесса как основа АСНИ

< Лекция 3 || Лекция 4: 12 || Лекция 5 >
Аннотация: Рассматривается методика построения "сквозных" математических моделей на основе статистических методов. Такие модели должны разрабатываться в самом начале проектирования, на этапе научно-исследовательских работ (НИР), и должны составлять основу автоматизированной системы научных исследований (АСНИ).
Ключевые слова: технологическая модель, автоматизированный технологический комплекс, математическая модель, средства автоматизации, очередь, структурная схема, технологический процесс, контролируемые параметры, контрольные точки, Построение математической модели, транзистор, показатели надежности, вектор, операции, информация, регистрация, количество информации, факторный анализ, обобщенные факторы, информативные параметры, математическая модель системы, регрессионный анализ, анализ, техническое задание, внешнее проектирование, статистические методы, грубые погрешности, погрешность, нормальное распределение, определение, гипотеза, эмпирическая формула, метода наименьших квадратов, оценка математического ожидания, выборочной средней, дисперсия, стандартное отклонение, коэффициент вариации, парная корреляция, доверительная вероятность, неравенство, Абсолютной величиной, распределение информации, активный эксперимент, объем выборки, нормальный закон, гипотеза нормального распределения исходной информации, обратное преобразование, преобразование данных, Гистограмма, нормализация, внутреннее проектирование

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

4.1. Исходные данные для построения математических моделей

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

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

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

На рис. 4.1 и 4.2 представлены такие схемы для технологических процессов изготовления видеоконтрольного устройства и транзисторов.

Структурная схема технологического процесса изготовления видеоконтрольного устройства (ВКУ)

Рис. 4.1. Структурная схема технологического процесса изготовления видеоконтрольного устройства (ВКУ)

Здесь: \overline{X1} - совокупность параметров комплектующих элементов;

\overline{X2} - совокупность контролируемых параметров модулей ВКУ и показателей надежности собранного прибора до регулировки;

\overline{X3} - совокупность показателей надежности готового ВКУ после регулировки;

\overline{P} - вектор выходных показателей надежности после испытаний;

\overline{Z1}, \overline{Z2}, \overline{Z3} - векторы неконтролируемых возмущений, воздействующих на технологические операции, соответственно, сборки, регулировки и испытаний, например, неконтролируемые параметры технологического оборудования, внешней среды;

\overline{U1}, \overline{U2}, \overline{U3} - векторы контролируемых воздействий на технологические операции сборки, регулировки испытаний, например, подаваемые напряжения, токи.

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

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

Структурная схема технологического процесса изготовления транзистора

увеличить изображение
Рис. 4.2. Структурная схема технологического процесса изготовления транзистора

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

Основными рекомендациями, выдаваемыми в результате проведения НИР и ОКР, должны быть перечень наиболее информативных ( контролируемых и регулируемых ) параметров, математические модели и алгоритмы управления, эскиз проектируемой системы.

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

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

Такой анализ исходной информации целесообразно проводить на самом первом этапе проектирования - формировании технического задания.

Итак, первый этап проектирования системы - внешнее проектирование - содержит в себе синтез, анализ, оптимизацию. Из ряда элементных баз выбирается та, которая бы удовлетворяла поставленным требованиям; к началу проектирования выбирается или синтезируется решение, отвечающее поставленным целям.

Так как цель проектирования - разработка не только допустимого варианта проекта, но и наиболее выгодного в техническом и экономическом отношении, то далее решается задача оптимизации проекта.

4.2. Предварительная обработка исходной информации: сущность и необходимость

Предварительная обработка результатов измерений или наблюдений необходима для того, чтобы в дальнейшем с наибольшей эффективностью, а главное - корректно, использовать для построения эмпирических зависимостей статистические методы.

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

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

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

Можно встретить указание, что аномальные значения измеряемой величины получаются в результате изменения условий эксперимента, однако это определение неполное.

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

Используя введенные обозначения (рис. 3.2), статистические данные о технологическом процессе и качестве продукции можно сгруппировать в виде двух матриц наблюдений \overline{Ч} и \overline{C}:

\overline{Ч}= (х_{\eta i});\,\, 
\overline{C}= (х_{\eta j}) ( 4.1)

где

  • х_{\eta i} - значение i -го параметра в \eta  - цикле;
  • \eta - номер цикла;
  • p_{\eta j} - значение j -го признака качества в \eta - цикле;
  • \eta  = 1, 2,…, n, n - число наблюдений технологических циклов;
  • i = 1, 2,…, k, k - количество параметров процесса;
  • j = 1, 2,…, l, l - количество признаков готового прибора.

Первый этап математического описания - статистические оценки характеристик генеральной совокупности (всех технологических циклов) по выборке.

Оценка математического ожидания - выборочное среднее арифметическое

\hat{x_i} = \cfrac{1}{n}\sum\limits_{\eta=1}^{n}{x_{\eta i}};\;\;\;
\hat{p_i} = \cfrac{1}{n}\sum\limits_{\eta=1}^{n}{p_{\eta i}} ( 4.2)

Оценка дисперсии - выборочная дисперсия

\hat{Д}_{xi} = \hat{S}_{xi}^2= \cfrac{1}{n-1}\sum\limits_{\eta=1}^{n}{(x_{\eta i}-\hat{ x}_{i})^2};\\
\hat{Д}_{pj} = \hat{S}_{pj}^2= \cfrac{1}{n-1}\sum\limits_{\eta=1}^{n}{(p_{\eta j}-\hat{p}_{j})^2} ( 4.3)

Оценка стандартного отклонения - выборочное среднеквадратиче-ское отклонение

\hat{S}_{xi}=+\sqrt{\hat{Д}_{xi}};\;\;\;
\hat{S}_{pj}=+\sqrt{\hat{Д}_{pj}} ( 4.4)

Оценка коэффициента вариации - выборочный коэффициент вариации

\hat{V}_{xi} =\cfrac{\hat{S}_{xi}}{\hat{x}_i};\;\;\;
\hat{V}_{pj} =\cfrac{\hat{S}_{pj}}{\hat{p}_j};\;\;\; ( 4.5)

Оценка коэффициентов парной корреляции - выборочный коэффициент парной корреляции

\hat{\rho}(x_i,x_\xi) =\cfrac{\sum\limits_{\eta=1}^{n}{(x_{\eta i} -\hat{x}_i)}\cdot(x_{\eta\xi} -\hat{x}_\xi)}
{(n-1)\hat{S}_{xi}\hat{S}_{x\xi}}; \\
\hat{\rho}(x_i,p_j) =\cfrac{\sum\limits_{\eta=1}^{n}{(x_{\eta i} -\hat{x}_i)}\cdot(p_{\eta j} -\hat{p}_j)}
{(n-1)\hat{S}_{xi}\hat{S}_{pj}}; \\
\hat{\rho}(p_\nu,p_j) =\cfrac{\sum\limits_{\eta=1}^{n}{(p_{\eta \nu} -\hat{p}_\nu)}\cdot(p_{\eta j} -\hat{p}_j)}
{(n-1)\hat{S}_{p\nu}\hat{S}_{pj}};\\ ( 4.5)

где {\xi} = 1, 2, …, k; v = 1, 2, …, l.

Отсев грубых погрешностей

Можно встретить большое количество различных рекомендаций для проведения отсева грубых погрешностей наблюдения (аномальных значений). Рассмотрим наиболее простые методы отсева грубых погрешностей. Если в распоряжении экспериментатора имеется выборка небольшого объема n \le 25, то можно воспользоваться методом вычисления максимального относительного отклонения:

|x_i - \overline{x}| / \overline{S} \le \tau_{1-p} ( 4.7)

где x_{i} - крайний (наибольший или наименьший) элемент выборки, по которой подсчитывались \overline{x} и \overline{S}, \tau_{1-p} - табличное значение статистики \tau, вычисленной при доверительной вероятности q = 1 - p.

Таким образом, для выделения аномального значения вычисляют

\tau = | x_{i}-\overline{x} | /S ( 4.8)

которое затем сравнивают с табличным значением \tau_{1-p}.

\tau\le\tau_{1-p} ( 4.9)

Если это неравенство соблюдается, то наблюдение не отсеивают; если не соблюдается, то наблюдение исключают. После исключения того или иного наблюдения или нескольких наблюдений характеристики эмпирического распределения должны быть перечитаны по данным сокращенной выборки.

Квантили распределения статистики t_\beta берут при уровнях значимости р = 0,10, р = 0,05, р=0,025, р = 0,01 или доверительной вероятности 1 - р = q = 0,90, 0,95, 0,975, 0,99. На практике обычно используют уровень значимости р = 0,05 (результат получается с доверительной вероятностью 95%).

Процедуру отсева можно повторить и для следующего по абсолютной величине максимального относительного отклонения, но предварительно необходимо пересчитать \overline{x} и \overline{S} для выборки нового объема n - 1.

Рассмотрим другой метод отсева грубых погрешностей для малой выборки. В этом случае вычисляют

\tau' = \cfrac{|x_i-\overline{x}|}
{\sqrt{(n-1)/n}\cdot \overline{S}} ( 4.10)

и полученный результат сравнивают с критическим значением, взятым из таблицы при соответствующих n и 1 - р.

Отсев грубых погрешностей можно провести и для больших выборок. Для практических целей лучше всего использовать таблицы распределения Стьюдента. Этот метод исключения аномальных значений для выборок большого объема отличается простотой, а таблицы распределения Стьюдента имеются практически в любой книге по математической статистике. Распределение Стьюдента относится к категории распределений, связанных с нормальным распределением. Подробно эти распределения рассмотрены в учебниках по математической статистике.

Известно, что критическое значение \tau_p ( р - процентная точка нормированного выборочного отклонения) выражается через критическое значение распределения Стьюдента t_{p,n-2}

\tau_(p,n) = \cfrac{ t_{(p,n-2)}\sqrt{n-1}}{\sqrt{n-2+[ t_{(p,n-2)}]^2}} ( 4.11)

На практике пользуются достаточно удобным методом отсева грубых погрешностей. Если рассматривать кривую нормального распределения (рис. 4.3), то это те значения, которые выходят за пределы \sigma, 2\sigma, 3\sigma в зависимости от объема выборки и доверительной вероятности. Они просто отбрасываются по "правилу 6 сигма" (на рисунке представлена половина графика).

Отсев грубых ошибок по правилу о, 2а, 3а.

Рис. 4.3. Отсев грубых ошибок по правилу о, 2а, 3а.

Например, для объема выборки n>50 измерений их отсев будет осуществляться при выходе за пределы:

-3\sigma  \text{ для }\alpha= 0,99; \\
-2\sigma  \text{ для }\alpha= 0,95; \\
-\sigma\text{ для }\alpha= 0,67.
< Лекция 3 || Лекция 4: 12 || Лекция 5 >