资源描述：
《数学建模实验报告.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、桂林电子科技大学2017-2018学年第1学期数学建模实验名称实验四、插值与拟合检查情况：成绩教师签名：院系计算机与信息安全学院专业计算机科学与技术学号姓名同作者实验日期2017年12月24日一、实验目的应用matlab处理数据并分析，主要学会并熟练掌握数据拟合和插值。二、实验内容题目1在化工生产中常常需要知道丙烷在各种温度T和压力P下的导热系数K。下面是实验得到的一组数据：T/°C68688787106106140140KP/103kPa9.798113.3249.007813.3559.791814.2779.656312.4630696K0.08480.08970.07620.08
2、070.06960.07530.06110.06510753试求T=99/°C和P=10.3×103kPa下的K。p2=[9.0078,13.355];k2=[0.0762,0.0807];%T=87℃p3=[9.7918,14.277];k3=[0.0696,0.0753];%T=106℃a2=polyfit(p2,k2,1);a3=polyfit(p3,k3,1);x1=polyval(a2,10.3);x2=polyval(a3,10.3);%x1,x2分别是P=10.3（103kPa）下87℃和106℃时的k值x=[87,106];y=[x1,x2];a=polyfit(x,y,
3、1);%求解多项式的系数z=polyval(a,99)%求解多项式z在x=99下对应的值。实验结果：z=0.0729题目2下表给出了某一海域以码为单位的直角坐标Oxy上一点（x，y）（水面一点）以英尺为单位的水深z，水深数据是在低潮时测得的，船的吃水深为5英尺，问在矩形区域（75，200）x（-50，150）里那些地方船要避免进入。x129140103.588185.5195105.5157.5107.57781162162117.5y7.5141.52314722.5137.585.5-6.5-81356.5-66.584-33.5z48686889988949插值法求解：代码：x=[
4、129140103.588185.5195105.5157.5107.57781162162117.5];y=[7.5141.52314722.5137.585.5-6.5-81356.5-66.584-33.5];z=[-4-8-6-8-6-8-8-9-9-8-8-9-4-9];cx=75:0.5:200;cy=-50:0.5:150;[cx,cy]=meshgrid(cx,cy);cz=griddata(x,y,z,cx,cy,'cubic');%等高线图红色区域为危险区域，所以船只要避免进入figure(2)contour(cx,cy,cz,[-5,-5],'r')holdonpl
5、ot(x,y,'*')xlabel('X'),ylabel('Y')实验结果：船的吃水深度为5英尺时，红色区域为危险区域，船要避免进入。题目3用给定的多项式，如y=x3-6x2+5x-3，产生一组数据(xi,yi，i=1,2,…,n),再在yi上添加随机干扰(可用rand产生(0,1)均匀分布随机数,或用rands产生N(0,1)分布随机数)，然后用xi和添加了随机干扰的yi作的3次多项式拟合，与原系数比较。如果作2或4次多项式拟合，结果如何？线性最小二乘拟合程序如下：x=1:0.5:10;y=x.^3-6*x.^2+5*x-3;y0=y+rand;f1=polyfit(x,y0,1)%
6、输出多项式系数y1=polyval(f1,x);%计算各x点的拟合值plot(x,y,'+',x,y1)gridontitle('一次拟合曲线');figure(2);f2=polyfit(x,y0,2)%2次多项式拟合y2=polyval(f2,x);plot(x,y,'+',x,y2);gridontitle('二次拟合曲线');figure(3);f4=polyfit(x,y0,4)%4次多项式拟合y3=polyval(f4,x);plot(x,y,'+',x,y3)gridontitle('四次拟合曲线');figure(4);f6=polyfit(x,y0,6)%6次多项式拟合
7、y4=polyval(f6,x);plot(x,y,'+',x,y4)gridontitle('六次拟合曲线');运行结果如下：依次为各个拟合曲线的系数（按降幂排列）f1=43.2000-149.0663f2=10.5000-72.300089.8087f4=0.00001.0000-6.00005.0000-2.5913f6=0.0000-0.00000.00001.0000-6.00005.0000-2.4199运行后，比较拟合
