《数学建模实验报告(共5页).docx》由会员分享,可在线阅读,更多相关《数学建模实验报告(共5页).docx(5页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上桂林电子科技大学2017-2018学年 第1学期 数学建模 实验名称实验四、插值与拟合检查情况:成绩 教师签名: 院 系计算机与信息安全学院专业计算机科学与技术学 号姓名同 作 者实验日期2017年12月24日一、 实验目的应用matlab处理数据并分析,主要学会并熟练掌握数据拟合和插值。二、 实验内容题目1在化工生产中常常需要知道丙烷在各种温度T和压力P下的导热系数K。下面是实验得到的一组数据:T/C68688787106106140140KP/103kPa9.798113.3249.007813.3559.791814.2779.656312.4630696K0.
2、08480.08970.07620.08070.06960.07530.06110.06510753试求 T=99/C和 P=10.3103 kPa下的K。p2=9.0078,13.355; k2=0.0762,0.0807; %T=87p3=9.7918,14.277; k3=0.0696,0.0753; %T=106a2=polyfit(p2,k2,1); a3=polyfit(p3,k3,1);x1=polyval(a2,10.3); x2=polyval(a3,10.3); %x1,x2 分别是 P=10.3(103 kPa) 下87和 106时的 k 值x=87,106;y=x1,x
3、2;a=polyfit(x,y,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.5z4868688998
4、8949插值法求解:代码:x=129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5;y=7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -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);%等高线图 红色区域为危险区域,所以船只要避免进
5、入figure(2)contour(cx,cy,cz,-5,-5,r)hold onplot(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-
6、3;y0=y+rand;f1=polyfit(x,y0,1)% 输出多项式系数y1=polyval(f1,x);% 计算各 x点的拟合值plot(x,y,+,x,y1)grid ontitle( 一次拟合曲线 );figure(2);f2=polyfit(x,y0,2)%2 次多项式拟合y2=polyval(f2,x);plot(x,y,+,x,y2);grid ontitle( 二次拟合曲线 );figure(3);f4=polyfit(x,y0,4)%4 次多项式拟合y3=polyval(f4,x);plot(x,y,+,x,y3)grid ontitle( 四次拟合曲线 );figure
7、(4);f6=polyfit(x,y0,6)%6 次多项式拟合y4=polyval(f6,x);plot(x,y,+,x,y4)grid ontitle( 六次拟合曲线 );运行结果如下:依次为各个拟合曲线的系数(按降幂排列)f1 =43.2000 -149.0663f2 = 10.5000 -72.3000 89.8087f4 =0.0000 1.0000 -6.0000 5.0000 -2.5913f6 = 0.0000 -0.0000 0.0000 1.0000 -6.0000 5.0000-2.4199运行后, 比较拟合后多项式和原式的系数, 发现四次多项式系数与原系数比较接近, 四次多项式的四次项系数很小。 作图后, 发现一次和二次多项式的图形与原函数的差别比较大, 属于欠拟合的情况, 而四次多项式和六次多项式符合得比较好。作图如下:三、 实验总结本次实验,主要学习和实践了插值和拟合这两个在实际应用中十分广泛的建模方法,同时,它们也是我们学习数学建模的基础。插值和拟合,基础但不简单,需要大量的训练才能很好的进行掌握。专心-专注-专业