当前位置: 首页 > 计算模拟 >matlab非线性拟合模型参数,求助各位大神了

matlab非线性拟合模型参数,求助各位大神了

作者 我是小璇璇
来源: 小木虫 500 10 举报帖子
+关注

本人想用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 返回小木虫查看更多

今日热帖
  • 精华评论
  • 独孤神宇

    引用回帖:
    4楼: Originally posted by 我是小璇璇 at 2019-11-22 14:31:11
    大神能不能给个详细的代码输入,感谢了...

    这个你照着写啊。

    不收敛 就换一下初始值 试试。

  • 独孤神宇

    引用回帖:
    6楼: Originally posted by 我是小璇璇 at 2019-11-23 12:40:32
    >> 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)))) ...

    你先将 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)

    括号什么的,你自己检查一下

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓