rwise if x>=s(k-1)&x<=s(k); l=(x-s(k-1))./(s(k)-s(k-1)); else if x>=s(k)&x<=s(k+1); l=(x-s(k+1))./(s(k)-s(k+1)); else l=0; end end end ff=ff+Rf(s(k))*l; %% 求插值函数值 end m=m+1; f(m)=ff; end %%% 作出曲线 x=-1:hh:1; plot(x,f,'r'); grid on hold on 三次样条插值源程序: (采用第一边界条件) clear . 10. n=input(' 将区间分为的等份数输入: \n'); %%% 插值区间 a=-1; b=1; hh=0.001 ;% 画图的步长 s=[a+(b-a)/n*[0:n]]; %%% 给定的定点, Rf 为给定的函数%%%% 第一边界条件 Rf"(-1),Rf"(1) v=5000*1/(1+25*a*a)^3-50/(1+25*a*a)^4; for k=1:n ;% 取出节点间距 h(k)=s(k+1)-s(k); end for k=1:n-1; % 求出系数向量 lamuda,miu la(k)=h(k+1)/(h(k+1)+h(k)); miu(k)=1-la(k); end %%%% 赋值系数矩阵 A for k=1:n-1; for p=1:n-1; switch p case k A(k,p)=2; case k-1 A(k,p)=miu(p+1); case k+1 A(k,p)=la(p-1); otherwise A(k,p)=0; end end end %%%% 求出 d阵 for k=1:n-1; switch k case 1 d(k)=6*f2c([s(k) s(k+1) s(k+2)])-miu(k)*v;
全文预览
