Пример построения полинома Лагранжа в пакете MATLAB.
function yy=lagrange(x,y,xx)
% число узлов интерполяции
N=length(x);
% число узлов, в которых высчитывается значение интерполяционного полинома
N_res=length(xx);
% создание нулевого массива значений интерполяционного полинома
yy=zeros(size(xx));
for k=1:N
% вычисление функции Лагранжа Li(X)
Li=ones(size(xx));
for j=[1:k-1, k+1:N]
for i=[1:N_res]
Li(i)=Li(i).*(xx(i)-x(j))/(x(k)-x(j));
end
end
% накопление суммы
yy = yy + y(k)*Li;
yy
end
Пример использования узлов Чебышева.
Заданная функция f(x) табулируется на интервале [-1,1] двумя способами: с шагом 0.2 и в узлах Чебышева. На двух полученных сетках строятся полиномы Лагранжа. Результаты визуализируются с помощью графиков. На графике, приведенном ниже, видно значительное отклонение полинома Лагранжа, построенного на сетке с шагом 0.2, от исходной функции и более приемлемое приближение функции для полинома, построенного на узлах Чебышева:
Узлы Чебышева:
Пример вычисления погрешности интерполирования.
Анализ погрешности замены исходной функции интерполяционным многочленом для таблиц с постоянным шагом:
Варианты лабораторных работ
| Номер
варианта
| Исходные данные
| |
| x
f(x)
| 1,4
0,3365
| 1,8
0.5878
| 2,3
0.8329
| 2,9
1.0647
| |
| x
f(x)
| 2,0
0,6931
| 2,5
0.9163
| 2,8
1.029б
| 3,3
1,1939
| |
| x
f(x)
| 4,0
1,3863
| 4,5
1,5041
| 4,9
1.5892
| 5,4
1.6864
| |
| x
f(x)
| 1,2
0,1823
| 1,6
0,4700
| 2,1
0.7419
| 2,6
1,6864
| |
| x
f(x)
| 2,2
0,7885
| 2,7
0,9933
| 3,1
1.1314
| 3,6
1,2809
| | | | | | | | |
Варианты лабораторных работ (продолжение)
| Номер
варианта
| Исходные данные
| | |
| x
f(x)
| 3,2
1,1632
| 3,6
1.2809
| 4,1
1,4110
| 4,6
1.5261
| | |
| x
f(x)
| 3,4
1,2238
| 3,9
1,3610
| 4.3
1,4586
| 4,9
1.5861
| | |
| x
f(x)
| 1,6
0,4700
| 2,1
0,7419
| 2,7
0,9933
| 8,2
1,1632
| | |
| x
f(x)
| 2,8
1,0296
| 3,1
1.1314
| 3,7
1,3083
| 4,2
1,4351
| | |
| x
f(x)
| 3,1
1,1314
| 3,6
1,2809
| 4,0
1,З863
| 4,6
1,5261
| | |
| x
f(x)
| 1,9
0,6419
| 2,5
0,9163
| 2,9
1,0647
| 3.4
1.2238
| | |
| x
f(x)
| 1,7
0,5306
| 2,2
0,7885
| 2,8
1.0296
| 3,2
1.1632
| | |
| x
f(x)
| 3,6
1.2809
| 4,2
1.4351
| 4,5
1,5041
| 5.2
1.6094
| | |
| x
f(x)
| 2,5
0,9163
| 2,9
1,0647
| 3,6
1.2809
| 4,1
1.4110
| | |
| x
f(x)
| 3,3
1,1939
| 3,9
1,3610
| 4,4
1,4816
| 5,0
1,6094
| |
| x
f(x)
| 1,1
0,0953
| 1,7
0,5306
| 2,4
0.6755
| 2,8
1,0296
| |
| x
f(x)
| 2,1
0.7419
| 2.5
0.9163
| 3,0
1,0986
| 3,5
1,2528
| |
| x
f(x)
| 3,2
1,1632
| 3,7
1.3083
| 4,3
1,4586
| 4,9
1,5892
| |
| x
f(x)
| 2,7
0.9933
| 3,3
1.1939
| 3,8
1.3350
| 4,6
1,5261
| |
| x
f(x)
| 1,0
0,0000
| 1,5
0.4055
| 2,1
0,7419
| 2,7
0.9933
| |
| x
f(x)
| 1,4
0,3365
| 1.9
0,6419
| 2,6
0,9555
| 3,0
1,0986
| |
| x
f(x)
| 3,1
1.1314
| 3.7
1,3083
| 4,2
1.4351
| 4.8
1.5686
| |
| x
f(x)
| 2,6
0.9555
| 3,2
1.1632
| 4.0
1.3863
| 4,5
1.5041
| |
| x
f(x)
| 1,6
0.4700
| 2,2
0,7885
| 2,7
.9933
| 8,4
1,2238
| |
| x
f(x)
| 2,1
0,7419
| 2,7
0,9933
| 3,3
1.1939
| 3,8
1.3350
| | | | | | | | | | | | |
Варианты лабораторных работ (окончание)
| Номер
варианта
| Исходные данные
| | |
| x
f(x)
| 2,6
0,9555
| 3.0
1,0986
|
1,3610
| 4,5
1,5041
| |
| x
f(x)
| 4,5
1,5041
| 4,9
1.5892
| 5,5
1,7047
| 6,0
1.7916
| |
| x
f(x)
| 3,5
1,2528
| 3,8
1,3350
| 4,5
1,5041
| 5,1
1.6292
| |
| x
f(x)
| 2,6
1,0296
| 3,3
1.1939
| 3,9
1,3610
| 4,6
1,5261
| |
| x
f(x)
| 4,1
1,4110
| 4,6
1,5261
| 5,2
1,6487
| 6,0
1,7918
| | | | | | | | | |
Аальтернативная стоимость. Кривая производственных возможностей В экономике Буридании есть 100 ед. труда с производительностью 4 м ткани или 2 кг мяса...
|
Вычисление основной дактилоскопической формулы Вычислением основной дактоформулы обычно занимается следователь. Для этого все десять пальцев разбиваются на пять пар...
|
Расчетные и графические задания Равновесный объем - это объем, определяемый равенством спроса и предложения...
|
Кардиналистский и ординалистский подходы Кардиналистский (количественный подход) к анализу полезности основан на представлении о возможности измерения различных благ в условных единицах полезности...
|
ЛЕЧЕБНО-ПРОФИЛАКТИЧЕСКОЙ ПОМОЩИ НАСЕЛЕНИЮ В УСЛОВИЯХ ОМС 001. Основными путями развития поликлинической помощи взрослому населению в новых экономических условиях являются все...
МЕТОДИКА ИЗУЧЕНИЯ МОРФЕМНОГО СОСТАВА СЛОВА В НАЧАЛЬНЫХ КЛАССАХ В практике речевого общения широко известен следующий факт: как взрослые...
СИНТАКСИЧЕСКАЯ РАБОТА В СИСТЕМЕ РАЗВИТИЯ РЕЧИ УЧАЩИХСЯ В языке различаются уровни — уровень слова (лексический), уровень словосочетания и предложения (синтаксический) и уровень
Словосочетание в этом смысле может рассматриваться как переходное звено от лексического уровня к синтаксическому...
|
Стресс-лимитирующие факторы Поскольку в каждом реализующем факторе общего адаптационного синдрома при бесконтрольном его развитии заложена потенциальная опасность появления патогенных преобразований...
ТЕОРИЯ ЗАЩИТНЫХ МЕХАНИЗМОВ ЛИЧНОСТИ В современной психологической литературе встречаются различные термины, касающиеся феноменов защиты...
Этические проблемы проведения экспериментов на человеке и животных В настоящее время четко определены новые подходы и требования к биомедицинским исследованиям...
|
|
