gistfile1.matlab

f = xlsread('C:\Documents and Settings\Administrator\桌面\恒定流-非恒定流关系图.xls')
fun = @(x, xdata)xdata(:, 3).* xdata(:, 4).^ x(1).* xdata(:, 5).^x(2).*xdata(:, 6).^ x(3)* x(4);
lsqcurvefit(fun, rand(4, 1), f, f(:, 2))

%% output
-1.99891623025666
0.627151427841047
-0.0645792624367950
8.20564570528408

f = xlsread('C:\Documents and Settings\Administrator\桌面\恒定流-非恒定流关系图.xls');
logf = log(f);
fun2 = @(x, xdata)xdata(:, 3)+ x(1)* xdata(:, 4)+ x(2)* xdata(:, 5)+ x(3)* xdata(:, 6)+ x(4);
c = [ff(:, 4:6), ones(size(ff, 1), 1)];
d = ff(:, 2) - ff(:, 3);
x3 = lsqlin(c, d, [], []);

% output
-1.69934227572503
1.00672073832549
-0.103933155909792
3.41567959152986