| 查看: 517 | 回复: 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
» 猜你喜欢
国社科通讯评审
已经有6人回复
-
面上的水很深
已经有21人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有15人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
-
生生命学部会评
已经有4人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68161.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.2小时
- 虫号: 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