Построение кубического сплайна в пакете MATLAB.
Function res=G(xx,x,f,M)
% исходные данные – x, f; М – значения вторых производных
n=size(x);
n=n(2);
h=zeros(1,n-1);
% вычисление шага
for i=1:n-1
h(i)=x(i+1) - x(i);
end;
for j=1:length(xx)
X=xx(j);
% определение номера интервала
for i=1:n-1
if x(i)<=X && X<x(i+1)
k=i;
break
end;
end;
% вычисление значения сплайна в промежуточной точке
yy(j)=M(k)*((x(k+1)-X)^3)/(6*h(k)) + M(k+1)*((X-x(k))^3)/(6*h(k)) +...
(f(k)-M(k)*(h(k)^2)/6)*(x(k+1)-X)/h(k)+...
(f(k+1)-M(k+1)*(h(k)^2)/6)*(X-x(k))/h(k);
end;
res=yy;
return
% вызов процедур и построение графика
>> x = [-4 -3.5 -2 0 2.1 2.5 5];
>> f=[0.5 0.2 -0.7 0 1.1 0.8 -1.6];
>> xx=linspace(-4,5,1000)
>> M = Work(x,f) % функция нахождения вторых производных
>> yy=G(xx,x,f,M);
>> plot(xx,yy,'r')
Построение кубического сплайна в пакете MATHCAD.
Задав исходные данные и определив шаг для каждого интервала, можно определить функции для формирования матриц A и H и функцию для построения кубического сплайна g(X):
Варианты лабораторных работ
Номер варианта
| Исходные данные
|
| x
y
| 1,4
0,3365
| 1,8
0,5878
| 2,3
0,8329
| 2,9
1,0647
| 3,2
1,1632
| 3,6
1,2809
|
| x
y
| 2,0
0,6931
| 2,5
0,9163
| 2,8
1,0296
| 3,3
1,1939
| 3,6
1,2809
| 4.0
1,3863
|
| x
y
| 4,0
1,3863
| 4,5
1,5041
| 4,9
1,5892
| 5,4
1,6864
| 5,7
1,7405
| 6,0
1,7918
|
| x
y
| 1,2
0,1823
| 1,6
0,4700
| 2,1
0,7419
| 2,6
0,9555
| 3,0
1,0986
| 3,3
1,1939
|
| x
y
| 2,2
0,7885
| 2,7
0,9933
| 3,1
1,1314
| 3,6
1,2809
| 4,0
1,3863
| 4,3
1,4586
|
| x
y
| 3,2
1,1632
| 3,6
1,2809
| 4,1
1,4110
| 4,6
1,5261
| 4,9
1,5892
| 5.4
1,6864
|
| x
y
| 3,4
1,2238
| 3,9
1,3610
| 4,3
1,4586
| 4,9
1,5892
| 5,2
1,6487
| 5,6
1,7228
|
| x
y
| 1,6
0,4700
| 2,1
0,7419
| 2,7
0,9933
| 3,2
1,1632
| 3,6
1,2809
| 4,1
1,4110
|
| x
y
| 2,8
1,0296
| 3,1
1,1314
| 3,7
1,3083
| 4,2
1,4351
| 4,6
1,5261
| 5,0
1,6094
|
| x
y
| 3,1
1,1314
| 3,6
1,2809
| 4,0
1,3863
| 4,6
1,5261
| 4,9
1,5892
| 5,3
1,6677
|
| x
y
| 1,9
0,6419
| 2,5
0,9163
| 2,9
1,0647
| 3,4
1,2238
| 3,6
1,2809
| 4,0
1,3863
|
| x
y
| 1,7
0,5306
| 2,2
0,7865
| 2,8
1,0296
| 3,2
1,1632
| 3,5
1,2528
| 4,0
1,3863
|
| x
y
| 3,6
1,2809
| 4,2
1,4351
| 4,5
1,5041
| 5,2
1,6487
| 5,5
1,7047
| 5,9
1,7750
|
Варианты лабораторных работ (окончание)
Номер варианта
| Исходные данные
|
| x
y
| 3,3
1,1939
| 3,9
1,3610
| 4,4
1,4816
| 5,0
1,6094
| 5,4
1,6864
| 5,9
1,7750
|
| x
y
| 1,1
0,0953
| 1,7
0,5306
| 2,4
0,8755
| 2,8
1,0296
| 3,3
1,1939
| 3,6
1,2809
|
| x
y
| 2,1
0,4718
| 2,5
0,9163
| 3,0
1,0986
| 3,5
1,2528
| 3,8
1,3350
| 4,2
1,4351
|
| x
y
| 3,2
1,1632
| 3,7
1,3083
| 4,3
1,4586
| 4,9
1,5892
| 5,2
1,6487
| 5,6
1,7228
|
| x
y
| 2,7
0.9933
| 3,3
1,1939
| 3.8
1,3350
| 4.6
1,5261
| 5,0
1,6094
| 5,5
1,7047
|
| x
y
| 1,0
0,0000
| 1,5
0,4055
| 2,1
0,7419
| 2,7
0,9933
| 3,0
1,0966
| 3,4
1,2238
|
| x
y
| 1,4
0,3365
| 1,9
0,6419
| 2,6
0,9555
| 3,0
1,0986
| 3,3
1,1939
| 3,6
1,2809
|
| x
y
| 3,1
1,1314
| 3,7
1,3083
| 4,2
1,4351
| 4.8
1,5686
| 5,2
1,6487
| 5,5
1,7047
|
| x
y
| 2,6
0,9555
| 3,2
1,1632
| 4,0
1,3863
| 4,5
1,5041
| 4,9
1,5692
| 5,4
1,6864
|
| x
y
| 1,6
0,4700
| 2,2
0,7885
| 2,7
0,9933
| 3,4
1,2238
| 3,6
1,2809
| 4,0
1,3836
|
| x
y
| 2,1
0,7419
| 2,7
0,9933
| 3,33
1,1939
| 3,8
1,3350
| 4,0
1,3863
| 4,4
1,4816
|
| x
y
| 2,6
0,9555
| 3,0
1,0986
| 3,9
1,3610
| 4,5
1,5041
| 4,8
1,5686
| 5,3
1,6677
|
| x
y
| 4,5
1,5041
| 4,9
1,5892
| 5,5
1,7047
| 6,0
1,7918
| 6,2
1,8245
| 6,5
1,8718
|
| x
y
| 3,5
1,2528
| 3,8
1,3350
| 4,5
1,5041
| 5,1
1,6292
| 5,4
1,6864
| 5,8
1,7579
|
| x
y
| 2,8
1,0296
| 3,3
1,1939
| 3,9
1,3610
| 4,6
1,5261
| 5,0
1,6094
| 5,5
1,7047
|
| x
y
| 4,1
1,4110
| 4,6
1,5261
| 5,2
1,6487
| 6,0
1,7918
| 6,2
1,8245
| 6,5
1,8718
|
ВСТРОЕННЫЕ ФУНКЦИИ ИНТЕРПОЛИРОВАНИЯ
Обзор компонентов Multisim Компоненты – это основа любой схемы, это все элементы, из которых она состоит. Multisim оперирует с двумя категориями...
|
Композиция из абстрактных геометрических фигур Данная композиция состоит из линий, штриховки, абстрактных геометрических форм...
|
Важнейшие способы обработки и анализа рядов динамики Не во всех случаях эмпирические данные рядов динамики позволяют определить тенденцию изменения явления во времени...
|
ТЕОРЕТИЧЕСКАЯ МЕХАНИКА Статика является частью теоретической механики, изучающей условия, при которых тело находится под действием заданной системы сил...
|
|
Седалищно-прямокишечная ямка Седалищно-прямокишечная (анальная) ямка, fossa ischiorectalis (ischioanalis) – это парное углубление в области промежности, находящееся по бокам от конечного отдела прямой кишки и седалищных бугров, заполненное жировой клетчаткой, сосудами, нервами и...
Основные структурные физиотерапевтические подразделения Физиотерапевтическое подразделение является одним из структурных подразделений лечебно-профилактического учреждения, которое предназначено для оказания физиотерапевтической помощи...
Почему важны муниципальные выборы? Туристическая фирма оставляет за собой право, в случае причин непреодолимого характера, вносить некоторые изменения в программу тура без уменьшения общего объема и качества услуг, в том числе предоставлять замену отеля на равнозначный...
|
|
Тема 5. Организационная структура управления гостиницей 1. Виды организационно – управленческих структур. 2. Организационно – управленческая структура современного ТГК...
Методы прогнозирования национальной экономики, их особенности, классификация В настоящее время по оценке специалистов насчитывается свыше 150 различных методов прогнозирования, но на практике, в качестве основных используется около 20 методов...
Методы анализа финансово-хозяйственной деятельности предприятия
Содержанием анализа финансово-хозяйственной деятельности предприятия является глубокое и всестороннее изучение экономической информации о функционировании анализируемого субъекта хозяйствования с целью принятия оптимальных управленческих...
|
|