f1=fittype('a*(1-b*y1.^2).*y2+c*y1+d/2*sin(e*t+f)','problem',{'y1','y2'},'coefficients',{'a','b','c','d','e','f'},'independent','t'); %定义拟合的方程
[FO,G,O]=fit(xdata,ydata,f1,'problem',{y1,y2},'StartPoint',[1,1,1,1,1,1]); %最小二乘拟合,也可以用'Lower'和'Upper'设置参数的范围
para=coeffvalues(FO); %获取参数a,b,c,d,e,f
a=para(1); %b,c,d,e,f依此类推
