Question

求大神，对函数1/y=a/(x+3)+b求最小二乘拟合的matlab程序，其中a、b为系数，x、y为已知列向量。谢谢啦！

Answer 1

clear;

x = [1  2  3  4  5  6  7  8  9 10];

y = [1.5377  1.0377  0.8710  0.7877  0.7377  0.7043  0.6805  0.6627  0.6488  0.6377];

fun = @(a,x)a(1)./(x+3)+a(2) ; % 拟合函数的形式

a0=[1,1]; % 猜测的初始参数

a = nlinfit(x,1./y,fun,a0)

plot(x,y,'b*',x,1./(a(1)./(x+3)+a(2)),'r')

a =

   -5.2213    1.9927