三次样条插值作业题Word下载.docx
《三次样条插值作业题Word下载.docx》由会员分享，可在线阅读，更多相关《三次样条插值作业题Word下载.docx（14页珍藏版）》请在冰豆网上搜索。
%自变量x与因变量y，两个边界条件的取值
IndVar=[0,1,2,3];
DepVar=[0,,2,];
LeftBoun=;
RightBoun=-1;
%区间长度向量，其各元素为自变量各段的长度
h=zeros（1,length（IndVar）-1）;
fori=1:
length（IndVar）-1
h（i）=IndVar（i+1）-IndVar（i）;
end
%为向量μ赋值
mu=zeros（1,length（h））;
length（mu）-1
mu（i）=h（i）/（h（i）+h（i+1））;
mu（i+1）=1;
%为向量λ赋值
lambda=zeros（1,length（h））;
lambda
（1）=1;
fori=2:
length（lambda）
lambda（i）=h（i）/（h（i-1）+h（i））;
%为向量d赋值
d=zeros（1,length（h）+1）;
d
（1）=6*（（DepVar
（2）-DepVar
（1））/（IndVar
（2）-IndVar
（1））-LeftBoun）/h
（1）;
length（h）
a=（DepVar（i）-DepVar（i-1））/（IndVar（i）-IndVar（i-1））;
b=（DepVar（i+1）-DepVar（i））/（IndVar（i+1）-IndVar（i））;
c=（b-a）/（IndVar（i+1）-IndVar（i-1））;
d（i）=6*c;
d（i+1）=6*（RightBoun-（DepVar（i+1）-DepVar（i））/（IndVar（i+1）-IndVar（i）））/h（i）;
%为矩阵A赋值
%将主对角线上的元素全部置为2
A=zeros（length（d）,length（d））;
length（d）
A（i,i）=2;
%将向量λ的各元素赋给主对角线右侧第一条对角线
length（d）-1
A（i,i+1）=lambda（i）;
%将向量d的各元素赋给主对角线左侧第一条对角线
A（i+1,i）=mu（i）;
%求解向量M
M=A\d'
;
%求解每一段曲线的函数表达式
Coefs_1=M（i）/（6*h（i））;
Part_1=conv（Coefs_1,...
conv（[-1,IndVar（i+1）],...
conv（[-1,IndVar（i+1）],[-1,IndVar（i+1）]）））;
S_1=polyval（Part_1,[IndVar（i）:
:
IndVar（i+1）]）;
Coefs_2=M（i+1）/（6*h（i））;
Part_2=conv（Coefs_2,...
conv（[1,-IndVar（i）],...
conv（[1,-IndVar（i）],[1,-IndVar（i）]）））;
S_2=polyval（Part_2,[IndVar（i）:
Coefs_3=（DepVar（i）-M（i）*h（i）^2/6）/h（i）;
Part_3=conv（Coefs_3,[-1,IndVar（i+1）]）;
S_3=polyval（Part_3,[IndVar（i）:
Coefs_4=（DepVar（i+1）-M（i+1）*h（i）^2/6）/h（i）;
Part_4=conv（Coefs_4,[1,-IndVar（i）]）;
S_4=polyval（Part_4,[IndVar（i）:
S=S_1+S_2+S_3+S_4;
plot（[IndVar（i）:
IndVar（i+1）],S,'
LineWidth'
%在样条插值曲线的相应位置标注该段曲线的函数表达式
text（i-1,polyval（Part_1,3）,...
['
\itS'
num2str（i）,'
（x）='
num2str（Coefs_1）,'
（'
num2str（IndVar（i+1））,'
-x）^{3}+'
...
num2str（Coefs_2）,'
（x-'
num2str（IndVar（i））,'
）^{3}+'
num2str（Coefs_3）,...
'
-x）+'
num2str（Coefs_4）,'
）'
],...
FontName'
'
TimesNewRoman'
FontSize'
14）
holdon
%过x=1和x=2两个横轴点作垂线%
line（[1,1],[,],'
LineStyle'
--'
）;
line（[2,2],[,],'
%为x轴和y轴添加标注
xlabel（'
\itx'
14,'
FontWeight'
bold'
ylabel（'
\its（x）'
'
Rotation'
0,'
最终，三次样条插值函数s（x）表达式为：
曲线的图像如图所示：
例2已知函数值表：
4
5
试求在区间[1,5]上满足上述函数表所给出的插值条件的三次自然样条插值函数
%本段代码解决作业题的例2
%自变量x与因变量y的取值
IndVar=[1,2,4,5];
DepVar=[1,3,4,2];
mu（i+1）=0;
lambda
（1）=0;
d
（1）=0;
d（i+1）=0;
text（i,polyval（Part_1,5）,...
%过x=2和x=4两个横轴点作垂线%
line（[4,4],[,],'
最终，三次自然样条插值函数s（x）表达式为：
例3课后习题与思考题第7题
试求在区间[,]上满足上述函数表所给出的插值条件的三次自然样条插值函数s（x）
求解出的三次样条插值函数将写成三弯矩方程的形式：
本题采用和例2基本相同的Matlab代码，只改变初始条件。
例4课后习题与思考题第6题
求
的二次插值式
，使：
并计算
的近似值并估计误差。
本题采用拉格朗日二次插值法进行计算：
%本段代码解决课本第2章习题与思考题第6题
IndVar=[100,121,144];
DepVar=[10,11,12];
%构造拉格朗日插值函数
Coefs_1=DepVar
（1）/...
（（IndVar
（1）-IndVar
（2））*（IndVar
（1）-IndVar（3）））;
Part_1=conv（Coefs_1,...
conv（[1,-IndVar
（2）],[1,-IndVar（3）]））;
f_1=polyval（Part_1,[IndVar
（1）:
IndVar（3）]）;
Coefs_2=DepVar
（2）/...
（（IndVar
（2）-IndVar
（1））*（IndVar
（2）-IndVar（3）））;
Part_2=conv（Coefs_2,...
conv（[1,-IndVar
（1）],[1,-IndVar（3）]））;
f_2=polyval（Part_2,[IndVar
（1）:
Coefs_3=DepVar（3）/...
（（IndVar（3）-IndVar
（1））*（IndVar（3）-IndVar
（2）））;
Part_3=conv（Coefs_3,...
conv（[1,-IndVar
（1）],[1,-IndVar
（2）]））;
f_3=polyval（Part_3,[IndVar
（1）:
f=f_1+f_2+f_3;
plot（[IndVar
（1）:
IndVar（3）],f,'
;
%在样条插值曲线的相应位置标注该段曲线的函数表达式
text（110,polyval（Part_1,110）+polyval（Part_2,110）+polyval（Part_3,110）,...
['
\itP_2（x）='
num2str（IndVar
（2））,...
）（x-'
num2str（IndVar（3））,'
）+'
num2str（Coefs_2）,'
num2str（IndVar
（1））,...
num2str（Coefs_3）,'
num2str（IndVar
（2））,'
\itP_2（x）'
%求115平方根的估计值
fun115=polyval（Part_1,115）+polyval（Part_2,115）+polyval（Part_3,115）
%求115平方根的准确值
sqrt（115）
%求误差，单位为%
est=100*（fun115-sqrt（115））/fun115
最终，拉格朗日插值函数
的表达式为：
从计算结果可以看出，
的准确值为，而通过拉格朗日插值法求出的近似值为，误差为%
的图像为：
升级会员 