ex3-9-3_2.m

来自「matlab程序」· M 代码 · 共 30 行

% 用6阶多项式对区间[0,2.5]上的误差函数进行最小二乘拟合

%误差函数:
% 	y = erf(x).
% 	y = 2/sqrt(pi) * integral from 0 to x of exp(-t^2) dt.

x=0:0.1:2.5;

y=erf(x);

P=polyfit(x,y,6)  %求拟合多项式

Px=poly2sym(round(P))


%有效拟合的区间性图示

x=0:0.1:5;

y=erf(x);

f=polyval(P,x);     %利用拟合多项式求值

plot(x,y,'bo',x,f,'r-')

axis([0,5,0,2])

legend('原始数据','拟合曲线')