matlab非线性拟合模型参数,求助各位大神了
本人想用matlab进行非线性拟合以求未知解
输入数据请见上传文件
拟合方程为y=(1.22/A)*exp(-x/A*B)+(1336896.616*C+43559.7084*1.635328*10^(-5)*x^(1/2))*(1-exp(-1.155*10^(-8)*D*(x-(1-exp(-Ex)/E))))
未知参量为 A;B;C;D;E
麻烦大神帮忙写下详细代码,万分感谢
y x
-7.03148E-4 0
-6.85533E-4 0.05
-6.68564E-4 0.1
-6.54465E-4 0.15
-6.42952E-4 0.2
-6.33536E-4 0.25
-6.2581E-4 0.3
-6.19483E-4 0.35
-6.143E-4 0.4
-6.10075E-4 0.45
-6.06671E-4 0.5
-6.03992E-4 0.55
-6.01877E-4 0.6
-6.00284E-4 0.65
-5.99166E-4 0.7
-5.98402E-4 0.75
-5.97889E-4 0.8
-5.97629E-4 0.85
-5.97607E-4 0.9
-5.97723E-4 0.95
-5.97965E-4 1
-5.98326E-4 1.05
-5.9874E-4 1.1
-5.99241E-4 1.15
-5.99756E-4 1.2
-6.00359E-4 1.25
-6.00946E-4 1.3
-6.01532E-4 1.35
-6.02059E-4 1.4
-6.02614E-4 1.45
-6.0312E-4 1.5
-6.0361E-4 1.55
-6.04079E-4 1.6
-6.04518E-4 1.65
-6.04945E-4 1.7
-6.05336E-4 1.75
-6.05706E-4 1.8
-6.06037E-4 1.85
-6.06324E-4 1.9
-6.06625E-4 1.95
-6.06839E-4 2
-6.07046E-4 2.05
-6.07218E-4 2.1
-6.07361E-4 2.15
-6.07486E-4 2.2
-6.07558E-4 2.25
-6.0764E-4 2.3
-6.07713E-4 2.35
-6.07735E-4 2.4
-6.07769E-4 2.45
-6.07738E-4 2.5
-6.07718E-4 2.55
-6.0765E-4 2.6
-6.07607E-4 2.65
-6.07543E-4 2.7
-6.07474E-4 2.75
-6.0738E-4 2.8
-6.07262E-4 2.85
-6.07169E-4 2.9
-6.07063E-4 2.95
-6.0691E-4 3
-6.06798E-4 3.05
-6.06708E-4 3.1
-6.06552E-4 3.15
-6.06399E-4 3.2
-6.06277E-4 3.25
-6.06041E-4 3.3
-6.05741E-4 3.35
-6.0555E-4 3.4
-6.05372E-4 3.45
-6.05216E-4 3.5
-6.05038E-4 3.55
-6.04912E-4 3.6
-6.04708E-4 3.65
-6.0454E-4 3.7
-6.04402E-4 3.75
-6.04243E-4 3.8
-6.04057E-4 3.85
-6.03891E-4 3.9
-6.03692E-4 3.95
-6.03496E-4 4
-6.03327E-4 4.05
-6.03149E-4 4.1
-6.02946E-4 4.15
-6.02724E-4 4.2
-6.02552E-4 4.25
-6.02329E-4 4.3
-6.0215E-4 4.35
-6.0193E-4 4.4
-6.01792E-4 4.45
-6.01675E-4 4.5
-6.01474E-4 4.55
-6.01285E-4 4.6
-6.01006E-4 4.65
-6.0077E-4 4.7
-6.00546E-4 4.75
-6.0027E-4 4.8
-6.00033E-4 4.85
-5.99815E-4 4.9
-5.99604E-4 4.95
-5.99493E-4 5
-5.99278E-4 5.05
-5.9909E-4 5.1
-5.98923E-4 5.15
-5.98686E-4 5.2
-5.98524E-4 5.25
-5.98536E-4 5.3
-5.98327E-4 5.35
-5.9811E-4 5.4
-5.97914E-4 5.45
-5.97673E-4 5.5
-5.97495E-4 5.55
-5.97312E-4 5.6
-5.97149E-4 5.65
-5.97178E-4 5.7
-5.9698E-4 5.75
-5.96758E-4 5.8
-5.96506E-4 5.85
-5.96354E-4 5.9
-5.96131E-4 5.95
-5.9611E-4 6
-5.96012E-4 6.05
-5.9575E-4 6.1
-5.95553E-4 6.15
-5.95317E-4 6.2
-5.95253E-4 6.25
-5.95225E-4 6.3
-5.95161E-4 6.35
-5.95047E-4 6.4
-5.9491E-4 6.45
-5.94786E-4 6.5
-5.94683E-4 6.55
-5.94479E-4 6.6
-5.94359E-4 6.65
-5.94221E-4 6.7
-5.94055E-4 6.75
-5.93973E-4 6.8
-5.93788E-4 6.85
-5.93597E-4 6.9
-5.93423E-4 6.95
-5.93331E-4 7
-5.93182E-4 7.05
-5.92998E-4 7.1
-5.92746E-4 7.15
-5.92643E-4 7.2
-5.92523E-4 7.25
-5.92329E-4 7.3
-5.92135E-4 7.35
-5.92066E-4 7.4
-5.91821E-4 7.45
-5.91553E-4 7.5
-5.91448E-4 7.55
-5.91436E-4 7.6
-5.91176E-4 7.65
-5.91187E-4 7.7
-5.91225E-4 7.75
-5.91123E-4 7.8
-5.90912E-4 7.85
-5.90902E-4 7.9
-5.90841E-4 7.95
-5.9048E-4 8
-5.90285E-4 8.05
-5.90194E-4 8.1
-5.90031E-4 8.15
-5.89733E-4 8.2
-5.8967E-4 8.25
-5.89613E-4 8.3
-5.89314E-4 8.35
-5.89116E-4 8.4
-5.89341E-4 8.45
-5.89127E-4 8.5
-5.88863E-4 8.55
-5.88947E-4 8.6
-5.88796E-4 8.65
-5.88622E-4 8.7
-5.8842E-4 8.75
-5.88219E-4 8.8
-5.88092E-4 8.85
-5.87835E-4 8.9
-5.8773E-4 8.95
-5.8771E-4 9
-5.87594E-4 9.05
-5.87336E-4 9.1
-5.87196E-4 9.15
-5.87199E-4 9.2
-5.87127E-4 9.25
-5.86974E-4 9.3
-5.86895E-4 9.35
-5.86766E-4 9.4
-5.86825E-4 9.45
-5.86926E-4 9.5
-5.86831E-4 9.55
-5.86626E-4 9.6
-5.8653E-4 9.65
-5.86548E-4 9.7
-5.8638E-4 9.75
-5.8628E-4 9.8
-5.86209E-4 9.85
-5.86099E-4 9.9
-5.86076E-4 9.95 返回小木虫查看更多
看看这个例子
http://muchong.com/t-12537865-1
这个你照着写啊。
不收敛 就换一下初始值 试试。
你先将 x y数据 输入。
x=
y=
fun=@(p,x)(1.22/p(1).*exp(-x./(p(1).*p(2)))+(1336896.616*p(3)+0.4602*x.^(1/2).*(1-exp(-4.822*10.^(-9).*p(4).*(x-(1-exp(-p(5).*x)./p(5)))))
p0 = [1;1;1;1;1]; %% 初始值
p = lsqcurvefit(fun,p0,x,y)
括号什么的,你自己检查一下
,