| 查看: 446 | 回复: 2 | |||
wangjunlitylg金虫 (小有名气)
|
[求助]
参数拟合：能运行但初值不进行迭代，输入什么值结果还是什么值
|
|
能运行的程序如下： function wenxian1 global k0 T s1 data=xlsread('D:文献.xlsx');T = data(:,1)';v= data(:,2); E0=232000;k0=8e12;kam=7130;s1=10/6;v0=0; beta0=[E0 kam]; [beta,resnorm,resid,exitflag,output,lambda] = lsqnonlin(@OptObjFunc,beta0,[],[],[],T,v0,v) vg=KineticsEqs(beta); figure,plot(T,v,'b-',T,vg,'k-') legend('Exp','Model') function f = OptObjFunc(beta,T,v0,v) vg = KineticsEqs(beta); f =v - vg'; function vg = KineticsEqs(beta) global a b c T s1 k0 d a=2*beta(2)^2; b=-k0*8.314/s1; c = 1/(beta(2)*sqrt(2*pi)); d=beta(1) ss=beta(1)+200*beta(2); xx=0; E=xx:ss; vg= quadv(@func,xx,ss); function y=func(E) global k0 T s1 a b c d y=k0./s1.*exp(-E./(8.314.*T)+b.*T.^2./E.*exp(-E./(8.314.*T))).*(c.*exp(-(E-d).^2./a)); 运行结果：还是原来附的初值，且模拟和实验差别很大，即好像不进行拟合， 我很努力了但就解决不了，谢谢大家，如果解决很感谢！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 文献.xlsx
2016-04-02 11:20:32, 10.81 K
» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2016-04-02 11:32:26
wangjunlitylg
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1048
- 红花: 2
- 帖子: 91
- 在线: 47.7小时
- 虫号: 337417
- 注册: 2007-04-02
- 性别: MM
- 专业: 传递过程
3楼2016-04-03 16:44:30