24小时热门版块排行榜    

查看: 1302  |  回复: 1

爱小燕

铜虫 (初入文坛)

[求助] [诚心求问!]利用matlab做动力学参数估计时遇到问题,求高手赐教!

利用Matlab做一个动力学参数的估计,指前因子和活化能的模拟,共9个独立反应18个反应参数需要拟合。优化计算用的是非线性回归中的最小二乘法进行模拟,解微分方程用的是ode45;程序运行很久,自动停止运行,出现以下错误:
Solver stopped prematurely.

lsqnonlin stopped because it exceeded the function evaluation limit,options.MaxFunEvals = 1800 (the default value).


ss =

   1.0e+07 *

  Columns 1 through 4

    0.1181    0.0000    2.1227    0.0001

  Columns 5 through 8

    0.5408    0.0770    0.0000    0.0075

  Columns 9 through 12

    0.7030    0.0000    0.0653    0.0000


rho2 =

   1.0e+06 *

  Columns 1 through 4

   -0.0157   -0.0000   -0.0465   -0.0000

  Columns 5 through 8

   -0.1278   -0.0528    0.0000   -0.0610

  Columns 9 through 12

   -3.0388    0.0000   -0.3056   -0.0000


F =

  Columns 1 through 4

   -1.3332   -0.4096   -1.3333   -0.8762

  Columns 5 through 8

   -1.3333   -1.3333    0.1523   -1.3333

  Columns 9 through 12

   -1.3333    0.5413   -1.3333   -1.2623



使用函数lsqnonlin()估计得到的参数值为:
Undefined function 'ci' for input arguments
of type 'double'.

Error in KINETIC0812>rho2_F (line 61)
fprintf('\tk1 = %.4f ±
%.4f\n',k(1),ci(1,2)-k(1))

Error in KINETIC0812 (line 30)
[rho2, F] = rho2_F(k,r,ss,Ne,Np);

新人金币较少,望见谅!
压缩文件中上传了部分文件,请指教!非常感谢!
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 相关文件.rar
  • 2014-08-16 13:45:56, 17.14 K

» 收录本帖的淘帖专辑推荐

Origin

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asedcft2010

新虫 (初入文坛)

你是用matlab求拟合参数吗?
2楼2015-03-01 16:32:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 爱小燕 的主题更新
信息提示
请填处理意见