内容发布更新时间 : 2024/11/15 4:34:04星期一 下面是文章的全部内容请认真阅读。

第七章

9.已知一组数据：

xi yi 试用解：

-1 10 0 9 1 7 2 5 3 4 4 3 5 0 6 -1 y?ax?b来拟合这组数据（写出matlab程）。

x=-1:6;

y=[10,9,7,5,4,3,0，-1]; p=polyfit(x,y,1); a=p(1);c=exp(p(2));

axy?ce10.对于以下数据表，试写出用来拟合这组数据的matlab程序。

xi 0 1.5 1 2.5 2 3.5 3 5.0 4 7.5 yi 解：

lny?lnc?ax

Y?lny,b?lnc,X?x

x=0:4; y=[1.5,2.5,3.5,5.0,7.5];Y=logy;X=x；

令

p=polyfit(X,Y,1);a=p(1);c=exp(p(2));

i xi yi 0 0 1.0000 1 0.25 1.2840 2 0.5 1.6487 给定数据表 3 0.75 2.1170 4 1.00 2.7183 5 1.25 3.4903 6 1.5 4.4817 试求出3次，4次多项式的曲线拟合，画出计算曲线。 具体步骤： 三次拟合： x=0:0.25:1.5；

y=[1.0000, 1.2840,1.6487 ,2.1170 , 2.7183, 3.4903, 4.4817 ]； plot(x,y,’r*’) p=polyfit(x,y,3) x1=0:0.1:1.5; y1=polyval(p,x1); hold on

plot(x1,y1,’g’)

过程中得出的三次拟合多项式的系数矩阵为p=[ 0.3659 0.2891 1.0642 0.9982] 曲线拟合图：

4.543.532.521.510.500.511.5

四次拟合： x=0:0.25:1.5；

y=[1.0000, 1.2840,1.6487 ,2.1170 , 2.7183, 3.4903, 4.4817 ]； plot(x,y,’r*’) p=polyfit(x,y,4) x1=0:0.1:1.5; y1=polyval(p,x1); hold on

plot(x1,y1,’g’)

曲线拟合图：

4.543.532.521.5100.511.5

过程中得出的三次拟合多项式的系数矩阵为p=[0.0912 0.0924 0.5422 0.9922 1.0001]

2.由实验给出如下实验表： X 0.0 0.1 Y 1.0 0.41 0.2 0.50 0.3 0.61 0.5 0.91 0.8 2.02 1.0 2.46 试求出3次，4次多项式的曲线拟合，画出计算曲线。 具体步骤： 三次拟合：

x=[0.0 ,0.1 ,0.2 ,0.3 ,0.5, 0.8 ,1.0]

y=[1.0 ,0.41 ,0.50 ,0.61 ,0.91, 2.02 ,2.46] plot(x,y,’r*’) p=polyfit(x,y,3) x1=0:0.01:1; y1=polyval(p,x1); hold on plot(x1,y1,’g’)

过程中得出的三次拟合多项式的系数矩阵为p=[-6.6221 12.8147 -4.6591 0.9266] 曲线拟合图：