Progress-servis55.ru

Новости из мира ПК
8 просмотров
Рейтинг статьи
1 звезда2 звезды3 звезды4 звезды5 звезд
Загрузка...

Fsolve matlab пример

Решение уравнений в matlab

Вы можете решать уравнения, содержащие переменные, с помощью команд
solve и fzero.

Разберем подробнее Matlab решение нелинейных уравнений, к примеру квадратного уравнения х 2 — 2х — 4 = 0, введите следующее:

syms x; solve (‘x^2 — 2*x -4=0’)

Разберем подробнее matlab решение линейных уравнений, к примеру вот такое уравнение х — 4 = 0, введите следующее:

syms x; solve (‘x -4=0’)

Обратите внимание, что уравнение, которое требуется решить, задано как
строка, то есть взято в одинарные кавычки. Ответ представляет собой точное
(символьное) решение 1+корень(5). Для получения числовых решений введите double
(ans) или vpa (ans), чтобы отобразить больше знаков. Ввод с командой solve
может также быть символьным выражением, но в этом случае программа MATLAB
потребует, чтобы правая часть выражения была заключена в скобки, и
фактически синтаксис решения уравнения х 2 — Зх = -7 будет выглядеть так:

syms x; solve (x^2 — 3*x + 7)

Ответ представляет собой точное (символьное) решение (3 + корень(19i))/2
(сложные числа, где буква i в ответе ставится для мнимой единицы V-1). Для
получения числовых решений введите double (ans) или vpa (ans), чтобы
отобразить больше знаков.
С помощью команды solve можно решать высокоуровневые полиномиальные
(многочленные) уравнения, равно как и многие другие типы уравнений. Можно
также решать уравнения, содержащие более чем одну переменную. Если
уравнений меньше, чем переменных, вам следует определить (как строки), какую
переменную (переменные) требуется вычислить. Например, введите solve ( ‘2*х — log (у) = 1’, ‘у’), чтобы решить уравнение 2х — log у = 1 для
переменной у при условии х. Подобным образом вы можете определить более чем
одно уравнение. Например:

[x, y] = solve (‘x^2 — y = 2’, ‘y — 2*x =5’)

Эта система уравнений имеет два решения. Программа MATLAB выдает решение,
выводя два значения х и два значения у для этих решений. Таким образом,
первое решение состоит из первого значения х и первого значения у. Вы можете
извлечь эти значения, введя в командную строку х (1) и у (1):

Второе решение можно извлечь, введя х (2) и у (2).
Обратите внимание, что в предыдущей команде solve мы назначили вывод в
векторной форме [х, у]. Если вы используете команду solve в системе
уравнений, не задавая вывод в векторной форме, в этом случае программа MATLAB не
отображает автоматически значения решения:

sol = solve (‘x^2 — y = 2’, ‘y — 2*x = 5’)

sol =
х: [2×1 sym]
у: [2×1 sym]

Чтобы увидеть векторы значений х и у, введите sol.x и sol.у. Чтобы увидеть
отдельные значения, введите sol.х (1) и sol.у (1), и т.п.

  • В этом примере вывод результата выполнения команды solve представляет собой структурный массив. Чтобы более подробно познакомиться с этим классом данных

Некоторые уравнения нельзя решить символически, и в таких случаях команда
solve пытается найти числовой ответ. Например:

solve (‘ sin (x) = 2 — x’)

Иногда бывает более одного решения, и вы можете не получить того, что
ожидаете, например:

solve (‘exp (-x) = sin (x) ‘ )

Ответ представляет собой комплексное число. Хотя оно является правильным
решением уравнения, существуют также решения, представленные
вещественными числами. Графики функций ехр (-х) и sin (x) показаны на Рис. 2.3;
каждая точка пересечения двух кривых представляет собой решение уравнения е -х = sin (х).

Вы можете в числовой форме найти (приблизительно) решения, показанные на
графике, с помощью команды fzero, которая ищет нулевое значение данной
функции в пределах заданного значения х. Решение уравнения е -х = sin (x)
равно нулю в функции е -х — sin (x), поэтому, чтобы найти приблизительное
решение при х = 0.5, введите следующее:

h = @(x) exp(-x) — sin(x);
fzero (h, 0.5)

Замените значение 0.5 на 3 и найдите следующее решение, и так далее.

Рис. 2.3. Две пересекающиеся кривые

Поэтому из выше всего сказанного можно сделать вывод, что вам необходимо просмотреть много дополнительной информации и альтернатив!

Fsolve matlab пример

предназначена функция fsolve.

Здесь x — вектор и F(x) — функция, которая возвращает значение вектора.

Наиболее полная информация о функции fsolve приведена в справочной системе MATLAB. Здесь приводится только краткое описание.

x = fsolve(fun,x0)
x = fsolve(fun,x0,options)
[x,fval,exitflag] = fsolve(. )
[x,fval,exitflag,output] = fsolve(. )
[x,fval,exitflag,output,jacobian] = fsolve(. )

fsolve находит корни (нули) системы нелинейных уравнений.

Таблица 4-1, Входные аргументы, содержит общее описание аргументов, передаваемых в fsolve. Данный подраздел приводит функционально-специфические детали для fun и options:

Подлежащая решению система уравнений.

fun есть такая функция, которая принимает вектор х и возвращает вектор F, нелинейные уравнения х. Функция fun может задаваться с помощью описателя функций

x = fsolve(@myfun,x0)
где myfun есть такая функция MATLAB, что

function F = myfun(x)
F = . % Расчет значений функции от x
fun так же может быть внутренним объектом.

Если к тому же может быть рассчитан Якобиан и установленная с помощью options = optimset(‘Jacobian’,’on’) опция options.Jacobian равна ‘on’, то функция fun во втором выходном аргументе должна возвращать значение Якобиана J, как матрицы от х.

Опции обеспечивают учет специфических деталей функции виде параметров options.

Таблица ниже содержат общее описание возвращаемых fsolve. аргументов. В этом разделе приводятся общие специфические детали для величин exitflag и output:

Описывает выходные условия.

  • > 0 Данная функция сходится к решению по х.
  • 0 Максимальное число оценки функции или итерации было превышено
  • LargeScale , поскольку она устанавливает преимущественное право, какой алгоритм использовать. Это только преимущественное право, поскольку должны быть выполнены определенные условия для того, что бы можно было использовать крупно-масштабный алгоритм. Для fsolve, поскольку система нелинейных уравнений не может быть недоопределенной, число уравнений (или иначе число возвращаемых fun элементовF) должно быть, по крайней мере, длины х, в противном следует использоваться средне-масштабный алгоритм.

    LargeScale. В случае установки ‘on’ используется, если это возможно, крупно-масштабный алгоритм. Для использования средне-масштабного алгоритма устанавливается значение ‘off’.

    Читать еще:  Matlab plot параметры

    Medium-Scale и Large-Scale Algorithms . Эти параметры используются как для средне-масштабного, так и крупно-масштабного алгоритмов:

    Проводится печать диагностической информации о минимизируемой функции.

    Уровень отображения. ‘off’ отображение не производится, ‘iter’ отображение проводится на каждой итерации, ‘final’ (принимается по умолчанию) отображение только конечной информации.

    При установке ‘on’, fsolve использует заданный пользователем Якобиан (определенный в fun), или информацию об Якобиане (при использовании JacobMult) для целевой функции. При установке ‘off’, fsolve аппроксимирует Якобиан с помощью конечных разностей.

    Максимально число допустимых расчетов функции.

    Максимальное число допустимых итераций.

    Конечное допустимое отклонение по значению функции.

    Конечное допустимое отклонение по значению х.

    Large-Scale Algorithm Only . Эти параметры используются только для крупно-масштабного алгоритма

    Указатель функции для функции множителей Якобиана.

    В случае крупно-масштабной задачи данная функция вычисляет произведение матриц Якобиана J*Y, J’*Y или J’*(J*Y) без действительного формирования J. Такая функция имеет форму

    где Jinfo и дополнительные параметры p1,p2. включают в себя используемые для расчета J*Y (или J’*Y или J’*(J*Y)) матрицы. Первый аргумент Jinfo должен быть таким же, как и возвращаемый целевой функцией fun второй аргумент.

    Переменные p1,p2. есть дополнительные параметры, передаваемые в fsolve (и в fun).

    fsolve (fun. options,p1,p2. )

    Y есть матрица с тем же самым числом строк, что и размерность данной задачи.

    flag определяет какое произведение нужно вычислять.

    If flag == 0 then W = J’*(J*Y).
    If flag > 0 then W = J*Y.
    If flag

    Примечание. ‘Jacobian’ должен быть установлен как ‘on’ для того, чтобы передать Jinfo из fun в jmfun.

    В каждом случае J явно не формируется.

    fsolve использует Jinfo для расчета предварительных данных.

    В качестве примера смотри Нелинейную Оптимизацию с Компактной, но Структурированной Матрицей Гессе и Ограничениями Типа Равенств.

    Разреженные шаблоны Якобиана для конечного дифференцирования. В случае, если в fun неразумно вычислять матрицу Якобиана J, то lsqnonlin может аппроксимировать J через заданные разреженные конечные разности и структура J — т.е. расположение не нулей – обеспечиваются как значения для JacobPattern. В наихудшем случае, когда эта структура неизвестна, можно установить JacobPattern, что бы плотная матрица и полные конечно-разностные аппроксимации вычислялись на каждой итерации (что принимается по умолчанию в случае неустановки JacobPattern). Это может быть чрезвычайно затратным для больших задач, поэтому обычно заслуживает внимания усилие по определению разреженной структуры.

    Максимальное число PCG (предварительно сопряженный градиент) итераций (Смотри ниже раздел Алгоритм).

    Верхняя полоса предварительной обработки для PCG. По умолчанию используется диагональная начальная подготовка (верхняя полоса из 0). Для некоторых задач увеличение полосы снижает число итераций PCG.

    Конечное допустимое число итераций PCG.

    Типичные значения х.

    Medium-Scale Algorithm Only. Эти параметры используются только для средне-масштабного алгоритма.

    Сравнение вводимых пользователем производных (Якобиана) с конечноразностными производными.

    Максимальное изменение в переменных для конечных-разностей.

    Минимальное изменение в переменных для конечных-разностей.

    Выбор алгоритма Левенберга-Макуарда вместо Гаусса-Ньютона.

    Выбор алгоритма линейного поиска.

    Пример 1. В данном примере находятся нули для системы из двух уравнений и двух неизвестных

    Таким образом, необходимо решить следующую систему уравнений от х

    Начнем с точки x0 = [-5 -5].

    Сперва запишем М-файл для расчета F или значений уравнений от х

    function F = myfun(x)
    F = [2*x(1) — x(2) — exp(-x(1));
    -x(1) + 2*x(2) — exp(-x(2))];

    Далее вызовем подпрограмму оптимизации

    x0 = [-5; -5]; % примем начальное приближение за решение
    options=optimset(‘Display’,’iter’); %Опция выходного отображения
    [x,fval] = fsolve(@myfun,x0,options) % вызов оптимизатора

    после 28 обращений к функциям нули будут найдены.

    Оптимальность первого порядка

    Оптимизация завершена успешно:

    Относительное изменение значений функции меньше, чем в OPTIONS.TolFun

    fval =
    1.0e-008 *
    -0.5320
    -0.5320

    Оптимизация завершена успешно:

    Пример 2. Найти матрицу x, которая удовлетворяет уравнению

    начнем с точки x= [1,1; 1,1].

    Сперва запишем М-файл, необходимый для расчета решаемых уравнений.

    function F = myfun(x)
    F = x*x*x-[1,2;3,4];

    Далее запустим подпрограмму оптимизации

    x0 = ones(2,2); % примем начальное приближение за решение
    options = optimset(‘Display’,’off’); % Выключим отображение
    [x,Fval,exitflag] = fsolve(@myfun,x0,options)

    x =
    -0.1291 0.8602
    1.2903 1.1612

    Fval =
    1.0e-03 *
    0.1541 -0.1163
    0.0109 -0.0243

    и остаток будет близок к нулю.

    sum(sum(Fval.*Fval))
    ans =
    3.7974e-008

    Данный метод основан на алгоритме метода нелинейных среднеквадратичных отклонений, используемого в lsqnonlin. Достоинство используемого метода среднеквадратичных отклонений состоит в том, что данная система уравнений не равна нулю вследствие малых погрешностей и поэтому алгоритм приходит в точку с малым остатком. Однако, если Якобиан системы является вырожденным, то алгоритм может сходиться в точку, не являющуюся решением системы уравнений (Смотри ниже Ограничения и Диагностика).

    Крупно-масштабная оптимизация. По умолчанию fsolve выберет крупно-масштабный алгоритм. Данный алгоритм является реализацией метода доверительных подпространств и основан на методе внутренних отражений Ньютона, описанного в [1], [2]. Каждая итерация включает в себя приближенное решение крупной линейной системы с помощью метода предварительно сопряженных градиентов(PCG).

    Средне-масштабная оптимизация. fsolve при установке options.LargeScale как ‘off’ использует метод Гаусса-Ньютона [3] с линейным поиском. В качестве альтернативы может быть выбран метод Левенберга-Макуарда [4], [5], [6] с линейным поиском. Выбор алгоритма проводится при помощи установки опции options.LevenbergMarquardt. Установка options.LevenbergMarquardt как ‘on’ (и опции options.LargeScale как ‘off’) выбирается метод Левенберга-Макуарда.

    Принимаемый по умолчанию алгоритм линейного поиска, т.е. опция options.LineSearchType установлена как ‘quadcubic’, обеспечивается методом смешанной квадратичной и кубической полиноминальной интерполяции и экстраполяции. Защищенный кубической полиномиальный метод может быть выбран установкой опции LineSearchType как ‘cubicpoly’. В общем случае, данный метод требует меньшего расчета функций, но большего обращения к расчету градиента. Таким образом, если градиенты приведены и могут быть вычислены без больших затрат, то метод с кубическим полиномиальным линейным поиском является предпочтительным.

    fsolve может сходиться к ненулевой точке и давать следующие сообщения:

    Оптимизатор затормозился в точке, которая не является минимумом.

    Повторите расчет с новой начальной точки.

    В этом случае снова выполните fsolve но с другой начальной точки.

    1. Разрешаемая функция должна быть непрерывной. В случае успеха fsolve дает только один корень. fsolve может сходиться к ненулевой точке, то в этом случае необходимо пытаться стартовать с другой точки.
    2. fsolve оперирует только с реальными переменными. Если х имеет комплексные значения, то эти переменные должны быть разделены на мнимую и реальные части.
    Читать еще:  Средства защиты конфиденциальной информации

    Иллюстрированный самоучитель по MatLab

    Вычисление нулей функции одной переменной

    Ряд функций системы MATLAB предназначен для работы с функциями. По аналогии с дескрипторами графических объектов могут использоваться объекты класса дескрипторов функций, задаваемых с помощью символа @, например:

    Примечание
    Подфункциями понимаются как встроенные функции, например sin(x) или ехр(х),так и функции пользователя, например f(x), задаваемые как т-файлы-функции
    .

    Численные значения таких функций, заданных дескрипторами, вычисляются с помощью функции feval:

    Для совместимости с прежними версиями можно записывать функции в символьном виде в апострофах, использование функции eval для их вычисления может быть более наглядно, не нужно создавать m-файл, но в учебном курсе мы будем стараться использовать новую нотацию, с использованием дескрипторов функций и feval, так как при этом программирование становится «более объектно-ориентированным», повышается скорость, точность и надежность численных методов. Поэтому, хотя везде в нижеследующем тексте вместо @fun можно подставить и символьное значение функции в апострофах, мы будем использовать нотацию @fun в дидактических целях. Все же иногда в интерактивном режиме можно использовать старую запись, чтобы не создавать m-файл функции.

    Довольно часто возникает задача решения нелинейного уравнения вида f(x) = 0 или /, (г) =/ 2 (дг). Последнее, однако, можно свести к виду f(x) =f 1 (х) – f 2 (х) = 0. Таким образом, данная задача сводится к нахождению значений аргумента х функции f(x) одной переменной, при котором значение функции равно нулю. Соответствующая функция MATLAB, решающая данную задачу, приведена ниже:

    • fzero(@fun,x) – возвращает уточненное значение х, при котором достигается нуль функции fun, представленной в символьном виде, при начальном значении аргумента х. Возвращенное значение близко к точке, где функция меняет знак, или равно NaN, если такая точка не найдена;
    • fzero(@fun,[x1 x2]) – возвращает значение х, при котором fun(x)=0 с заданием интервала поиска с помощью вектора x=[xl х2], такого, что знак fun(x(D) отличается от знака fun(x(2)). Если это не так, выдается сообщение об ошибке. Вызов функции fzero с интервалом гарантирует, что fzero возвратит значение, близкое к точке, где fun изменяет знак;
    • fzero(@fun,x.tol) – возвращает результат с заданной погрешностью tol;
    • fzero(@fun,x.tol.trace) – выдает на экран информацию о каждой итерации;
    • fzero(@fun,х.tol.trace,Р1.Р2,…) – предусматривает дополнительные аргументы, передаваемые в функцию fun(x.Pl,P2,…). При задании пустой матрицы для tol или trace используются значения по умолчанию. Пример:

    Для функции fzero ноль рассматривается как точка, где график функции fun пересекает ось х, а не касается ее. В зависимости от формы задания функции fzero реализуются следующие хорошо известные численные методы поиска нуля функции: деления отрезка пополам, секущей и обратной квадратичной интерполяции. Приведенный ниже пример показывает приближенное вычисление р/2 из решения уравнения cos(x)=0 с представлением косинуса дескриптором:

    В более сложных случаях настоятельно рекомендуется строить график функции f(x) для приближенного определения корней и интервалов, в пределах которых они находятся. Ниже дан пример такого рода (следующий листинг представляет собой содержимое m-файла fun1.m):

    Из рисунка нетрудно заметить, что значения корней заключены в интервалах [0.5 1], [2 3] и [5 6]. Найдем их, используя функцию fzero:

    Обратите внимание на то, что корень х3 найден двумя способами и что его значения в третьем знаке после десятичной точки отличаются в пределах заданной погрешности tol =0.001. К сожалению, сразу найти все корни функция fzero не в состоянии. Решим эту же систему при помощи функции fsolve из пакета Optimization Toolbox, которая решает систему нелинейных уравнений вида f(x)=0 методом наименьших квадратов, ищет не только точки пересечения, но и точки касания, fsolve имеет почти те же параметры (дополнительный параметр – задание якобиана) и почти ту же запись, что и функция lsqnonneg, подробно рассмотренная ранее. Пример:

    Для решения систем нелинейных уравнений следует также использовать функцию solve из пакета Symbolic Math Toolbox. Эта функция способна выдавать результат в символьной форме, а если такого нет, то она позволяет получить решение в численном виде. Пример:

    [Пакеты расширения Symbolic Math ToolBox и Extended Symbolic Math Toolbox MATLAB 6.0 используют ядро Maple V Release 5 [30-35] и являются поэтому исключением: они пока не поддерживают новую нотацию с использованием дескрипторов функций. – Примеч. ред.]

    Решение систем нелинейных уравнений в Matlab

    Доброго времени суток! В этой статье мы поговорим о решении систем нелинейных алгебраических уравнений в Matlab. Вслед за решением нелинейных уравнений, переходим к их системам, рассмотрим несколько методов реализации в Matlab.

    Общая информация

    Итак, в прошлой статье мы рассмотрели нелинейные уравнения и теперь необходимо решить системы таких уравнений. Система представляет собой набор нелинейных уравнений (их может быть два или более), для которых иногда возможно найти решение, которое будет подходить ко всем уравнениям в системе.

    В стандартном виде, количество неизвестных переменных равно количеству уравнений в системе. Необходимо найти набор неизвестных переменных, которые при подставлении в уравнения будут приближать значение уравнения к 0. Иногда таких наборов может быть несколько, даже бесконечно много, а иногда решений не существует.

    Чтобы решить СНАУ, необходимо воспользоваться итеративными методами. Это методы, которые за определенное количество шагов получают решение с определенной точностью. Также очень важно при решении задать достаточно близкое начальное приближение, то есть такой набор переменных, которые близки к решению. Если решается система из 2 уравнений, то приближение находится с помощью построение графика двух функций.

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

    Читать еще:  Система уравнений matlab

    Оператор Matlab для решения СНАУ

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

    Решить систему нелинейных уравнений с точность 10 -2 :
    cos(x-1) + y = 0.5
    x-cos(y) = 3

    Нам дана система из 2 нелинейных уравнений и сначала лучше всего построить график. Воспользуемся командой ezplot в Matlab, только не забудем преобразовать уравнения к стандартному виду, где правая часть равна 0:

    Функция ezplot строит график, принимая символьную запись уравнения, а для задания цвета и толщины линии воспользуемся функцией set. Посмотрим на вывод:

    Как видно из графика, есть одно пересечение функций — то есть одно единственное решение данной системы нелинейных уравнений. И, как было сказано, по графику найдем приближение. Возьмем его как (3.0, 1.0). Теперь найдем решение с его помощью:

    Создадим функцию m-файлом fun.m и поместим туда следующий код:

    Заметьте, что эта функция принимает вектор приближений и возвращает вектор значений функции. То есть, вместо x здесь x(1), а вместо y — x(2). Это необходимо, потому что fsolve работает с векторами, а не с отдельными переменными.

    И наконец, допишем функцию fsolve к коду построения графика таким образом:

    Таким образом у нас образуется два m-файла. Первый строит график и вызывает функцию fsolve, а второй необходим для расчета самих значений функций. Если вы что-то не поняли, то в конце статьи будут исходники.

    И в конце, приведем результаты:

    xr (это вектор решений) =
    3.3559 1.2069

    fr (это значения функций при таких xr, они должны быть близки к 0) =
    1.0e-09 *
    0.5420 0.6829

    ex (параметр сходимости, если он равен 1, то все сошлось) =
    1

    И, как же без графика с ответом:

    Метод простых итераций в Matlab для решения СНАУ

    Теперь переходим к методам, которые запрограммируем сами. Первый из них — метод простых итераций. Он заключается в том, что итеративно приближается к решению, конечно же, с заданной точностью. Алгоритм метода достаточно прост:

    1. Если возможно, строим график.
    2. Из каждого уравнения выражаем неизвестную переменную след. образом: из 1 уравнения выражаем x1, из второго — x2, и т.д.
    3. Выбираем начальное приближение X0, например (3.0 1.0)
    4. Рассчитываем значение x1, x2. xn, которые получили на шаге 2, подставив значения из приближения X0.
    5. Проверяем условие сходимости, (X-X0) должно быть меньше точности
    6. Если 5 пункт не выполнился, то повторяем 4 пункт.

    И перейдем к практике, тут станет все понятнее.
    Решить систему нелинейных уравнений методом простых итераций с точность 10 -2 :
    cos(x-1) + y = 0.5
    x-cos(y) = 3

    График мы уже строили в предыдущем пункте, поэтому переходим к преобразованию. Увидим, что x из первого уравнения выразить сложно, поэтому поменяем местами уравнения, это не повлияет на решение:

    x-cos(y) = 3
    cos(x-1) + y = 0.5

    Далее приведем код в Matlab:

    В этой части мы выразили x1 и x2 (у нас это ‘x’ и ‘y’) и задали точность.

    В этой части в цикле выполняются пункты 4-6. То есть итеративно меняются значения x и y, пока отличия от предыдущего значения не станет меньше заданной точности.

    k = 10
    x = 3.3587
    y = 1.2088

    Как видно, результаты немного отличаются от предыдущего пункта. Это связано с заданной точностью, можете попробовать поменять точность и увидите, что результаты станут такими же, как и при решении стандартным методом Matlab.

    Метод Ньютона в Matlab для решения СНАУ

    Решение систем нелинейных уравнений в Matlab методом Ньютона является более эффективным, чем использование метода простых итераций. Сразу же представим алгоритм, а затем перейдем к реализации.

    1. Если возможно, строим график.
    2. Выбираем начальное приближение X0, например (3.0 1.0)
    3. Рассчитываем матрицу Якоби w, это матрица частных производных каждого уравнения, считаем ее определитель для X0.
    4. Находим вектор приращений, который рассчитывается как dx = -w -1 * f(X0)
    5. Находим вектор решения X = X0 + dx
    6. Проверяем условие сходимости, (X-X0) должно быть меньше точности

    Далее, решим тот же пример, что и в предыдущих пунктах. Его график мы уже строили и начальное приближение останется таким же.
    Решить систему нелинейных уравнений методом Ньютона с точность 10 -2 :
    cos(x-1) + y = 0.5
    x-cos(y) = 3

    Перейдем к коду:

    Сначала зададим начальное приближение. Затем необходимо просчитать матрицу Якоби, то есть частные производные по всем переменным. Воспользуемся символьным дифференцированием в Matlab, а именно командой diff с использованием символьных переменных.

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

    В этой части кода выполняем первую итерацию, чтобы получить вектор решения и сравнивать его с вектором начального приближения. Отметим, чтобы посчитать значение символьной функции в Matlab, необходимо воспользоваться функцией subs. Эта функция заменяет переменную на числовое значение. Затем функция double рассчитает это числовое значение.

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

    В этой части кода выполняется цикл по расчету решения с заданной точностью. Еще раз отметим, что если в первой итерации до цикла были использованы функции diff, double и subs для вычисления производной в Matlab, то в самом цикле матрица якоби была явно задана этими частными производными. Это сделано, чтобы показать возможности среды Matlab.

    За 3 итерации достигнут правильный результат. Также важно сказать, что иногда такие методы могут зацикливаться и не закончить расчеты. Чтобы такого не было, мы прописали проверку на количество итераций и запретили выполнение более 100 итераций.

    Заключение

    В этой статье мы познакомились с основными понятиями систем нелинейных алгебраических уравнений в Matlab. Рассмотрели несколько вариантов их решения, как стандартными операторами Matlab, так и запрограммированными методами простых итераций и Ньютона.

  • Ссылка на основную публикацию
    Adblock
    detector