我用lsqcurvefit函数和自定义的多元参数的函数拟合实验数据，但是得到的拟合曲线总是不太理想，请大侠帮偶看看是咋回事呢。

实验数据见上传的附件。

目的是拟合实验数据中第一列定义为Xdata，第六列定义为Ydata，要利用自定义的函数拟合该数据。

命令行如下：
clear
load Kopie_von_B282A070830_temp.TXT; % load experimental data
DRS=Kopie_von_B282A070830_temp;

xdata=DRS(:,1); % frequency
ydata=DRS(:,6); % eps''

% fit function
% eps2=-imag(x(1)./(1.0+(i*2*pi*xdata*x(4)).^x(2)).^x(3)+x(6)./
%(1.0+(i*2*pi*xdata*x(9)).^x(7)).^x(8))+x(5)+x(10)*1e-14./(8.8542e-12*2*pi*xdata);
eps2=inline('-imag(x(1)./(1.0+(i*2*pi*xdata*x(4)).^x(2)).^x(3)+x(6)./(1.0+(i*2*pi*xdata*x(9)).^x(7)).^x(8))+x(5)+x(10)*1e-14./(8.8542e-12*2*pi*xdata)','x','xdata');

% define parameters
deps1 = 6;
alpha1 = 0.65;
beta1 = 0.6;
tau1 = 2.8e-7;
eps2off = 0.1;

deps2 = 2.8;
alpha2 = 0.5;
beta2 = 0.5;
tau2 = 1;
sigmadc = 1.5e2;

x0=[deps1,alpha1,beta1,tau1,eps2off,deps2,alpha2,beta2,tau2,sigmadc];
% matrix of the initial value of parameters

opt=optimset('TolFun',1e-5);
[x,resnorm,residual]=lsqcurvefit(eps2,x0,xdata,ydata,[],[],opt);

eps2=-imag(x(1)./(1.0+(i*2*pi*xdata*x(4)).^x(2)).^x(3)+x(6)./...
(1.0+(i*2*pi*xdata*x(9)).^x(7)).^x(8))+x(5)+x(10)*1e-14./(8.8542e-12*2*pi*xdata);

loglog(xdata,ydata,'o',xdata,eps2,'or');

hold on

说明下，上述命令中用的初始值x0，别人用gnuplot拟合比我的效果好多了，不解ing

kCHtswNw.txt (13.45 KB) [求助]我的拟合为什么不好呢？

[此贴子已经被作者于2007-9-22 0:18:59编辑过]

谢谢版主的解答。


我所用的拟合方程是一个模型的方程，实验数据应该是符合这个方程所描述的模型的。另外我感到不解的是，为什么别人用gnuplot这个软件拟合的效果比我的好多了（你和出来的曲线几乎和实验曲线完全重合，且他所用的拟合方程和我用的完全一样）。