Exponenta Banner
exponenta event banner

Создать код для fsolve

В этом примере показано, как генерировать код C для решения систем нелинейных уравнений с помощью fsolve.

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

Система нелинейных уравнений для решения

e − e − (x1 + x2) = x2 (1 + x12) x1cos (x2) + x2sin (x1) = 12.

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

e − e − (x1 + x2) − x2 (1 + x12) = 0x1cos (x2) + x2sin (x1) − 12 = 0.

Шаги создания кода

  1. Запишите функцию, которая вычисляет левую часть двух уравнений. Для создания кода программа должна выделять все массивы при их создании и не должна изменять их размеры после создания.

    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

  2. Запись функции, устанавливающей проблему и вызовы 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

  3. Создайте конфигурацию для создания кода. В этом случае используйте 'mex'.

    cfg = coder.config('mex');

  4. Создание кода для solveroot функция.

    codegen -config cfg solveroot

  5. Проверьте созданный код, запустив созданный файл с именем solveroot_mex.mexw64 или аналогичные.

    [x,fval] = solveroot_mex
    x =

    0.3532
    0.6061


fval =

   1.0e-14 *

   -0.1998
   -0.1887

См. также

| | (кодер MATLAB)

Связанные темы

Документация по панели инструментов оптимизации

Поддержка