Учебное пособие 1324
.pdfМИНИСТЕРСТВО НАУКИ И ВЫСШЕГО ОБРАЗОВАНИЯ РОССИЙСКОЙ ФЕДЕРАЦИИ
Федеральное государственное бюджетное образовательное учреждение высшего образования
«Воронежский государственный технический университет»
Кафедра конструирования и производства радиоаппаратуры
СИНТЕЗ МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ
МЕТОДИЧЕСКИЕ УКАЗАНИЯ
к выполнению курсовой работы по дисциплине «Основы САПР»
для студентов направления 11.03.03 «Конструирование
и технология электронных средств» (профиль «Проектирование и технология радиоэлектронных средств») всех форм обучения
Воронеж 2021
УДК 519.6:519.193 ББК 22.193:22.311я7
Составитель ст. преп. О. Н. Чирков
Синтез математических моделей: методические указания к выполнению курсовой работы по дисциплине «Основы САПР» для студентов направления 11.03.03 «Конструирование и технология электронных средств» профиль («Проектирование и технология радиоэлектронных средств») всех форм обучения / ФГБОУ ВО «Воронежский государственный технический университет»; cост.: О. Н. Чирков. Воронеж: Изд-во ВГТУ, 2021. 41 с.
Основной целью является приобретение студентами знаний, умений и практических навыков по синтезу адекватных математических моделей на основе уравнения регрессии в виде полноквадратичной зависимости.
Предназначены для выполнения курсовой работы по дисциплине «Основы САПР» для студентов 3 курса.
Методические указания подготовлены в электронном виде и содержатся в файле KURS_SAPR1.pdf.
Библиогр.: 7 назв.
УДК 519.6:519.193 ББК 22.193:22.311я7
Рецензент - О. Ю. Макаров, д-р техн. наук, проф. кафедры конструирования и производства радиоаппаратуры ВГТУ
Издается по решению редакционно-издательского совета Воронежского государственного технического университета
ПОЛУЧЕНИЕ МОДЕЛИ МЕТОДОМ ОЦКП
Исходными данными являются базовые значения факторов (число факторов k=n=3) и шаги варьирования. Задана матрица планирования эксперимента и результаты трёх дублирующих эскпериментов (для каждого эксперимента проведено 3 дублирующих опыта, n=3 – количество факторов, m=3 – количество дублирующих опытов). Общее количество экспериментов в методе ортогонального центрального композиционного планирования
N = 2n + 2n +1
Обозначим L – порядковый номер эксперимента, L = 1,…,N . В случае трёхфакторного эксперимента N=15 (15 экспериментов). Результаты всех опытов запишем в виде матрицы размерности 15х3, обозначим её элементыY lj, где l- номер эксперимента, а j-номер дублирующего опыта.
Результаты опытов
Y11Y12Y13
Y21Y22Y23
Y15,1Y15,2Y15,3
Находим среднее значение в каждой серии опытов: Базовые значения факторов (от 0 до 10)
YL = 1 ∑3 YLJ ,L =1,15
3 j =1
Y1 = 13 (Y11 +Y12 +Y13 ) = 13 (34,1+34+34,1)
Y2 = 13 (Y21 +Y22 +Y23) = 34,06 ≈ 34,1
Шаги варьирования факторов (от 0 до 1)
3
Исходные данные для ОЦКП |
|
|
|
|||||||
|
|
d |
= 0.1 |
X61 |
=1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
d2 = 0.2 |
X62 |
= 2 |
|
|
|
|
|
|
|
|
d3 |
= 0.3 |
X63 |
= 3 |
|
|
|
|
|
L |
|
X1 |
|
X2 |
|
X3 |
|
Yl1 |
Yl2 |
Yl3 |
|
|
|
|
|
|
|
|
|
|
|
l1 |
|
|
|
|
|
|
|
|
|
|
l2 |
|
-1.000 |
-1.000 |
|
-1.000 |
|
34.1 |
34.0 |
34.1 |
|
l3 |
|
|
|
|
|
|
|
|
|
|
l4 |
|
1.000 |
-1.000 |
|
-1.000 |
|
49.4 |
34.0 |
34.0 |
|
l5 |
|
|
|
|
|
|
|
|
49.4 |
34.0 |
l6 |
|
-1.000 |
1.000 |
|
-1.000 |
|
34.0 |
34.0 |
49.4 |
|
l7 |
|
|
|
|
|
|
|
|
49.4 |
34.0 |
l8 |
|
-1.000 |
-1.000 |
|
1.000 |
|
34.1 |
49.4 |
49.4 |
|
l9 |
|
|
|
|
|
|
|
49.4 |
34.0 |
49.4 |
l10 |
|
1.000 |
1.000 |
|
-1.000 |
|
34.1 |
49.4 |
49.4 |
|
l11 |
|
|
|
|
|
|
|
49.4 |
41.3 |
41.4 |
l12 |
|
-1.000 |
1.000 |
|
1.000 |
|
49.4 |
41.3 |
41.3 |
|
l13 |
|
|
|
|
|
|
|
41.4 |
41.4 |
41.3 |
l14 |
|
1.000 |
-1.000 |
|
1.000 |
|
51.2 |
51.2 |
41.3 |
|
l15 |
|
|
|
|
|
|
|
32.6 |
32.6 |
41.3 |
|
|
1.000 |
1.000 |
|
1.000 |
|
41.3 |
41.4 |
51.2 |
|
|
|
|
|
|
|
|
|
41.4 |
41.3 |
32.6 |
|
|
0.000 |
0.000 |
|
0.000 |
|
41.4 |
|
|
|
|
|
1.215 |
0.000 |
|
0.000 |
|
41.3 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
-1.215 |
0.000 |
|
0.000 |
|
|
|
|
|
|
|
0.000 |
1.215 |
|
0.000 |
|
|
|
|
|
|
|
0.000 |
-1.215 |
|
0.000 |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
0.000 |
0.000 |
|
1.215 |
|
|
|
|
|
|
|
0.000 |
0.000 |
|
-1.215 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Матрица планирования Z
Общий вид матрицы планирования эксперимента ОЦКП:
4
|
|
Z11 |
Z12 |
Z13 |
Z = {Zlj |
}= |
Z21 |
Z22 |
Z23 |
|
|
|
|
|
|
|
Z15,1 |
Z15,2 |
Z15,3 |
где l-номер опыта, от 1 до15, j-номер фактора, от 1 до 3. Найдём дисперсию в каждой серии опытов по
формуле:
DYl =Dl = |
1 ∑(Yl −Yl )2 |
= |
|
1[(Yl1 −Yl )2 +(Yl2 −Yl )2 +(Yl3 −Yl )2], |
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
j=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
m−1 |
J |
|
2 |
|
|
|
|
|
|
|
||||||||
l = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D1 =51 |
[(Y1 − |
|
)2 +(Y1 − |
|
)2 +(Y1 − |
|
)2]= |
|
|
|
|
|
|
|
|||||||||
Y1 |
Y1 |
Y1 |
|
|
|
|
|
|
|
||||||||||||||
2 |
1 |
|
|
|
|
|
2 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|||
=1(−0.1)2 |
=1 |
0.0 =0.00 |
|
|
|
|
|
|
|
|
|
|
|||||||||||
2 |
|
|
|
2 |
1 |
5 |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
Аналогично рассчитывают остальные дисперсии. |
|
|||||||||||||||||||||
Найдём коэффициенты нормированной модели |
|
|
|
||||||||||||||||||||
Y=β0 +β1X1+β2X2+β3X3 +β1 X1X2+β1 X1X3 +β2 X2X3 +β1 X12+β2 X22+β3 X32, |
|||||||||||||||||||||||
|
|
|
|
x1−x1δ |
|
|
|
|
|
|
2 |
3 |
3 |
1 |
2 |
3 |
|||||||
где = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
∆x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
X |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X2 =x2∆−x2x2δ
X3 =x3∆−x3x3δ
где X1, X2, X3-нормированные значения факторов.
ЛИНЕЙНЫЕ КОЭФФИЦИЕНТЫ β1, β2, β3 находят по формуле:
5
15
βj = C1 ∑Zel ×Yl , j =1,2,3.
l=1
Zlj - элементы матрицы планирования экспериментов, при этом 1-й столбец матрицы планирования скалярно умножается на столбец средних значений.
Коэффициенты β2,β3 рассчитываются аналогично, но вместо первого столбца берутся соответственно второй и третий.
β1 =C1[Z11 ×Y1 +Z21 ×Y2 +... +Z15,1 ×Y15]=C1[(−1)×Y1 +1×Y2 +(−1)Y3 +(−1)×Y4 +... +0×Y15]
СМЕШАННЫЕ КОЭФФИЦИЕНТЫ β12, β13, β23 находят по формуле:
|
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
= C2[Z1i ×Z1 j × |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
] |
|
||||||||
βij |
= C2∑Zli ×Zlj |
× |
|
Yl |
|
Y1 |
|
+Z2i ×Z2 j × |
Y2 |
+... +Z15,i ×Z15, j |
Y15 |
|
|||||||||||||||||||||||||||||||
|
|
l =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
при этом перемножаются два столбца (i-ый и j-ый) из матрицы |
|||||||||||||||||||||||||||||||||||||||||||
планирования и столбец средних значений. |
|||||||||||||||||||||||||||||||||||||||||||
Например, требуется найти |
|||||||||||||||||||||||||||||||||||||||||||
β12 |
=? |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i =1, j =2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
β12 |
=C2 |
[Z11 ×Z12 |
× |
|
|
|
|
+ Z21 ×Z22 × |
|
|
+ Z31 ×Z32 × |
|
+... + Z15,1 ×Z15,2 × |
|
]= |
||||||||||||||||||||||||||||
Y1 |
Y2 |
Y3 |
Y15 |
||||||||||||||||||||||||||||||||||||||||
=C1 [(−1) ×(−1) |
|
|
|
+(+1)(−1) |
|
+(−1)(1) |
|
+... +0 ×0 × |
|
]=... |
|
|
|||||||||||||||||||||||||||||||
Y1 |
Y2 |
Y3 |
Y15 |
||||||||||||||||||||||||||||||||||||||||
Аналогично можно найти β13 и β23. |
|||||||||||||||||||||||||||||||||||||||||||
КВАДРАТИЧНЫЕ КОЭФФИЦИЕНТЫ β11, β22, β33 находят |
|||||||||||||||||||||||||||||||||||||||||||
по формуле: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
βjj =C3∑15 |
(Zlj2 − ji) |
|
,j = |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
Yl |
1,3 |
||||||||||||||||||||||||||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
γ =(2n−p +2α2) N =[23 +2(α)2] 15 |
|||||||||||||||||||||||||||||||||||||||||||
β11 |
=?J =1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
β11 |
=C3[(Z112 −γ) |
|
+(Z212 −γ) |
|
+...+(Z15,12 −γ) |
|
]=C3[((−1)2 −γ) |
|
+((12) −γ) |
|
+...+(0 −γ) |
|
]=... |
||||||||||||||||||||||||||||||
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
Y1 |
Y2 |
Y15 |
Y1 |
Y2 |
Y15 |
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Аналогично находим β22 (j=2), β33(j=3).
СВОБОДНЫЙ ЧЛЕН рассчитываем по формуле:
|
15 |
3 |
|
15 |
|
|
|
||
β0 = |
1 |
∑ |
|
−γ∑βjj = |
1 |
∑ |
|
−γ(β11 + β22 + β33 ) |
|
Yl |
Yl |
|
|||||||
|
|
|
|||||||
|
15 l =1 |
|
j =1 |
15 l =1 |
|
||||
I. |
|
|
|
Проверим |
|
|
значимость |
найденных |
|
коэффициентов по t-критерию Стьюдента. |
|
Дисперсии коэффициентов находят по следующим формулам. |
|||||||||||||||||||
Для линейных коэффициентов: |
|
|
|
|
|
||||||||||||||
Dв=(D1+D2+…+D15)/15 – дисперсия воспроизводимости |
|||||||||||||||||||
опытов, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D(βj) = |
|
C1 |
|
DB |
= C1 DB , |
|
|
||||||||||||
|
|
|
|
|
|
|
m |
|
|
|
|
3 |
|
|
|
|
|
|
|
D(β1) = D(β2) = D(β3 ) = |
C1 |
DB |
|
||||||||||||||||
Для проверки значимости |
3 |
|
|
|
|
||||||||||||||
смешанных коэффициентов |
|||||||||||||||||||
используются расчётные формулы: |
|
|
|||||||||||||||||
D(βj) = |
C2 |
|
DB = |
C2 |
|
|
DB |
|
|
|
|
|
|||||||
|
m |
3 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
C2 |
|
|
|
|
||||||
D(β12 ) = D(β13 ) = D(β23 ) = |
|
|
DB |
|
|||||||||||||||
|
|
||||||||||||||||||
Для |
|
|
|
проверки |
|
|
3 |
|
|
|
квадратичных |
||||||||
|
|
|
значимости |
||||||||||||||||
коэффициентов используются расчётные формулы: |
|||||||||||||||||||
D(βj) |
= |
C3 |
|
DB = |
|
C3 |
|
DB |
|
|
|
|
|
||||||
m |
|
3 |
|
|
C3 |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
D(β11 ) = D(β22 ) = D(β33 ) = |
|
DB |
|
||||||||||||||||
3 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
C0 |
|
|
|
3 |
|
|
|
|
|
|
|
|
|
||||
D(β0 ) = |
|
DB −γ 2 ∑D(β jj ) |
|
|
|
|
|
||||||||||||
|
|
m |
|
|
|
|
j =1 |
|
|
|
|
|
|
|
|||||
D(β0 ) = |
|
C0 |
DB −γ 2 (D(β11 ) + D(β22 ) + D(β33 )) |
||||||||||||||||
3 |
|||||||||||||||||||
Для |
|
|
|
|
|
|
значимости |
свободного члена |
|||||||||||
|
проверки |
|
|
|
используются расчётные формулы:
7
Значения табличных коэффициентов приведены в [9], с. 87. Для проверки значимости коэффициентов находим расчётные значения t-критерия Стьюдента по формуле:
t расч (β) = |
|
|
β |
|
|
|
|
|
|
|
|||
|
|
|
|
|
||
|
|
|
|
|
|
|
D(β) |
|
|||||
|
|
|
Для всех коэффициентов достаточно найти одно значение t крит по таблице критических значений для t-критерия Стьюдента при f=N(m-1)=15(3-1)=30
Если t расч(β)<t крит, то коэффициент β незначимый, его исключают из модели, приравнивая к нулю. Например, если t
расч(β12)<t крит., то β12=0
Если t расч (β)>t крит, то коэффициент значимый и его оставляют в модели.
После того, как в модели остались только значимые коэффициенты, нужно проверить адекватность полученной математической модели по критерию Фишера, то есть сравнить значения Y, полученные при расчёте по нормированной модели (1) с средними значениями по каждой сери и опытов. При расчете по нормированной модели в качестве значений X1, X2, X3 выбирают L-ую строку матрицы планирования и находят при
L =1,15
Y = β0 + β1 ×Z 1 + β2 ×Z 2 + β3 ×Z 3 + β12 ×Z 1 ×Zl 2 + β13 ×Z 1 ×Z 3 + + β23 ×Zl 2 ×Zl3 + β11 ×Z 12 + β22 ×Z 22 + β33 ×Z 32
Например, если все коэффициенты значимые
Y1 =β0 +β1 (−1) +β2(−1) +β3 (−1) +β12(−1)(−1) +β13 (−1)(−1) +β23(−1)(−1) +β11 (−1)2 +
+β22(−1)2 +β33 (−1)2 __ и_ т.д.__Y2,Y3 ,...,Y15
Сравнить между собой значения результатов эксперимента с значениями, рассчитанными по найденной математической модели позволяет дисперсия адекватности
8
DA = N 1−d ∑N (Y −Y )2 ,_ где_ N =15,
=1
а d - количество незначимых коэффициентов, которые мы исключили из модели(приравняли к нулю).
Например, если d=0
DA = 151 [(Y1 −Y1 )2 +(Y2 −Y2 )2 +... +(Y15 −Y15 )2 ]
Расчётное значение критерия Фишера
Fрасч = DA
DB
где Dв -дисперсия воспроизводимости, которую мы использовали при проверке значимости коэффициентов. Табличное (критическое) значение критерия Фишера
F крит находят по таблице.
Степени свободы f1 и f2 выбирают по правилу.
1)Если DА<Dв, то f1=N-d=15-d f2=N(m-1)=15(3-1)=30
2)Если Dв<DA, то f1=N(m-1)=30 f2=N-d=15-d
Если F расч. < F крит., то получена адекватная нормированная модель (1)
Если F расч. > F крит., то модель неадекватна, её использовать нельзя. Для получения адекватной модели рекомендуется
уменьшить шаги варьирования (∆X1= d1, ∆X2= d2, ∆X3= d3). II. Если нормированная модель (1) адекватна, то нужно перейти к реальным физическим величинам.
Для этого в модель (1) с учётом того, что незначимые коэффициенты =0 нужно подставить
X1 = x1∆−x1x1δ ,X2 = x2∆−x2x2δ ,X3 = x3∆−x3x3δ
9
МЕТОДИЧЕСКИЕ УКАЗАНИЯ ПО КУРСОВОМУ ПРОЕКТИРОВАНИЮ
1. ЦЕЛЬ КУРСОВОГО ПРОЕКТА Целью курсового проекта является приобретение студентами
знаний, умений и практических навыков по синтезу адекватных математических моделей на основе уравнения регрессии в виде полноквадратичной зависимости.
2. ЗАДАНИЕ НА КУРСОВУЮ РАБОТУ Для данных эксперимента, соответствующих
индивидуальному варианту студента (определяемого по двум последним цифрам зачётки) провести проверку выполнения условия воспроизводимости опытов, получить экономичную математическую модель и оценить её работоспособность.
3. ОБЪЕМ ЗАДАНИЯ Для выполнения курсового проекта необходимо
осуществить следующие операции:
-провести обоснование выбора метода синтеза экспериментально-статистической модели;
-проверить условия применимости регрессионного анализа;
-расчёт коэффициентов регрессии и оценка значимости;
-проверить адекватность модели и получить модель в
реальных физических величинах.
Оформление пояснительной записки к курсовой работе проводится в соответствии с стандартом ВГТУ. Пояснительная записка должна включать:
-титульный лист; -задание на курсовой проект;
-лист «Замечания руководителя»; -содержание; -введение; -основные разделы; -заключение; -список литературы; -приложения.
10