资源描述:
《北京科技大学数学实验第五次教学内容.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、…………………………………………………………最新精品资料推荐……………………………………………………《数学实验》报告实验名称 Matlab拟合与插值2013年12月…………………………………………………………最新精品资料推荐……………………………………………………11…………………………………………………………最新精品资料推荐……………………………………………………一、【实验目的】1.学习Matlab的一些基础知识,主要多项式及其相关计算等;2.熟悉Matlab中多项式的拟合,编写一些相关的Matlab命令等;3.熟悉
2、Matlab中多项式的插值,并编写一些相关的Matlab命令等;4.完成相关的练习题。二、【实验任务】1.在钢线碳含量对于电阻的效应的研究中,得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形,计算当x=0.45时的电阻值.碳含量x0.100.300.400.550.700.800.95电阻y1518192122.623.8262.在某种添加剂的不同浓度之下对铝合金进行抗拉强度试验,得到数据如下,现分别使用不同的插值方法,对其中没有测量的浓度进行推测,并估算出浓度X=18及26时的抗压强度Y的值
3、.浓度X1015202530抗压强度Y25.229.831.231.729.43.用不同方法对在(-3,3)上的二维插值效果进行比较.…………………………………………………………最新精品资料推荐……………………………………………………11…………………………………………………………最新精品资料推荐……………………………………………………三、【实验程序】1.在钢线碳含量对于电阻的效应的研究中,得到以下数据.分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形,计算当x=0.45时的电阻值.碳含量x0.100.30
4、0.400.550.700.800.95电阻y1518192122.623.826M文件clc;clf;x=[0.10.30.40.550.70.80.95];y=[1518192122.623.826];p1=polyfit(x,y,1);p3=polyfit(x,y,3);p5=polyfit(x,y,5);x1=0.1:0.05:1;y1=polyval(p1,x1);y3=polyval(p3,x1);y5=polyval(p5,x1);plot(x,y,'rp',x1,y1,'b-',x1,y3,'g-.',
5、x1,y5,'m--');legend('拟合点','一次拟合','三次拟合','五次拟合');disp('以下为当x=0.45时的电阻值:')disp('一阶拟合函数值'),g1=polyval(p1,0.45)disp('三阶拟合函数值'),g3=polyval(p3,0.45)disp('五阶拟合函数值'),g5=polyval(p5,0.45)…………………………………………………………最新精品资料推荐……………………………………………………11…………………………………………………………最新精品资料推荐………
6、……………………………………………2.在某种添加剂的不同浓度之下对铝合金进行抗拉强度试验,得到数据如下,现分别使用不同的插值方法,对其中没有测量的浓度进行推测,并估算出浓度X=18及26时的抗压强度Y的值.浓度X1015202530抗压强度Y25.229.831.231.729.4M文件clc;clf;x=[1015202530];y=[25.229.831.231.729.4];xi=10:0.05:30;yi1=interp1(x,y,xi,'*nearest');yi2=interp1(x,y,xi,'*line
7、ar');yi3=interp1(x,y,xi,'*spline');yi4=interp1(x,y,xi,'*cubic');plot(x,y,'rp',xi,yi1,'g-',xi,yi2,'c-.',xi,yi3,'m--',xi,yi4,'b:')gridon;legend('原始数据','最近点插值','线性插值','样条插值','立方差值')…………………………………………………………最新精品资料推荐……………………………………………………11…………………………………………………………最新精品资料推荐……
8、………………………………………………3.用不同方法对在(-3,3)上的二维插值效果进行比较.M文件clf;t=-3:1:3;[x,y]=meshgrid(t);z=x.^2/16-y.^2/9;t1=-3:0.5:3;[x1,y1]=meshgrid(t1);z1=x1.^2/16-y1.^2/9;[xi,yi]=meshgrid