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