КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Общие положения теории планирования эксперимента
В ТПЭ исследуемый объект (реальный объект, модель объекта) рассматривается как "черный ящик", имеющий входы v (управляемые независимые параметры, или внутренние параметры) и выходы y [3, 6].По традиции входные управляемые переменные (внутренние параметры) обозначим: v = x1, x2,..., хn. Переменные v принято называть факторами. Теория ПЭ изучает только активный тип экспериментов, когда имеется возможность независимо и целенаправленно менять значения факторов v во всем требуемом диапазоне. Факторы в эксперименте бывают качественными и количественными. Качественные факторы можно квантифицировать или приписать им числовые обозначения, тем самым перейти к количественным значениям. В дальнейшем будем считать, что все факторы являются количественными и представлены непрерывными величинами (если другое не оговорено особо). Переменным v можно сопоставить геометрическое понятие факторного пространства – пространства, координатные оси которого соответствуют значениям факторов. Совокупность конкретных значений всех факторов образует точку в многомерном факторном пространстве. Примерами факторов являются: интенсивность потока запросов к базе данных, скорость передачи данных по каналу, объем запоминающего устройств. Кроме того, на объект воздействуют возмущающие факторы, они являются случайными и не поддаются управлению. Область планирования задается интервалами возможного изменения факторов xi,min < xi < xi,max для i =1, 2, …, n, где n =k – количество факторов. В теории ПЭ часто используют нормализацию факторов, т.е. преобразование натуральных значений факторов в безразмерные (нормированные) величины. Переход к безразмерным значениям xi задается преобразованием
где хi – натуральное (реальное) значение фактора, хiб = хi 0 – базовое значение основного уровня фактора, соответствующее нулю в безразмерной шкале, Δ хi – интервал (шаг) варьирования фактора хi. Выходные параметры объекта проектирования у1, у2,..., уm называются переменными состояния, или функцией отклика. Различают экономические и технологические переменные состояния. В качестве экономических используют производительность, себестоимость и другие показатели. Технологическими переменными служат качество продукта, выход целевого продукта, надежность получаемых изделий и др. Объект исследования может иметь несколько переменных состояния, их количество следует сократить до минимума. Опыт показывает, что в большинстве случаев удается ограничиться одной переменной состояния, и тогда вектор У превращается в скалярную величину у. Если переменных состояния несколько, то эксперимент проводится по каждой из них, а затем решается компромиссная задача. При выборе переменной состояния необходимо учитывать следующие требования: 1) переменная состояния должна иметь количественную характеристику, т. е. измеряться; 2) переменная состояния должна однозначно измерять эффективность объекта исследования; это требование эквивалентно корректной постановке задачи; 3) переменная состояния должна быть статистически эффективной, т. е. обладать возможно меньшей дисперсией при проведении опытов; это позволяет хорошо различать опыты. Правильный выбор переменной состояния объекта исследования повышает шансы экспериментатора на успех. При выборе факторов нужно выполнять следующие требования: 1) фактор должен быть регулируемым, т. е. с помощью определенного регулирующего устройства фактор можно изменять от значения x1 до значения x2 (системотехники называют это операционной определенностью). 2) точность измерения и управления факторов должна быть известна и достаточно высока (хотя бы на порядок выше точности измерения выходной переменной); очевидно, что низкая точность измерения факторов уменьшает возможности воспроизведения эксперимента. К факторам и переменным состояния одновременно также предъявляется ряд требований: 1) факторы и переменные состояния должны иметь области определения, заданные технологическими или принципиальными ограничениями; области определения факторов должны быть таковы, чтобы при различных их комбинациях переменные состояния не выходили за свои ограничения; 2) между факторами и переменными состояниями должно существовать однозначное соответствие; оно позволит в основном эксперименте построить математическую модель объекта исследования и решить поставленную задачу эксперимента. Метод наименьших квадратов позволяет получить описание объекта по любым данным, лишь бы матрица системы нормальных уравнений была невырожденной. Поэтому с появлением ЭВМ возникла идея — получать математические описания технологических процессов, пользуясь в качестве исходных данных результатами нормальной эксплуатации процесса. В реальных условиях технологический процесс все время испытывает случайные колебания режима. Сегодня значения контролируемых факторов несколько иные, чем вчера, а завтра будут еще немного другими. Нельзя ли каждое изменение режима рассматривать как эксперимент, и, обработав совокупность таких «экспериментов» методом наименьших квадратов, получить описание процесса, а затем использовать это описание для управления и оптимизации? Такой подход получил название пассивного эксперимента. Достоинство пассивного эксперимента—отсутствие затрат на опыты: данные получаются «сами собой». Но надежды, возлагавшиеся на этот метод, в большинстве случаев не оправдались. Анализ неудач пассивного эксперимента выявил несколько их причин. Во-первых, в нормальных условиях колебания режима малы, опытные точки находятся близко одна к другой. Хорошо известно, что чем ближе опытные точки, тем сильнее влияют на описание случайные ошибки. Действительно, различия в получаемых значениях отклика при этом малы, и эти малые различия плохо выделяются на фоне шума — случайных ошибок. Поэтому значения коэффициентов регрессии оцениваются со значительными ошибками. Во-вторых, в пассивном эксперименте факторы сильно коррелированны. Это делает крайне ненадежным анализ влияния отдельных факторов—всегда может оказаться, что влияет не данный фактор, а другой, с ним коррелированный. В-третьих, сами значения факторов в производственных условиях часто измеряются с заметными ошибками; поэтому применение метода наименьших квадратов в его обычном варианте становится некорректным. В связи с этим в теории эксперимента любой эксперимент, при планировании которого не учтено влияние плана эксперимента на статистические свойства получаемых оценок, часто называют пассивным. Ему противопоставляют активный эксперимент, в основе которого лежит планирование эксперимента (то есть управление условиями его проведения) Наиболее широко в настоящее время используются планы так называемого экстремального эксперимента, разработанные для определения оптимальных условий протекания процессов в объектах исследования. Оптимум определяется по математической модели объекта исследования, которую ищут в виде полиномиального уравнения (уравнения регрессии). Логику появления полинома как математической модели объекта исследования можно объяснить следующим образом. Исследователь полагает, что математическую модель объекта принципиально можно представить дифферен-циальными уравнениями. В общем виде искомое решение можно представить функцией: y= F (X,b), (1.2) где у— переменная состояния объекта исследования; Х— матрица факторов; b — матрица коэффициентов. Коэффициенты b полинома можно интерпретировать как коэффициенты ряда Тейлора, в который «удается» разложить решение в окрестностях некоторой точки. Пользуясь статистическими методами и учитывая конечность экспериментальных данных, можно получить оценки коэффициентов регрессии b в уравнении (1). Уравнение (1) называют уравнением регрессии и широко используют для получения математической модели объекта исследования. Прежде чем приступить к планированию эксперимента, необходимо убедиться в том, что опыты воспроизводимы, Для этой цели проводят N серий из k параллельных опытов в рассматриваемой области изменения влияющих факторов. Результаты опытов заносятся в таблицу: Таблица 1.1
Определяется среднее арифметическое значение функции отклика для любой серии опытов j (j=1,…,N):
Оценивается дисперсия для каждой серии параллельных опытов:
Определяется расчётное значение критерия Кохрена: .
По специальным таблицам определяют табличное значение Кохрена - GT. Они зависят от доверительной вероятности P(как правило P=0.95),от N и от f=k-1.
Если выполняется условие Gp GT, то опыты считаются воспроизводимым Определяется погрешность эксперимента. Оценка дисперсии воспроизводимости рассчитывается по формуле: , с ней связано число степеней свободы N*(k-1) В последовательности реализации планов можно выделить следующие этапы: 1) оценка априорной информации и отсеивание факторов, несущественных для конкретного объекта исследования; 2) получение математической модели объекта в виде линейной функции отклика; 3) поиск оптимальной области объекта по линейной функции отклика; 4) получение математической модели объекта исследования области оптимума в виде нелинейной функции отклика; 5) поиск оптимальной координаты факторного пространства в области оптимума. Экстремальный эксперимент заканчивается определением оптимальной координаты, соответствующей оптимальным условиям протекания процессов объекта исследования. В подавляющем большинстве РЭС и технологические процессы их изготовления являются сложными, на процесс влияет не один, а ряд факторов. Возможны два подхода к исследованию таких многофакторных систем. Первый можно описать формулой: «Изменяй факторы по одному». Исследование системы разбивается на серии, в пределах каждой из которых изменяется (варьируется) лишь один фактор, а остальные неизменны. В следующей серии изменяется второй фактор и т. д. Идея другого подхода— построить план эксперимента, предусматривающий изменение всех влияющих факторов, с тем, чтобы этот план обеспечивал максимум точности, минимум корреляции и другие хорошие статистические свойства. Такой эксперимент называют многофакторным. Долгое время в науке господствовал первый подход. Его главное преимущество—наглядность: данные каждой серии легко поддаются интерпретации. Один из важных результатов теории планирования эксперимента заключается в том, что второй подход—значительно эффективнее первого. При том же объеме эксперимента и той же точности опытов получается большая точность результатов. Геометрическим образом совокупности независимых переменных х и зависимой переменной у является пространство n+1 измерения, где n — число независимых переменных; (n+1)-е измерение относится к у. В этом пространстве зависимости у от всех х соответствует n-мерная поверхность, которую обычно называют поверхностью отклика (результат опыта рассматривается как отклик системы на опыт—заданную совокупность независимых переменных, или входов). План эксперимента указывает расположение опытных точек в n-мерном пространстве независимых переменных (факторном пространстве), или иными словами, условия всех опытов, которые следует провести. Чаще всего план эксперимента задается в виде матрицы планирования — прямоугольной таблицы, каждая строка которой отвечает условиям определенного опыта, а каждый столбец — значениям какой-то из независимых переменных в разных опытах. Полным факторным экспериментом называется система опытов, содержащая все возможные неповторяющиеся комбинации уровней варьирования факторов.Для удобства вычислений коэффициентов регрессии все факторы в ходе полного факторного эксперимента варьируют на двух уровнях: нижнем -1 и верхнем +1, соответствующих значениям нормированных переменных X1,X2,……Xn.Таким образом, полным факторным экспериментом называется система опытов, содержащая все возможные неповторяющиеся комбинации уровней варьирования факторов. В табл. 1.2 приведены условия опытов полного трехфакторного эксперимента. Эти опыты соответствуют в факторном пространстве вершинам куба с центром в начале координат. Таблица 1.2
Каждый фактор принимает лишь два значения — варьируется на двух уровнях, верхнем и нижнем. Поэтому общее число экспериментов N=2 n. Пример. В четырех опытах исследуется влияние 3 факторов:температуры T,К, давления р, МПа, и времени t, с, на выход продукта. Здесь в любом из опытов температура—либо 1000 К (нижний уровень), либо 1200 К (верхний уровень); аналогично варьируются р и t. Таблица 1.3
Матрица планирования эксперимента для этого случая может иметь вид: Таблица 1.4
Из таблицы видны основные принципы построения матриц планирования полного факторного эксперимента: Матрица планирования ПФЭ обладает следующими свойствами:
где j- номер опыта; i— номер фактора,(l¹m) Свойство, выраженное последним уравнением, называется ортогональностью матрицы. Оно позволяет вычислять коэффициенты регрессии по простым формулам независимо друг от друга Расширенная матрица—это матрица, дополненная столбцами, учитывающими взаимодействия факторов. На практике, как правило, ограничиваются парными взаимодействиями. Коэффициенты регрессии рассчитываются методом наименьших квадратов. Основное условие метода формулируется следующим образом: коэффициенты регрессии определяются на основании минимизации суммы квадратов отклонений между экспериментальными уэ, и рассчитанными по уравнению регрессии yр значениями функции отклика:
После определения коэффициентов регрессии определяем значимость этих коэффициентов. Все коэффициенты подразделяются на значимые и незначимые. Для определения значимости коэффициентов регрессии сравнивается погрешность вычисления коэффициента с погрешностью экспериментальных данных - . Вычисляется доверительный интервал: . Здесь tT - табличное значение критерия Стьюдента, которое находится по числу степеней свободы и доверительной вероятности. Тогда значимость оценивают, сравнивая абсолютные значения коэффициента и доверительного интервала. .Незначимые коэффициенты отбрасываются из уравнения регрессии, после чего записывается окончательный вид уравнения регрессии. Это уравнение проверяется на адекватность. Для этого вычисляется оценка дисперсии адекватности:
Здесь B- число значимых коэффициентов регрессии. Вычисляют расчётное значение критерия Фишера: По таблице находят табличное значение критерия Фишера. Оно зависит от доверительной вероятности P, числа степеней свободы fад= N-B и f= N*(k-1). На основании этого делается вывод об адекватности или неадекватности уравнения регрессии. Уравнение регрессии считается адекватным, если выполняется условие: Fр FT. Чем меньше B, тем больше N-B—в этом одна из главных целей, достигаемых при исключении незначимых членов. Если уравнение неадекватно, переходят к более сложной модели (например, повышают степень многочлена), для чего обычно требуется постановка добавочных опытов. Иногда можно обойтись без дополнительного эксперимента, если соответствующим образом преобразовать переменные у или х. Совокупность основных уровней всех факторов представляет собой точку в пространстве параметров, называемую центральной точкой плана или центром эксперимента. С геометрической точки зрения нормализация факторов равноценна линейному преобразованию пространства факторов, при котором проводятся две операции: перенос начала координат в точку, соответствующую значениям основных уровней факторов; сжатие – растяжение пространства в направлении координатных осей. Активный эксперимент включает: систему воздействий, при которых воспроизводится функционирование объекта; регистрацию отклика объекта. План эксперимента задает совокупность данных, определяющих количество, условия и порядок реализации опытов. Опыт составляет элементарную часть эксперимента и предусматривает воспроизведение исследуемого явления в конкретных условиях с последующей регистрацией результата. В условиях случайности в одних и тех же условиях проводятся параллельные (повторные) опыты в интересах получения статистически устойчивых результатов. Опыт u предполагает задание конкретных значений факторам v u = v 1 u , v 2 u , …, vku, а совокупность значений факторов во всех N точках плана эксперимента образует матрицу плана
Строки матрицы соответствуют опытам, столбцы – факторам, элемент матрицы viz задает значение z- го фактора в i- м опыте. Вектор y называется откликом. В ТПЭ обычно изучается ситуация, в которой вектор отклика y состоит из одного элемента y. При наличии нескольких составляющих вектора y, каждую из них можно исследовать отдельно. Зависимость отклика от факторов носит название функции отклика, а геометрическое представление функции отклика – поверхности отклика. Функция отклика рассматривается как показатель качества или эффективности объекта. Этот показатель является функцией от параметров – факторов. На практике широкое распространение получили простые функции вида М{ y' } = bf (v), где b= (b 0, b 1, …, bh ) – вектор неизвестных параметров модели размерности h +1, f (v)=(f 0(v), f 1(v), …, fh (v)) – вектор заданных базисных функций, М{ y' } – математическое ожидание функции отклика. Такое представление функции отклика соответствует линейной по параметрам модели регрессионного анализа, т.е. функция отклика есть линейная комбинация базисных функций от факторов. Вследствие влияния на результаты экспериментов случайных воздействий истинные значения коэффициентов можно определить только приближенно. Оценку b = (b0, b1, …, b h) вектора неизвестных параметров b находят по результатам экспериментов, в ходе которых получают значения yu при заданных значениях факторов v u. Эти оценки обычно рассчитываются с помощью метода наименьших квадратов (МНК) на основе выборок значений факторов и откликов системы на воздействия [8]. В качестве оценки b вектора b выбирается такое значение, которое минимизирует , где y'u – вычисленное на модели значение функции отклика в u -й точке факторного пространства. Приравнивая нулю частные производные от данной квадратичной формы, взятые по переменным b0, b1, …, b h, можно получить систему уравнений вида , где i = 0, 1, 2, …, h. Значение b находят путем решения этой системы уравнений. Решение системы возможно при линейной независимости базисных функций. Если не принимать специальных мер, то оценки коэффициентов b станут взаимозависимыми, и полученное выражение для функции отклика можно рассматривать только как интерполяционную формулу, что затрудняет ее физическую интерпретацию и последующие расчеты. Однако, формируя специальным образом матрицу плана, можно получить независимые значения b. И эти величины будут характеризовать вклад каждого фактора в значение функции отклика. Итак, задача заключается в определении общей формы записи функции отклика y'. В большинстве случаев вид этой функции, получаемый из теоретических соображений, является сложным для практического применения, а при неполном знании объекта вообще неизвестен. По данным причинам функцию целесообразно представить в универсальном, удобном для практического применения виде, чему соответствует представление в виде полинома. Тогда системой базисных функций является совокупность степенных функций с целыми неотрицательными значениями показателей степени. Полиномиальная форма представления функции отклика примет вид
где e – случайная величина, характеризующая ошибку опыта. Такая функция отклика линейна относительно неизвестных коэффициентов и будет полностью определена, если заданы степень полинома и коэффициенты. Степень полинома обычно задается исследователем априорно и уточняется в ходе исследования. На практике наибольшее распространение получили полиномы первого и второго порядка, соответственно линейные и квадратичные модели. Коэффициенты полинома принято называть эффектами факторов. Иногда функцию отклика целесообразно представить в другом виде, например, в виде степенной функции, так как достижение заданной точности требует применения полинома высокого порядка. Однако использование функций, нелинейных относительно неизвестных параметров, усложняет вычисления, затрудняет оценку их свойств. В некоторых случаях задачу можно упростить путем искусственного преобразования нелинейной функции в линейную. При этом требуется соответствующее преобразование и результатов экспериментов. Применение ТПЭ основано на ряде допущений, а именно [2, 6]: функция отклика содержит в своем составе неслучайную и случайную составляющую. Многие показатели качества автоматизированных систем обработки информации носят случайный характер. Это требует многократного повторения опытов в одних и тех же условиях в целях получения статистически устойчивых результатов, а получаемые оценки показателей должны обладать свойствами состоятельности, эффективности, несмещенности и достаточности. Оценки типовых показателей формируются путем усреднения результатов наблюдений. Поэтому при достаточно большом количестве наблюдений можно считать, что случайная составляющая e распределена по нормальному закону с нулевым математическим ожиданием, что позволяет получить несмещенную оценку математического ожидания функции отклика в конкретной точке плана. Будем также считать, что величина e имеет дисперсию, не зависящую от значений факторов. Иначе говоря, результаты, полученные путем усреднения повторных опытов в каждой точке плана, представляют собой независимые, нормально распределенные случайные величины; факторы v 1, v 2, …, vk измеряются с пренебрежимо малой ошибкой по сравнению с ошибкой в определении величины y (учет помех в задании факторов приводит к трудно разрешимым проблемам в оценке коэффициентов функции отклика). Ошибка в определении значения функции отклика объясняется не столько погрешностью измерений, сколько влиянием на результат работы системы неучтенных или случайных факторов, например различиями в формируемой последовательности случайных чисел при статистическом моделировании; дисперсии среднего значения функции отклика в различных точках равны друг другу (выборочные оценки дисперсии однородны). Это означает, что при многократных повторных наблюдениях над величиной yu при некотором наборе значений v 1 u , v 2 u , …, vku, получаемая оценка дисперсии среднего значения не будет отличаться от оценки дисперсии, полученной при многократных наблюдениях для любого другого набора значений независимых переменных v 1 s , v 2s, …, vks. Указанные допущения позволяют использовать для расчетов коэффициентов полинома МНК, который дает эффективные и несмещенные оценки коэффициентов и обеспечивает простоту проведения самих расчетов. Применение МНК, вообще говоря, не требует соблюдения нормального распределения результатов наблюдения. Этот метод в любом случае дает решение, минимизирующее сумму квадратов отклонений результатов наблюдения от значений функции отклика. Допущение о нормальном распределении используется при проведении различного рода проверок, например, при проверке адекватности функции отклика и экспериментальных данных. Естественно, что точность оценок коэффициентов функции отклика повышается с увеличением числа опытов, по которым вычисляются коэффициенты.
Дата добавления: 2014-12-27; Просмотров: 1182; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |