Мет.указания по лабораторным работам. Лабораторная работа 1. 3 Основы теории погрешностей 3 Лабораторная работа 2. 11 Основы теории погрешностей 11
Скачать 0.72 Mb.
|
Лабораторная работа №12.«Аппроксимация, построение аппроксимирующих кривых»Цель: научиться подбирать аппроксимирующий полином первой и второй степени и строить соответствующие аппроксимирующие кривые. Задание: Подобрать аппроксимирующий полином первой и второй степени, используя метод наименьших квадратов. Используя среду MatLab подобрать аппроксимирующий полином первой, второй и третьей (если необходимо, то четвертой и т.д.) степени и построить соответствующие аппроксимирующие кривые. Найти значение функции в точке х. Варианты заданий:
Теоретический материал: Аппроксимация функций с помощью метода наименьших квадратов Метод наименьших квадратов применяется при обработке результатов эксперимента для аппроксимации (приближения) экспериментальных данных аналитической формулой. Конкретный вид формулы выбирается, как правило, из физических соображений. Такими формулами могут быть: Y=ax+b, Y=ax2+bx+c, Y=anxn+an-1xn-1+….+a1x+a0 Y=aedx+c, y=a/x+b и другие. Пример решения в MatLab: Введем исходные данные >> x=[1 2.2 2.4 2.7 3.1 3.5 4.5 5]; >> y=[9.054 15.077 15.754 18.3 17.984 15.852 1.772 -13.042]; Вычислим вектор коэффициентов полинома Y=ax+b >> a=polyfit(x,y,1) a = -5.4276 26.6479 Вычислим значения полинома на интервале от 10 до 20 >> x1=0.9:0.5:5.1; >> y1=polyval(a,x1); Построим график полинома и экспериментальных точек в одной графической области >> plot(x1,y1,'-k',x,y,'ok') По графику определяем что данная кривая неподходит (см. Рисунок 5). Поэтому далее вычисляем вектор коэффициентов полинома Y=ax2+bx+c >> a=polyfit(x,y,2) a = -4.7042 23.6773 -11.5522 Вычислим значения полинома на интервале от 10 до 20 >> x1=0.9:0.5:5.1; >> y1=polyval(a,x1); Построим график полинома и экспериментальных точек в одной графической области >> plot(x1,y1,'-k',x,y,'ok') Рисунок 5. График функции Имеем график (см. Рисунок 6) Найдем еще вектор коэффициентов полинома Y=ax3+bx2+cх+d >> a=polyfit(x,y,3) a = -1.3534 7.2257 -6.9270 10.0393 Вычислим значения полинома на интервале от 10 до 20 >> x1=0.9:0.5:5.1; >> y1=polyval(a,x1); Рисунок 6. График функции Построим график полинома и экспериментальных точек в одной графической области >> plot(x1,y1,'-k',x,y,'ok') Из графика (см. Рисунок 7) можно сделать вывод, что полином вида Y=ax3+bx2+cх+d является аппроксимурующей кривой для исходных данных Рисунок 7. График функции Содержание отчета: 1. Титульный лист. 2. Цель лабораторной работы. 3. Исходные данные, указываемые в задании и необходимые для достижения поставленной цели. 4. Расчетная часть: описание выполнения задания. 5. Выводы и анализ полученных результатов. Контрольные вопросы: В каком случае возникает необходимость постановки задачи аппроксимации? Сформулируйте задачу аппроксимации функции? Какая кривая называется линией регрессии? В чем заключается задача аппроксимации? Для чего используется аппроксимация функции на практике? В каком виде Вы искали линию регрессии? Для чего используется метод наименьших квадратов? Что называется аппроксимирующей кривой? |