matlab拟合曲线出现虚数

fangzhensgz

clc;clear
x=[1        2        3        4        5        6        7        8        9        10        11        12        13        14        15        16        17        18        19        20        21        22        23        24        25        26        27        28        29        30        31        32        33        34        35        36        37        38        39        40        41        42        43        44        45        46        47        48        49        50        51        52        53        54        55        56        57        58        59        60
]
y=[1.9749E-05        1.8547E-05        1.99946E-05        2.45418E-05        2.57757E-05        2.73701E-05        2.95021E-05        3.10504E-05        3.24112E-05        3.37821E-05        3.48578E-05        3.67454E-05        3.70243E-05        3.74074E-05        3.77932E-05        3.78787E-05        3.80667E-05        3.81555E-05        3.88504E-05        3.91431E-05        3.95376E-05        3.97312E-05        3.96231E-05        3.96164E-05        3.96102E-05        3.98062E-05        4.00025E-05        3.98968E-05        4.01947E-05        4.05936E-05        4.06905E-05        4.06867E-05        4.09857E-05        4.09823E-05        4.108E-05        4.12787E-05        4.13767E-05        4.13741E-05        4.14724E-05        4.147E-05        4.15686E-05        4.1768E-05        4.1766E-05        4.16632E-05        4.17621E-05        4.17604E-05        4.19602E-05        4.20594E-05        4.22594E-05        4.22579E-05        4.25589E-05        4.25575E-05        4.28586E-05        4.30589E-05        4.32592E-05        4.33588E-05        4.30553E-05        4.35582E-05        4.35571E-05        4.39593E-05
]
fx =inline('(1./(63000.*(x./(x+0.47)).^0.5))+1./(a(1).*(x/(x+a(2)).^0.5)).*(1-exp((-(a(1).*(x./(x+a(2)).^0.5).*x./(a(3).*x)./(x+a(4))).^0.5)))','a','x')
a=nlinfit(x,y,fx,[1 1 1 1])

fire051

拟合过程中导致平方根里面出现了负数。可能是拟合方程与数据不太匹配。