Учебное пособие 1491
.pdf157
4.4.4. Лабораторное задание и методические указания по его выполнению
Предположим, что распределение функции поля U=U( x , y ) в конструкции РЭС описывается дифференциальным уравнением
где = [a,b]x[c,d] прямоугольная область.
Рис. 4.15
Линейка меню состоит из следующих компонент:
158
Рис. 4.16
А. Меню ―Файл‖, в котором можно загрузить и списать параметры дифференциального уравнения и разностной сетки:
Б. Меню ―Опции‖, в котором можно сохранить результаты расчета значений дифференциального уравнения на разностной сетке, создать журнал записей и изменить направление интерполяции.
2U( x ,y ) |
|
2U( x ,y ) |
U( x ,y ) |
U( x ,y ) |
U( x ,y ) |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.14) |
а |
a |
22 |
|
a |
a |
a |
a U( x ,y ) =b |
1 |
||||||
11 |
+ |
|
+ 12 |
|
+ 1 |
|
+ 2 |
|
+0 |
|
||||
|
x 2 |
|
|
y 2 |
|
x y |
x |
|
y |
|
|
Рис. 4.17 |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с граничными условиями |
|
|
|
|
|
|
|
|||||||
|
|
|
U( x ,y ) |
= b2, , |
|
|
|
|
|
|
(4.1 |
159
После выбора в меню пункта ―Направление интерполяции‖ на экране появится окно следующего вида:
Рис. 4.18
160
Влевой части окна выбирается направление интерполяции по оси Х и по оси Y соответственно, а в правой части окна выводится визуальное отображение направления интерполяции.
Впанели дифференциального уравнения непосредственно
задаются коэффициенты при |
соответствующем порядке |
производной функции. |
|
Рис. 4.19
Вводимые коэффициенты представляют собой конечную суперпозицию элементарных функций двух переменных х и у. В программе реализованы следующие функции и операции:
операции сложения, вычитания, умножения, деления;
тригонометрические функции - sin, cos, tg, ctg, arcsin, arccos, arctg, arcctg;
натуральный логарифм - ln, экспонента – exp; модуль - abs, знак числа – sgm;
степенная и показательные функции реализуются знаком ^, например, x - будет записываться в виде x^ .
Замечание: эта операция реализована на x, R , поэтому для
реализации степенной функции с целой степенью рекомендуется использовать знак умножения. Панель задания разностной граничных условий имеет вид:
161
Окно разностной сетки обеспечивает пользователю следующие возможности. Значения функции после вычисления можно узнать посредством нажатия мышью на соответствующую точку в разностной сетке, а само значение будет отображено в нижней части панели задания граничных условий. Там же будет отображена точность вычислений.
Рис. 4.20
Задание значений функции на границе разностной сетки производится аналогично заданию коэффициентов при производных функции. В редакторах min x, max x, min y, max у производится задание границ разностной сетки. Сохранение значений границ и значения шага сетки производится после нажатия кнопки ENTER в соответствующем редакторе. После задания дифференциального уравнения значений на границе разностной сетки и шага разностной сетки можно произвести вычисление значений функции во внутренних точках разностной сетки посредством нажатия кнопки CALCULATE.
162
Рис. 4.21
Результатом работы с программным комплексом является протокол пользователя - файл имя.log , где ―имя‖ задается студентом при сохранении данных в режиме включенного журнала записей. По результатам расчета следует оценить эффективность метода конечных разностей и качество полученных решений.
4.4.5. Контрольные вопросы
1. Какова цель лабораторной работы?
163
2. В чем заключается лабораторное задание? Пояснить ход его выполнения.
3. Какие данные являлись исходными для Вашего варианта?
4.Какую разностную сетку Вы построили? Поясните выбор шага сетки.
5.Какие разностную схемы Вы выбрали и почему?
6.Какие разностные схемы не обеспечивают требуемой точности для Вашего варианта?
7.Какие шаблоны имеют выбранные Вами разностные схемы?
164
11. Сформулируйте выводы по данной лабораторной работе.
4.4.6. Указания по оформлению отчета
Отчет по каждой лабораторной работе должен содержать наименование и цель работы, краткие теоретические сведения в виде ответов на вопросы домашнего задания, ход и результаты выполнения лабораторного задания, где приводятся исходные данные, промежуточные расчеты и результаты работы с программным комплексом с необходимыми пояснениями. Отчет
8.Рассчитайте суммарную погрешность каждой схемы. завершается кратким перечнем приобретенных навыков и выводами
орезультатах проделанной работы. Оформление отчета выполняется
9. Проведите анализ машинного решения. |
в соответствии с ГОСТ. |
10. Перечислите приобретенные при выполнении работы знания и навыки.
Исходные данные вариантов к лабораторной работе N 8
|
Вариант 9 |
|
|
Вариант 10 |
|
Вариант 11 |
|
Вариант 12 |
|
||||
|
X1 |
|
X2 |
Y |
X1 |
X2 |
Y |
X1 |
X2 |
Y |
X1 |
X2 |
Y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.1 |
|
0.2 |
9.3 |
0.4 |
0.2 |
10.9 |
0.5 |
0.3 |
11.4 |
0.2 |
0.4 |
12.6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.3 |
|
0.6 |
9.9 |
1.2 |
0.5 |
10.7 |
0.2 |
0.9 |
11.2 |
0.1 |
0.2 |
12.3 |
|
0.4 |
|
0.8 |
9.2 |
1.6 |
0.4 |
10.5 |
0.9 |
1.2 |
11.6 |
0.3 |
0.6 |
12.9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.5 |
|
1.0 |
9.5 |
2.0 |
0.7 |
10.4 |
0.3 |
1.5 |
11.0 |
0.4 |
0.8 |
12.2 |
|
0.6 |
|
1.2 |
9.8 |
2.4 |
0.2 |
10.8 |
0.5 |
1.8 |
11.4 |
0.5 |
1.0 |
12.5 |
|
0.7 |
|
1.4 |
9.1 |
2.8 |
0.5 |
10.7 |
0.2 |
2.1 |
11.8 |
0.6 |
1.2 |
12.8 |
|
0.8 |
|
1.6 |
9.4 |
3.2 |
0.4 |
10.5 |
0.9 |
2.4 |
11.2 |
0.7 |
1.4 |
12.1 |
|
0.9 |
|
1.8 |
9.7 |
3.6 |
0.8 |
10.4 |
0.3 |
2.7 |
11.6 |
0.8 |
1.6 |
12.4 |
|
0.25 |
|
0.5 |
9.75 |
2.0 |
0.2 |
10.9 |
0.5 |
1.75 |
11.0 |
0.9 |
1.8 |
12.7 |
|
0.15 |
|
0.3 |
9.45 |
0.6 |
0.5 |
10.7 |
0.2 |
0.45 |
11.6 |
0.25 |
0.5 |
12.5 |
|
0.35 |
|
0.7 |
9.05 |
1.4 |
0.4 |
10.5 |
0.9 |
1.05 |
11.4 |
0.15 |
0.3 |
12.4 |
|
0.45 |
|
0.9 |
9.35 |
1.8 |
0.8 |
10.4 |
0.3 |
1.35 |
11.8 |
0.35 |
0.7 |
12.5 |
|
0.55 |
|
1.1 |
9.65 |
2.2 |
0.2 |
10.9 |
0.5 |
1.65 |
11.2 |
0.45 |
0.9 |
12.5 |
|
0.65 |
|
1.3 |
9.95 |
2.6 |
0.5 |
10.7 |
0.2 |
1.95 |
11.6 |
0.55 |
1.1 |
12.5 |
|
0.75 |
|
1.5 |
9.25 |
3.0 |
0.4 |
10.5 |
0.9 |
2.25 |
11.0 |
0.65 |
1.3 |
12.9 |
|
0.85 |
|
1.7 |
9.55 |
3.4 |
0.8 |
10.4 |
0.3 |
2.55 |
11.4 |
0.75 |
1.5 |
12.2 |
|
0.95 |
|
1.9 |
9.85 |
3.8 |
0.2 |
10.9 |
0.5 |
2.85 |
11.8 |
0.85 |
1.7 |
12.5 |
|
1.25 |
|
1.5 |
9.75 |
2.0 |
0.5 |
10.7 |
0.2 |
1.75 |
11.0 |
0.95 |
1.9 |
12.8 |
|
1.2 |
|
1.3 |
9.4 |
1.5 |
0.4 |
10.5 |
0.9 |
1.4 |
11.5 |
1.25 |
1.5 |
12.7 |
|
1.3 |
|
1.6 |
9.9 |
2.2 |
0.8 |
10.4 |
0.3 |
1.9 |
11.2 |
1.2 |
1.3 |
12.4 |
|
1.4 |
|
1.8 |
9.2 |
2.6 |
0.2 |
10.9 |
0.5 |
2.2 |
11.6 |
1.3 |
1.6 |
12.9 |
|
1.5 |
|
2.0 |
9.5 |
3.0 |
0.5 |
10.7 |
0.2 |
2.5 |
11.0 |
1.4 |
1.8 |
12.2 |
|
1.6 |
|
2.2 |
9.8 |
3.4 |
0.4 |
10.5 |
0.9 |
2.8 |
11.4 |
1.5 |
2.0 |
12.5 |
|
0.25 |
|
0.5 |
9.75 |
2.0 |
0.2 |
10.9 |
0.5 |
1.75 |
2.0 |
1.6 |
2.2 |
2.8 |
5.Курсовая работа на тему Получение математической модели методом ОЦКП
Задание на курсовую работу Получить модель технологического процесса
по экспериментальным данным варианта, проверить модель на адекватность и оценить значимость факторов, влияющих на качество технологического процесса.
165
Исходные данные для ОЦКП
|
|
|
|
|
|
|
|
|
|
|
|
X1 |
X2 |
X3 |
|
|
Yl1 |
Yl2 |
Yl3 |
|
|
|
|
|
|
|
|
|
|
l=1 |
|
-1.000 |
-1.000 |
-1.000 |
|
|
34.1 |
34.0 |
34.1 |
l=2 |
|
1.000 |
-1.000 |
-1.000 |
|
|
49.4 |
34.0 |
34.0 |
l=3 |
|
-1.000 |
1.000 |
-1.000 |
|
|
34.0 |
49.4 |
34.0 |
l=4 |
|
-1.000 |
-1.000 |
1.000 |
|
|
34.1 |
34.0 |
49.4 |
l=5 |
|
1.000 |
1.000 |
-1.000 |
|
|
49.4 |
49.4 |
34.0 |
l=6 |
|
-1.000 |
1.000 |
1.000 |
|
|
34.1 |
49.4 |
49.4 |
l=7 |
|
1.000 |
-1.000 |
1.000 |
|
|
49.4 |
34.0 |
49.4 |
l=8 |
|
1.000 |
1.000 |
1.000 |
|
|
49.4 |
49.4 |
49.4 |
l=9 |
|
0.000 |
0.000 |
0.000 |
|
|
41.4 |
41.3 |
41.4 |
l=10 |
|
1.215 |
0.000 |
0.000 |
|
|
51.2 |
41.3 |
41.3 |
l=11 |
|
-1.215 |
0.000 |
0.000 |
|
|
32.6 |
41.4 |
41.3 |
l=12 |
|
0.000 |
1.215 |
0.000 |
|
|
41.3 |
51.2 |
41.3 |
l=13 |
|
0.000 |
-1.215 |
0.000 |
|
|
41.4 |
32.6 |
41.3 |
l=14 |
|
0.000 |
0.000 |
1.215 |
|
|
41.4 |
41.4 |
51.2 |
l=15 |
|
0.000 |
0.000 |
-1.215 |
|
|
41.3 |
41.3 |
32.6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l-номер опыта, от 1 до15
X1 X 2 X 3
d1 0.1 d2 0.2 d3 0.3
X1 X 2 X3
Задана матрица планирования эксперимента и результаты эскпериментов. Для каждого эксперимента проведено 3 дублирующих опыта. Обозначим
n=3 – количество факторов,
m=3 – количество дублирующих опытов, N 2n 2n 1
общее количество экспериментов.
L – порядковый номер эксперимента, L = 1,…,N .
N=15 (15 экспериментов). Результаты всех опытов запишем в виде матрицы размерности 15х3, обозначим еѐ элементы
Y lj, l-номер эксперимента, j-номер дублирующего опыта. Результаты опытов
166 |
167 |
|
Y11Y12Y13
Y21Y22Y23
Y15,1Y15,2Y15,3
Находим среднее значение в каждой серии опытов
Найдѐм дисперсию в каждой серии опытов.
|
|
|
|
|
|
|
|
1 |
|
|
m |
|
|
|
|
|
|
|
|
1 |
(Yl1 |
|
|
|
|
|
|
, |
|||||
DYl |
|
Dl |
|
|
(YlJ |
Yl )2 |
тюкюm 3 |
Yl |
)2 (Yl 2 |
Yl |
)2 (Yl 3 |
Yl |
)2 |
||||||||||||||||||||
m |
|
2 |
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
1 j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
l |
1.15 длякаждогоэксперимента |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
1 |
(Y11 |
|
|
|
|
|
|
|
(Y13 |
|
)2 |
1 |
(34.1 34.1)2 |
|
(34.1 34.1)2 |
|||||||||||||||
D1 |
Y1 )2 (Y12 |
Y1 )2 |
Y1 |
(34 34.1)2 |
|||||||||||||||||||||||||||||
|
2 |
||||||||||||||||||||||||||||||||
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
1 |
|
( 0.1)2 |
|
1 |
0.01 0.005 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
АналогичнонайтиD2 , D3 ,..., D15
I. |
|
Проверим выполнение |
условия воспроизводимости |
||||||
G |
|
|
max D1, D2 ,..., D15 |
|
|
max D1, D2 ,...D15 |
|
... |
|
расч |
|
|
|||||||
|
15 |
|
|
D1 D2 ... D15 |
|
||||
|
|
|
Dj |
|
|
||||
|
|
|
|
|
|
|
|
||
|
|
|
l 1 |
|
|
|
|
|
|
|
|
опытов (однородность |
дисперсий) |
по G-критерию |
Кокрена.
Gкрит по табл. 1.6 на с.140 книги Львович, Фролов при
168
f1 m 1 3 1 2 f2 N 15
Если Gрасч.<Gкрит., то условие воспроизводимости опытов выполнено и можно применять метод ОЦКП, иначе расчѐт окончен.
II.Найдѐм коэффициенты нормированной модели
Y 0 1 X1 2 X 2 3 X 3 12 X1 X 2 13 X1 X 3 23 X 2 X 3 11 X12 22 X 22 33 X 32 , (1)
гдеX1 |
|
x1 x1 |
||||
x1 |
||||||
|
|
|
|
|||
X 2 |
|
x2 x2 |
|
|||
x2 |
||||||
|
|
|
||||
X 3 |
|
|
x3 x3 |
|
||
|
x3 |
|||||
|
|
|
X1, X2, X3-нормированные значения факторов, x1, x2 , x3 реальныезначенияфакторов
по формулам:
(C1,C2,C3 из таблицы 7.13 с. 87 Львович,Фролов)
ЛИНЕЙНЫЕ КОЭФФИЦИЕНТЫ
15
j C1 Zel Yl , j 1,2,3.
l 1
1 C1 Z11 Y1 Z21 Y2 ... Z15,1 Y15 C1 ( 1) Y1 1 Y2 ( 1)Y3 ( 1) Y4 ... 0 Y15
Zlj - элементы матрицы планирования экспериментов
Т.е. 1-й столбец матрицы планирования скалярно умножается на столбец средних значений.
Коэффициенты 2, 3 рассчитываются аналогично, но вместо первого столбца берутся соответственно второй и третий.
СМЕШАННЫЕ КОЭФФИЦИЕНТЫ12, 13, 23 находят по формуле:
Т.е. перемножаются 2 столбца i-ый и j-ый из матрицы планирования и столбец средних значений.
169
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.
КВАДРАТИЧНЫЕ КОЭФФИЦИЕНТЫ
|
15 |
|
C2 Z1i Z1 j |
|
|
|
|
|
|
ij |
C2 Zli Zlj |
Yl |
Y1 |
Z2i Z2 j |
Y2 |
... Z15,i Z15, j |
Y15 |
||
|
l 1 |
|
|
|
|
|
|
|
находят по формуле: Аналогично 22, 33
(j=2) (j=3)
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
jj |
C3 (Zlj |
2 ji)Yl , j 1,3 |
|
|
|
|
|
|
|
|
|
|
||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(2n p 2 |
2 ) |
N |
|
23 2( )2 |
|
15 |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
11 |
? J 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
C3 (Z112 ) |
|
|
|
|
... (Z15,12 ) |
|
C3 (( 1)2 ) |
|
((12 ) ) |
|
... (0 ) |
|
... |
||||||
11 |
Y1 |
(Z212 )Y2 |
Y15 |
Y1 |
Y2 |
Y15 |
III.Проверим значимость найденных коэффициентов по t-
критерию Стьюдента.
Дисперсии коэффициентов находят по формулам: для линейных коэффициентов:
Dв =( D1 + D2 + … + D15)/15 – дисперсия воспроизводимости, была рассчитана в пункте I.
|
|
|
|
|
|
170 |
|
|
|
D( j) |
C1 D |
|
C1 D , гдеD |
|
_ _ |
_ |
|||
|
m 3 |
B |
|
3 |
B |
B |
|
|
|
|
|
|
|
|
l |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
l 1 |
|
|
D( 1 ) D( 2 ) D( 3 )
|
t расч |
( ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
D( ) |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
t |
расч |
( |
0 ) |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
D( 0 ) |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
t |
расч |
( |
2 ) |
|
|
|
|
|
2 |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
D( 2 ) |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
t расч |
(1 ) |
|
|
|
|
1 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
D( 1 ) |
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
C31 DB
D(j) Cm3 DB C33 DB
D(11) D(22 ) D(33 ) C33 DB
D(j) Cm2 DB C32 DB
D(12 ) D(13 ) D(23 ) C32 DB
171
для смешанных коэффициентов
t расч |
(12 ) |
|
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
D(12 ) |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
t |
расч |
( |
13 |
) |
|
|
|
|
13 |
|
|
|
|
|
|
|
|
t расч |
(33 ) |
|
|
33 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
D(13 ) |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D(33 ) |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
t расч ( |
) |
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
D(3 ) |
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для всех коэффициентов находят одно значение t крит. По таблице П1.4 с. 139 (Львович, Фролов) при
f=N(m-1)=15(3-1)=30
Если t расч( )<t крит, то коэффициент незначимый, его исключают
t расч |
( |
22 ) |
|
|
|
22 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|||
D( 22 ) |
|
|||||||
|
|
|
|
|
из модели, приравнивая к нулю.
Например, если t расч(12)<t крит., то 12=0
172
|
|
|
|
|
C0 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
||
D(0 ) |
|
|
DB 2 D( jj ) |
|
|
|
|
||||||||||||||||
|
|
m |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
D( |
|
) |
|
C0 |
D 2 (D( |
|
|
) D( |
|
) D( |
|
)) |
|||||||||||
0 |
|
|
11 |
22 |
33 |
||||||||||||||||||
|
|
|
3 |
|
|
B |
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
11 |
|
|
|
|
|
|
|
|
|
|
|||
|
t расч |
(11 ) |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
D(11 ) |
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Для квадратичных коэффициентов |
|
|
|||||||||||||||||||||
|
|
|
|
15 |
|
|
|
|
|
3 |
|
|
|
15 |
|
|
|
|
|
||||
0 |
|
1 |
|
|
|
|
jj |
1 |
|
|
(11 22 33 ) |
||||||||||||
|
|
Yl |
Yl |
||||||||||||||||||||
|
|
|
|
||||||||||||||||||||
|
|
15 l 1 |
|
|
|
|
j 1 |
|
|
15 l 1 |
|
|
|
|
Для свободного члена
Если t расч ()>t крит, то коэффициент значимый и его оставляют в модели.
IV. После того, как в модели остались только значимые коэффициенты, нужно проверить адекватность полученной математической модели по критерию Фишера, то есть
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 Zl 3 11 Z 12 22 Z 2 2 33 Z 3 2
173
V.сравнить значения Y, полученные при расчѐте по нормированной модели (1) с средними результатами эксперимента (значениями Y).
При расчете по нормированной модели в качестве значений X1, X2, X3 выбирают L-ую строку матрицы
планирования и находят при L==1,…,15Сравнить между собой значения
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y , Y ,..., Y |
_ и _ Y , Y ,...Y |
|||||||||||||||||||||
|
1 |
|
2 |
|
|
|
15 |
|
1 |
2 |
|
15 |
|
|||||||||
|
|
|
позволяет дисперсия адекватности |
|
|
|
||||||||||||||||
|
|
|
|
1 |
|
|
N |
2 |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
( |
Y |
Y ) |
, _ где _ N 15, |
|
|
|
|
|
|
|||||||
DA |
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
N |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
а |
|
|
|
d-количество незначимых коэффициентов, которые мы |
|||||||||||||||
|
|
|
исключили из модели(приравняли к нулю). |
|||||||||||||||||||
|
|
|
Например, если d=0 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
1 |
( |
|
Y )2 |
|
|
|
|
|
|
|
|
||||||||
D |
A |
|
Y |
(Y |
Y )2 |
... (Y |
Y )2 |
|||||||||||||||
|
||||||||||||||||||||||
|
15 |
1 |
|
1 |
|
2 |
|
2 |
15 |
15 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчѐтное значение критерия Фишера Dв-дисперсия воспроизводимости, которую мыиспользовали при проверке значимости коэффициентов.
174
Табличное (критическое) значение критерия Фишера F крит находят по таблице П1.3 на с.138 при числах Степеней свободы 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).
Повторный расчѐт не проводить
VI. Если нормированная модель (1) адекватна, то нужно перейти к реальным физическим величинам.
Для этого в модель (1) с учѐтом того, что незначимые
|
|
|
x x |
|
|
x x |
|
|
|
x x |
|
X |
|
1 1 |
, X |
|
2 2 |
, X |
|
3 3 |
|||
1 |
|
2 |
|
3 |
|
||||||
|
|
x1 |
|
x2 |
|
|
x3 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
176 |
|
|
|
175 |
|
|
|
|
|
ВПИ, 1990. 92 с. |
|
|
|
|
|
|
|
|
|
12. Русак И.М., Луговский В.П. Технические средства ПЭВМ: |
|
|
Список литературы |
|
|
|
Справочник / Под ред. И.М. Русака.- Мн.:Высш шк., 1996.- 304 с.: ил. |
||||
1. |
Деньдобренько В.Н., Малика А.С. Автоматизация |
|
|
|
|||||
конструирования РЭА: Учебник для вузов. М.: Высшая школа, 1980, |
|
|
|
||||||
384 с. |
|
|
|
|
|
|
|
ОГЛАВЛЕНИЕ |
|
2. |
Морозов К.К., Одиноков В.Г., Курейчик В.М. |
|
|
|
|||||
Автоматизированное проектирование конструкций радиоэлектронной |
ВВЕДЕНИЕ |
|
|||||||
аппаратуры: Учеб. пособие для вузов. М.: Радио и связь, 1983. 280 с. |
|
|
|
|
|||||
3. |
Норенков |
И.П., |
Маничев |
В.В. |
Системы |
1. |
ПРОГРАММА ДИСЦИПЛИНЫ |
2 |
|
автоматизированного |
проектирования |
электронной |
и |
|
|
|
|||
вычислительной аппаратуры: Учеб. пособие для вузов. М.: Высшая |
2. КОНТРОЛЬНЫЕ РАБОТЫ |
39 |
|||||||
школа, 1983, 289 с. |
|
|
|
|
|
|
|
|
|
4. |
Керимов В.Г., Багиров С.А. Автоматизированное |
3. |
КУРС ЛЕКЦИЙ |
48 |
|||||
проектирование конструкций. М.: Машиностроение, 1985, 224 с. |
|
|
|
|
|||||
5. |
Корячко В.П., Курейчик В.М., Норенков И.П. |
4. |
ЛАБОРАТОРНЫЙ ПРАКТИКУМ |
117 |
|||||
Теоретические основы САПР: Учебник для вузов. М.: |
|
|
|
||||||
Энергоатомиздат, 1987.400 с. |
|
|
|
|
5. |
КУРСОВАЯ РАБОТА |
164 |
||
6. |
Автоматизация схемотехнического проектирования / Под |
|
|
|
|||||
ред. В.Н. Ильина. М.: Радио и связь, 1987. 386 с. |
|
|
|
СПИСОК ЛИТЕРАТУРЫ |
175 |
||||
|
|
|
|
|
|
|
7. Кофанов Ю.Н. Теоретические основы конструирования, технологии и надежности радиоэлектронных средств: Учебник для вузов. М.: Радио и связь, 1991. 360 с.
8.Львович Я.Е., Фролов В.Н. Теоретические основы конструирования, технологии и надежности РЭА: Учеб. пособие для вузов. М.: Радио и связь, 1986. 192 с.
9.Яншин А.А. Теоретические основы конструирования, технологии и надежности ЭВА: Учеб. пособие для вузов. М.: Радио и связь, 1983, 312 с.
10.Норенков И.П. Введение в автоматизированное проектирование технических устройств и систем: Учеб. пособие для вузов. М.: Высшая школа, 1980. 311 с.
11.Техническое обеспечение САПР : Учеб.пособие /В.А.Дыбой, В.Е.Межов, В.М. Питолин, И.Г. Проценко; Воронеж :