curvefit04.m

数值分析与计算的PPT幻灯片与MAT源码

%curvefit04.m     已知经验公式的拟合方法

x0=[1.00 1.25 1.50 1.75 2.00];
y0=[1.629 1.756 1.876 2.008 2.135];
plot(x0,y0,'or')                 %画出离散点
pause
n=1;
a=polyfit(x0,y0,n);              %拟合多项式
p=poly2sym(a)                    %求出多项式的表达式
xi=1:0.01:2;                     %细分区间
yi=polyval(a,xi);                %求多项式的值
hold on
plot(xi,yi,'-b','MarkerSize',20) %绘制拟合曲线图
xlabel('x') 
legend('离散点','拟合曲线')
B=a(1),A=a(2)                    %求出A,B的值           
a=exp(A),b=B                     %转化成a、b的原值
a=vpa(a),b=vpa(b)
x=sym('x');                      %定义符号函数
p1=a*exp(b*x)                    %写出经验公式