Методика выполнения задания №3 «Создание М-файла»
Для решения системы уравнений создадим два М-файла – файл-функцию для решения системы уравнений методом Крамера и файл-сценарий, где будут определяться исходные матрицы и вычисляться вектор корней системы.
3.1.Создайте файл-функцию. Для этого М-файла на панели инструментов нажмите Новый М-файл (New M-file), либо выполните команды File – New—M-file.
Открывается окно М-файла. Сюда введите операторы, соответствующие решению системы уравнений методом Крамера. Назовите функцию Kramer. Входными аргументами являются матрица системы А и вектор свободных членов b, выходными – вектор Х. Название функции, набор входных и выходных параметров определяются первым оператором файла-функции, имеющим следующий вид: function X = Kramer (A, b)
Далее рассчитывается определитель матрицы системы А, вводятся новые матрицы А1, А2, А3 и А4, где соответственно 1-й, 2-й, 3-й и 4-й столбцы заменяются столбцом свободных членов b, рассчитываются определители новых матриц и элементы вектора корней системы Х. Файл-функция имеет вид, представленный на рисунке 1.3.
Рисунок 1.3 – Файл-функция решения системы уравнений методом Крамера
Вектор Х представляет из себя вектор-столбец. На это указывает нумерация элементов вектора, где первая цифра номер строки, вторая – номер столбца. Если пронумеровать элементы вектора Х, как Х(1), Х(2), Х(3) и Х(4), то вектор Х будет представлять из себя вектор-строку.
3.2. Сохраните файл-функцию под именем Kramer в папке work каталога MATLAB. Важно – имя файла должно совпадать с именем функции (Kramer).Для этого нажмите в строке меню Файл (File) – Сохранить как (Save As), либо на Панели инструментов нажмите Сохранить. Появляется окно Сохранить как ((Save As)с содержимым папки work (рисунок 1.4). Нажать Сохранить (Save).
Рисунок 1.4 – Сохранение М-файла
3.3.Создайте файл-сценарий, где будут указаны исходные матрицы, произведен вызов созданного ранее файла-функции Kramer и решена система уравнений матричным методом. Блок-схема алгоритма решения задачи представлена на рисунке 1.5. Предопределенный процесс «Решение системы уравнений методом Крамера» представляет из себя файл-функцию Kramer, введенную ранее.
Аналогично пункту 3.1 создайте новый М-файл. Запишите в него соответствующий приведенному алгоритму следующий программный код:
% Решение системы уравнений
A=[1 1 1 0; 1 -1 3 1; -1 2 0 3; 3 1 1 -1]
b=[3; 7; -2; 7]
% Решение методом Крамера
Xk=Kramer(A,b)
% Решение матричным методом
Xm=inv(A)*b
3.4. Запустите на выполнение файл-сценарий. Для этого в строке меню выберите Отладка (Debug) – Сохранить и запустить (Save and Run).Следует помнить, что на решение запускается именно файл-сценарий, а не файл-функция. Происходит выполнение программы. Результаты отображаются в основном окне MATLAB.
Рисунок 1.5 Блок-схема алгоритма
решения задачи
В Окне команд отображаются результаты выполнения каждой строки программы, где в конце строки не стоит знак «точка с запятой» (;). Это удобно при отладке программ. Если вы не хотите вывода промежуточных результатов, то в конце соответствующей строки программы поставьте знак «точка с запятой». Ошибки отображаются в Окне команд красным цветом.
3.5. Проанализируйте результаты расчетов. В Рабочей области (Workspace) окна MATLAB отображаются переменные и массивы, введенные в программу и рассчитанные в ней. В данном случае там будут отображаться матрица А, вектор-столбец b, векторы-столбцы результатов решения системы, методом Крамера – Xk, матричным методом Xm. Щелкнув дважды левой кнопкой мыши по любому массиву или переменной из Рабочей области, можно раскрыть их значения в специальном окне Array Editor: (рисунок 1.6).
Рисунок 1.6 – Результаты решения задачи
3.6.Завершите работу MATLAB/