Программирование в MATLAB

2



Кафедра: Информационные технологии





















ПРОГРАММИРОВАНИЕ В MATLAB






Операторы MATLAB



Операторы цикла



Цикл for



Синтаксис

for count=start:step:final

команды MATLAB

end



Описание

count – переменная цикла,

start – ее начальное значение,

final – ее конечное значение,

step – шаг, на который увеличивается count при каждом следующем заходе в цикл

цикл заканчивается, как только значение count становится больше final.



Пример

Пусть требуется вывести семейство кривых для x€ [0, 2π], которое задано функцией, зависящей от параметра

y (x, a) = e-ax sin x,

для значений параметра а от -0.1 до 0.1. Ниже приведен листинг файл-программы для вывода семейства кривых.



Листинг программы

Figure

x = [0:pi/30:2*pi];

for a = -0.1:0.02:0.1

y = exp (-a*x).*sin(x);

hold on

plot (x, y)

end



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



Цикл while



Синтаксис

while условие цикла

команды MATLAB

end



Описание

Цикл работает, пока выполняется (истинно) условие цикла. Для задания условия выполнения цикла допустимы следующие операции отношения:



Обозначение

Операция отношения

==

Равенство

<

Меньше

<=

Меньше или равно

>=

Больше или равно

>

Больше

~=

Не равно



Задание более сложных условий производится с применением логических операторов. Логические операторы приведены в следующей таблице




Оператор

Условие

Эквивалентная запись

Логическое «и»

x<3 и k=4

(x<3) & (k==4)

Логическое «или»

x=1 или x=2

(x==1) | (x==2)

Отрицание «не»

a≠1.9

~(a==1.9)



Пример



Операторы ветвления



Условный оператор if



Синтаксис

if условие

команды MATLAB

end



Описание



Если условие верно, то выполняются команды MATLAB, размещенные между if и end, а если условие не верно, то происходит переход к командам, расположенных после end.



Пример



Условный оператор elseif



Синтаксис

if условие1

команды1

elseif условие2

команды2

………………………

elseif условиеn

командыn

else

команды

end



Описание

В зависимости от выполнения того или иного условия работает соответствующая ветвь программы, если все условия неверны, то выполняются команды, размещенные после else.



Пример



Оператор switch



Синтаксис



switch переменная

case значение1

команды1

case значение2

команды2

……………………

case значениеn

командыn

otherwise

команды

end


Каждая ветвь определяется оператором case, переход в нее выполняется тогда, когда переменная оператора switch принимает значение, указанное после case, или одно из значение из списка case. После выполнения какой-либо из ветвей происходит выход из switch, при этом значения, заданные в других case, уже не проверяются. Если подходящих значений для переменной не нашлось, то выполняется ветвь программы, соответствующая otherwise.



Пример



Прерывания цикла. Исключительные ситуации.



Оператор break

Синтаксис



break



Оператор break используется при организации циклических вычислений: for…end, while…end. При выполнении условия

if условие

break

end

оператор break заканчивает цикл (for или while) и происходит выполнение операторов, которые расположены в строках, следующих за end. В случае вложенных циклов break осуществляет выход из внутреннего цикла.



Обработка исключительных ситуаций, оператор try…catch



Синтаксис

try

операторы, выполнение которых

может привести к ошибке

catch

операторы, которые следует выполнить

при возникновении ошибки в блоке

между try и catch

end



Описание

Конструкция try…catch позволяет обойти исключительные ситуации (ошибки, приводящие к окончанию работы программы, например, обращение к несуществующему файлу) и предпринять некоторые действия в случае их возникновения.



Пример



Сервисные функции



disp – осуществляет вывод текста или значения переменной в командное окно

input – осуществляет запрос на ввод с клавиатуры. Используется при создании приложений с интерфейсом из командной строки.

eval – выполняет содержимое строки или строковой переменной, как команды MATLAB

clear – удаляет переменные рабочей среды.

сlc – производит очистку командного окна



Более подробную информацию об этих и других функциях можно узнать, выполнив в командной строке

help имя_функции


Задания на выполнение лабораторной работы



Номер конкретного варианта задания определяется преподавателем.



Задание №1

Данное задание подразумевает нахождение для некоторой совокупности данных алгебраического интерполяционного многочлена степени n: Pn(x).

Цель работы:

Необходимо составить программу вычисления коэффициентов алгебраического интерполяционного многочлена Pn(x)=a0+a1x+ … +anxn.

Методические указания:

Пусть, например, имеется следующая совокупность данных:



i

0

1

2

3

хi

1,2

1,4

1,6

1,8

yi

8,3893

8,6251

8,9286

8,9703



Коэффициенты a0, a1, …, an определяются из решения системы уравнений:







Здесь n – порядок интерполяционного многочлена,

n+1 – количество заданных пар точек (x, y),

a0, a1,… an– искомые коэффициенты многочлена Pn(x)=a0+a1x+ … +anxn).

Требования к программе

При выполнении данной работы необходимо:

		Задать границы отрезка [x0, xn], на котором строится интерполяционный многочлен P(x)

		Задать n – количество отрезков интерполяции (или, что то же самое, степень многочлена)

		Примечание: x0, xn, n вводятся с клавиатуры.

	

		Для получения исходных данных (x, y) (количество пар точек (xi, yi), по которым строится интерполяционный многочлен P(x) – n1=n+1) предусмотреть:

Ввод произвольно расположенных узлов xi, i=0, n с клавиатуры

Расчет узлов xi, i=0, n, соответствующих равномерному расположению аргумента x на отрезке [x0, xn]

В пп. 1,2 значения yi, i=0, n либо вводятся с клавиатуры (если исходная функция неизвестна), либо вычисляются по заданной функции f(x). Выражение, определяющее функцию, вводится с клавиатуры и должно соответствовать правилам записи выражений в MATLAB

Ввод данных (xi,yi, i=0, n) из файла

		Решить систему уравнений для определения коэффициентов многочлена P(x)

		Построить графики исходной таблично заданной функции и многочлена P(x)

		Если исходные данные заданы в виде функции f(x), построить график погрешности интерполяции /f(x) – P(x)/. Рассчитать максимальное по модулю значение погрешности интерполяции на заданном промежутке.

При выполнении последнего пункта на отрезке [x0, xn] взять не менее 500 точек для проведения расчетов




Задание №2

Интерполяция сплайнами

Цель работы:

Необходимо составить программу вычисления коэффициентов и построения сплайн-функции S(x), «склеенной» из кусков многочленов 3го порядка Si(x), которые имеют специальную форму записи:



,



функция Si(x) определена на отрезке [xi-1, xi],









Требования к программе

При выполнении данной работы необходимо:

		Задать границы отрезка [x0, xn], на котором строится сплайн-функция S(x)

		Задать n – количество отрезков интерполяции, на каждом из которых строится кубический многочлен Si(x).

		Примечание: x0, xn, n вводятся с клавиатуры.

		Организовать ввод исходных данных (x, y) (количество пар точек (xi, yi), по которым строится сплайн-функция S(x), n1=n+1), предусмотрев:

Ввод произвольно расположенных узлов xi, i=0, n с клавиатуры

Расчет узлов xi, i=0, n, соответствующих равномерному расположению аргумента x на отрезке [x0, xn]

В пп. 1,2 значения yi, i=0, n либо вводятся с клавиатуры (если исходная функция неизвестна), либо вычисляются по заданной функции f(x). Выражение, определяющее функцию, вводится с клавиатуры и должно соответствовать правилам записи выражений в MATLAB

Ввод данных (xi, yi, i=0, n) из файла

		Для краевых условий 1–3 рассчитать коэффициенты i кубических сплайнов

	S1''(x0)=0, S3''(x3)=0

	S1'(x0)=f '(x0), S3'(x3)=f '(x3)

	S1''(x0)=f «(x0), S3''(x0)=f «(x3)

		Для определения коэффициентов естественного кубического сплайна (краевые условия 1) необходимо решить следующую систему уравнений:

		





Коэффициенты σ0 =0,σn =0

		Построить графики исходной функции и сплайн-функций для всех трех типов краевых условий.

		Построить графики функций погрешности сплайн-интерполяции f(x) – S(x) для всех трех типов краевых условий.

		Рассчитать максимальные по модулю значения погрешности интерполяции на заданном промежутке [x0, xn] для каждого вида сплайна.

Примечание:

В пакете MATLAB индексы одномерных и двумерных массивов начинаются с 1, а не с 0. Учесть это при составлении программы.


Задание №3

Аппроксимация функции методом наименьших квадратов (МНК).

Данное задание подразумевает нахождение для некоторой совокупности данных аппроксимирующей функции (многочлена степени m), построенной методом наименьших квадратов (МНК).

Цель работы:

Необходимо составить программу нахождения коэффициентов многочлена φ(x)=a0+a1*x+… an * xm методом наименьших квадратов.

Рекомендации по выполнению работы:

Пусть, например, имеется следующая совокупность данных:



хi

1,2

1,4

1,6

1,8

2,0

2,2

2,4

2,6

2,8

3,0

yi

8,3893

8,6251

8,9286

8,9703

9,1731

9,1784

8,8424

8,7145

8,3077

7,9611



Поиск необходимых коэффициентов осуществляется следующим образом:







где 	n – количество точек (x, y),

m – степень искомого многочлена,

a0, a1, …, am – искомые коэффициенты (φ(x)=a0+a1x+ … +amxm).

Требования к программе

При выполнении данной работы необходимо:

		Задать границы отрезка [x1, xn], на котором строится аппроксимирующая функция φ(x)=a0+a1*x+… an * xm

		Задать m – степень многочлена

		Примечание: x1, xn, m вводятся с клавиатуры.

		Для получения исходных данных (x, y), по которым строится аппроксимирующая функция φ(x)=a0+a1*x+… an* xm предусмотреть:

Ввод произвольно расположенных узлов xi, i=1, n с клавиатуры

Расчет узлов xi, i=1, n, соответствующих равномерному расположению аргумента x на отрезке [x1, xn]

В пп. 1,2 значения yi, i=1, n либо вводятся с клавиатуры (если исходная функция неизвестна), либо вычисляются по заданной функции f(x). Выражение, определяющее функцию, вводится с клавиатуры и должно соответствовать правилам записи выражений в MATLAB

Ввод данных (xi, yi, i=1, n) из файла

		Решить систему уравнений для определения коэффициентов многочлена φ(x)

		Построить графики исходной таблично заданной функции и многочлена φ(x)

		Если исходные данные заданы в виде функции f(x), построить график погрешности интерполяции /f(x) – φ(x)/. Рассчитать максимальное по модулю значение погрешности интерполяции на заданном промежутке.

При выполнении последнего пункта на отрезке [x1, xn] взять не менее 500 точек для проведения расчетов




Требования к оформлению лабораторной работы



Отчет должен содержать:

Постановку задачи

Текст программы

Результаты тестирования

Примечание тексты программ должны быть снабжены комментариями.






Список литературы



			Ануфриев И.Е. Самоучитель Matlab 5.3/6.x – СПб.: БХВ-Петербург, 2003. – 736 с.: ил.

			В.П. Дьяконов MATLAB 6.5 SPI/7 + Simulink 5/6 в математике и моделировании. Серия «Библиотека профессионала». – М.: СОЛОН-Пресс, 2005. – 576 с.: ил.

			Ануфриев И.Е., Смирнов А.Б., Смирнова Е.Н. 	MathLab 7. – СПб.: БХВ-Петербург, 2005. – 1104 с.: ил.