| 查看: 474 | 回复: 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
» 猜你喜欢
一志愿重庆大学085700资源与环境总分308求调剂
已经有6人回复
-
268求调剂
已经有7人回复
-
302求调剂
已经有11人回复
-
306求调剂
已经有4人回复
-
材料学学硕080502 337求调剂-一志愿华中科技大学
已经有4人回复
-
求调剂
已经有6人回复
-
313求调剂
已经有3人回复
-
南昌大学材料专硕311分求调剂
已经有6人回复
-
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.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
