三次样条差值-matlab通用程序

103 阅读 0 评论 68 点赞

我是靠谱客的博主无限板凳，这篇文章主要介绍三次样条差值-matlab通用程序，现在分享给大家，希望可以做个参考。

数值分析课程的大作业，教材《数值分析》李乃成.梅立泉

clear
x=input('请按照格式[x1,x2,x3...]格式输入y=f(x)函数已知点的横坐标xi='); %三次样条差值函数
y=input('请按照格式[y1,y2,y3...]格式输入y=f(x)函数已知点对应的纵坐标yi=');
n=size(x,2); %特别注意，matlab中的矩阵编号是从1开始的，而教材上的矩阵编号是从0开始的
for k=2:n %计算h(i)
h(k)=x(k)-x(k-1);
end
for k=2:(n-1) %计算μ和λ
mu(k)=h(k)/(h(k)+h(k+1));
lumbda(k)=1-mu(k);
end
for k=2:(n-1)
d(k)=6*((y(k+1)-y(k))/h(k+1)-(y(k)-y(k-1))/h(k))/(h(k)+h(k+1)); %计算di
end
fprintf('边界条件类型选择:n1.已知f(a)和f(b)的二阶导数n2.