MATLAB第五次实验

1、一、实验目的学习利用matlab作拟合和插值二、实验任务P130 8.在钢线碳含量对于电阻的效应的研究中，得到以下数据，分别用一次、三次、五次多项式拟合曲线来拟合这组数据并画出图形。碳含量x0.100.300.400.550.700.800.95电阻y1518192122.623.826P130 10.在某种添加剂的不同浓度之下对铝合金进行抗拉强度试验，得到数据如下，现分别使用不同的插值方法，对其中间没有测量的浓度进行推测，并估算出浓度X=18与26时的抗压强度丫的值。浓度X1015202530抗压强度Y25.229.831.231.729.42 2P130 12.用不同方法对z 上的二维插值

2、效果进行比拟。三、实验程序P130 8.x=0.10 0.30 0.40 0.55 0.70 0.80 0.95;y=15 18 19 21 22.6 23.8 26;p仁polyfit(x,y,1);p3=polyfit(x,y,3);p5=polyfit(x,y,5);disp(' 一阶拟合函数'),f1 = poly2str(p1, 'x') disp('三阶拟合函数'),f3 = poly2str(p3, 'x') disp('五阶拟合函数'),f5 = poly2str(p5, 'x')

3、x1=0:0.01:0.95;y1=polyval(p1,x1);y3=polyval(p3,x1);五次拟合')y5=polyval(p5,x1); plot(x,y,'rp',x1,y1,'-',x1,y3,'k-',x1,y5,'g-') legend('拟合点','一次拟合','三次拟合',P130 10.x = 10 :5 : 30;%浓度 xy =25.2, 29.8, 31.2, 31.7, 29.4;%抗压强度p4=polyfit(x,y,4);x1=10 :

4、 0.1 : 30;y1=in terp1(x, y, x1, '*n earest');%最近点插值y2=in terp1(x, y, x1, '*li near');%线性插值y3=in terp1(x, y, x1, '*spli ne');%样条插值y4=in terp1(x, y, x1, '*cubic');%立方插值plot(x, y, 'ro', x1, y1, '-', x1, y2, '-', x1, y3, 'k-.', x1, y4, '

5、;m:')legend('原始数据','最近点插值','线性插值','样条插值','立方插值')P130 12.给出数据点x, y = meshgrid(-3 : 0.5 : 3); z = x.A2./16 - y.A2./9;%x1, y1 = meshgrid(-3 : 0.1 : 3); z1 = x1.A2./16 - y1.A2./9; figure(1)subplot(1,2, 1), mesh(x, y, z), title('数据点')subplot(1,2, 2), me

6、sh(x1, y1, z1), title('数据图像')%确定插值点%最近点插值%线性插值%样条插值%立方插值绘制使用4种方法得到的图形 最近点插值') 线性插值') 样条插值') 立方插值')xi, yi = meshgrid(-3 : 0.15 : 3); zi1 = in terp2(x, y, z, xi, yi, '*n earest'); zi2 = in terp2(x, y, z, xi, yi, '*li near'); zi3 = in terp2(x, y, z, xi, yi, '*spli ne'); zi4 = in terp2(x, y, z, xi, yi, '*cubic'); figure(2) %翻开另一个图形窗口，subplot(2, 2, 1), mesh(xi, yi, zi1), title(' subplot(2, 2, 2), mesh(xi, yi, zi2), title(' subplot(2, 2, 3), mes