matlab求解方程系数

x1=[1;4];
y1=[0.01;0.8];
p=fittype(\'(1+a1*(x1-b1)^-2)^-1\', \'independent\',\'x1\');%定义函数定义自变量
opt=fitoptions(p);
opt.StartPoint=[0.1 0.1];%设置初始参数
f1=fit(x1,y1,p,opt)