Lab5-7 completed

Author: Throttle

lab05_07/report/parts/20-analysis.tex

@@ -127,6 +127,67 @@ \section{Часть 2}
 
 Устранение разрыва функции производится путём замены аргумента, при выполнении условия $|x_i| < \varepsilon$, на константное значение $\varepsilon$, где $\varepsilon = 1e-1$. При $x_1 \times x_2 < 0$ значение функции устанавливается равным $45$.
 
+В таблице \ref{tb:tab21} представлены результаты работы рассмотренных методов для точности $0.01$.
+
+\begin{table}[!ht]
+\caption{Сводная таблица результатов работы методов}
+\begin{tabular}{|p{0.03\textwidth}|p{0.39\textwidth}|p{0.20\textwidth}|p{0.11\textwidth}|p{0.13\textwidth}|}
+\hline
+\No & Метод & Количество вычислений функции & $X$ & $f(X)$\\
+\hline
+1 & правильный симплекс & 25 & [2.794894, 0.143244] & 2.384087 \\
+\hline
+2 & деформируемый симплекс & 24 & 0.[2.465024, 0.166927] & 2.38653 \\
+\hline
+3 & случайного поиска с возвратом &  114 & [1.615551, 0.235544] & 2.395212 \\
+\hline
+5 & fminsearch & 82 & [3.970522, 0.100013] & 2.382102 \\
+\hline
+\end{tabular}
+\label{tb:tab21}
+\end{table}
+
+В таблице \ref{tb:tab22} представлены результаты работы рассмотренных методов для точности $0.0001$.
+
+\begin{table}[!ht]
+\caption{Сводная таблица результатов работы методов}
+\begin{tabular}{|p{0.03\textwidth}|p{0.39\textwidth}|p{0.20\textwidth}|p{0.11\textwidth}|p{0.13\textwidth}|}
+\hline
+\No & Метод & Количество вычислений функции & $X$ & $f(X)$\\
+\hline
+1 & правильный симплекс & 37 & [2.794894, 0.143244] & 2.384087 \\
+\hline
+2 & деформируемый симплекс & 30 & 0.[2.465024, 0.166927] & 2.386538 \\
+\hline
+3 & случайного поиска с возвратом &  559 & [1.877914, 0.206477] & 2.390224 \\
+\hline
+5 & fminsearch & 143 & [3.968544, 0.100000] & 2.382102 \\
+\hline
+\end{tabular}
+\label{tb:tab22}
+\end{table}
+
+
+В таблице \ref{tb:tab23} представлены результаты работы рассмотренных методов для точности $0.000001$.
+
+\begin{table}[!ht]
+\caption{Сводная таблица результатов работы методов}
+\begin{tabular}{|p{0.03\textwidth}|p{0.39\textwidth}|p{0.20\textwidth}|p{0.11\textwidth}|p{0.13\textwidth}|}
+\hline
+\No & Метод & Количество вычислений функции & $X$ & $f(X)$\\
+\hline
+1 & правильный симплекс & 76 & [2.801648, 0.140676] & 2.383914 \\
+\hline
+2 & деформируемый симплекс & 110 & 0.[[3.962550, 0.100045] & 2.382104 \\
+\hline
+3 & случайного поиска с возвратом &  1704 & [1.390272, 0.270747] & 2.402598 \\
+\hline
+5 & fminsearch & 207 & [3.968503, 0.100000] & 2.382102 \\
+\hline
+\end{tabular}
+\label{tb:tab23}
+\end{table}
+
 
 %
 % % В начале раздела  можно напомнить его цель

lab05_07/report/rpz.pdf

@@ -15,28 +15,28 @@
 %% Первая часть
 %optimizer.draw(-6, 2, 0, 9, 100);
 % использование метода минизации по правильному симплексу
-%optimizer.simple_simplex(x0, a, eps3);
+%optimizer.simple_simplex(x0, a, eps2);
 
 % использование метода минизации по деформируемому симплексу
 %optimizer.downhill_simplex(x0, a, eps3);
 
 % использование метода минимизациии случайного поиска в возвратом
-%optimizer.returnrandsearch(x0, a, eps3);
+%optimizer.rand_search(x0, a, eps3);
 
 %optimizer.fminsearch(x0, eps3);
 
 %% Вторая часть
 optimizer.draw(-2, 5, -2, 5, 100);
 
 % использование метода минизации по правильному симплексу
-%optimizer.simple_simplex(x02, a2, eps2);
+%optimizer.simple_simplex(x02, a2, eps3);
 
 % использование метода минизации по деформируемому симплексу
-%optimizer.downhill_simplex(x02, a2, eps2);
+%optimizer.downhill_simplex(x02, a2, eps3);
 
 % использование метода минимизациии случайного поиска в возвратом
-%optimizer.returnrandsearch(x02, a2, eps2);
-
-optimizer.fminsearch(x02, eps2);
+optimizer.rand_search(x02, a2, eps3);
+%
+%optimizer.fminsearch(x02, eps3);
 
 hold off;

lab05_07/source/launcher.m

@@ -26,17 +26,30 @@
             %y = 4*x(1)*x(2) + 7*x(1)*x(1) + 4*x(2)*x(2) + 6*sqrt(5)*x(1) - 12*sqrt(5)*x(2) + 51;
             %return;
             % 2àÿ ôóíêöèÿ
-            if abs(x(1)) < this.eps
-                x(1) = x(1) + sign(x(1)) * this.eps;
-            end
 
-            if abs(x(2)) < this.eps
-                x(2) = x(2) + sign(x(2)) * this.eps;
-            end
 
-            if x(1) < 0 || x(2) < 0
+            if x(1) <= 0 || x(2) <= 0
                 y = 500;
             else
+                if abs(x(1)) < this.eps
+                    s1 = sign(x(1));
+                    if abs(s1) > 0
+                        x(1) = x(1) + sign(x(1)) * this.eps;
+                    else
+                        x(1) = this.eps;
+                    end
+                
+                end
+
+                if abs(x(2)) < this.eps
+                    s2 = sign(x(2));
+                    if abs(s2)>0
+                        x(2) = x(2) + sign(x(2)) * this.eps;
+                    else
+                        x(2) = this.eps;
+                    end
+
+                end
                 y = x(2)*x(2)*x(2) + 2 * x(1) * x(2) + 1 / sqrt(x(1) * x(2));
             end
         end
@@ -61,8 +74,8 @@ function draw(this, a1, b1, a2, b2, n)
                 end
             end
 
-            figure(2);
-            mesh(x2,x1,f);
+            %figure(2);
+            %mesh(x2,x1,f);
 
             figure(1);
             levels = -30:5:50;