Аппроксимация функций обобщенным рядом Фурье. Методические указания к лабораторному практикуму и самостоятельной работе студентов
.pdf
Рисунок 4.7 – Результаты аппроксимации
Как видим из анализа графиков эмпирической и аппроксимированной плотностей распределения вероятностей, представленных на рисунке 4.7, расхождения между линиями графиков минимальны, а среднеквадратическая ошибка составила 0.01, что является приемлемым для практических целей результатом и говорит о достаточно высоком качестве проведенной аппроксимации.
4) Сохранение результатов. На заключительном шаге этапа аппроксимации пользователю предлагается:
–сохранить полученные результаты в файл, либо
–сохранить их в базе данных для последующего анализа и сравнения с результатами, полученными другими методами аппроксимации или для других функций.
При сохранении результатов в файл пользователь самостоятельно назначает имя файла и место его сохранения. Программой предусмотрен следующий формат хранения результатов аппроксимации (рисунок 4.8):
21
Рисунок 4.8 – Формат хранения результатов аппроксимации
В первой строке результирующего файла указывается тип ортогональных полиномов использовавшихся для аппроксимации. Во второй и третьей строках указываются минимальное и максимальное значения исходной выборки соответственно. В четвертой строке указывается оптимальное количество использованных для аппроксимации полиномов. В четвертой и последующих строках перечисляются коэффициенты, начиная с нулевого, стоящие при многочленах соответствующих порядков. Используя полученные коэффициенты, можно записать аналитический вид получившейся аппроксимирующей функции.
При сохранении результатов в базу данных пользователю предлагается самостоятельно ввести название сохраняемых результатов. Рекомендуется генерировать название, придерживаясь нижеследующих правил. Название должно содержать (рисунок 4.9):
–фамилию исследователя;
–аппроксимируемый параметр (процесс, плотность, функция распределения).
Рисунок 4.9 – Генерация названия для сохранения результатов в базу данных
После сохранения результатов в базу данных можно изменить параметры аппроксимации и попытаться получить приближение анализируемой характеристики с меньшей среднеквадратической ошибкой, или же просто подобрать другую систему аппроксимирующих полиномов для дальнейшего анализа и сравнения результатов, сохраненных в базе данных.
22
4.3 Анализ результатов
Программное средство «АЗР» позволяет сохранять полученные на предыдущих этапах результаты в базу данных с целью их дальнейшего анализа. Общий вид программы на этапе анализа результатов, хранящихся в базе данных, представлен ан рисунке 4.10.
Рисунок 4.10 – База данных о результатах аппроксимаций
Как видно из рисунка 4.10, для управления базой данных предназначены три кнопки:
1)кнопка «Удалить запись» – удаляет выбранную запись из базы данных;
2)кнопка «Добавить в график» – добавляет результаты аппроксимации, соответствующие выбранной записи, на график;
3)кнопка «Очистить график» – удаляет все построения на графике.
С помощью данного блока программы можно сравнивать полученные различными методами результаты как по числовым критериям – среднеквадратической ошибке и ее квадрату, так и визуально, анализируя построенные графики аппроксимаций.
23
5 Содержание письменного отчета
Отчет должен быть оформлен на листах формата А4 в соответствии с действующими правилами оформления и содержать:
1)титульный лист;
2)постановку задачи с исходными данными для анализа;
3)краткое изложение теории;
4)результаты компьютерной обработки данных;
5)анализ полученных результатов;
6)содержательная интерпретация полученных результатов.
6 Вопросы к защите лабораторной работы
1)Сформулировать общую постановку задачи аппроксимации функции.
2)Опишите алгоритм аппроксимации функции обобщенным рядом Фурье.
3)Какие системы ортогональных и ортонормированных полиномов, используемых для аппроксимации функций, вам известны?
4)Для аппроксимации каких функций удобно использовать системы полиномов Лежандра и Чебышева?
5)Для аппроксимации каких функций удобно использовать систему полиномов Лагерра?
6)Для аппроксимации каких функций удобно использовать систему полиномов Эрмита?
7)Что означают термины «Серединная аппроксимация», «Левая аппроксимация», «Правая аппроксимация»?
8)Каким образом определяется оптимальное число членов аппроксимирующего ряда Фурье?
9)Какие функции называются ортогональными?
10)Какие функции называются ортонормированными?
11)Как рассчитывается среднеквадратическая погрешность аппроксимации?
12)Какие действия необходимо предпринять в случае, когда аппроксимируемая функция и система ортогональных многочленов определены на различных множествах?
13)Для чего нужна предварительная обработка эмпирических данных?
24
Список использованных источников
1Фихтенгольц, Г.М. Курс дифференциального и интегрального исчисления, т. III. – М.: Наука, 1966.
2Романенко, А. Ф. Аппроксимативные методы анализа случайных процессов / А. Ф. Романенко, Г. А. Сергеев . - М. : Энергия, 1974. - 177 с. : черт. -
Библиогр.: с. 175-176.
3Прохоров, С.А. Аппроксимативный анализ случайных процессов. – 2-е
изд., перераб. и доп. / СНЦ РАН, 2001. – 380с., ил.
4Вентцель, А. Д. Курс теории случайных процессов: учеб. пособие для вузов / А. Д. Вентцель . - М. : Наука, 1975. - 320 с.
5Марон, И. А. Дифференциальное и интегральное исчисление в примерах и задачах (функции одной переменной): учеб. пособие для втузов.- 2-е
изд., стер. - М. : Наука, 1973. - 400 с.
6Шевцов, Г. С. Численные методы линейной алгебры: учеб. пособие для математических направлений и специальностей / Г. С. Шевцов, О. Г. Крюкова,
Б. И. Мызникова . - Москва : Финансы и статистика, 2008. - 480 с. : ил.. - Библиогр.:
с. 473. - Предм. указ.: с. 475.
7Чураков, Е.П. Математические методы обработки экспериментальных данных в экономике: Учеб. пособие / Е.П. Чураков. - М.: Финансы и статистика, 2004. - 240 с.
8Треногин, В.А. Функциональный анализ: учеб. для вузов / В.А.
Треногин.- 3-е изд., испр. - М.: Физматлит, 2002. - 488 с.
9 Геворкян, П.С. Высшая математика. Интегралы, ряды, ТФКП,
дифференциальные уравнения: учебное пособие / П.С. Геворкян. – М.: Физматлит, 2007.
10 Каханер, Д. Численные методы и программное обеспечение = Numerical Methods and Software / Д. Каханер, К. Моулер, С. Нэш ; пер. с англ. под ред. Х. Д. Икрамова .- 2-е изд., стер. - М. : Мир, 2001. - 575 с.
25
Приложение А
(обязательное)
Исходные данные
Таблица А.1 – Таблица значений функции, измеренной в конечном числе точек
№ |
х |
f(x) |
|
№ |
х |
f(x) |
|
№ |
х |
f(x) |
1 |
0.000 |
1.109 |
|
41 |
0.157 |
4.302 |
|
81 |
0.314 |
1.520 |
2 |
0.004 |
11.350 |
|
42 |
0.161 |
50.645 |
|
82 |
0.318 |
22.914 |
3 |
0.008 |
2.127 |
|
43 |
0.165 |
1.669 |
|
83 |
0.322 |
1.427 |
4 |
0.012 |
3.953 |
|
44 |
0.169 |
1.004 |
|
84 |
0.325 |
1.052 |
5 |
0.016 |
2.870 |
|
45 |
0.173 |
45.447 |
|
85 |
0.329 |
19.775 |
6 |
0.020 |
1.861 |
|
46 |
0.176 |
1.074 |
|
86 |
0.333 |
6.086 |
7 |
0.024 |
1.605 |
|
47 |
0.180 |
1.197 |
|
87 |
0.337 |
5.234 |
8 |
0.027 |
1.484 |
|
48 |
0.184 |
1.060 |
|
88 |
0.341 |
1.477 |
9 |
0.031 |
3.260 |
|
49 |
0.188 |
1.982 |
|
89 |
0.345 |
116.609 |
10 |
0.035 |
2.334 |
|
50 |
0.192 |
1.603 |
|
90 |
0.349 |
1.473 |
11 |
0.039 |
2.708 |
|
51 |
0.196 |
1.010 |
|
91 |
0.353 |
13.710 |
12 |
0.043 |
39.242 |
|
52 |
0.200 |
1.217 |
|
92 |
0.357 |
1.022 |
13 |
0.047 |
1.075 |
|
53 |
0.204 |
1.357 |
|
93 |
0.361 |
1.701 |
14 |
0.051 |
5.883 |
|
54 |
0.208 |
1.518 |
|
94 |
0.365 |
1.965 |
15 |
0.055 |
2.692 |
|
55 |
0.212 |
1.451 |
|
95 |
0.369 |
1.228 |
16 |
0.059 |
6.353 |
|
56 |
0.216 |
3.149 |
|
96 |
0.373 |
3.919 |
17 |
0.063 |
1.506 |
|
57 |
0.220 |
2.221 |
|
97 |
0.376 |
1.476 |
18 |
0.067 |
1.060 |
|
58 |
0.224 |
1.032 |
|
98 |
0.380 |
4.412 |
19 |
0.071 |
37.448 |
|
59 |
0.227 |
4.018 |
|
99 |
0.384 |
2.068 |
20 |
0.075 |
1.822 |
|
60 |
0.231 |
1.560 |
|
100 |
0.388 |
2.264 |
21 |
0.078 |
1.796 |
|
61 |
0.235 |
21.772 |
|
101 |
0.392 |
1.820 |
22 |
0.082 |
1.542 |
|
62 |
0.239 |
1.328 |
|
102 |
0.396 |
3.650 |
23 |
0.086 |
86.456 |
|
63 |
0.243 |
2.720 |
|
103 |
0.400 |
2.117 |
24 |
0.090 |
1.810 |
|
64 |
0.247 |
3.184 |
|
104 |
0.404 |
1.401 |
25 |
0.094 |
1.108 |
|
65 |
0.251 |
1.640 |
|
105 |
0.408 |
8.029 |
26 |
0.098 |
1.188 |
|
66 |
0.255 |
3.290 |
|
106 |
0.412 |
19.516 |
27 |
0.102 |
1.203 |
|
67 |
0.259 |
60.013 |
|
107 |
0.416 |
1.666 |
28 |
0.106 |
1.610 |
|
68 |
0.263 |
2.669 |
|
108 |
0.420 |
20.102 |
29 |
0.110 |
1.064 |
|
69 |
0.267 |
1.496 |
|
109 |
0.424 |
1.520 |
30 |
0.114 |
1.914 |
|
70 |
0.271 |
1.624 |
|
110 |
0.427 |
1.380 |
31 |
0.118 |
1.184 |
|
71 |
0.275 |
1.008 |
|
111 |
0.431 |
14.409 |
32 |
0.122 |
9.259 |
|
72 |
0.278 |
1.670 |
|
112 |
0.435 |
4.129 |
33 |
0.125 |
1.220 |
|
73 |
0.282 |
1.107 |
|
113 |
0.439 |
13.008 |
34 |
0.129 |
9.822 |
|
74 |
0.286 |
8.687 |
|
114 |
0.443 |
10.597 |
35 |
0.133 |
1.369 |
|
75 |
0.290 |
10.382 |
|
115 |
0.447 |
4.249 |
36 |
0.137 |
1.349 |
|
76 |
0.294 |
1.209 |
|
116 |
0.451 |
2.187 |
37 |
0.141 |
1.774 |
|
77 |
0.298 |
1.836 |
|
117 |
0.455 |
1.468 |
38 |
0.145 |
2.498 |
|
78 |
0.302 |
2.783 |
|
118 |
0.459 |
2.738 |
39 |
0.149 |
1.152 |
|
79 |
0.306 |
2.370 |
|
119 |
0.463 |
1.497 |
40 |
0.153 |
2.799 |
|
80 |
0.310 |
1.006 |
|
120 |
0.467 |
1.074 |
26
Продолжение таблицы А.1
№ |
х |
f(x) |
|
№ |
х |
f(x) |
|
№ |
х |
f(x) |
121 |
0.471 |
2.330 |
|
166 |
0.647 |
2.663 |
|
211 |
0.824 |
5.042 |
122 |
0.475 |
5.967 |
|
167 |
0.651 |
7.413 |
|
212 |
0.827 |
2.003 |
123 |
0.478 |
2.149 |
|
168 |
0.655 |
1.496 |
|
213 |
0.831 |
4.627 |
124 |
0.482 |
1.113 |
|
169 |
0.659 |
1.108 |
|
214 |
0.835 |
2.380 |
125 |
0.486 |
1.911 |
|
170 |
0.663 |
2.589 |
|
215 |
0.839 |
1.656 |
126 |
0.490 |
1.495 |
|
171 |
0.667 |
1.237 |
|
216 |
0.843 |
1.710 |
127 |
0.494 |
3.381 |
|
172 |
0.671 |
7.599 |
|
217 |
0.847 |
1.276 |
128 |
0.498 |
2.874 |
|
173 |
0.675 |
3.781 |
|
218 |
0.851 |
21.347 |
129 |
0.502 |
1.916 |
|
174 |
0.678 |
2.261 |
|
219 |
0.855 |
3.046 |
130 |
0.506 |
8.300 |
|
175 |
0.682 |
42.610 |
|
220 |
0.859 |
1.215 |
131 |
0.510 |
7.774 |
|
176 |
0.686 |
1.107 |
|
221 |
0.863 |
6.158 |
132 |
0.514 |
1.432 |
|
177 |
0.690 |
1.008 |
|
222 |
0.867 |
1.531 |
133 |
0.518 |
1.456 |
|
178 |
0.694 |
4.033 |
|
223 |
0.871 |
4.644 |
134 |
0.522 |
1.461 |
|
179 |
0.698 |
3.345 |
|
224 |
0.875 |
4.242 |
135 |
0.525 |
7.613 |
|
180 |
0.702 |
26.640 |
|
225 |
0.878 |
8.161 |
136 |
0.529 |
1.014 |
|
181 |
0.706 |
1.677 |
|
226 |
0.882 |
1.090 |
137 |
0.533 |
1.260 |
|
182 |
0.710 |
1.340 |
|
227 |
0.886 |
1.303 |
138 |
0.537 |
1.132 |
|
183 |
0.714 |
1.156 |
|
228 |
0.890 |
4.859 |
139 |
0.541 |
1.438 |
|
184 |
0.718 |
6.117 |
|
229 |
0.894 |
2.889 |
140 |
0.545 |
1.511 |
|
185 |
0.722 |
1.368 |
|
230 |
0.898 |
1.243 |
141 |
0.549 |
8.666 |
|
186 |
0.725 |
1.165 |
|
231 |
0.902 |
1.795 |
142 |
0.553 |
1.936 |
|
187 |
0.729 |
1.217 |
|
232 |
0.906 |
8.509 |
143 |
0.557 |
16.392 |
|
188 |
0.733 |
1.131 |
|
233 |
0.910 |
1.141 |
144 |
0.561 |
1.699 |
|
189 |
0.737 |
2.708 |
|
234 |
0.914 |
10.757 |
145 |
0.565 |
1.238 |
|
190 |
0.741 |
1.228 |
|
235 |
0.918 |
55.916 |
146 |
0.569 |
2.589 |
|
191 |
0.745 |
1.902 |
|
236 |
0.922 |
1.167 |
147 |
0.573 |
1.947 |
|
192 |
0.749 |
2.001 |
|
237 |
0.925 |
1.004 |
148 |
0.576 |
3.002 |
|
193 |
0.753 |
96.658 |
|
238 |
0.929 |
1.610 |
149 |
0.580 |
2.461 |
|
194 |
0.757 |
2.390 |
|
239 |
0.933 |
11.550 |
150 |
0.584 |
1.471 |
|
195 |
0.761 |
6.118 |
|
240 |
0.937 |
1.314 |
151 |
0.588 |
1.631 |
|
196 |
0.765 |
6.574 |
|
241 |
0.941 |
1.338 |
152 |
0.592 |
2.878 |
|
197 |
0.769 |
1.611 |
|
242 |
0.945 |
1.694 |
153 |
0.596 |
1.522 |
|
198 |
0.773 |
10.680 |
|
243 |
0.949 |
1.509 |
154 |
0.600 |
2.860 |
|
199 |
0.776 |
1.008 |
|
244 |
0.953 |
8.671 |
155 |
0.604 |
2.315 |
|
200 |
0.780 |
1.018 |
|
245 |
0.957 |
1.214 |
156 |
0.608 |
7.534 |
|
201 |
0.784 |
1.862 |
|
246 |
0.961 |
1.288 |
157 |
0.612 |
2.486 |
|
202 |
0.788 |
19.196 |
|
247 |
0.965 |
2.809 |
158 |
0.616 |
3.014 |
|
203 |
0.792 |
4.743 |
|
248 |
0.969 |
1.272 |
159 |
0.620 |
3.749 |
|
204 |
0.796 |
2.533 |
|
249 |
0.973 |
24.308 |
160 |
0.624 |
1.022 |
|
205 |
0.800 |
21.543 |
|
250 |
0.976 |
1.015 |
161 |
0.627 |
1.357 |
|
206 |
0.804 |
2.333 |
|
251 |
0.980 |
2.048 |
162 |
0.631 |
1.047 |
|
207 |
0.808 |
1.078 |
|
252 |
0.984 |
11.268 |
163 |
0.635 |
1.260 |
|
208 |
0.812 |
39.960 |
|
253 |
0.988 |
8.830 |
164 |
0.639 |
2.698 |
|
209 |
0.816 |
1.160 |
|
254 |
0.992 |
1.659 |
165 |
0.643 |
6.462 |
|
210 |
0.820 |
1.375 |
|
255 |
0.996 |
1.217 |
27
