测控系统原理与设计实验报告--测试数据的拟合与插值.docx

《测控系统原理与设计实验报告--测试数据的拟合与插值.docx》由会员分享，可在线阅读，更多相关《测控系统原理与设计实验报告--测试数据的拟合与插值.docx（5页珍藏版）》请在第一文库网上搜索。

1、实验报告（理工类）开课学院及实验室：实验时间：年月日学生姓名学号成绩学生所在学院年级/专业/班12级/自动化/西华理T课程名称测控系统原理与设计课程代码实验项目名称测试数据的拟合与插值项目代码指导教师项目学分-实验目的学习利用MAT1AB编程实现拟合法以及求误差。二实验内容1拟合法：（工）po1yfit函数T1AB的po1yfit函数用于多项式拟合，其语法为：P=po1fit（xzPk）;其中，X,y分别是横纵坐标向量，它们不仅元素个数相同，而且同为行向量或同为列向量。k为非负整数，是待拟合的多项式的最高次数。P是输出项，为待拟合的多项式的系数向量（由高次到低次排列）。在MAT1AB中输入以下

2、代码：X=10,11,12,13;y=8.0000,7.84317.6923,7.5471;P=po1yfit（x,y,1）；运行结果：P=-0.15099.5065所以拟合得的函数就是：y=-0.1509X+9.5065.在进行多项式拟合时，必须注意的是，拟合的精度是有限的，一般而言，需要满足以下条件:记m为不重复的横坐标的数目，则拟合次数kv=m1,在此前提下尽量使用低次多项式进行拟合。(2)po1yva1函数po1yva1,顾名思义就是“多项式的值”，该函数的功能是将已知数据代入拟合得的多项式求值。语法格式：y=poyva1(pzx);其中，P是已经拟合的多项式(比如说()中的P),X是

3、自变量组成的向量，y是所求值组成的向量。输入以下代码：X=10,11,12;y=8.0000,7.8431,7.6923;P=po1yfit(xy,1)t=1,2,3,4,5;s=po1yva1(p,t)运行结果：-0.15389.53758.76829.38369.22989.07598.9221(3)计算多项式拟合的方差已知原始数据X和y,拟合得到多项式p,判断拟合效果好坏的一个重要指标是方差，方差的计算方法是e=sum(y-po1yva1(p,x).2).PoIyVaI(P,X)得到拟合值向量，y是真实值向量，两者相减得到真实值和拟合值的差值向量，“人2”表示对矩阵中的每一个元素进行平方

4、运算，于是得到差值向量中每一个元素的平方，SUm是求和函数，显然就是求差值向量元素的平方和，而这就是方差。输入以下代码：X=10,11,12;y=8.0000,7.84317.6923;P=po1yfit(x,y,1)e=sum(y-po1yva1(pAx).a2)运行结果：-0.15389.5375e=6.2017e-006可见，方差能0,表明拟合效果很好。2插值法：根据教材P149之表451的数据，取前十个数据可得。输入以下代码：t=10:1:19;p=8.0007.84317.69237.54717.40747.27277.14287.01746.89656.7796;x=10:19;y

5、1=iterp1(t3p5x);%默认为1ineary2=iterp1(t,p5x5nearest);y3=iterp1(t,p5x5sp1ine);y4=iterp1(t,p5xjpchip);y5=iterp1(t,p5xjcubic);subp1ot(253,1);Pot(5y1515p,ro);tit1e(1iear);subp1ot(253,2);pot(x5y2,t5p5ro);tit1e(earest);subp1ot(253,3);pot(x5y3,t5p5ro);tit1e(sp1ie);subp1ot(253,4);pot(x5y4,t5p5ro);tit1e(pchip);subp1ot(253,5);pot(x5y5,t5p5ro);tit1e(cubic);运行结果:三实验总结