Сгенерируйте код для fsolve

Сгенерируйте код для `fsolve`

Этот пример показывает, как сгенерировать код С для решения систем нелинейных уравнений с fsolve.

Уравнение для решения

Система нелинейных уравнений, которые нужно решить, является

$\begin{matrix} e^{- e^{- (x_{1} + x_{2})}} = x_{2} (1 + x_{1}^{2}) \\ x_{1} \cos (x_{2}) + x_{2} \sin (x_{1}) = \frac{1}{2} . \end{matrix}$

Преобразуйте уравнения в вид F (x) = 0.

$\begin{matrix} e^{- e^{- (x_{1} + x_{2})}} - x_{2} (1 + x_{1}^{2}) = 0 \\ x_{1} \cos (x_{2}) + x_{2} \sin (x_{1}) - \frac{1}{2} = 0. \end{matrix}$

Шаги генерации кода

Напишите функцию, которая вычисляет левую часть двух уравнений. Для генерации кода ваша программа должна выделять все массивы, когда они созданы, и не должна изменять их размеры после создания.
```
function F = root2d(x)
F = zeros(2,1); % Allocate return array
F(1) = exp(-exp(-(x(1)+x(2)))) - x(2)*(1+x(1)^2);
F(2) = x(1)*cos(x(2)) + x(2)*sin(x(1)) - 0.5;
end
```
Написание функции, которая настраивает задачу и вызовы fsolve. Функция должна ссылаться на root2d как указатель на функцию, а не как имя.
```
function [x,fval] = solveroot
options = optimoptions('fsolve','Algorithm','levenberg-marquardt','Display','off');
fun = @root2d;
rng default
x0 = rand(2,1);
[x,fval] = fsolve(fun,x0,options);
end
```
Создайте строение для генерации кода. В этом случае используйте 'mex'.
```
cfg = coder.config('mex');
```
Сгенерируйте код для solveroot функция.
```
codegen -config cfg solveroot
```
Протестируйте сгенерированный код, запустив сгенерированный файл с именем solveroot_mex.mexw64 или аналогичный.
```
[x,fval] = solveroot_mex
```
```
x =

    0.3532
    0.6061


fval =

   1.0e-14 *

   -0.1998
   -0.1887
```

См. также

fsolve | optimoptions | codegen (MATLAB CODER)

Документация по Optimization Toolbox

Поддержка

Памятка переводчика

1. Если смысл перевода понятен, то лучше оставьте как есть и не придирайтесь к словам, синонимам и тому подобному. О вкусах не спорим.

2. Не дополняйте перевод комментариями “от себя”. В исправлении не должно появляться дополнительных смыслов и комментариев, отсутствующих в оригинале. Такие правки не получится интегрировать в алгоритме автоматического перевода.

3. Сохраняйте структуру оригинального текста - например, не разбивайте одно предложение на два.

4. Не имеет смысла однотипное исправление перевода какого-то термина во всех предложениях. Исправляйте только в одном месте. Когда Вашу правку одобрят, это исправление будет алгоритмически распространено и на другие части документации.

5. По иным вопросам, например если надо исправить заблокированное для перевода слово, обратитесь к редакторам через форму технической поддержки.

Документация