matlab拟合参数的代码,matlab 对方程参数的拟合

Constant c =2*10^(-5);

ConstStr H= (sqrt(-(4*k1^3*k2^2*L^4+(-32*c*k1^2*k2^3+(8*k1^2-12*c*k1^3)*k2^2+8*k1^3*k2)*L^3+(48*c^2*k1^2*k2^3+(13*c^2*k1^3-100*c*k1^2+4*k1)*k2^2+(20*k1^2-22*c*k1^3)*k2+k1^3)*L^2+(-24*c^3*k1^2*k2^3+(-6*c^3*k1^3+112*c^2*k1^2-40*c*k1)*k2^2+(20*c^2*k1^3-78*c*k1^2+16*k1)*k2-2*c*k1^3+2*k1^2)*L+4*c^4*k1^2*k2^3+(c^4*k1^3-32*c^3*k1^2-8*c^2*k1)*k2^2+(-6*c^3*k1^3+46*c^2*k1^2-32*c*k1+4)*k2+c^2*k1^3-6*c*k1^2+k1)/k1)/(2*3^(3/2)*k1*k2^2)+(k1*(-16*k2^3*L^3-42*k2^2*L^2-21*k2*L-2)+c*(k1*(24*k2^3*L^2+51*k2^2*L+15*k2)-18*k2^2)+c^2*k1*(-12*k2^3*L-15*k2^2)-18*k2^2*L+2*c^3*k1*k2^3-9*k2)/(54*k1*k2^3))^(1/3)+(k1*(4*k2^2*L^2+7*k2*L+1)+c*k1*(-4*k2^2*L-5*k2)+c^2*k1*k2^2+3*k2)/(9*k1*k2^2*(sqrt(-(4*k1^3*k2^2*L^4+(-32*c*k1^2*k2^3+(8*k1^2-12*c*k1^3)*k2^2+8*k1^3*k2)*L^3+(48*c^2*k1^2*k2^3+(13*c^2*k1^3-100*c*k1^2+4*k1)*k2^2+(20*k1^2-22*c*k1^3)*k2+k1^3)*L^2+(-24*c^3*k1^2*k2^3+(-6*c^3*k1^3+112*c^2*k1^2-40*c*k1)*k2^2+(20*c^2*k1^3-78*c*k1^2+16*k1)*k2-2*c*k1^3+2*k1^2)*L+4*c^4*k1^2*k2^3+(c^4*k1^3-32*c^3*k1^2-8*c^2*k1)*k2^2+(-6*c^3*k1^3+46*c^2*k1^2-32*c*k1+4)*k2+c^2*k1^3-6*c*k1^2+k1)/k1)/(2*3^(3/2)*k1*k2^2)+(k1*(-16*k2^3*L^3-42*k2^2*L^2-21*k2*L-2)+c*(k1*(24*k2^3*L^2+51*k2^2*L+15*k2)-18*k2^2)+c^2*k1*(-12*k2^3*L-15*k2^2)-18*k2^2*L+2*c^3*k1*k2^3-9*k2)/(54*k1*k2^3))^(1/3))+(-2*k2*L+c*k2-1)/(3*k2);

Variable L, A;

Function A=(L*(e1*k1*H+e2*k1*k2*H^2))/(1+k1*H+k1*k2*H^2);

Data;

//% L            A

5.00E-07        0.01

1.00E-06        0.02

1.50E-06        0.032

2.00E-06        0.042

2.50E-06        0.05

3.00E-06        0.061

3.50E-06        0.07

4.00E-06        0.08

4.50E-06        0.09

5.00E-06        0.1

5.50E-06        0.11

6.00E-06        0.119

6.50E-06        0.128

7.00E-06        0.136

7.50E-06        0.148

8.00E-06        0.157

8.50E-06        0.166

9.00E-06        0.174

9.50E-06        0.178

1.00E-05        0.188

1.05E-05        0.195

1.10E-05        0.2

1.15E-05        0.2

1.20E-05        0.203

1.25E-05        0.204

1.30E-05        0.206

1.35E-05        0.205

1.40E-05        0.206

1.45E-05        0.208

1.50E-05        0.21

1.55E-05        0.211

1.60E-05        0.212

1.65E-05        0.213

1.70E-05        0.213

1.75E-05        0.214

1.80E-05        0.214

1.85E-05        0.215