4楼: Originally posted by yanshima at 2012-12-01 08:39:01
这个往下再咋做呢？我用lsqcurvefit(ydata,p,xdata,y)，可是说我是错误的。能麻烦您写详细些嘛。谢谢...

Variable sf,x,y;
Function y=a/b*sqrt((1-exp(-x/a))^2-(1-exp(-sf/a))^2);
Data;
sf           x               y
0        0.209735        11.0164
0        0.473382        29.377
0        1.10326        64.5246
0.15993        1.25993        70.8197
0.26454        1.36454        75.541
0.36882        1.46882        79.2131
0.52481        1.62481        83.4098
0.68029        1.78029        86.0328
0.78457        1.88457        89.7049
0.94022        2.04022        92.8525
1.14776        2.24776        97.0492
1.3553        2.4553        101.246
1.56284        2.66284        105.443
1.7702        2.8702        109.115
2.02895        3.12895        112.262
2.23632        3.33632        115.934
2.44369        3.54369        119.607
2.70226        3.80226        122.23
2.90946        4.00946        125.377
3.11633        4.21633        127.475
3.32336        4.42336        130.098
3.63348        4.73348        132.721
3.94343        5.04343        134.82
4.20218        5.30218        137.967
4.6154        5.7154        140.59
5.02827        6.12827        142.164
5.4927        6.5927        143.738
5.85369        6.95369        144.262

p=[a,b];
xdata=[x,sf]
ydata=p(1)/p(2)*sqrt((1-exp(-xdata(:,1)./p(1))).^2-(1-exp(-xdata(:,2)./p(1))).^2)