fsolve
В этом примере показано, как сгенерировать код С для решения систем нелинейных уравнений с fsolve
.
Система нелинейных уравнений, чтобы решить
Преобразуйте уравнения в форму F (x) = 0.
Запишите функцию, которая вычисляет левую сторону этих двух уравнений. Для генерации кода ваша программа должна выделить все массивы, когда они создаются и не должны изменять свои размеры после создания.
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)