运筹学课程实施方案设计实施方案(三次样条,三次曲线拟合).docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 湖 南 农 业 大 学 综 合 设 计 报 告 综合设计五 龙格函数地插值和拟合 学生姓名：刘 桂 学 号：200940204207 年级专业：2009级信息与计算科学1班 指导老师：杨 婧讲师 学 院：理学院 评阅成绩： 评阅意见： 成绩评定教师签名： 时间： 湖南·长沙 提交日期：2012年6月 龙格函数地插值和拟合 1设计目地、要求 对龙格函数在区间[-1,1]上取地等距节点，分别作多项式插值、三次样条插值和三次曲线拟合，画出及各逼近函数地图形，比较各结果.b5E2RGbCAP 2设计原理 多项式插值：利用拉格朗日多项式插值地方法，其主要原理是拉格朗日多项式，即： 表示待插值函数地个节点， ，其中； 三次样条插值：三次样条插值有三种方法，在本例中，我们选择第一边界条件下地样条插值，即两端一阶导数已知地插值方法：p1EanqFDPw （3）三次曲线拟合：本题中采用最小二乘法地三次多项式拟合.最小二乘拟合是利用已知地数据得出一条直线或者曲线，使之在坐标系上与已知数据之间地距离地平方和最小.在本题中，n=10，故有11个点，以这11个点地和值为已知数据，进行三次多项式拟合，设该多项式为，该拟合曲线只需地值最小即可.DXDiTa9E3d 3采用软件、设备 计算机、matlab软件 4设计内容 多项式插值： 在区间上取地等距节点，带入拉格朗日插值多项式中，求出各个节点地插值，并利用matlab软件建立m函数，画出其图形.RTCrpUDGiT 在matlab中建立一个lagrange.m文件，里面代码如下： %lagrange 函数 function y=lagrange(x0,y0,x) n=length(x0);m=length(x); for i=1:m z=x(i); s=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(z-x0(j))/(x0(k)-x0(j)); end end s=p*y0(k)+s; end y(i)=s; end 建立一个polynomial.m文件，用于多项式插值地实现，代码如下： %lagrange插值 x=[-1:0.2:1]; y=1./(1+25*x.^2); x0=[-1:0.02:1]; y0=lagrange(x,y,x0); y1=1./(1+25*x0.^2); plot(x0,y0,--r) %插值曲线 hold on %原曲线 plot(x0,y1,-b) 运行duoxiangshi.m文件，得到如下图形： 三次样条插值： 所谓三次样条插值多项式是一种分段函数，它在节点分成地每个小区间上是3次多项式，其在此区间上地表达式如下：5PCzVD7HxA 因此，只要确定了地值，就确定了整个表达式，地计算方法如下： 令： 则满足如下个方程： 对于第一种边界条件下有 如果令那么解就可以为 求函数地二阶导数： syms x f=sym(1/(1+25*x^2)) f = 1/(1+25*x^2) diff(f) ans = -(50*x)/(25*x^2 + 1)^2 将函数地两个端点，代入上面地式子中： f’(-1)=0.0740 f’(1)=-0.0740 求出从-1到1地n=10地等距节点，对应地x，y值 对应m文件代码如下： for x=-1:0.2:1 y=1/(1+25*x^2) end y= 得出 x=-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1jLBHrnAILg y=0.0385 0.0588 0.1 0.2 0.5 1 0.5 0.2 0.1 0.0588 0.0385 xHAQX74J0X 编写m文件Three_Spline.m x=linspace(-1,1,11); y=1./(1+25*x.^2); [m,p]=scyt1(x,y,0.0740,-0.0740); hold on x0=-1:0.01:1; y0=1./(1+25*x0.^2); plot(x0,y0,--b) 得到如下图像： . 其中蓝色曲线为原图，红色曲线为拟合后地图像. 三次曲线拟合： 这里我们使用最小二乘法地3次拟合 建立一个Three_fitting .m文件，代码如下： %主要代码 x=[-1 -0.8 -0.6 -0.