Question

matlab做最小二乘拟合

基函数已经给定了，m(1) ,m(2),m(3)为系数，
x=[15,30,45,60,75,90,105,120,135,150,165,180,195,210,225,240,255,270];
y=[7,7.9,8.8,10,11,12.6,14,14.6,15.2,15.7,16.2,16.4,16.7,17,17.3,17.5,17.7,18.0];
其中a,b为一个通过人为给定的数，例如令a=3，b=0.05

f=m(1)*(1-exp(a*x).*erfc(sqrt(a*x)))+m(2)*(1-exp(x*b).*erfc(sqrt(x*b)))+m(3)*2*exp(a*x).*sqrt(a*x).*(1/sqrt(pi)*exp(-a*x)-sqrt(a*x).*erfc(sqrt(a*x)))

怎么写程序啊

Answer 1

程序为：

x=[15,30,45,60,75,90,105,120,135,150,165,180,195,210,225,240,255,270];

y=[7,7.9,8.8,10,11,12.6,14,14.6,15.2,15.7,16.2,16.4,16.7,17,17.3,17.5,17.7,18.0];

p=polyfit(x,y,2)

x2=0:300;

y2=polyval(p,x2);

plot(x,y,'o',x2,y2)

结果为：

p =

  -0.000193441119138   0.098295493636051   5.181372549019607

p所得三个值就是对应你要求的m(1) ,m(2),m(3)三个系数

追问

polyfit是多项式拟合啊，这个有规定的拟合函数，怎么解决啊

追答

看你要拟合成什么函数了。用这个语句cftool(x,y)