3-36.m

来自「《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计」· M 代码 · 共 16 行

x = (0: 0.1: 2.5)';
% 给出一组数据，为误差函数的一个区间
y = erf(x); 
p = polyfit(x,y,6)
% 计算该区间内6阶拟合多项式的向量

x1 = (0: 0.1: 5)';
% 将区间增长一倍
y1 = erf(x1);
% 计算误差函数在新区间内的函数值
f = polyval(p,x1);
% 计算6阶拟合曲线在新区间内的取值
plot(x1,y1,'o',x1,f,'-')
% 绘图，比较前后区间内的曲线拟合效果，如图3-7
axis([0  5  0  2])