24小时热门版块排行榜

返回列表

piao2004ha

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 98.8
散金: 32
帖子: 143
在线: 69.1小时
虫号: 398719
注册: 2007-06-11
性别: GG
专业: 材料

[交流] 【求助】紧急求助非线性拟合问题已有3人参与

数据是：tdata ydata
2.3 -1062.197934
2.32 -1063.983696
2.34 -1065.73183
2.36 -1067.430216
2.38 -1069.07733
2.4 -1070.65617
2.42 -1072.203151
2.44 -1073.66999
2.46 -1075.087911
2.48 -1076.431205
2.5 -1077.712286
2.52 -1078.929926
2.54 -1080.080415
2.56 -1081.154955
2.58 -1082.14071
2.6 -1083.0537
2.62 -1083.877276
2.64 -1084.614006
2.66 -1085.291156
2.68 -1085.890718
2.7 -1086.431199
2.72 -1086.904206
2.74 -1087.306713
2.76 -1087.655476
2.77 -1087.805687
2.78 -1087.942434
2.79 -1088.066537
2.8 -1088.167443
2.81 -1088.264643
2.82 -1088.349063
2.83 -1088.416231
2.84 -1088.481484
2.85 -1088.530527
2.86 -1088.577998
2.87 -1088.605638
2.88 -1088.62841
2.89 -1088.638199
2.9 -1088.640171
2.91 -1088.630768
2.92 -1088.621723
2.93 -1088.602651
2.94071 -1088.570683
2.95 -1088.540892
2.96 -1088.50636
2.97 -1088.463236
2.98 -1088.41405
2.99 -1088.367545
3 -1088.31302
3.01 -1088.25003
3.02 -1088.18769
3.03 -1088.13481
3.04 -1088.06474
3.05 -1088.005696
3.06 -1087.935479
3.07 -1087.866723
3.08 -1087.804594
3.1 -1087.669303
3.12 -1087.529295
3.14 -1087.38405
3.16 -1087.240991
3.18 -1087.092895
3.2 -1086.955574
3.22 -1086.805261
3.24 -1086.655028
3.26 -1086.518322
3.28 -1086.36468
3.3 -1086.213376
3.32 -1086.061328
3.35 -1085.8127
3.38 -1085.55651
3.41 -1085.282171
3.44 -1084.783413
3.47 -1084.676517
3.5 -1084.355237
3.53 -1084.02303
3.56 -1083.676324
3.59 -1083.333467
3.62 -1082.992051
3.66 -1082.553911
3.7 -1082.112386
3.74 -1081.672979
3.78 -1081.237767
3.82 -1080.808863
3.86 -1080.377062
3.9 -1079.974347
3.94 -1079.600722
3.98 -1079.250442
4.02 -1078.929443
4.06 -1078.63086
4.1 -1078.343271

欲拟合的方程：f=c(1)*(exp(-2*c(2)*(tdata-c(3)))-2*exp(-c(2)*(tdata-c(3))))
c(1) ，c(2)， c(3)为想得到的参数；

以下是自己写的，结果和实际相差太大了：

M文件：
function f=fitfun1(c,tdata,ydata)
f=c(1)*(exp(-2*c(2)*(tdata-c(3)))-2*exp(-c(2)*(tdata-c(3))));

代码:
c0=[1 1 1];
c=lsqcurvefit('fitfun1',c0,tdata,ydata);
f=fitfun1(c,tdata);
plot(tdata,ydata,'o',tdata,ydata,tdata,f)
title('E-r curve fitting')
xlabel r(A)
ylabel E(eV)
求助其非线性拟合的程序代码，谢谢，这个困扰我很长时间了！！十分感谢！！

回复此楼