РУХУ ПО градієнт )% -------------
format bank = b (2:4) ';% = 0.004; = 6; j = 1: ng (j, 1:3) = X0 + (j-1) * gamma * b. * dX;
Елементи матриці G говорять про те, що ціна 1-го товару при русі по градієнту вийде за кордон його закупівельної ціни на 6-му кроці. Тому при проведенні дослідів, які забезпечують рух по градієнту, обмежимося п'ятьма кроками:
disp ('РУХ ПО градієнт )% ---------- = profit3norm (G (2: ng-1, :))
[[Ym, kYm] = max (Y);
Враховуючи, що на основному рівні обраного дизайну отримано значення прибутку 4445.15 руб., можна припустити, що реалізовано сходження по градієнту. Далі треба переконатися, що на третьому кроці сходження було отримано значення прибутку, яке статистично значимо відрізняється від значення на основному рівні. З цією метою почнемо продавати товари за новими цінами:
disp ('ПОРІВНЯННЯ СЕРЕДНІХ )%
N = 4; = G (ones (1, N), :) = profit3norm (D0) ('Середня прибуток на основному рівні матриці дизайну:') mean = sum (Y0)/N = G ([kYm kYm kYm kYm], :) = profit3norm (D) ('Середня прибуток на 3-му кроці руху по градієнту:') = sum (Ymgrad)/N = [Y0 Ymgrad]; ('Оцінка p-value відмінності середніх ') short = anoval (gradmean)
Результати порівняння середніх, представлені на рис. 1 і 2, говорять про статистично значущому їх відмінності і дають підставу перейти до побудови моделі 2-го порядку з метою більш точного визначення координат екстремальних продажів. p align="justify"> disp ('ЦЕНТРАЛЬНИЙ КОМПОЗИЦІЙНИЙ ДИЗАЙН')%
format
% xm = [249.5 419.5
% alpha = 0.2;
[b, bint, r, rint, stats] = regress (Y, X, alpha); = [b bint] short = stats (3)
На підставі отриманої моделі
В
можна знайти оцінки координат екстремуму
format bank = [2 * b (2) 0 0, 0 2 * b (3) 0, 0 0 2 * b (4)],
В = [-b (5);-b (6);-b <7)]
хm = А ^ -1 * B = X0G '+ xm. * dX' = b (1) + sum (b (2:4). * xm (1:3)) + sum (b (5:7). * xm (1:3). ^ 2)
Графічна ілюстрація результатів дана на рис. 3 і 4. На цих малюнках ромбом позначені істинні координати екстремуму, задані в нашому випадку імітаційної моделлю. br/>В
Рис. 3
Задача 2
Перевірити гіпотезу про нормальність розподілу прибутку з продажів, що генеруються імітаційної моделлю profit3norm (D), де D - матриця дизайну експерименту.
Для проведення дослідів візьмемо точку факторного простору з координатами. bank = [300 490 580]; = 200; = ones (N...
