syms f_1 f_2 r_1 r_2
G=(f_1*f_2+r_2*f_2^2)/(r_1*f_1^2+2*f_1*f_2+r_2*f_2^2);
Gdr_1=f_1^2*(f_1*f_2 + r_2*f_2^2)/(2*f_1*f_2 + f_1^2*r_1 + f_2^2*r_2)^2;
Gdr_2=f_2^2*(f_1*f_2 + r_2*f_1^2)/(2*f_1*f_2 + f_1^2*r_1 + f_2^2*r_2)^2;
%定义符号表达式
f_11=[0.925 0.846 0.763 0.674 0.579 0.479 0.256 0.133];
f_22=[0.075 0.154 0.237 0.326 0.421 0.521 0.744 0.867];
F_22=[0.441 0.699 0.753 0.828 0.864 0.900 0.968 0.983];
r_11=0.09;r_22=10;
%准备给符号表达式赋值
Gdr_1=subs(Gdr_1,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
Gdr_2=subs(Gdr_2,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
G=subs(G,{'f_1','f_2','r_1','r_2'},{f_11,f_22,r_11,r_22});
d=F_22-G;
myfun=@(x,xdata)xdata(1,:)*x(1)+xdata(2,:)*x(2);
xdata=[Gdr_1;Gdr_2];
x0=[0.1 10];
x = lsqcurvefit(@(x,xdata) myfun(x,xdata),x0,xdata,d)