На рис. 1.3.5-4 приведен script-файл с именем zadasha2, где введена матрица, а, состоящая из трех строк и четырех столбцов. Организован цикл по количеству формируемых элементов массива b путем обращения к функции sred_ar(). В функцию передается массив а, номер строки (i) и количество элементов в строке (m). Вывод элементов массива b предусмотрен в столбец.

Рис. 1.3.5-4. Программа формирования массива b

Функция sred_ar() (рис. 1.3.5-5) предназначена для формирования i-го элемента массива b, равного среднему арифметическому элементов строки массива а.

Рис. 1.3.5-5. Функция sred_ar(), вычисляющая среднее арифметическое
элементов строки массива a

В результате запуска script-файле с именем zadasha2 в окно Command Window выводится столбец элементов массива b

Рис. 1.3.5-6. Запуск script-файла zadasha2

Пример. 1.3.5-3. Задать действительные числа a,b, натуральное n (a<b) и вычислить выражение , где , если

Решение задачи требует разработки двух функций: fab(a, h, i), предназначенной для вычисления i-го слагаемого (рис. 1.3.5-7) и sumf(a, h, n), предназначенной для вычисления заданного выражения (рис. 1.3.5-8).

Рис. 1.3.5-7. Функция fab(), вычисляющая значение i-го слагаемого

Рис. 1.3.5-8. Функция sumf(), вычисляющая заданное выражение

Запуск на выполнение осуществляется из командного окна к функции sumf(). Предварительно переменным a,b и n присваиваются числовые значения. Проверка правильности ввода исходных данных предусмотрена в функции sumf(). Вычисления выполняются, и результат выводится на экран только в случае, если b>a, иначе в командной строке появляется сообщение «Ошибка в исходных данных b<a» (рис. 1.3.5-9).

Рис. 1.3.5-9. Запуск функции sumf() на выполнение

1.3.6. Лабораторная работа по теме

«Средства алгоритмизации и программирования

в Matlab»

1. Вопросы, подлежащие изучению

1. Виды m-файлов.

2. Создание и сохранение новых, и открытие ранее созданных m-файлов.

3. Особенности script-файлов и m-функций.

4. Запуск на выполнение script-файла из текстового редактора.

5. Запуск на выполнение script-файла из командного окна.

6. Обращения к script-файлам и m-функциям.

7. Средства языка программирования в системе Matlab.

8. Основные операторы m-языка их назначение и форматы.

2. Общее задание

1. Изучите материал Темы 1.3 (п.п. 1.3.1 – 1.3.5).

2. Выберите индивидуальное задание из табл. 1.3.6-1.

3. Разработайте m-функции для реализации стандартных алгоритмов: вычисления конечных сумм, разветвлений, поиска минимума и максимума в последовательности данных и т.п.

4. Введите и сохраните m-функции на внешнем носителе.

5. Создайте новый script-файл, в который введите код программы, описывающий логику решения поставленной задачи.

6. Сохраните script-файл в текущем каталоге.

7. Произведите отладку script-файла, запуская его на выполнение из текстового редактора командой Run.

8. Подготовьте и введите исходные данные для решения поставленной задачи;

9. Выполните script-файл из командной строки окна Command Window.

10. Сохраните текст рабочего окна на внешнем носителе.

11. Предоставьте результаты работы преподавателю, ответьте на поставленные вопросы.

12. Выполните команду clear all для очистки Рабочей среды.

13. Оформите отчет по выполненной работе.

3. Варианты индивидуальных заданий

Таблица 1.3.6-1

4. Содержание отчета

1. В форме комментариев:

• Название лабораторной работы

• ФИО студента, номер группы

• № варианта

• Индивидуальное задание

2. Протокол вычислений (сессии) в окне Command Window, снабженный необходимыми комментариями.

1.3.7. Контрольные вопросы по теме

1. Что такое script-файл и каковы его особенности?

2. Каким образом script-файл запускается на выполнение?

3. Что такое m-функция?

4. В чем отличие script-файла от m-функции?

5. Может ли m-функция иметь несколько выходных параметров?

6. Обращение к m-функции.

7. Формат оператора input().

8. Как с использованием оператора if…end реализовать стандартное, усеченное и вложенное разветвление?

9. Формат оператора множественного разветвления switch.

10. Формат оператора регулярного цикла for…end, особенности задания значений переменной цикла.

11. Назначение операторов continue и brek.

12. Оператор итеративного цикла while…end и его структура.

Раздел 2. Технология решения
вычислительных задач средствами MatLab

Тема 2.1. Решение нелинейных уравнений

2.1.1. Численное решение нелинейных уравнений

2.1.2. Лабораторная работа по теме

2.1.3. Контрольные вопросы по теме

2.1.1. Численное решение нелинейных уравнений

В общем виде уравнение можно представить, как f(x) = 0. В зависимости от вида функции f(x) различают алгебраические и трансцендентные уравнения.

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

Решением уравнения f(x)=0 называется совокупность корней , при которых уравнение обращается в тождество . Однако, точные значения корней могут быть найдены аналитически только для некоторых типов уравнений. Еще меньше возможностей при получении точного решения трансцендентных уравнений. Следует отметить, что задача нахождения точных значений корней не всегда корректна. Так, если коэффициенты уравнения являются приближенными числами, точность вычисленных значений корней заведомо не может превышать точности исходных данных. Эти обстоятельства заставляют рассматривать возможность отыскания корней уравнения с ограниченной точностью (приближенных корней).

З адача нахождения корня уравнения с заданной точностью  (>0) считается решенной, если найдено приближенное значение , которое отличается от точного значения корня не более чем на значение :

Для отделения корней нелинейных уравнений применяются графический и аналитический методы. Для уточнения корней с заданной степенью точности существует множество численных методов. Самыми распространенными из них являются: метод Ньютона, метод хорд, метод итераций и метод половинного деления [3].

Процесс нахождения приближенного корня нелинейного уравнения состоит из двух этапов:

1. отделение корня уравнения (локализация корня на отрезке);

2. уточнение корня с заданной точностью.

Пример 2.1.1-1. Отделить и уточнить корень уравнения .

Согласно теореме о существовании и единственности корня на отрезке, найдем отрезок, на концах которого функция имеет разные знаки, а первая производная непрерывна и знакопостоянна (рис. 2.1.1-1).

Рис. 2.1.1-1. Отделение корня нелинейного уравнения

Условия существования и единственности корня на отрезке [0;1] выполняются. Команды построения графика функции f(x) на отрезке [0;1] и график функции приведены на рис. 2.1.1-2.

Рис. 2.1.1-2. Построение графика функции f(x)

Для решения нелинейных уравнений в среде Matlab используются встроенные функции: solve(), fzero(), root().

Функция fzero() используется для нахождения вещественных корней уравнений вида f(х)=0. Алго­ритм, реализованный этой функцией, представляет собой комбинацию метода дихотомии (деления пополам), метода секущих и метода обрат­ной квадратичной интерполяции. В простейшем варианте обращения кроме указателя на функцию, корень которой необходимо найти, задается окрестность х0, с которой начинается поиск: х = fzero(f, x0).

Аргумент f может быть задан одним из следующих способов:

• формулой с неизвестным х, заключенной в апострофы;

• именем m-файла (в апострофах и без расширения m);

• указателем на функцию (например, @f_name).

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

Аргумент х0 может быть задан одним из двух способов:

• вектором [a;b], на концах которого функция f(x) меняет знак, что гарантирует нахождение, по крайней мере, одного корня на этом интервале;

• скалярным значением, в окрестности которого предполагается нахождение корня (в этом случае функция fzero() сама пытается найти отрезок с центром в заданной точке х0, на концах которого функция f(x) меняет знак).

Из графика (рис. 2.1.1-2) следует, что корень находится на интервале [0;1]. Используем полученную информацию и обратимся к функ­ции Matlab fzero(). При этом если мы хотим получить не только значение корня, но и узнать значение функции в найденной точке, то к функции Matlab fzero() можно обратиться с двумя выходными параметрами. На рис. 2.1.1-3 рассмотрено обращение к функции fzero() как с одним, так и с двумя выходными параметрами.

Рис. 2.1.1-3. Уточнение корня уравнения с использованием функции fzero()

Вместо явного задания выражения для функции f(x) можно создать соответствующую функцию, запомнив ее в виде m-файла (рис. 2.1.1-4).

Рис. 2.1.1-4. Описание левой части в виде m-файла

Для символьного (аналитического) решения уравнений в MatLab используется функция solve(), которая может иметь следующие форматы:solve('f(x)',x),solve('f(x)'), где 'f(x)'– левая часть решаемого уравнения, заключенная в апострофах; x – искомая символьная неизвестная (описанная как syms x).

Рассмотрим технологию определения корня с помощью функции solve() на следующем примере (рис. 2.1.1-5).

Пример 2.1.1-2. Решить уравнение 2^x–4∙x+3=0 аналитически.

Рис. 2.1.1-5. Уточнение корня уравнения с использованием функции solve()

Функция solve() в ряде случаев позволяет определить все корни уравнения f(x)=0. Она не требует информации о начальном значении корня или области его изоляции. Поэтому в случае решения ряда нелинейных уравнений, она не находит всех корней уравнения.

Функция Matlab roots() используется для вычисления корней полинома вида . Перед ее использованием создается список коэффициентов (даже нулевых), а затем вводится сама функция roots(), у которой в качестве аргумента указывается имя списка коэффициентов (рис. 2.1.1-6).

Рис. 2.1.1-6. Уточнение корня уравнения с использованием функции roots()

Вычисление корней полинома с использованием команды roots(а) показало, что у полинома есть четыре действительных корня: х₁=-1, x₂=1, x₃=2, x₄=8 и два мнимых корня x₅=i и x₆=-i.

2.1.2. Лабораторная работа по теме