Создание функции и находит свою производную

Создайте вектор $$x$$ и матрица $$A$$ как переменные символьной матрицы. Задайте функцию $\mathit{f}\left(\mathit{x}\right)=\sin \left({\mathit{x}}^{\mathit{T}}\mathit{A}\mathit{x}\right)$.

syms x [2 1] matrix
syms A [2 2] matrix
f = sin(x.'*A*x)

f = $$\sin \left(x^{\mathrm{T}}\hspace{0.17em}A\hspace{0.17em}x\right)$$

Вычислите производный D из функции $$f(x)$$ относительно вектора $$x$$. Производный D отображен в компактном матричном обозначении в терминах $$x$$ и $$A$$.

D = diff(f,x)

D = $$\cos \left(x^{\mathrm{T}}\hspace{0.17em}A\hspace{0.17em}x\right)\hspace{0.17em}\left(x^{\mathrm{T}}\hspace{0.17em}A+x^{\mathrm{T}}\hspace{0.17em}{A}^{\mathrm{T}}\right)$$

Преобразуйте symmatrix Объекты к sym Объекты

Переменные x символьной матрицыAF, и D symmatrix объекты. Эти объекты представляют матрицы, векторы и скаляры в компактном матричном обозначении. Чтобы показать компоненты этих переменных, преобразуйте symmatrix объекты к sym объекты с помощью symmatrix2sym.

xsym = symmatrix2sym(x)

xsym = 
$$\left(\begin{array}{c}{x}_1\\ x_2\end{array}\right)$$

Asym = symmatrix2sym(A)

Asym = 
$$\left(\begin{array}{cc}{A}_{1,1}& A_{1,2}\\ A_{2,1}& A_{2,2}\end{array}\right)$$

fsym = symmatrix2sym(f)

fsym = $$\sin \left(x_1\hspace{0.17em}\left(A_{1,1}\hspace{0.17em}{x}_1+A_{1,2}\hspace{0.17em}{x}_2\right)+x_2\hspace{0.17em}\left(A_{2,1}\hspace{0.17em}{x}_1+A_{2,2}\hspace{0.17em}{x}_2\right)\right)$$

Dsym = symmatrix2sym(D)

Dsym = $$\left(\begin{array}{cc}\cos \left(x_1\hspace{0.17em}\left(A_{1,1}\hspace{0.17em}{x}_1+A_{1,2}\hspace{0.17em}{x}_2\right)+x_2\hspace{0.17em}\left(A_{2,1}\hspace{0.17em}{x}_1+A_{2,2}\hspace{0.17em}{x}_2\right)\right)\hspace{0.17em}\left(2\hspace{0.17em}{A}_{1,1}\hspace{0.17em}{x}_1+A_{1,2}\hspace{0.17em}{x}_2+A_{2,1}\hspace{0.17em}{x}_2\right)& \cos \left(x_1\hspace{0.17em}\left(A_{1,1}\hspace{0.17em}{x}_1+A_{1,2}\hspace{0.17em}{x}_2\right)+x_2\hspace{0.17em}\left(A_{2,1}\hspace{0.17em}{x}_1+A_{2,2}\hspace{0.17em}{x}_2\right)\right)\hspace{0.17em}\left(A_{1,2}\hspace{0.17em}{x}_1+A_{2,1}\hspace{0.17em}{x}_1+2\hspace{0.17em}{A}_{2,2}\hspace{0.17em}{x}_2\right)\end{array}\right)$$

Замените числовыми значениями и найдите минимум

Предположим, что вы интересуетесь случаем где значение $$A$$ [2 -1; 0 3]. Замените этим значением в функциональный fsym.

fsym = subs(fsym,Asym,[2 -1; 0 3])

fsym = $$\sin \left(3\hspace{0.17em}{x_2}^2+x_1\hspace{0.17em}\left(2\hspace{0.17em}{x}_1-x_2\right)\right)$$

Замените значением $$A$$ в производный Dsym

Dsym = subs(Dsym,Asym,[2 -1; 0 3])

Dsym = $$\left(\begin{array}{cc}\cos \left(3\hspace{0.17em}{x_2}^2+x_1\hspace{0.17em}\left(2\hspace{0.17em}{x}_1-x_2\right)\right)\hspace{0.17em}\left(4\hspace{0.17em}{x}_1-x_2\right)& -\cos \left(3\hspace{0.17em}{x_2}^2+x_1\hspace{0.17em}\left(2\hspace{0.17em}{x}_1-x_2\right)\right)\hspace{0.17em}\left(x_1-6\hspace{0.17em}{x}_2\right)\end{array}\right)$$

Затем примените символьный функциональный solve получить корень производной.

[xmin,ymin] = solve(Dsym,xsym,'PrincipalValue',true);
x0 = [xmin; ymin]

x0 = 
$$\left(\begin{array}{c}0\\ 0\end{array}\right)$$

Постройте функцию $$f(x)$$ вместе с решением для экстремума $$x_0$$. Установите интервал графика на $$-1<x_1<1$$ и $$-1<x_2<1$$ в качестве второго аргумента fsurf. Используйте fplot3 построить координаты решения для экстремума.

fsurf(fsym,[-1 1 -1 1])
hold on
fplot3(xmin,ymin,subs(fsym,xsym,x0),'ro')
view([-68 13])

Аппроксимированная функция около ее минимума

Можно аппроксимировать многомерную функцию вокруг точки $$x_0$$ с многочленом с помощью Разложения Тейлора.

Здесь, термин $$\nabla f(x_0)$$ вектор градиента, и $$H(f(x_0))$$ матрица Гессиана многомерной функции $\mathit{f}\left(\mathit{x}\right)$ вычисленный в $$x_0$$.

Найдите матрицу Гессиана и возвратите результат как переменную символьной матрицы.

H = diff(f,x,x.')

H = $$-\sin \left(x^{\mathrm{T}}\hspace{0.17em}A\hspace{0.17em}x\right)\hspace{0.17em}\left(A^{\mathrm{T}}\hspace{0.17em}x+A\hspace{0.17em}x\right)\hspace{0.17em}\left(x^{\mathrm{T}}\hspace{0.17em}A+x^{\mathrm{T}}\hspace{0.17em}{A}^{\mathrm{T}}\right)+\cos \left(x^{\mathrm{T}}\hspace{0.17em}A\hspace{0.17em}x\right)\hspace{0.17em}\left(A^{\mathrm{T}}+A\right)$$

Преобразуйте матрицу Гессиана $$H(f(x_0))$$ к sym тип данных, который представляет матрицу в ее форме компонента. Используйте subs оценивать матрицу Гессиана для $$A$$= [2 -1; 0 3] в минимальной точке $$x_0$$.

Hsym = symmatrix2sym(H);
Hsym = subs(Hsym,Asym,[2 -1; 0 3]);
H0 = subs(Hsym,xsym,x0)

H0 = 
$$\left(\begin{array}{cc}4& -1\\ -1& 6\end{array}\right)$$

Оцените вектор градиента $$\nabla f(x_0)$$ в $$x_0$$.

D0 = subs(Dsym,xsym,x0)

D0 = $$\left(\begin{array}{cc}0& 0\end{array}\right)$$

Вычислите приближение Тейлора к функции около ее минимума.

fapprox = subs(fsym,xsym,x0) + D0*(xsym-x0) + 1/2*(xsym-x0).'*H0*(xsym-x0)

fapprox = 
$$x_1\hspace{0.17em}\left(2\hspace{0.17em}{x}_1-\frac{x_2}{2}\right)-x_2\hspace{0.17em}\left(\frac{x_1}{2}-3\hspace{0.17em}{x}_2\right)$$

Постройте приближение функций на том же графике, который показывает $$f(x)$$ и $$x_0$$.

hold on
fsurf(fapprox,[-1 1 -1 1])
zlim([-1 3])