- •1.1 Цель работы
- •1.2 Постановка задачи
- •1.4 Порядок выполнения лабораторной работы
- •1.5 Программа аппроксимации экспериментальных данных
- •1.6 Содержание отчета
- •1.7 Контрольные вопросы
- •2.4 Программа для интерполяции кубическими сплайнами
- •2.6 Содержание отчета
- •2.7 Варианты
- •2.8 Контрольные вопросы
- •3.2.3 Проверка адекватности модели
- •4.1 Цель работы
- •4.3.1 Метод сканирования
- •4.3.1.1 Сущность метода
- •4.3.1.2 Программа поиска методом сканирования
- •4.3.2 Метод половинного деления (метод дихотомии)
- •4.3.2.1. Сущность метода
- •4.3.2.2. Программа поиска методом половинного деления
- •4.3.3 Метод золотого сечения
- •4.3.3.1 Сущность метода
- •4.3.3.2 Программа поиска методом золотого сечения
- •4.5 Порядок выполнения работы
- •4.6 Содержание отчета
- •4.7 Контрольные вопросы
- •5.1 Цель работы
- •5.2 Теоретические основы
- •5.3 Постановка задачи
- •5.4 Порядок выполнения работы
- •5.5 Содержание отчета
- •5.6 Варианты
- •5.7 Контрольные вопросы
- •6.1 Цель работы
- •6.2 Теоретические основы
- •6.2.1 Математические основы решения задачи (при К=2)
- •6.2.2 Программа поиска симплексным методом
- •6.3 Постановка задачи
- •6.4 Порядок выполнения работы
- •6.5 Содержание отчета
- •6.6 Варианты
- •6.7 Контрольные вопросы
8
4)Нажмите левой кнопкой мышки на кнопку «Параметры». Откроется диалоговое окно «Параметры поиска решения». Установите опции, как указано на рисунке 1.3. Закройте окно параметров нажав на кнопку
«ОК».
5)Для начала поиска решения нажмите на кнопку «Выполнить».
6)Обработка экспериментальных данных производится всеми десятью уравнениями, приведенными в таблице 1.2, и в качестве окончательного уравнения для описания выбирается уравнение, для которого получалось наименьшее значение Φ.
7)Качество аппроксимации визуально можно наблюдать на графике, приведенном на листе «Graf».
8)После окончательного выбора модели, запишите данные обработки экспериментальных данных.
9)После выполнения лабораторной работы необходимо составить отчет о проделанной работе.
1.5Программа аппроксимации экспериментальных данных
dim Y 200 ,X 200 ,S 10,200 ,A 50
'В правойчасти уравнения 2 напишите аппроксимирующее
'выражение, например y a * x^b * exp c*x . |
1 |
||
deffnz a,b,c,x a*x^b*exp c*x |
|
' 2 |
|
'Введите число экспериментальныхданных |
|||
N1 |
13 |
|
|
'Ввод исходных данных |
|
|
|
For i 1 To N1 |
4, 3 |
.Value |
|
x i Worksheets "Lab1" .Cells i |
|||
y i Worksheets "Lab1" .Cells i |
4, 4 |
.Value |
|
Next i |
|
|
'Введите число неизвестных параметров в аппроксимирующем урав нении
N2 2
N2 N2 1
'Зануление коэффициентов алгебраического уравнения
For I 0 To N2: C I 0:For J 0 To N2: S J,I 0:Next J:Next I 'Формирование коэффициентов нормальных уравнений
For I 1 To N1
'Ввод исходных уравнений
'Приведите правую часть 1 |
к линейному виду относительно параметров |
||
а,в,с. |
lna |
b*lnx |
c*x, тогда для этого уравнения |
Например, lny |
|||
f 0 1;f 1 |
lnx;f 2 |
x; fc |
lny. |
'Коэффициентыаппроксимирующегоуравнения
'A 0 *f 0 |
A 1 *f 1 |
... A i *f i |
9 |
|
... A m *f m fc. |
||||
|
'Ниже вводите выражениядля f i и fc |
|||
|
f 0 |
1 |
|
|
|
f 1 |
log x i |
|
|
|
f 2 |
x i |
|
|
|
fc |
log y i |
|
|
|
'Формирование коэффициентов нормальныхуравнений |
|||
For V 0 To N2:For J |
0 To N2: S V,J |
F J *F V S V,J :Next J |
||
C V |
C V |
F V *Fc: Next V:Next I |
|
|
|
'Решение системы уравнений методом Гаусса |
|||
For J 0 To N2 1 |
|
|
||
R S J,J |
|
|
|
|
ForI JToN2 |
|
|
||
S J,I |
S J,I /R |
|
|
|
Next I |
|
|
|
|
C J C J /R |
|
|
||
For I |
J 1 To N2 |
|
|
|
R S I,J |
|
|
|
|
For V J To N2 |
|
|
||
S I,V |
S I,V R*S J,V |
|
|
Next V
C I C I C J *R Next I
Next J
For J N2 To 0 Step 1 A J C J /S J,J
for I J 1 to 0 step 1 C I C I S I,J *A J next I
next J
'Определение параметров уравнения 'Введите формулы для расчета параметров уравнения 1 через коэффици енты A I .
aEXP A 0
b A 1
cA 2
sym 0
For I 1 To N1 zz fnz a, b, c, x i G ZZ Y I
sym sym G^2
Выводданныхнапечать:x i ;zz; y i
Next I
10
Выводзначениядисперсииаппроксимациинапечать:sym/N1;a, b, c
End
1.6 Содержание отчета
1Теоретические основы метода наименьших квадратов.
2Экспериментальные данные, которые обрабатывались.
3Значения функции Φ, соответствующие приведенным в таблице
1.2уравнениям.
4Уравнение, которое подобрано для аппроксимации экспериментальных данных.
5Параметры уравнения, определенные МНК.
1.7Контрольные вопросы
1Объясните сущность метода наименьших квадратов.
2Как находится минимум функции невязки f(x,a) в методе наименьших квадратов?
3Какая система уравнений решается для нахождения коэффициентов аппроксимирующего уравнения?
4Каким методом решается система алгебраических уравнений?
5Каким свойством должно обладать аппроксимирующее уравнение относительно своих коэффициентов?
6В чем сущность решения системы алгебраических уравнений методом Гаусса?
7Что означает «идентификация математической модели»?
8Что означает « параметрическая идентификация математической модели»?
9Приведите уравнение y/x2 =a+bx+cln(y) к линейному виду относительно своих коэффициентов.
11
ЛАБОРАТОРНАЯ РАБОТА № 2 ОБРАБОТКА ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ
КУБИЧЕСКИМИ СПЛАЙНАМИ
2.1 Цель работы
Освоение метода и алгоритма интерполяции экспериментальных данных кубическими сплайнами.
2.2 Область применения сплайнов
Сплайны могут использоваться для сглаживания экспериментальных данных, искаженных влиянием помех. Затем сглаженные данные используются для построения модели объекта.
Кроме того, сплайны могут использоваться при численном интегрировании функций, что позволяет более точно найти значение интеграла. Интеграл от сплайна легко вычисляется, так как сплайн является кусочнополиномиальной функцией.
Эффективно также использование сплайнов при численном дифференцировании функции. В этом случае по известным значениям функции вначале строится интерполяционный или сглаживающий сплайн, а затем находится производная от него.
2.3Интерполяция функции одного переменного кубическим сплайном
Постановка задачи. На отрезке [a, b] оси x выделен ряд точек
x0 , x1 ,...,xi ,..., xn , причем x0 =a, xn =b, x0 <x1 ...<xi <...xn . В каждой точ-
ке xi, которая называется узлом сетки, задано значение fi функции f(x).
Требуется найти функцию |
, такую, что: |
|
|
|
|
|||
1) |
|
и ее производные первого и второго порядков непрерывны на |
||||||
отрезке [a, b], и в том числе в узлах сетки, т.е. |
0 |
" |
0 ; |
|
||||
|
0 |
0 , |
0 |
0 , " |
(2.1) |
|||
2) |
на каждом отрезке [xi-1, xi] |
является кубическим многочле- |
||||||
ном вида |
|
|
|
, |
|
1,2,..., ; |
|
|
3) |
в узлах сетки xi |
= fi , i=0,1,2,...,n; |
|
(2.2) |
||||
|
|
|
|
|||||
4) |
на границах отрезка вторые производные функции равны нулю: |
|||||||
|
|
" |
" |
0 |
|
|
|
(2.3) |
Доказано, что эта задача имеет единственное решение, и это решение определяет наиболее гладкую функцию, проходящую через заданные точки.