MATLAB三次样条插值之三弯矩法

MATLAB三次样条插值的三弯矩法首先，据说这个程序并不完美。为了解决一般(1，2，n)格式并调用追赶程序，对于节点数量少于三个的情况没有分类讨论。我希望一些朋友能给我更好的方法。首先，通过函数sanwanj，即追赶法方程的四个向量参数，得到方程的系数矩阵。然后调用(intersanwj函数中的)追赶法来获得三次样条分段函数的系数因子。然后进行多项式组合，得到分段函数的解析表达式。最后，程序通过判断输入值的间隔自动选择相应的分段函数，并计算校正点的值。附录：追踪程序% % % % % % % % % % % % % % % %函数newv，w，newu，newd=sanwj(x，y，x0，y0，y1a，y1b)% 3弯矩样条插值%插值点将输入两次，x0 y0将单独输入。建议用y2a和y2b代替%边界条件a的二阶导数y1a和b的二阶导数y1b，以区别于三角码。n=长度(x)；m=长度(y)；如果m=n错误(x或y输入错误，再次出现)；目标v=1(n-1，1)；u=1(n-1，1)；d=零(n-1，1)；w=2 * 1(n1)；h0=x(1)-x0；h=零(n-1，1)；对于k=1:n-1h(k)=x(k1)-x(k)；目标v(1)=h0/(h0 h(1)；u(1)=1-v(1)；d(1)=6 *(y(2)-y(1)/h(1)-(y(1)-y0)/h0)/(h0 h(1)；%对于k=2:n-1v(k)=h(k-1)/(h(k-1)h(k)；u(k)=1-v(k)；d(k)=6 *(y(k-1)-y(k)/h(k)-(y(k)-y(k-1)/h(k-1)/(h(k-1)h(k)；目标五世；1；newu=1；u；d0=6 *(y(1)-y0)/h0-y1a)/h0；d(n)=6 *(y1b-(y(n)-y(n-1)/h(n-1)/h(n-1)；newd=d0；d；% % % % % % % % % % % % %函数intersanwj(x，y，x0，y0，y1a，y1b)% 3弯矩样条插值%第一部分n=长度(x)；m=长度(y)；如果m=n错误(x或y输入错误，再次出现)；目标%重定义hh=零(n，1)；h(1)=x(1)-x0；对于k=2:nh(k)=x(k)-x(k-1)；目标% sptep1调用三弯矩函数a，b，c，d=sanwj(x，y，x0，y0，y1a，y1b)；%三对角方程M=chase(a，b，c，d)。%查找插值函数三次样条函数表达式(三个弯矩)插值 n；符号十；fprintf(S0-1: n)；S(1)=收集(1/6)*M(2)*(X-x0)。3-(1/6)*M(1)*(X-x(1).3 (y(1)-(M(2)*h(1)。2)/6)*(X-x0)-(y0-(M(1)*h(1).2)/6)*(x-x(1)/h(1)；对于k=2:nfprintf(S % d-u 360 n，k-1，k)；s(k)=收集(1/6)*M(k 1)*(X-x(k-1)。3-(1/6)*M(k)*(X-x(k).3 (y(k)-(M(k 1)*h(k)。2)/6)*(x-x(k-1)-(y(k-1)-(m(k)*h(k).2)/6)*(x-x(k)/h(k)；目标S=S显示(S)；Fprintf(上面是样条函数(三个弯矩)的解析表达式，以手写形式显示如下： n)；漂亮；%第二部分%您想继续运行程序吗Myloop=输入