| 查看: 393 | 回复: 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) |
回复此楼» 猜你喜欢
网上报道青年教师午睡中猝死、熬夜猝死的越来越多，主要哪些原因引起的？
已经有5人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有6人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有16人回复
-
版面费该交吗
已经有13人回复
-
面上可以超过30页吧？
已经有6人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
什么是人一生最重要的？
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】fminsearchi求解参数问题
已经有2人回复
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
wlx0914
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1187.3
- 散金: 76
- 红花: 4
- 帖子: 180
- 在线: 38.6小时
- 虫号: 903494
- 注册: 2009-11-15
- 性别: GG
- 专业: 化工热力学和基础数据
4楼2010-11-08 14:35:43