МегаПредмет

ПОЗНАВАТЕЛЬНОЕ

Сила воли ведет к действию, а позитивные действия формируют позитивное отношение


Как определить диапазон голоса - ваш вокал


Игровые автоматы с быстрым выводом


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


Целительная привычка


Как самому избавиться от обидчивости


Противоречивые взгляды на качества, присущие мужчинам


Тренинг уверенности в себе


Вкуснейший "Салат из свеклы с чесноком"


Натюрморт и его изобразительные возможности


Применение, как принимать мумие? Мумие для волос, лица, при переломах, при кровотечении и т.д.


Как научиться брать на себя ответственность


Зачем нужны границы в отношениях с детьми?


Световозвращающие элементы на детской одежде


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


Как слышать голос Бога


Классификация ожирения по ИМТ (ВОЗ)


Глава 3. Завет мужчины с женщиной


Оси и плоскости тела человека


Оси и плоскости тела человека - Тело человека состоит из определенных топографических частей и участков, в которых расположены органы, мышцы, сосуды, нервы и т.д.


Отёска стен и прирубка косяков Отёска стен и прирубка косяков - Когда на доме не достаёт окон и дверей, красивое высокое крыльцо ещё только в воображении, приходится подниматься с улицы в дом по трапу.


Дифференциальные уравнения второго порядка (модель рынка с прогнозируемыми ценами) Дифференциальные уравнения второго порядка (модель рынка с прогнозируемыми ценами) - В простых моделях рынка спрос и предложение обычно полагают зависящими только от текущей цены на товар.

Нахождение линейной и квадратичной аппроксимирующих функций





Курсовая работа по предмету

Информатика

Тема: Численные методы и оптимизация

Выполнил: Студент 2 курса

Специальность 210406

Студ. билет №7СС10014

Цветков А. В.

Преподаватель: Саладаев Е. Н.

Нижний Новгород

Год


 

Задание

1. Получить линейную и квадратичную аппроксимирующие функции заданной функции y(x) методом наименьших квадратов для степенного базиса:

- выполнить расчет параметров линейной аппроксимирующей функции;

- получить систему нормальных уравнений и выполнить расчет ее параметров для проведения квадратичной аппроксимирующей функции;

- проверить полученные результаты;

- вычислить значения аппроксимирующих функций в узлах аппроксимации;

- построить график заданной функции (множество заданных точек) и графики функций линейной и квадратичной аппроксимации;

- оценить качество аппроксимации.

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

- отделить корни уравнения;

- проверить (аналитически) условия сходимости применяемых методов решения уравнений. В случае необходимости привести уравнение к виду, обеспечивающему сходимость процесса приближения к корню;

- выбрать начальные приближения;

- записать рекуррентную формулу для уточнения корня и произвести по ней расчеты;

- оценить погрешности.

3. Вычислить при разбиении отрезка интегрирования на n1=10 и n2=20 интервалов, где x1, x2 - корни уравнения :

- оценить погрешность.

4. Определить точку экстремума функции методами одномерной оптимизации (с точностью ):

- проверить условие унимодальности функции и выбрать начальный отрезок оптимизации;

- провести расчеты по сокращению отрезка оптимизации и проверить условие окончания поиска минимума (максимума) функции.


 

Нахождение линейной и квадратичной аппроксимирующих функций

 

Исходная функция задана таблицей.

i
xi
y(x) 0.476 1.496 1.15 -0.196 -1.176 -1.376

 

Линейная аппроксимирующая функция имеет вид

Найдем параметры линейной аппроксимации:

 

Проверим полученные данные средствами MathCad:

или

Искомая линейная аппроксимирующая функция примет вид:

 

Квадратичная аппроксимирующая функция имеет вид . Для определения параметров многочлена второй степени составим и решим систему нормальных уравнений.

Запишем в следующую таблицу элементы матрицы Грама и столбец свободных членов:

i x0 x1 x2 x3 x4 y xy x2y
0.467 5.137 56.507
1.496 19.448 252.824
1.15 17.25 258.75
-0.196 -3.332 -56.644
-1.176 -22.344 -424.536
-1.376 -28.896 -606.816
0.365 -12.737 -519.915

 

Система нормальных уравнений примет вид:

 

 

Найдем решение этой системы уравнений:

Проверим полученные данные средствами MathCad:

Искомая квадратичная аппроксимирующая функция примет вид:

Значения исходной функции y(x), а также аппроксимирующих функций Flin(x) и Fsq(x) в узлах аппроксимации приведены в таблице 3:



x
y(x) 0.467 1.496 1.15 -0.196 -1.176 -1.376
Flin(x) 1.392 0.862 0.332 -0.198 -0.728 -1.258
Fsq(x) 0.954 1.078 0.898 0.414 -0.374 -1.466

 

Графики функций линейной и квадратичной аппроксимации показаны на рисунке

Оценим качество аппроксимации:

Для линейной функции: ρlin=0.654

Для квадратичной функции: ρsq=0.547

Получили, что ρ21, из чего можно сделать вывод, что квадратичная аппроксимация более качественная.

 

2. Решение уравнения Fsq(x)=0 с точностью E = 10-3.

 

Решим уравнение .

Найдем производные функции Fsq(x):

Для отделения корней уравнения составим таблицу знаков функции Fsq(x) и ее производной Fsq’(x):

x
Fsq(x) -0.206 0.526 0.414 -0.374
Sign Fsq(x) - + + -
Fsq’(x) 0.442 0.29 -0.318 -0.47
Sign Fsq’(x) + + - -

 

На отрезках [7;9] и [17;19] функция Fsq(x) меняет знак, т.е. существует, по крайней мере, по одному корню, принадлежащему каждому из этих отрезков. Поскольку знак первой (и второй) производной на выбранных отрезках остается постоянным, то можно сказать, что функция на этих отрезках монотонна. Следовательно, выбранные отрезки содержат по одному корню.

 

1) Уточним первый корень уравнения Fsq(x)=0 на отрезке [7;9] методом Ньютона.

Проверим выполнение достаточных условий сходимости метода Ньютона:

Fsq(x) непрерывна на отрезке [7;9]

Fsq(7)∙Fsq(9)<0

Fsq’(x) и Fsq”(x) отличны от нуля и сохраняют знаки при .

Выберем начальное приближение x0, удовлетворяющее условию

Возьмем x0=8, поскольку Fsq”(x)=-0.076<0 и Fsq(8)<0.

Рекуррентная формула для уточнения корня уравнения по методу Ньютона:

Уточним корень уравнения на отрезке [7;9].

Условие окончания поиска корня:

Результат запишем в таблицу:

i xi Fsq(xi) Fsq’(xi) xi-1 xi-xi-1
0.198 0.366    
7.45902 -0.01112 0.40711 -0.54098
7.48633 -2.83566∙10-5 0.40504 7.45902 0.02732
7.4864 -1.8625∙10-10 0.40503 7.48633 7.00095∙10-5

 

Оценим погрешность вычисления 1-го корня после 3 итераций, используя формулу , где , .

Получим .

Проверим найденный корень уравнения средствами MathCad:

Таким образом, x1=7.4864.

 

2) Уточним второй корень уравнения Fsq(x)=0 на отрезке [17;19] методом итераций.

Fsq(x) дифференцируема и имеет разные знаки на отрезке [17;19].

Найдем итерирующую функцию. Итерирующая функция обеспечивает выполнение условия сходимости для . Выберем параметр λ по правилу:

В нашем случае r=0.47 и Fsq’(x)>0, следовательно, . Пусть λ=-2.

Запишем рекуррентную формулу для вычисления приближений к корню:

Выберем начальное приближение к корню, например, x0=2.

Условие окончания поиска корня:

Результат запишем в таблицу:

i xi Fsq(xi) xi-1 xi-xi-1
0.058    
17.826 0.12541 -0.174
17.44978 0.26328 17.826 -0.37622
16.65994 0.51774 17.44978 -0.78984
15.10671 0.87985 16.65994 -1.55323

 

Оценим погрешность вычисления 2-го корня после 4 итераций, используя формулу , где .

 

Получим проверим найденный корень уравнения средствами MathCad:

 

Таким образом, x2=18.14518.

 

3. Вычисление определенного интеграла методами Симпсона, трапеций и средних прямоугольников, полагая n1=10 и n2=20.

 

Вычислим интеграл

 

 

Метод средних прямоугольников:

Метод трапеций:

 

Метод Симпсона:

Результаты занесем в таблицу:

n Метод средних прямоугольников Метод трапеций Метод Симпсона
Ip = 7.416 It = 7.387 Is = 7.786
Ip = 7.614 It = 7.618 Is = 7.708

 

Оценка погрешности по правилу Рунге:

 

Для методов средних прямоугольников и трапеций при k=2

Rср.пр. = 0.065

Rтрап. = 0.077

Для метода Симпсона при k=4

RСимп. = 5.167∙10-3

 

Вычислим значение средствами MathCad:

 

 

 

 

4. Определение точки экстремума функции Fsq(x) методами одномерной оптимизации с точностью E=10-2.

 

Для нахождения точки экстремума функции применим методы дихотомии и золотого сечения.

Выясним, унимодальна ли функция Fsq(x), и выберем отрезок неопределенности. Рассмотрим отрезок [12.5;13.5].

 

Получили, что функция Fsq(x) дважды дифференцируема на рассматриваемом отрезке, ее первая производная Fsq’(x) не убывает на этом отрезке, а вторая производная Fsq’’(x)>0. Таким образом, на отрезке [12.5;13.5] функция Fsq(x) унимодальна, и, следовательно, этот отрезок может быть выбран в качестве начального отрезка неопределенности. Найдем точку экстремума методами дихотомии и золотого сечения.

Условия прекращения поиска:

 

Метод дихотомии:

i a b x1 x2 Fsq(x1) Fsq(x2) b-a
12.5 13.5 12.998 13.002 13.91808 13.92592
12.998 13.5 13.247 13.251 14.40893 14.41686 0.502
13.247 13.5 13.3715 13.3755 14.65613 14.66409 0.253
13.247 13.3755 13.30925 13.31325 14.53238 14.54033 0.1285
13.247 13.31325 13.27812 13.28212 14.47062 14.47855 0.06625
13.247 13.28212 13.26256 13.26656 14.43977 14.4477 0.03513
13.247 13.26656 13.25478 13.25878 14.42435 14.43227 0.01956
13.25478 13.26656 13.25867 13.26267 14.43206 14.43998 0.01178
13.25478 13.26267 13.25673 13.26073 14.4282 14.43613 7.89063∙10-3

 

После 8 итераций получили xmin=13.25873 и Fsq(xmin)=14.43217.

 

Метод золотого сечения:

 

i a b x1 x2 Fsq(x1) Fsq(x2) b-a
12.5 13.5 12.882 13.118 13.69101 14.15405
12.882 13.5 13.118 13.26392 14.15405 14.44247 0.618
13.118 13.5 13.26392 13.35408 14.44247 14.62146 0.382
13.118 13.35408 13.20818 13.26392 14.3321 14.44247 0.23608
13.20818 13.35408 13.26392 13.29834 14.44247 14.51073 0.14589
13.20818 13.29834 13.24262 13.26392 14.40026 14.44247 0.09016
13.24262 13.29834 13.26392 13.27706 14.44247 14.46851 0.05572
13.24262 13.27706 13.25578 13.26392 14.42632 14.44247 0.03444
13.25578 13.27706 13.26392 13.26893 14.44247 14.45239 0.02128
13.25578 13.26893 13.2608 13.26392 14.43628 14.44247 0.01315
13.25578 13.26392 13.25889 14.43249 14.43249 14.43628 8.14633∙10-3

 

 

 

После 10 итераций получили xmin=13.25985 и Fsq(xmin)=14.43439.

 

Выполним задачу оптимизации средствами MathCad:

 


 

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

1. Конспект лекций по информатике

2. Методические указания и контрольные задания по дисциплине Информатика. Численные методы и оптимизация. Г.К.Сосновиков, И.Б.Юскова, Н.П.Муравьев. 2008 г.

3. Практическое решение инженерных и научных задач на ПК с использованием математических пакетов. В.Н.Шакин, Т.И.Семенова, О.М.Кравченко. 2010 г.





©2015 www.megapredmet.ru Все права принадлежат авторам размещенных материалов.