| 查看: 390 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
wlx0914金虫 (小有名气)
|
[交流]
【求助】fminsearchi求解参数
|
||
|
根据xx和T的数据，关联方程参数a,b,c ，方程为 ln xx= a+b/T+c*ln T 下面是我的程序，提示有错，请高手指教 function y=myfun() clear all clc ExpData =[292.15 0.00848105 0.200113301 296.09 0.009329695 0.200113301 300.09 0.010580997 0.200113301 304.09 0.011993794 0.200113301 308.29 0.013486483 0.200113301 312.49 0.015286746 0.200113301 316.9 0.017209812 0.2001133 321.39 0.019224455 0.200113301 325.59 0.021446333 0.2001133 329.89 0.023818241 0.200113301 334.13 0.026117718 0.200113301 338.08 0.028576812 0.200113301 341.98 0.031254491 0.200113301 345.86 0.033917449 0.200113301 ]; T1= ExpData(:,1); xx1= ExpData(:,2); T=T1';xx=xx1'; %用fminsearch()—求解非线性Nelder-Mead单纯形法（非线性数据拟合）问题 A0= [1.22 1 1]; A=fminsearch(@Eqs,A0,[],T,xx) function f =Eqs(A,T,xx) % 模型方程 x1=exp(A(1)+A(2)/T+A(3)*log(T)); f=sum((xx-x1).^2) |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
最失望的一年
已经有10人回复
-
求推荐英文EI期刊
已经有5人回复
-
存款400万可以在学校里躺平吗
已经有27人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】fminsearchi求解参数问题
已经有2人回复
wlx0914
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1187.3
- 散金: 76
- 红花: 4
- 帖子: 180
- 在线: 38.6小时
- 虫号: 903494
- 注册: 2009-11-15
- 性别: GG
- 专业: 化工热力学和基础数据
4楼2010-11-08 14:35:43
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
2楼2010-11-08 10:45:31
change0618
铁杆木虫 (著名写手)
方丈大师
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★
wlx0914(金币+15): 2010-11-08 14:36:07
robert2020(金币+3):辛苦了！ 2010-11-08 23:00:36
wlx0914(金币+15): 2010-11-08 14:36:07
robert2020(金币+3):辛苦了！ 2010-11-08 23:00:36
|
function y=myfun() clear all clc ExpData =[292.15 0.00848105 0.200113301 296.09 0.009329695 0.200113301 300.09 0.010580997 0.200113301 304.09 0.011993794 0.200113301 308.29 0.013486483 0.200113301 312.49 0.015286746 0.200113301 316.9 0.017209812 0.2001133 321.39 0.019224455 0.200113301 325.59 0.021446333 0.2001133 329.89 0.023818241 0.200113301 334.13 0.026117718 0.200113301 338.08 0.028576812 0.200113301 341.98 0.031254491 0.200113301 345.86 0.033917449 0.200113301 ]; T1= ExpData(:,1); xx1= ExpData(:,2); T=T1';xx=xx1'; %用fminsearch()—求解非线性Nelder-Mead单纯形法（非线性数据拟合）问题 A0= [1.22 1 1]; A1=fminsearch(@Eqs,A0,[],T,xx); fprintf('\n\tA1 = \t%.4f\t%.4f\t%.4f\n',A1) options=optimset('TolFun',1e-10,'display','off'); A2= nlinfit(T,xx,@model,A0,options); TT=290:346; XX1=model(A1,TT); XX2=model(A2,TT); figure(1) plot(T,xx,'ro',TT,XX1,'b') grid on; figure(2) plot(T,xx,'ro',TT,XX2,'b') grid on; fprintf('\n\tA2 = \t%.4f\t%.4f\t%.4f\n',A2) function y=model(A,T) y=exp(A(1)+A(2)./T+A(3)*log(T)); function f =Eqs(A,T,xx) % 模型方程 x1=model(A,T); f=sum((xx-x1).^2) ; |
3楼2010-11-08 10:52:21