Построение кубического сплайна в пакете 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
|
ВСТРОЕННЫЕ ФУНКЦИИ ИНТЕРПОЛИРОВАНИЯ
Важнейшие способы обработки и анализа рядов динамики Не во всех случаях эмпирические данные рядов динамики позволяют определить тенденцию изменения явления во времени...
|
ТЕОРЕТИЧЕСКАЯ МЕХАНИКА Статика является частью теоретической механики, изучающей условия, при которых тело находится под действием заданной системы сил...
|
Теория усилителей. Схема Основная масса современных аналоговых и аналого-цифровых электронных устройств выполняется на специализированных микросхемах...
|
Логические цифровые микросхемы Более сложные элементы цифровой схемотехники (триггеры, мультиплексоры, декодеры и т.д.) не имеют...
|
Тема: Изучение приспособленности организмов к среде обитания Цель:выяснить механизм образования приспособлений к среде обитания и их относительный характер, сделать вывод о том, что приспособленность – результат действия естественного отбора...
Тема: Изучение фенотипов местных сортов растений Цель: расширить знания о задачах современной селекции. Оборудование:пакетики семян различных сортов томатов...
Тема: Составление цепи питания Цель: расширить знания о биотических факторах среды. Оборудование:гербарные растения...
|
Психолого-педагогическая характеристика студенческой группы
Характеристика группы составляется по 407 группе очного отделения зооинженерного факультета, бакалавриата по направлению «Биология» РГАУ-МСХА имени К...
Общая и профессиональная культура педагога: сущность, специфика, взаимосвязь Педагогическая культура- часть общечеловеческих культуры, в которой запечатлил духовные и материальные ценности образования и воспитания, осуществляя образовательно-воспитательный процесс...
Устройство рабочих органов мясорубки Независимо от марки мясорубки и её технических характеристик, все они имеют принципиально одинаковые устройства...
|
|