Как было отмечено ранее, файл-функция (процедура) должна начинаться со строки заголовка
function [<ПКВ>] = <имя процедуры>(<ПВВ>).
Если перечень конечных (выходных) величин (ПКВ) содержит только один объект (в общем случае - матрицу), то файл-функция представляет собой обычную функцию (одной или нескольких переменных). Фактически даже в этом простейшем случае файл-функция является уже процедурой в обычном смысле других языков программирования, если выходная величина является вектором или матрицей. Первая строка в этом случае имеет вид:
function <имя переменной> = <имя процедуры>(<ПВВ>).
Если же в результате выполнения файл-функции должны быть определены (вычислены) несколько объектов (матриц), такая файл-функция представляет собой уже более сложный объект, который в программировании обычно называется процедурой, или подпрограммой. Общий вид первой строки в этом случае становится таким:
function [y1, y2,..., y] = <имя процедуры>(<ПВВ>), т. е. перечень выходных величин y1, y2,..., y должен быть представлен как вектор-строка с элементами y1, y2,..., y (все они могут быть матрицами).
В простейшем случае функции одной переменной заголовок приобретет вид:
function y = func (x),
где func - имя функции (М-файла).
В качестве примера рассмотрим составление М-файла для функции
Для этого следует активизировать меню File командного окна MatLAB и выбрать в нем сначала команду New, а затем команду M-file. На экране появится окно текстового редактора. В нем нужно набрать такой текст:
function y = F1(x,d)
% Процедура, которая вычисляет значение функции
% y = (d3)*ctg(x)*sqrt(sin(x)4-cos(x)4).
% Обращение y = F1(x,d).
y = (d^3)*cot(x). *sqrt(sin(x). ^4-cos(x). ^4);
После этого необходимо сохранить этот текст в файле под именем F1.m.
Необходимый М-файл создан. Теперь можно пользоваться этой функцией при расчетах. Так, если ввести команду
» y = F1(1, 0.1)
то получим результат
y = 4. 1421e-004.
Следует заметить, что аналогично можно получить сразу вектор всех значений указанной функции при разных значениях аргумента, если последние собрать в некоторый вектор. Так, если сформировать вектор
» zet= 0:0. 3:1. 8;
и обратиться в ту же процедуру
» my = F1(zet,1),
то получим:
Warning: Divide by zero
my =
Columns 1 through 4
Na + Infi 0 + 2. 9369i 0 + 0. 8799i 0. 3783
Columns 5 through 7
0. 3339 0. 0706 -0. 2209
Примечания.
1. Возможность использования сформированной процедуры как для отдельных чисел, так и для векторов и матриц обусловлена применением в записи соответствующего М-файла вместо обычных знаков арифметических действий их аналогов с предшествующей точкой.
2. Во избежание вывода на экран нежелательных промежуточных результатов, необходимо в тексте процедуры все вычислительные операторы завершать символом "; ".
3. Как показывают приведенные примеры, имена переменных, указанные в заголовке файл-функции могут быть любыми (совпадать или нет с именами, используемыми при обращении к этой файл-функции), т. е. носят формальный характер. Важно, чтобы структура обращения полностью соответствовала структуре заголовка в записи текста М-файла и чтобы переменные в этом обращении имели тот же тип и размер, как и в заголовке М-файла.
Чтобы получить информацию о созданной процедуре, достаточно набрать в командном окне команду:
» help f1,
и в командном окне появится
