这些数据如何用matlab进行非线性拟合

五十多分

x1=[29.31 25.84 25.66 25.49 25.4 25.35 25.3 22.34 20.97 18.79 18.67 18.56 18.51 18.48 18.44 14.81 14.04 13.33 12.05 11.47 25.72 23.03 22.77 22.53 22.41 22.35 22.28 18.71 16.98 16.82 16.66 16.58 16.54 16.5 22.72 16.78 19.09 14.39];
>> x2=[5.69 1.51 2.5 4.96 7.42 9.87 14.78 4.36 4.1 1.1 1.83 3.64 5.44 7.24 10.84 2.91 2.76 2.63 2.38 2.27 5.01 1.35 2.22 4.4 6.56 8.72 13.04 3.66 1 1.65 3.27 4.88 6.49 9.71 4.43 3.29 3.74 2.83];
>> x3=[2.24 1.23 1.64 2.24 2.6 2.87 3.23 2.24 2.24 1.23 1.64 2.24 2.6 2.87 3.23 2.24 2.24 2.24 2.24 2.24 3.72 2.08 2.75 3.72 4.33 4.76 5.35 3.72 2.08 2.75 3.72 4.33 4.76 5.35 4.73 7.44 7.44 7.44];
>> x4=[0.00262 0.0021 0.0021 0.0021 0.0021 0.0021 0.0021 0.00152 0.00116 0.00063 0.00063 0.00063 0.00063 0.00063 0.00063 0.00034 0.00034 0.00037 0.00046 0.00052 0.00262 0.0021 0.0021 0.0021 0.0021 0.0021 0.0021 0.00088 0.00051 0.00051 0.00051 0.00051 0.00051 0.00051 0.00262 0.00077 0.00193 0.00064];
>> y=[11.9 11.63 10.91 10.13 9.75 9.5 9.19 7.51 6.54 6.4 6 5.57 5.36 5.22 5.05 6.2 6.77 7.2 7.11 5.95 15.96 14.42 13.54 12.59 12.11 11.81 11.42 9.6 9.92 9.31 8.66 8.33 8.12 7.86 18.09 11.3 22.46 16.43];
>>y1=log(y);
>> f=inline('a(1)-a(2).*log(x1)+a(3).*log(x2)+a(4).*log(x3)+log(x4)','a','x1','x2','x3','x4');