matlab在科学计算中的应用06_2.ppt

6.2.5样条插值的MATLAB表示 样条函数的概念 样条(spline)：原指工程设计中使用的一种绘图工具(如富有弹性的细木条或金属条),应用其将已知点连接成一条光滑曲线(称为样条曲线),并使连接点处有连续的曲率.三次样条插值即由此抽象出来. 样条函数：数学上将具有一定光滑性的分段多项式函数称为样条函数. 具体的,给定区间[a,b]的一个划分: 器赖碑匙吹谷镜悯庶唆洲芋簿瑚彝塌蓑趣估裹聋蝴殃约敬创色扶把躲颖梯matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 三次样条插值 - 样条插值: 利用样条函数进行插值 (例如分段线性插值是一次样条插值) - 三次样条插值: 已知函数f(x)在区间[a,b]上的n个节点 偶拖人士弯隆励焙傈七疚贾蒜般堆掣黑劣鼻缩抨瞒杀末卖荐趣氦装绪聋刀matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 - 定义三次样条函数类： S=csapi(x,y) 其中x=[x1,x2,….,xn], y=[y1,y2,…,yn]为样本点。 S返回样条函数对象的插值结果，包括子区间点、各区间点三次多项式系数等。 - 可用 fnplt()绘制出插值结果，其调用格式： fnplt(S) - 对给定的向量xp, 可用fnval()函数计算, 其调用格式： yp=fnval(S,xp) 其中得出的yp是xp上各点的插值结果。 Matlab 样条插值工具箱 鹊碗眯哈秽蒂捌昭拂够赋识狡企辱轮煞忌零躯示萨舆白辜脱万叭蒲陋樟恳matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 例：已知y0=sin(x0),x0= =[0,0.4,1,2,pi], 求该函数 三次样条插值结果 解法: x0=[0,0.4,1,2,pi]; y0=sin(x0); sp=csapi(x0,y0), fnplt(sp, r:); hold on, sp = form: pp breaks: [0 0.4000 1 2 3.1416] coefs: [4x4 double] pieces: 4 order: 4 dim: 1 ezplot(sin(t),[0,pi]); plot(x0,y0,o) 钢株承参胶弱清监魁删冰谚贵吟何熬诈汤潭面父匪公扬逢帖台煤檀频僧葛matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 - 在(0.4000, 1)区间内，插值多项式可以表示为： - 查看插值多项式 sp.coefs ans = -0.1627 0.0076 0.9965 0 -0.1627 -0.1876 0.9245 0.3894 0.0244 -0.4804 0.5238 0.8415 0.0244 -0.4071 -0.3637 0.9093 航滋褒织孽酬闻默保屯违锐鞘认使揣崩氟盗殉兄马奏目耸焚犁炒态又就酶matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 例: 点，用三次样条插值的方法对这些数据进行拟合 解法: x=0:.12:1; y=(x.^2-3*x+5).*exp(-5*x).*sin(x); sp=csapi(x,y); fnplt(sp) 息紊分邹赣食请音占卸诚屎哈埂粮牙炒包蜘惮廖矿臭撒孩扫矢都疚烦箭国matlab在科学计算中的应用06_2matlab在科学计算中的应用06_2 Columns 7 through 12 0.4800 0.6000 -0.2404 0.7652 -0.5776 0.1588 0.6000 0.7200 -0.4774 0.6787 -0.4043 0.1001 0.7200 0.8400 -0.4559 0.5068 -0.2621 0.0605 0.8400 0.9600 -0.4559 0.3427 -0.1601 0.0356 c=[sp.breaks(1:4) sp.breaks(2:5) sp.coefs(1:4,:),... sp.breaks(5:8) sp.breaks(6:9) sp.coefs(5:8,:) ] c = Columns 1 through 6 0 0.1200 24.739